[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