changepd function

Update model fit statistics depending on calculation for pD