/[MITgcm]/MITgcm/model/src/correction_step.F
ViewVC logotype

Diff of /MITgcm/model/src/correction_step.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.14 by cnh, Sun Feb 4 14:38:46 2001 UTC revision 1.15 by jmc, Tue Feb 20 15:06:21 2001 UTC
# Line 56  C     Step forward zonal velocity Line 56  C     Step forward zonal velocity
56        DO j=jMin,jMax        DO j=jMin,jMax
57         DO i=iMin,iMax         DO i=iMin,iMax
58          uVel(i,j,k,bi,bj)=( gUNm1(i,j,k,bi,bj)          uVel(i,j,k,bi,bj)=( gUNm1(i,j,k,bi,bj)
59       &       -deltaTmom*hxFac*gBaro*etaSurfX(i,j)       &       -deltaTmom*hxFac*gBaro*implicSurfPress*etaSurfX(i,j)
60  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
61       &       -deltaTmom*hx3dFac*gravity*_recip_dxC(i,j,bi,bj)*       &       -deltaTmom*hx3dFac*gravity*_recip_dxC(i,j,bi,bj)*
62       &          (cg3d_x(i,j,k,bi,bj)-cg3d_x(i-1,j,k,bi,bj))       &          (cg3d_x(i,j,k,bi,bj)-cg3d_x(i-1,j,k,bi,bj))
# Line 70  C     Step forward meridional velocity Line 70  C     Step forward meridional velocity
70        DO j=jMin,jMax        DO j=jMin,jMax
71         DO i=iMin,iMax         DO i=iMin,iMax
72          vVel(i,j,k,bi,bj)=( gVNm1(i,j,k,bi,bj)          vVel(i,j,k,bi,bj)=( gVNm1(i,j,k,bi,bj)
73       &       -deltaTmom*hyFac*gBaro*etaSurfY(i,j)       &       -deltaTmom*hyFac*gBaro*implicSurfPress*etaSurfY(i,j)
74  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
75       &       -deltaTmom*hy3dFac*gravity*_recip_dyC(i,j,bi,bj)*       &       -deltaTmom*hy3dFac*gravity*_recip_dyC(i,j,bi,bj)*
76       &          (cg3d_x(i,j,k,bi,bj)-cg3d_x(i,j-1,k,bi,bj))       &          (cg3d_x(i,j,k,bi,bj)-cg3d_x(i,j-1,k,bi,bj))
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22