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

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

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

revision 1.76 by heimbach, Mon Aug 13 18:05:26 2001 UTC revision 1.79 by adcroft, Thu Aug 16 17:12:24 2001 UTC
# Line 58  C     myThid - Thread number for this in Line 58  C     myThid - Thread number for this in
58        INTEGER myThid        INTEGER myThid
59    
60  C     == Local variables  C     == Local variables
 C     maskUp                   o maskUp: land/water mask for W points  
61  C     fVer[STUV]               o fVer: Vertical flux term - note fVer  C     fVer[STUV]               o fVer: Vertical flux term - note fVer
62  C                                      is "pipelined" in the vertical  C                                      is "pipelined" in the vertical
63  C                                      so we need an fVer for each  C                                      so we need an fVer for each
# Line 77  C     bi, bj Line 76  C     bi, bj
76  C     k, kup,        - Index for layer above and below. kup and kDown  C     k, kup,        - Index for layer above and below. kup and kDown
77  C     kDown, km1       are switched with layer to be the appropriate  C     kDown, km1       are switched with layer to be the appropriate
78  C                      index into fVerTerm.  C                      index into fVerTerm.
 C     tauAB - Adams-Bashforth timestepping weight: 0=forward ; 1/2=Adams-Bashf.  
       _RS maskUp  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
79        _RL fVerU   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)        _RL fVerU   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)
80        _RL fVerV   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)        _RL fVerV   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)
81        _RL phiHyd  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL phiHyd  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
# Line 91  C     tauAB - Adams-Bashforth timesteppi Line 88  C     tauAB - Adams-Bashforth timesteppi
88        _RL sigmaX  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL sigmaX  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
89        _RL sigmaY  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL sigmaY  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
90        _RL sigmaR  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL sigmaR  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
       _RL tauAB  
91    
92  C This is currently used by IVDC and Diagnostics  C This is currently used by IVDC and Diagnostics
93        _RL ConvectCount (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL ConvectCount (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
# Line 100  C This is currently used by IVDC and Dia Line 96  C This is currently used by IVDC and Dia
96        INTEGER jMin, jMax        INTEGER jMin, jMax
97        INTEGER bi, bj        INTEGER bi, bj
98        INTEGER i, j        INTEGER i, j
99        INTEGER k, km1, kup, kDown        INTEGER k, km1, kp1, kup, kDown
100    
101  Cjmc : add for phiHyd output <- but not working if multi tile per CPU  Cjmc : add for phiHyd output <- but not working if multi tile per CPU
102  c     CHARACTER*(MAX_LEN_MBUF) suff  c     CHARACTER*(MAX_LEN_MBUF) suff
# Line 240  C       (note: this loop will be replace Line 236  C       (note: this loop will be replace
236       I         myThid )                               I         myThid )                        
237          ENDIF          ENDIF
238    
239    #ifdef  INCLUDE_CALC_DIFFUSIVITY_CALL
240    C--      Calculate the total vertical diffusivity
241            DO k=1,Nr
242             CALL CALC_VISCOSITY(
243         I        bi,bj,iMin,iMax,jMin,jMax,k,
244         O        KappaRU,KappaRV,
245         I        myThid)
246           ENDDO
247    #endif
248    
249  C--     Start of dynamics loop  C--     Start of dynamics loop
250          DO k=1,Nr          DO k=1,Nr
251    
# Line 248  C--       kup    Cycles through 1,2 to p Line 254  C--       kup    Cycles through 1,2 to p
254  C--       kDown  Cycles through 2,1 to point to current layer  C--       kDown  Cycles through 2,1 to point to current layer
255    
256            km1  = MAX(1,k-1)            km1  = MAX(1,k-1)
257              kp1  = MIN(k+1,Nr)
258            kup  = 1+MOD(k+1,2)            kup  = 1+MOD(k+1,2)
259            kDown= 1+MOD(k,2)            kDown= 1+MOD(k,2)
260    
# Line 272  C        distinguishe between Stagger an Line 279  C        distinguishe between Stagger an
279       I        myThid )       I        myThid )
280           ENDIF           ENDIF
281    
 #ifdef  INCLUDE_CALC_DIFFUSIVITY_CALL  
 C--      Calculate the total vertical diffusivity  
          CALL CALC_VISCOSITY(  
      I        bi,bj,iMin,iMax,jMin,jMax,k,  
      I        maskUp,  
      O        KappaRU,KappaRV,  
      I        myThid)  
 #endif  
   
282  C--      Calculate accelerations in the momentum equations (gU, gV, ...)  C--      Calculate accelerations in the momentum equations (gU, gV, ...)
283  C        and step forward storing the result in gUnm1, gVnm1, etc...  C        and step forward storing the result in gUnm1, gVnm1, etc...
284           IF ( momStepping ) THEN           IF ( momStepping ) THEN
285             CALL CALC_MOM_RHS(  #ifndef DISABLE_MOM_FLUXFORM
286               IF (.NOT. vectorInvariantMomentum) CALL MOM_FLUXFORM(
287       I         bi,bj,iMin,iMax,jMin,jMax,k,kup,kDown,       I         bi,bj,iMin,iMax,jMin,jMax,k,kup,kDown,
288       I         phiHyd,KappaRU,KappaRV,       I         phiHyd,KappaRU,KappaRV,
289       U         fVerU, fVerV,       U         fVerU, fVerV,
290       I         myTime, myThid)       I         myTime, myThid)
291    #endif
292    #ifndef DISABLE_MOM_VECINV
293               IF (vectorInvariantMomentum) CALL MOM_VECINV(
294         I         bi,bj,iMin,iMax,jMin,jMax,k,kup,kDown,
295         I         phiHyd,KappaRU,KappaRV,
296         U         fVerU, fVerV,
297         I         myTime, myThid)
298    #endif
299             CALL TIMESTEP(             CALL TIMESTEP(
300       I         bi,bj,iMin,iMax,jMin,jMax,k,       I         bi,bj,iMin,iMax,jMin,jMax,k,
301       I         phiHyd, phiSurfX, phiSurfY,       I         phiHyd, phiSurfX, phiSurfY,
# Line 395  Cjmc(end) Line 402  Cjmc(end)
402         ENDDO         ENDDO
403        ENDDO        ENDDO
404    
405  #ifndef EXCLUDE_DEBUGMODE  #ifndef DISABLE_DEBUGMODE
406        If (debugMode) THEN        If (debugMode) THEN
407         CALL DEBUG_STATS_RL(1,EtaN,'EtaN (DYNAMICS)',myThid)         CALL DEBUG_STATS_RL(1,EtaN,'EtaN (DYNAMICS)',myThid)
408         CALL DEBUG_STATS_RL(Nr,uVel,'Uvel (DYNAMICS)',myThid)         CALL DEBUG_STATS_RL(Nr,uVel,'Uvel (DYNAMICS)',myThid)
Legend:
Removed from v.1.76  
changed lines
  Added in v.1.79
  ViewVC Help
Powered by ViewVC 1.1.22