MLAdjust/input.AhStTn/data

021f4f33e7 Jean*0001 # ====================
                0002 # | Model parameters |
                0003 # ====================
                0004 #
                0005 # Continuous equation parameters
                0006  &PARM01
                0007  viscAr=5.E-5,
                0008 #- This is the background viscosity used; total=viscah+viscleith
                0009 # viscAh=1.E2,
                0010 # viscA4=1.E2,
                0011 #- This is the background viscosity used on vorticity points
                0012 # viscAhZ=1.E2,
                0013 # viscA4Z=1.E2,
                0014 #- This is the background viscosity used on divergence points
                0015 # viscAhD=1.E2,
                0016 # viscA4D=1.E2,
                0017 #- viscAhD and viscAhZ are set to viscah if they are unset initially
                0018 #
                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,
                0032 # This is harmonic Smagorinsky Coefficient
                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)
                0043  useFullLeith=.TRUE.,
                0044 #- This uses the Strain-Tension form for calculating viscous terms
                0045  useStrainTensionVisc=.TRUE.,
                0046 #no_slip_sides=.FALSE.,
                0047  no_slip_bottom=.FALSE.,
                0048  vectorInvariantMomentum=.TRUE.,
                0049  ivdc_kappa=10.,
                0050  implicitDiffusion=.TRUE.,
                0051  implicitViscosity=.TRUE.,
                0052  f0=7.29E-5,
                0053  beta=0.E-11,
                0054  useJamartWetPoints=.TRUE.,
                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,
                0061  rhoNil=1035.,
                0062  implicitFreeSurface=.TRUE.,
                0063  staggerTimeStep=.TRUE.,
                0064  tempAdvScheme=33,
                0065  saltAdvScheme=33,
                0066  readBinaryPrec=64,
                0067  writeBinaryPrec=64,
                0068  &
                0069
                0070 # Elliptic solver parameters
                0071  &PARM02
                0072  cg2dMaxIters=1000,
                0073  cg2dTargetResidual=1.E-11,
                0074  cg3dMaxIters=40,
                0075  cg3dTargetResidual=1.E-9,
                0076  &
                0077
                0078 # Time stepping parameters
                0079  &PARM03
                0080  nIter0=0,
                0081 #endTime=864000.,
                0082 #- for testing, only run for 12 iterations:
                0083  nTimeSteps=12,
                0084  deltaT=1200.,
                0085  abEps=0.01,
                0086 #- for testing, reduce amount of output by commenting out output frequencies:
                0087 # pchkptFreq=43200.,
                0088 #taveFreq=7200.,
                0089 #dumpFreq=7200.,
                0090 #diagFreq=7200.,
                0091  monitorSelect=2,
                0092  monitorFreq=1200.,
                0093  &
                0094
                0095 # Gridding parameters
                0096  &PARM04
                0097  usingCartesianGrid=.TRUE.,
                0098  delR= 40*5.0,
                0099  dxSpacing=1.E3,
                0100  dySpacing=1.E3,
                0101  &
                0102
                0103 # Input datasets
                0104  &PARM05
                0105  bathyFile='topo_sl.bin',
                0106  hydrogThetaFile='thetaInitial.bin',
                0107  hydrogSaltFile='spiceInitial.bin',
                0108  &