/[MITgcm]/MITgcm/pkg/mom_common/mom_calc_visc.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_common/mom_calc_visc.F

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

revision 1.46 by jmc, Thu Aug 1 20:13:38 2013 UTC revision 1.47 by jmc, Fri Apr 4 20:08:11 2014 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "MOM_COMMON_OPTIONS.h"  #include "MOM_COMMON_OPTIONS.h"
5    #ifdef ALLOW_AUTODIFF
6    # include "AUTODIFF_OPTIONS.h"
7    #endif
8    
9  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
10  CBOP  CBOP
# Line 78  C     == Global variables == Line 81  C     == Global variables ==
81  #include "EEPARAMS.h"  #include "EEPARAMS.h"
82  #include "PARAMS.h"  #include "PARAMS.h"
83  #include "MOM_VISC.h"  #include "MOM_VISC.h"
84  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF
85  #include "tamc.h"  #include "tamc.h"
86  #include "tamc_keys.h"  #include "tamc_keys.h"
87  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF */
88    
89  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
90  C     myThid               :: my thread Id number  C     myThid               :: my thread Id number
# Line 162  c     IF ( useVariableVisc ) THEN Line 165  c     IF ( useVariableVisc ) THEN
165  C---- variable viscosity :  C---- variable viscosity :
166    
167         recip_dt = 1. _d 0         recip_dt = 1. _d 0
168         IF ( deltaTmom.NE.0. ) recip_dt = 1. _d 0/deltaTmom         IF ( deltaTMom.NE.0. ) recip_dt = 1. _d 0/deltaTMom
169    
170         IF ( useHarmonicVisc .AND. viscAhReMax.NE.0. ) THEN         IF ( useHarmonicVisc .AND. viscAhReMax.NE.0. ) THEN
171          viscAhRe_max=SQRT(2. _d 0)/viscAhReMax          viscAhRe_max=SQRT(2. _d 0)/viscAhReMax
# Line 306  C--   end if calcLeith Line 309  C--   end if calcLeith
309          DO i=2-OLx,sNx+OLx-1          DO i=2-OLx,sNx+OLx-1
310  CCCCCCCCCCCCCCC Divergence Point CalculationsCCCCCCCCCCCCCCCCCCCC  CCCCCCCCCCCCCCC Divergence Point CalculationsCCCCCCCCCCCCCCCCCCCC
311    
312  #ifdef    ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
313  # ifndef AUTODIFF_DISABLE_LEITH  # ifndef AUTODIFF_DISABLE_LEITH
314             lockey_2 = i+olx + (sNx+2*olx)*(j+oly-1)             lockey_2 = i+olx + (sNx+2*olx)*(j+oly-1)
315       &                      + (sNx+2*olx)*(sNy+2*oly)*(lockey_1-1)       &                      + (sNx+2*olx)*(sNy+2*oly)*(lockey_1-1)
# Line 315  CADJ &     = comlev1_mom_ijk_loop , key= Line 318  CADJ &     = comlev1_mom_ijk_loop , key=
318  CADJ STORE viscAh_ZSmg(i,j)  CADJ STORE viscAh_ZSmg(i,j)
319  CADJ &     = comlev1_mom_ijk_loop , key=lockey_2, byte=isbyte  CADJ &     = comlev1_mom_ijk_loop , key=lockey_2, byte=isbyte
320  # endif  # endif
321  #endif    /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
322    
323  C These are (powers of) length scales  C These are (powers of) length scales
324           L2 = L2_D(i,j,bi,bj)           L2 = L2_D(i,j,bi,bj)
Legend:
Removed from v.1.46  
changed lines
  Added in v.1.47
  ViewVC Help
Powered by ViewVC 1.1.22