/[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.16 by jmc, Tue Mar 6 17:10:29 2001 UTC revision 1.17 by jmc, Thu Mar 8 20:35:39 2001 UTC
# Line 9  C     | o Corrects the horizontal flow f Line 9  C     | o Corrects the horizontal flow f
9  C     |   slope.                                                 |  C     |   slope.                                                 |
10  C     \==========================================================/  C     \==========================================================/
11        SUBROUTINE CORRECTION_STEP( bi, bj, iMin, iMax, jMin, jMax,        SUBROUTINE CORRECTION_STEP( bi, bj, iMin, iMax, jMin, jMax,
12       I                       K, etaSurfX, etaSurfY,       I                       K, phiSurfX, phiSurfY,
13       I                       myCurrentTime, myThid )       I                       myCurrentTime, myThid )
14        IMPLICIT NONE        IMPLICIT NONE
15    
# Line 23  C     == Global variables == Line 23  C     == Global variables ==
23  #include "CG3D.h"  #include "CG3D.h"
24  #endif  #endif
25  C     == Routine Arguments ==  C     == Routine Arguments ==
26  C     etaSurfX, etaSurfY           - Surface slope  C     phiSurfX, phiSurfY           - Surface Potential gradient
27  C     bi,bj,iMin,iMax,jMin,jMax, K - Loop counters  C     bi,bj,iMin,iMax,jMin,jMax, K - Loop counters
28  C     myThid                       - Instance number for  C     myThid                       - Instance number for
29  C                                    this call to S/R CORRECTION_STEP  C                                    this call to S/R CORRECTION_STEP
30  C     myCurrentTime                - Current simulation time for this instance.  C     myCurrentTime                - Current simulation time for this instance.
31        _RL  etaSurfX(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL  phiSurfX(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
32        _RL  etaSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL  phiSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
33        INTEGER bi,bj,iMin,iMax,jMin,jMax        INTEGER bi,bj,iMin,iMax,jMin,jMax
34        INTEGER K        INTEGER K
35        INTEGER myThid        INTEGER myThid
# 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*implicSurfPress*etaSurfX(i,j)       &       -deltaTmom*hxFac*implicSurfPress*phiSurfX(i,j)
59  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
60       &       -deltaTmom*hx3dFac*gravity*_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))       &          (cg3d_x(i,j,k,bi,bj)-cg3d_x(i-1,j,k,bi,bj))
62  #endif  #endif
63       &                    )*_maskW(i,j,k,bi,bj)       &                    )*_maskW(i,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*implicSurfPress*etaSurfY(i,j)       &       -deltaTmom*hyFac*implicSurfPress*phiSurfY(i,j)
73  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
74       &       -deltaTmom*hy3dFac*gravity*_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))       &          (cg3d_x(i,j,k,bi,bj)-cg3d_x(i,j-1,k,bi,bj))
76  #endif  #endif
77       &                    )*_maskS(i,j,k,bi,bj)       &                    )*_maskS(i,j,k,bi,bj)
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22