/[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.17 by jmc, Thu Mar 8 20:35:39 2001 UTC revision 1.18 by adcroft, Fri Jun 29 17:14:49 2001 UTC
# Line 20  C     == Global variables == Line 20  C     == Global variables ==
20  #include "PARAMS.h"  #include "PARAMS.h"
21  #include "GRID.h"  #include "GRID.h"
22  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
23  #include "CG3D.h"  #include "SOLVE_FOR_PRESSURE3D.h"
24  #endif  #endif
25  C     == Routine Arguments ==  C     == Routine Arguments ==
26  C     phiSurfX, phiSurfY           - Surface Potential gradient  C     phiSurfX, phiSurfY           - Surface Potential gradient
# Line 58  C     Step forward zonal velocity Line 58  C     Step forward zonal velocity
58       &       -deltaTmom*hxFac*implicSurfPress*phiSurfX(i,j)       &       -deltaTmom*hxFac*implicSurfPress*phiSurfX(i,j)
59  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
60       &       -deltaTmom*hx3dFac*_recip_dxC(i,j,bi,bj)*       &       -deltaTmom*hx3dFac*_recip_dxC(i,j,bi,bj)*
61       &          (cg3d_x(i,j,k,bi,bj)-cg3d_x(i-1,j,k,bi,bj))       &          (phi_nh(i,j,k,bi,bj)-phi_nh(i-1,j,k,bi,bj))
62  #endif  #endif
63       &                    )*_maskW(i,j,k,bi,bj)       &                    )*_maskW(i,j,k,bi,bj)
64          gUNm1(i,j,k,bi,bj)=gU(i,j,k,bi,bj)          gUNm1(i,j,k,bi,bj)=gU(i,j,k,bi,bj)
# Line 72  C     Step forward meridional velocity Line 72  C     Step forward meridional velocity
72       &       -deltaTmom*hyFac*implicSurfPress*phiSurfY(i,j)       &       -deltaTmom*hyFac*implicSurfPress*phiSurfY(i,j)
73  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
74       &       -deltaTmom*hy3dFac*_recip_dyC(i,j,bi,bj)*       &       -deltaTmom*hy3dFac*_recip_dyC(i,j,bi,bj)*
75       &          (cg3d_x(i,j,k,bi,bj)-cg3d_x(i,j-1,k,bi,bj))       &          (phi_nh(i,j,k,bi,bj)-phi_nh(i,j-1,k,bi,bj))
76  #endif  #endif
77       &                    )*_maskS(i,j,k,bi,bj)       &                    )*_maskS(i,j,k,bi,bj)
78          gVNm1(i,j,k,bi,bj)=gV(i,j,k,bi,bj)          gVNm1(i,j,k,bi,bj)=gV(i,j,k,bi,bj)
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22