# beta-plane model for the MITgcm
# adapted from the MITgcm rotating tank tutorial

# equation
 &PARM01
 tRef=486*22.0,
 sRef=486*0.0,
 viscAh=2.1E-6,
 viscAz=2.1E-6,
 no_slip_sides=.TRUE.,
 no_slip_bottom=.TRUE.,
 diffKhT=1.3E-7,
 diffKzT=1.3E-7,
 f0=-2.,
 eosType='LINEAR',
 sBeta =0.,
 saltStepping=.FALSE.,
 gravity=-9.81,
 rhoConst=1043.0,
 rhoNil=1043.0,
 tAlpha=3.3E-4,
 heatCapacity_Cp=3764.0,
 rigidLid=.TRUE.,
 implicitFreeSurface=.FALSE.,
 nonHydrostatic=.TRUE.,
 metricTerms=.TRUE.,
 useCoriolis=.TRUE.,
 use3dCoriolis=.TRUE.,
 useEnergyConservingCoriolis=.TRUE.,
 vectorInvariantMomentum=.TRUE.,
 tempAdvScheme=7,
 staggerTimeStep=.TRUE.,
 readBinaryPrec=64,
 &

# solver
 &PARM02
 cg2dMaxIters=1000,
 cg2dTargetResidual=1.E-7,
 cg3dMaxIters=10,
 cg3dTargetResidual=1.E-9,
 &

# time
 &PARM03
 nIter0=0,
 nTimeSteps=2048000,
 deltaT=0.0015625,
 abEps=0.1,
 pChkptFreq=160.0,
 chkptFreq=80.0,
 dumpFreq=320.0,
 monitorSelect=2,
 monitorFreq=2048000,
 &

# grid
 &PARM04
 usingCylindricalGrid=.TRUE.,
 dXspacing=3.,
 dYspacing=0.0003125,
 delZ=486*0.0003125,
 ygOrigin=0.02375,
 &

# input
 &PARM05
 hydrogThetaFile='theta.bin',
 bathyFile='bathy.bin', 
 tCylIn  = 18.,
 tCylOut = 22.,
 &
