/[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.10 by cnh, Fri Nov 6 22:44:45 1998 UTC revision 1.11 by adcroft, Mon Mar 22 15:54:03 1999 UTC
# Line 18  C     \================================= Line 18  C     \=================================
18  #include "PARAMS.h"  #include "PARAMS.h"
19  #include "GRID.h"  #include "GRID.h"
20  #include "CG2D.h"  #include "CG2D.h"
21    #ifdef ALLOW_NONHYDROSTATIC
22    #include "CG3D.h"
23    #endif
24  C     == Routine Arguments ==  C     == Routine Arguments ==
25  C     etaSurfX, etaSurfY           - Surface slope  C     etaSurfX, etaSurfY           - Surface slope
26  C     bi,bj,iMin,iMax,jMin,jMax, K - Loop counters  C     bi,bj,iMin,iMax,jMin,jMax, K - Loop counters
# Line 33  C     myCurrentTime                - Cur Line 36  C     myCurrentTime                - Cur
36    
37  C     == Local variables ==  C     == Local variables ==
38        INTEGER i,j        INTEGER i,j
39        _RL hxFac,hyFac,rRhoNil        _RL hxFac,hyFac
40          _RL hx3dFac,hy3dFac
41    
42  C     On/off scaling paramters  C     On/off scaling paramters
43        hxFac = pfFacMom        hxFac = pfFacMom
44        hyFac = pfFacMom        hyFac = pfFacMom
45          IF ( nonHydrostatic ) THEN
46            hx3dFac = pfFacMom
47            hy3dFac = pfFacMom
48          ELSE
49            hx3dFac = 0.
50            hy3dFac = 0.
51          ENDIF
52    
53  C     Step forward zonal velocity  C     Step forward zonal velocity
54        DO j=jMin,jMax        DO j=jMin,jMax
55         DO i=iMin,iMax         DO i=iMin,iMax
56          uVel(i,j,k,bi,bj)=( gUNm1(i,j,k,bi,bj)          uVel(i,j,k,bi,bj)=( gUNm1(i,j,k,bi,bj)
57       &       -deltaTmom*hxFac*gBaro*etaSurfX(i,j)       &       -deltaTmom*hxFac*gBaro*etaSurfX(i,j)
58    #ifdef ALLOW_NONHYDROSTATIC
59         &       -deltaTmom*hx3dFac*gravity*_recip_dxC(i,j,bi,bj)*
60         &          (cg3d_x(i,j,k,bi,bj)-cg3d_x(i-1,j,k,bi,bj))
61    #endif
62       &                    )*_maskW(i,j,k,bi,bj)       &                    )*_maskW(i,j,k,bi,bj)
63          gUNm1(i,j,k,bi,bj)=gU(i,j,k,bi,bj)          gUNm1(i,j,k,bi,bj)=gU(i,j,k,bi,bj)
64         ENDDO         ENDDO
# Line 54  C     Step forward meridional velocity Line 69  C     Step forward meridional velocity
69         DO i=iMin,iMax         DO i=iMin,iMax
70          vVel(i,j,k,bi,bj)=( gVNm1(i,j,k,bi,bj)          vVel(i,j,k,bi,bj)=( gVNm1(i,j,k,bi,bj)
71       &       -deltaTmom*hyFac*gBaro*etaSurfY(i,j)       &       -deltaTmom*hyFac*gBaro*etaSurfY(i,j)
72    #ifdef ALLOW_NONHYDROSTATIC
73         &       -deltaTmom*hy3dFac*gravity*_recip_dyC(i,j,bi,bj)*
74         &          (cg3d_x(i,j,k,bi,bj)-cg3d_x(i,j-1,k,bi,bj))
75    #endif
76       &                    )*_maskS(i,j,k,bi,bj)       &                    )*_maskS(i,j,k,bi,bj)
77          gVNm1(i,j,k,bi,bj)=gV(i,j,k,bi,bj)          gVNm1(i,j,k,bi,bj)=gV(i,j,k,bi,bj)
78         ENDDO         ENDDO
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22