<div>Hello, all,</div><div><br></div><div>I want to calculate the pressure gradient term in the momentum equation. In fact, I can directly diagnose this term using diagnostics pkg (which is called Um_dPdx & Vm_dPdy). But when I try to calculate this term, I found the result is not similar with the diagnostics output at all. </div><div><br></div><div>My detailed method is as follow:</div><div><br></div><div>dPdx=d(PH-PHI_SURF)dx, PH is the direct output from model (in my mind, this term is hydrostatic pressure anomal, the vertical integral of rho-rhoconst), PHI_SURF is a diagnostics term which mean surf pressure (like g*eta).</div><div><br></div><div>I am not sure where is wrong of this calculation. So could anyone help me? Thanks a lot.</div><div><br></div><div>Cheers,</div><div><br></div><div>Yang </div>