/[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.8 by cnh, Thu Aug 20 19:25:05 1998 UTC revision 1.9 by cnh, Thu Aug 20 20:05:00 1998 UTC
# Line 19  C     \================================= Line 19  C     \=================================
19  #include "GRID.h"  #include "GRID.h"
20  #include "CG2D.h"  #include "CG2D.h"
21  C     == Routine Arguments ==  C     == Routine Arguments ==
22  C     etaSurfX, etaSurfY - Surface slope  C     etaSurfX, etaSurfY           - Surface slope
23  C     bi,bj,iMin,iMax,jMin,jMax, K - Loop counters  C     bi,bj,iMin,iMax,jMin,jMax, K - Loop counters
24  C     myThid - Instance number for this call to S/R CORRECTION_STEP  C     myThid                       - Instance number for
25  C     myCurrentTime - Current simulation time for this instance.  C                                    this call to S/R CORRECTION_STEP
26    C     myCurrentTime                - Current simulation time for this instance.
27        _RL  etaSurfX(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL  etaSurfX(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
28        _RL  etaSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL  etaSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
29        INTEGER bi,bj,iMin,iMax,jMin,jMax        INTEGER bi,bj,iMin,iMax,jMin,jMax
# Line 42  C     Step forward zonal velocity Line 43  C     Step forward zonal velocity
43        DO j=jMin,jMax        DO j=jMin,jMax
44         DO i=iMin,iMax         DO i=iMin,iMax
45          uVel(i,j,k,bi,bj)=( gUNm1(i,j,k,bi,bj)          uVel(i,j,k,bi,bj)=( gUNm1(i,j,k,bi,bj)
46       &       -deltaTmom*hxFac*gBaro *etaSurfX(i,j)       &       -deltaTmom*hxFac*gBaro*etaSurfX(i,j)
47       &                    )*_maskW(i,j,k,bi,bj)       &                    )*_maskW(i,j,k,bi,bj)
48          gUNm1(i,j,k,bi,bj)=gU(i,j,k,bi,bj)          gUNm1(i,j,k,bi,bj)=gU(i,j,k,bi,bj)
49         ENDDO         ENDDO
# Line 52  C     Step forward meridional velocity Line 53  C     Step forward meridional velocity
53        DO j=jMin,jMax        DO j=jMin,jMax
54         DO i=iMin,iMax         DO i=iMin,iMax
55          vVel(i,j,k,bi,bj)=( gVNm1(i,j,k,bi,bj)          vVel(i,j,k,bi,bj)=( gVNm1(i,j,k,bi,bj)
56       &       -deltaTmom*hyFac*gBaro *etaSurfY(i,j)       &       -deltaTmom*hyFac*gBaro*etaSurfY(i,j)
57       &                    )*_maskS(i,j,k,bi,bj)       &                    )*_maskS(i,j,k,bi,bj)
58          gVNm1(i,j,k,bi,bj)=gV(i,j,k,bi,bj)          gVNm1(i,j,k,bi,bj)=gV(i,j,k,bi,bj)
59         ENDDO         ENDDO
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22