GetPdynPara_BC function

Compute P-dynamics parameters using the bias correction method from BRW (2012)