/[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.21 by edhill, Thu Oct 9 04:19:18 2003 UTC revision 1.22 by jmc, Tue Nov 8 02:14:10 2005 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
 #include "PACKAGES_CONFIG.h"  
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6  CBOP  CBOP
# Line 9  C     !ROUTINE: CORRECTION_STEP Line 8  C     !ROUTINE: CORRECTION_STEP
8  C     !INTERFACE:  C     !INTERFACE:
9        SUBROUTINE CORRECTION_STEP( bi, bj, iMin, iMax, jMin, jMax,        SUBROUTINE CORRECTION_STEP( bi, bj, iMin, iMax, jMin, jMax,
10       I                       K, phiSurfX, phiSurfY,       I                       K, phiSurfX, phiSurfY,
11       I                       myCurrentTime, myThid )       I                       myTime, myThid )
12  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
13  C     *==========================================================*  C     *==========================================================*
14  C     | S/R CORRECTION_STEP                                        C     | S/R CORRECTION_STEP                                      
# Line 22  C     !USES: Line 21  C     !USES:
21        IMPLICIT NONE        IMPLICIT NONE
22  C     == Global variables ==  C     == Global variables ==
23  #include "SIZE.h"  #include "SIZE.h"
 #include "DYNVARS.h"  
24  #include "EEPARAMS.h"  #include "EEPARAMS.h"
25  #include "PARAMS.h"  #include "PARAMS.h"
26  #include "GRID.h"  #include "GRID.h"
27    #include "DYNVARS.h"
28  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
29  #include "SOLVE_FOR_PRESSURE3D.h"  #include "NH_VARS.h"
30  #endif  #endif
31    
32  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
# Line 36  C     phiSurfX, phiSurfY           - Sur Line 35  C     phiSurfX, phiSurfY           - Sur
35  C     bi,bj,iMin,iMax,jMin,jMax, K - Loop counters  C     bi,bj,iMin,iMax,jMin,jMax, K - Loop counters
36  C     myThid                       - Instance number for  C     myThid                       - Instance number for
37  C                                    this call to S/R CORRECTION_STEP  C                                    this call to S/R CORRECTION_STEP
38  C     myCurrentTime                - Current simulation time for this instance.  C     myTime                       - Current simulation time for this instance.
39        _RL  phiSurfX(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL  phiSurfX(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
40        _RL  phiSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL  phiSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
41        INTEGER bi,bj,iMin,iMax,jMin,jMax        INTEGER bi,bj,iMin,iMax,jMin,jMax
42        INTEGER K        INTEGER K
43        INTEGER myThid        INTEGER myThid
44        _RL myCurrentTime        _RL myTime
45    
46  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
47  C     == Local variables ==  C     == Local variables ==
# Line 71  C     Step forward zonal velocity Line 70  C     Step forward zonal velocity
70          uVel(i,j,k,bi,bj)=( gU(i,j,k,bi,bj)          uVel(i,j,k,bi,bj)=( gU(i,j,k,bi,bj)
71       &       -deltaTmom*hxFac*implicSurfPress*phiSurfX(i,j)       &       -deltaTmom*hxFac*implicSurfPress*phiSurfX(i,j)
72  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
 ceh3 needs an IF ( useNONHYDROSTATIC ) THEN  
73       &       -deltaTmom*hx3dFac*_recip_dxC(i,j,bi,bj)*       &       -deltaTmom*hx3dFac*_recip_dxC(i,j,bi,bj)*
74       &          (phi_nh(i,j,k,bi,bj)-phi_nh(i-1,j,k,bi,bj))       &          (phi_nh(i,j,k,bi,bj)-phi_nh(i-1,j,k,bi,bj))
75  #endif  #endif
# Line 85  C     Step forward meridional velocity Line 83  C     Step forward meridional velocity
83          vVel(i,j,k,bi,bj)=( gV(i,j,k,bi,bj)          vVel(i,j,k,bi,bj)=( gV(i,j,k,bi,bj)
84       &       -deltaTmom*hyFac*implicSurfPress*phiSurfY(i,j)       &       -deltaTmom*hyFac*implicSurfPress*phiSurfY(i,j)
85  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
 ceh3 needs an IF ( useNONHYDROSTATIC ) THEN  
86       &       -deltaTmom*hy3dFac*_recip_dyC(i,j,bi,bj)*       &       -deltaTmom*hy3dFac*_recip_dyC(i,j,bi,bj)*
87       &          (phi_nh(i,j,k,bi,bj)-phi_nh(i,j-1,k,bi,bj))       &          (phi_nh(i,j,k,bi,bj)-phi_nh(i,j-1,k,bi,bj))
88  #endif  #endif
Legend:
Removed from v.1.21  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22