Po_surf = rHalf(Ks) + (rHalf(Ks+1)-rHalf(Ks))*
& *(Pfld(i,j,bi,bj)-rC(ks))/(rHalf(2*ks-1)-rHalf(2*ks))
& *(Pfld(i,j,bi,bj)-rC(ks))/(rHalf(2*ks+1)-rHalf(2*ks))
_RL hRef(2*Nr+1), rHalf(2*Nr+1)
rHalf(2*K-1) = rF(K)
rHalf(2*K) = rC(K)
rHalf(2*Nr+1) = rF(Nr+1)
ddPI=atm_Cp*( ((rHalf( k )/atm_Po)**atm_kappa)
& -((rHalf(k+1)/atm_Po)**atm_kappa) )
& ' K=',k*0.5,' ; r=',rHalf(k),' ; phiRef/g=',
_RL rHalf(2*Nr+1)
rHalf(2*k-1) = rF(k)
rHalf(2*k) = rC(k)