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 |
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 |
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 |
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) |
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) |