[MITgcm-support] correction_step.F geometry question
Haferman, Jeffrey (Jeff) (CIV)
jlhaferm at nps.edu
Thu Apr 19 12:48:24 EDT 2007
Hello all -
In routine "correction_step.F" consider the following snippet of code:
C Step forward meridional velocity
DO j=jMin,jMax
DO i=iMin,iMax
vVel(i,j,k,bi,bj)=( gV(i,j,k,bi,bj)
& -deltaTmom*hyFac*implicSurfPress*phiSurfY(i,j)
#ifdef ALLOW_NONHYDROSTATIC
& -deltaTmom*hy3dFac*_recip_dyC(i,j,bi,bj)*
& (phi_nh(i,j,k,bi,bj)-phi_nh(i,j-1,k,bi,bj))
#endif
& )*_maskS(i,j,k,bi,bj)
ENDDO
ENDDO
My understanding is that the non-hydrostatic pressure phi_nh is defined
at the cell-center, but that vVel is defined at the cell interfaces?
So, vVel at point j would be defined at a point between the reference
points for phi_nh(j) and phi_nh(j-1)... is this correct?
Jeff
More information about the MITgcm-support
mailing list