/[MITgcm]/MITgcm/model/src/timestep.F
ViewVC logotype

Diff of /MITgcm/model/src/timestep.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.33 by edhill, Thu Oct 9 04:19:18 2003 UTC revision 1.34 by edhill, Tue Oct 28 22:57:59 2003 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4    #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
6    
7  CBOP  CBOP
# Line 48  C     == Local variables == Line 49  C     == Local variables ==
49        _RL phxFac,phyFac, psFac        _RL phxFac,phyFac, psFac
50        _RL     gUtmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL     gUtmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
51        _RL     gVtmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL     gVtmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
52  #ifdef INCLUDE_CD_CODE  #ifdef ALLOW_CD_CODE
53        _RL     guCor(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL     guCor(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
54        _RL     gvCor(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL     gvCor(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
55  #endif  #endif
# Line 86  C-    Initialize local arrays (not reall Line 87  C-    Initialize local arrays (not reall
87         DO i=1-Olx,sNx+Olx         DO i=1-Olx,sNx+Olx
88          gUtmp(i,j) = 0. _d 0          gUtmp(i,j) = 0. _d 0
89          gVtmp(i,j) = 0. _d 0          gVtmp(i,j) = 0. _d 0
90  #ifdef INCLUDE_CD_CODE  #ifdef ALLOW_CD_CODE
91          guCor(i,j) = 0. _d 0          guCor(i,j) = 0. _d 0
92          gvCor(i,j) = 0. _d 0          gvCor(i,j) = 0. _d 0
93  #endif  #endif
# Line 154  C-     for CD-scheme, compute gU,Vtmp = Line 155  C-     for CD-scheme, compute gU,Vtmp =
155          ENDDO          ENDDO
156        ENDIF        ENDIF
157    
158  #ifdef INCLUDE_CD_CODE  #ifdef ALLOW_CD_CODE
 #ifdef ALLOW_MOM_FLUXFORM  
159        IF (useCDscheme) THEN        IF (useCDscheme) THEN
160  C-     Step forward D-grid velocity using C-grid gU,Vtmp = gU,V^n + forcing  C-     Step forward D-grid velocity using C-grid gU,Vtmp = gU,V^n + forcing
161  C      and return coriolis terms on C-grid (guCor,gvCor)  C      and return coriolis terms on C-grid (guCor,gvCor)
162          CALL MOM_CDSCHEME(          CALL CD_SCHEME(
163       I                    bi,bj,k, dPhiHydX,dPhiHydY, gUtmp,gVtmp,       I                  bi,bj,k, dPhiHydX,dPhiHydY, gUtmp,gVtmp,
164       O                    guCor,gvCor,       O                  guCor,gvCor,
165       I                    myTime, myIter, myThid)       I                  myTime, myIter, myThid)
166          DO j=jMin,jMax          DO j=jMin,jMax
167           DO i=iMin,iMax           DO i=iMin,iMax
168            gUtmp(i,j) = gU(i,j,k,bi,bj)            gUtmp(i,j) = gU(i,j,k,bi,bj)
# Line 172  C      and return coriolis terms on C-gr Line 172  C      and return coriolis terms on C-gr
172           ENDDO           ENDDO
173          ENDDO          ENDDO
174        ENDIF        ENDIF
175  #endif /* ALLOW_MOM_FLUXFORM */  #endif /* ALLOW_CD_CODE */
 #endif /* INCLUDE_CD_CODE */  
176    
177  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
178        IF (.NOT. vectorInvariantMomentum        IF (.NOT. vectorInvariantMomentum
Legend:
Removed from v.1.33  
changed lines
  Added in v.1.34
  ViewVC Help
Powered by ViewVC 1.1.22