/[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.13 by adcroft, Fri Feb 2 21:04:48 2001 UTC revision 1.16 by jmc, Tue Mar 6 17:10:29 2001 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
# Line 18  C     == Global variables == Line 19  C     == Global variables ==
19  #include "EEPARAMS.h"  #include "EEPARAMS.h"
20  #include "PARAMS.h"  #include "PARAMS.h"
21  #include "GRID.h"  #include "GRID.h"
 #include "CG2D.h"  
22  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
23  #include "CG3D.h"  #include "CG3D.h"
24  #endif  #endif
# Line 55  C     Step forward zonal velocity Line 55  C     Step forward zonal velocity
55        DO j=jMin,jMax        DO j=jMin,jMax
56         DO i=iMin,iMax         DO i=iMin,iMax
57          uVel(i,j,k,bi,bj)=( gUNm1(i,j,k,bi,bj)          uVel(i,j,k,bi,bj)=( gUNm1(i,j,k,bi,bj)
58       &       -deltaTmom*hxFac*gBaro*etaSurfX(i,j)       &       -deltaTmom*hxFac*gBaro*implicSurfPress*etaSurfX(i,j)
59  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
60       &       -deltaTmom*hx3dFac*gravity*_recip_dxC(i,j,bi,bj)*       &       -deltaTmom*hx3dFac*gravity*_recip_dxC(i,j,bi,bj)*
61       &          (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 69  C     Step forward meridional velocity Line 69  C     Step forward meridional velocity
69        DO j=jMin,jMax        DO j=jMin,jMax
70         DO i=iMin,iMax         DO i=iMin,iMax
71          vVel(i,j,k,bi,bj)=( gVNm1(i,j,k,bi,bj)          vVel(i,j,k,bi,bj)=( gVNm1(i,j,k,bi,bj)
72       &       -deltaTmom*hyFac*gBaro*etaSurfY(i,j)       &       -deltaTmom*hyFac*gBaro*implicSurfPress*etaSurfY(i,j)
73  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
74       &       -deltaTmom*hy3dFac*gravity*_recip_dyC(i,j,bi,bj)*       &       -deltaTmom*hy3dFac*gravity*_recip_dyC(i,j,bi,bj)*
75       &          (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.13  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22