/[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.78 by heimbach, Tue Aug 14 00:19:42 2001 UTC revision 1.79 by adcroft, Thu Aug 16 17:12:24 2001 UTC
# Line 76  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.  
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 89  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 284  C        distinguishe between Stagger an Line 282  C        distinguishe between Stagger an
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.78  
changed lines
  Added in v.1.79
  ViewVC Help
Powered by ViewVC 1.1.22