MITgcm

479fcf10d8 Bayl*0001 # ====================
                0002 # | Model parameters |
                0003 # ====================
                0004 #
                0005 # Continuous equation parameters
                0006  &PARM01
f8b1cc36a1 Bayl*0007  viscAr=5.E-5,
6d286f5656 Jean*0008 #- This is the background viscosity used; total=viscah+viscleith
f8b1cc36a1 Bayl*0009 # viscAh=1.E2,
                0010 # viscA4=1.E2,
6d286f5656 Jean*0011 #- This is the background viscosity used on vorticity points
f8b1cc36a1 Bayl*0012 # viscAhZ=1.E2,
                0013 # viscA4Z=1.E2,
6d286f5656 Jean*0014 #- This is the background viscosity used on divergence points
f8b1cc36a1 Bayl*0015 # viscAhD=1.E2,
                0016 # viscA4D=1.E2,
6d286f5656 Jean*0017 #- viscAhD and viscAhZ are set to viscah if they are unset initially
f8b1cc36a1 Bayl*0018 #
6d286f5656 Jean*0019 #- This is a maximum gridscale Reynolds number limiter
                0020 # viscAhReMax=10.,
                0021 # viscA4ReMax=10.,
                0022 #- This is a minimum autoscaling viscosity
                0023 # viscAhGridMin=.01,
                0024 # viscA4GridMin=.01,
                0025 #- This is the standard Leith Viscosity, prop to grad(curl v_h)
                0026 #viscC2leith = 1.85,
                0027  viscC4leith = 1.85,
                0028 # -This is an additional Leith Viscosity, prop to grad(div.v_h)
                0029 #  It traps instabilities that have no signal in grad(curl v_h)
                0030 #viscC2leithD = 1.85,
                0031  viscC4leithD = 1.85,
479fcf10d8 Bayl*0032 # This is harmonic Smagorinsky Coefficient
6d286f5656 Jean*0033 # 0.2-0.9 = value from 3-d turbulence theory; Smagorinsky (1993)
                0034 # 2.2-4   = stable oceanic values (Griffies and Hallberg, 2000)
                0035 # 1. smag is approximately the same as 1.85 leith+leithd for this grid spacing
                0036 #viscC2smag = 1.,
                0037  viscC4smag = 1.,
                0038 #- This is the max Viscosity used, be it Smag or Leith, and it scales
                0039 #  with grid size and timestep
                0040 #viscAhGridMax = 1.,
                0041  viscA4GridMax = 1.,
                0042 #- This uses the full grad(vort) and grad(div)
479fcf10d8 Bayl*0043  useFullLeith=.TRUE.,
6d286f5656 Jean*0044 #- This uses the Strain-Tension form for calculating viscous terms
                0045 #useStrainTensionVisc=.TRUE.,
479fcf10d8 Bayl*0046  no_slip_sides=.FALSE.,
                0047  no_slip_bottom=.FALSE.,
6d286f5656 Jean*0048  vectorInvariantMomentum=.TRUE.,
479fcf10d8 Bayl*0049  ivdc_kappa=10.,
                0050  implicitDiffusion=.TRUE.,
                0051  implicitViscosity=.TRUE.,
6d286f5656 Jean*0052  f0=7.29E-5,
479fcf10d8 Bayl*0053  beta=0.E-11,
6d286f5656 Jean*0054  useJamartWetPoints=.TRUE.,
479fcf10d8 Bayl*0055  eosType='LINEAR',
                0056 # Using Temperature as only Active Tracer
                0057  tAlpha=2.E-4,
                0058 # Using Salinity as Spice
                0059  sBeta =0.E-4,
                0060  gravity=9.81,
6d286f5656 Jean*0061  rhoNil=1035.,
479fcf10d8 Bayl*0062  implicitFreeSurface=.TRUE.,
                0063  staggerTimeStep=.TRUE.,
                0064  tempAdvScheme=33,
                0065  saltAdvScheme=33,
6d286f5656 Jean*0066  readBinaryPrec=64,
                0067  writeBinaryPrec=64,
7a2d61f89c Jean*0068  &
479fcf10d8 Bayl*0069
                0070 # Elliptic solver parameters
                0071  &PARM02
                0072  cg2dMaxIters=1000,
                0073  cg2dTargetResidual=1.E-11,
                0074  cg3dMaxIters=40,
                0075  cg3dTargetResidual=1.E-9,
7a2d61f89c Jean*0076  &
479fcf10d8 Bayl*0077
                0078 # Time stepping parameters
                0079  &PARM03
7a2d61f89c Jean*0080  nIter0=0,
846d98fdd8 Jean*0081 #endTime=864000.,
4bee903946 Jean*0082 #- for testing, only run for 12 iterations:
91be9bd3f4 Bayl*0083  nTimeSteps=12,
                0084  deltaT=1200.,
479fcf10d8 Bayl*0085  abEps=0.01,
4bee903946 Jean*0086 #- for testing, reduce amount of output by commenting out output frequencies:
f8b1cc36a1 Bayl*0087 # pchkptFreq=43200.,
4bee903946 Jean*0088 #taveFreq=7200.,
91be9bd3f4 Bayl*0089 #dumpFreq=7200.,
4bee903946 Jean*0090 #diagFreq=7200.,
2f6f180016 Jean*0091  monitorSelect=2,
91be9bd3f4 Bayl*0092  monitorFreq=1200.,
7a2d61f89c Jean*0093  &
479fcf10d8 Bayl*0094
                0095 # Gridding parameters
                0096  &PARM04
                0097  usingCartesianGrid=.TRUE.,
7a2d61f89c Jean*0098  delR= 40*5.0,
                0099  dxSpacing=1.E3,
                0100  dySpacing=1.E3,
                0101  &
479fcf10d8 Bayl*0102
                0103 # Input datasets
                0104  &PARM05
                0105  bathyFile='topo_sl.bin',
                0106  hydrogThetaFile='thetaInitial.bin',
                0107  hydrogSaltFile='spiceInitial.bin',
7a2d61f89c Jean*0108  &