/[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.161 by jmc, Mon Mar 5 18:21:12 2012 UTC revision 1.162 by jmc, Sun Mar 18 22:19:45 2012 UTC
# Line 168  C     guDissip   :: dissipation tendency Line 168  C     guDissip   :: dissipation tendency
168  C     gvDissip   :: dissipation tendency (all explicit terms), v component  C     gvDissip   :: dissipation tendency (all explicit terms), v component
169  C     KappaRU    :: vertical viscosity  C     KappaRU    :: vertical viscosity
170  C     KappaRV    :: vertical viscosity  C     KappaRV    :: vertical viscosity
171  C     iMin, iMax     - Ranges and sub-block indices on which calculations  C     iMin, iMax :: Ranges and sub-block indices on which calculations
172  C     jMin, jMax       are applied.  C     jMin, jMax    are applied.
173  C     bi, bj  C     bi, bj     :: tile indices
174  C     k, kup,        - Index for layer above and below. kup and kDown  C     k          :: current level index
175  C     kDown, km1       are switched with layer to be the appropriate  C     km1, kp1   :: index of level above (k-1) and below (k+1)
176  C                      index into fVerTerm.  C     kUp, kDown :: Index for interface above and below. kUp and kDown are
177    C                   are switched with k to be the appropriate index into fVerU,V
178        _RL fVerU   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)        _RL fVerU   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)
179        _RL fVerV   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)        _RL fVerV   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)
180        _RL phiHydF (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL phiHydF (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
# Line 191  C                      index into fVerTe Line 192  C                      index into fVerTe
192        INTEGER jMin, jMax        INTEGER jMin, jMax
193        INTEGER bi, bj        INTEGER bi, bj
194        INTEGER i, j        INTEGER i, j
195        INTEGER k, km1, kp1, kup, kDown        INTEGER k, km1, kp1, kUp, kDown
196    
197  #ifdef ALLOW_DIAGNOSTICS  #ifdef ALLOW_DIAGNOSTICS
198        LOGICAL dPhiHydDiagIsOn        LOGICAL dPhiHydDiagIsOn
# Line 521  C--      Calculate accelerations in the Line 522  C--      Calculate accelerations in the
522  C        and step forward storing the result in gU, gV, etc...  C        and step forward storing the result in gU, gV, etc...
523           IF ( momStepping ) THEN           IF ( momStepping ) THEN
524  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
525  # if (defined NONLIN_FRSURF) && (defined ALLOW_MOM_FLUXFORM)  # ifdef NONLIN_FRSURF
526  #  ifndef DISABLE_RSTAR_CODE  #  if (defined ALLOW_MOM_FLUXFORM) && !(defined DISABLE_RSTAR_CODE)
527  CADJ STORE dWtransC(:,:,bi,bj)  CADJ STORE dWtransC(:,:,bi,bj)
528  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte
529  CADJ STORE dWtransU(:,:,bi,bj)  CADJ STORE dWtransU(:,:,bi,bj)
# Line 530  CADJ &     = comlev1_bibj_k, key=kkey, b Line 531  CADJ &     = comlev1_bibj_k, key=kkey, b
531  CADJ STORE dWtransV(:,:,bi,bj)  CADJ STORE dWtransV(:,:,bi,bj)
532  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte
533  #  endif  #  endif
534  # endif  CADJ STORE fVerU(:,:,:)
535  #endif  CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte
536    CADJ STORE fVerV(:,:,:)
537    CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte
538    # endif /* NONLIN_FRSURF */
539    #endif /* ALLOW_AUTODIFF_TAMC */
540             IF (.NOT. vectorInvariantMomentum) THEN             IF (.NOT. vectorInvariantMomentum) THEN
541  #ifdef ALLOW_MOM_FLUXFORM  #ifdef ALLOW_MOM_FLUXFORM
 C  
542                CALL MOM_FLUXFORM(                CALL MOM_FLUXFORM(
543       I         bi,bj,iMin,iMax,jMin,jMax,k,kup,kDown,       I         bi,bj,k,iMin,iMax,jMin,jMax,
544       I         KappaRU, KappaRV,       I         KappaRU, KappaRV,
545       U         fVerU, fVerV,       U         fVerU(1-OLx,1-OLy,kUp),   fVerV(1-OLx,1-OLy,kUp),
546         O         fVerU(1-OLx,1-OLy,kDown), fVerV(1-OLx,1-OLy,kDown),
547       O         guDissip, gvDissip,       O         guDissip, gvDissip,
548       I         myTime, myIter, myThid)       I         myTime, myIter, myThid)
549  #endif  #endif
550             ELSE             ELSE
551  #ifdef ALLOW_MOM_VECINV  #ifdef ALLOW_MOM_VECINV
 C  
 # ifdef ALLOW_AUTODIFF_TAMC  
 #  ifdef NONLIN_FRSURF  
 CADJ STORE fVerU(:,:,:)  
 CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte  
 CADJ STORE fVerV(:,:,:)  
 CADJ &     = comlev1_bibj_k, key=kkey, byte=isbyte  
 #  endif  
 # endif /* ALLOW_AUTODIFF_TAMC */  
 C  
552               CALL MOM_VECINV(               CALL MOM_VECINV(
553       I         bi,bj,iMin,iMax,jMin,jMax,k,kup,kDown,       I         bi,bj,k,iMin,iMax,jMin,jMax,
554       I         KappaRU, KappaRV,       I         KappaRU, KappaRV,
555       U         fVerU, fVerV,       I         fVerU(1-OLx,1-OLy,kUp),   fVerV(1-OLx,1-OLy,kUp),
556         O         fVerU(1-OLx,1-OLy,kDown), fVerV(1-OLx,1-OLy,kDown),
557       O         guDissip, gvDissip,       O         guDissip, gvDissip,
558       I         myTime, myIter, myThid)       I         myTime, myIter, myThid)
559  #endif  #endif
Legend:
Removed from v.1.161  
changed lines
  Added in v.1.162
  ViewVC Help
Powered by ViewVC 1.1.22