dPIdp = (atm_Rd/atm_Po)*
Bo_surf(i,j,bi,bj) = dPIdp*thetaConst
Bo_surf(i,j,bi,bj) = dPIdp*tRef(kSurfC(i,j,bi,bj))
_RL dPIdp