Warning, /verification/1D_ocean_ice_column/input_ad/data is written in an unsupported language. File is not indexed.
view on githubraw file Latest commit 11c3150c on 2021-10-18 18:00:35 UTC
8f749e3979 Ian *0001 # ====================
0002 # | Model parameters |
0003 # ====================
0004 #
0005 # Continuous equation parameters
0006 #
0007 # tRef - Reference vertical potential temperature (deg C)
0008 # sRef - Reference vertical salinity (PSU)
0009 # viscAh - Horizontal eddy viscosity coefficient (m^2/s)
0010 # viscAz - Vertical eddy viscosity coefficient (m^2/s)
0011 # diffKhT - Horizontal temperature diffusivity (m^2/s)
0012 # diffKzT - Vertical temperature diffusivity (m^2/s)
0013 # diffKhS - Horizontal salt diffusivity (m^2/s)
0014 # diffKzS - Vertical salt diffusivity (m^2/s)
0015 # f0 - Reference coriolis parameter,
0016 # south edge of f on beta plane (1/s)
0017 # beta - df/dy (s^-1.m^-1)
0018 # tAlpha - Linear EOS thermal expansion coefficient (1/oC)
0019 # sBeta - Linear EOS haline contraction coefficient (1/ppt)
0020 # gravity - Acceleration due to gravity (m/s^2)
0021 # gBaro - Accel. due to gravity used in barotropic equation (m/s^2)
0022 # rigidLid - Set to true to use rigid lid
0023 # implicitFreeSurface - Set to true to use implicit free surface
0024 # eosType - Flag for linear or polynomial equation of state
0025 # momAdvection - On/Off flag for momentum self transport
0026 # momViscosity - On/Off flag for momentum mixing
0027 #
0028 &PARM01
973b42676b Ian *0029 # an's T,S profiles @ 23 levels
657cdb5f3f Ian *0030 tRef = -1.930,-1.931,-1.932,-1.2854,-0.6,-0.6,-0.9,-1.2,-1.2,-1.3,
8f749e3979 Ian *0031 -1.4,-1.4,-1.5,-1.5,-1.5,-1.5,-1.3,-0.9,-0.3, 0.2,
63584c545e Jean*0032 0.5, 0.5, 0.5,
5c05ca2ad9 Ian *0033 # tRef to use for 11,000 time step adjoint calculations
266294ca66 Ian *0034 # tRef = -1.8420,-1.8464,-1.8392,-1.2854,-0.6,-0.6,-0.9,-1.2,-1.2,-1.3,
0035 # -1.4,-1.4,-1.5,-1.5,-1.5,-1.5,-1.3,-0.9,-0.3, 0.2,
0036 # 0.5, 0.5, 0.5
8f749e3979 Ian *0037 sRef = 29.0079,29.0080,29.0086,29.0775,30.7,31.6,31.9,32.1,32.3,32.4,
0038 32.5,32.7,32.8,32.9,33.1,33.4,33.8,34.2,34.5,34.7,
63584c545e Jean*0039 34.8,34.8,34.8,
8f749e3979 Ian *0040 no_slip_sides=.FALSE.,
0041 no_slip_bottom=.TRUE.,
0042 viscAz=1.93e-5,
0043 viscAh=5.E4,
0044 diffKhT=0.0,
0045 diffKzT=1.46e-7,
0046 diffKhS=0.0,
0047 diffKzS=1.46e-7,
bdf2bf0470 Ian *0048 # beta=1.E-11,
05b7bb78c9 Jean*0049 # 1-D setups work best with an f-plane due to unequal Coriolis force on north
0050 # and south faces
0051 selectCoriMap=0,
8f749e3979 Ian *0052 rigidLid=.FALSE.,
0053 implicitFreeSurface=.TRUE.,
0054 eosType='JMD95Z',
0055 readBinaryPrec=32,
0056 writeBinaryPrec=32,
0057 saltStepping=.TRUE.,
0058 tempStepping=.TRUE.,
0059 momStepping=.TRUE.,
0060 implicitDiffusion=.TRUE.,
0061 implicitViscosity=.TRUE.,
0062 allowFreezing=.FALSE.,
05b7bb78c9 Jean*0063 #- set wrong celsius2K to reproduce old results:
0064 celsius2K=273.16,
0065 HeatCapacity_Cp = 3986.0,
8f749e3979 Ian *0066 gravity = 9.8156,
05b7bb78c9 Jean*0067 rhoConst = 1027.0,
0e64d12c66 Ian *0068 rhoConstFresh = 999.8,
8f749e3979 Ian *0069 useCDscheme=.FALSE.,
0070 staggerTimeStep=.TRUE.,
0071 multiDimAdvection=.false.,
0072 tempAdvScheme=30,
0073 saltAdvScheme=30,
0074 vectorInvariantMomentum=.TRUE.,
63584c545e Jean*0075 useRealFreshWaterFlux = .TRUE.,
11c3150c71 Mart*0076 useSingleCpuIO = .TRUE.,
973b42676b Ian *0077 # debugLevel=1,
8f749e3979 Ian *0078 &
0079
0080 # Elliptic solver parameters
0081 #
0082 # cg2dMaxIters - Maximum number of 2d solver iterations
0083 # cg2dTargetResidual - Solver target residual
0084 #
0085 &PARM02
0086 cg2dMaxIters=1000,
0087 cg2dTargetResidual=1.D-13,
0088 &
0089
0090 # Time stepping parameters
0091 #
0092 # startTime - Integration starting time (s)
0093 # endTime - Integration ending time (s)
0094 # tauCD - CD scheme coupling timescale (s)
0095 # deltaTMom - Timestep for momemtum equations (s)
0096 # deltaTtracer - Tracer timestep (s)
0097 # deltaTClock - Timestep used as model "clock" (s)
0098 # abEps - Adams-Bashforth stabilising factor
0099 # pChkPtFreq - Frequency of permanent check pointing (s)
0100 # chkPtFreq - Frequency of rolling check pointing (s)
0101 # dumpFreq - Frequency at which model state is stored (s)
0102 # tauThetaClimRelax - Relaxation to climatology time scale (s)
0103 # tauSaltClimRelax - Relaxation to climatology time scale (s)
0104 #
0105 &PARM03
0106 startTime=0.0,
657cdb5f3f Ian *0107 nTimeSteps= 10,
973b42676b Ian *0108 # nTimeSteps= 11000,
8f749e3979 Ian *0109 deltaTtracer=3600.0,
0110 deltaTClock =3600.0,
0111 cAdjFreq=0.,
0112 abEps=0.1,
e37fc501c5 Jean*0113 tracForcingOutAB=1,
8f749e3979 Ian *0114 pChkptFreq=0.,
0115 chkptFreq= 0.,
0116 dumpFreq = 0.,
e37fc501c5 Jean*0117 taveFreq = 0.,
973b42676b Ian *0118 # monitorFreq=86400.,
0119 # adjDumpFreq=86400.,
0120 # adjMonitorFreq = 86400.,
e37fc501c5 Jean*0121 monitorFreq=1.,
0122 adjDumpFreq=1.,
0123 adjMonitorFreq = 1.,
8f749e3979 Ian *0124 &
0125
0126 # Gridding parameters
0127 #
0128 # usingSphericalPolarGrid - On/Off flag for spherical polar coordinates
0129 # delX - Zonal grid spacing (degrees)
0130 # delY - Meridional grid spacing (degrees)
0131 # delZ - Vertical grid spacing (m)
0132 # ygOrigin - Southern boundary latitude (degrees)
0133 #
0134 &PARM04
bdf2bf0470 Ian *0135 usingCartesianGrid=.TRUE.,
0136 dXspacing=5000.,
0137 dYspacing=5000.,
0138 # usingSphericalPolarGrid=.TRUE.,
0139 # delX=1*2.E0,
0140 # delY=1*2.E0,
8f749e3979 Ian *0141 delZ = 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.01,
63584c545e Jean*0142 10.03, 10.11, 10.32, 10.80, 11.76, 13.42, 16.04 , 19.82, 24.85,
0143 31.10, 38.42, 46.50, 55.00, 63.50, 71.58,
bdf2bf0470 Ian *0144 # ygOrigin=65.,
0145 # xgOrigin=280.,
0146 # rSphere = 6371.D3,
8f749e3979 Ian *0147 &
0148
0149 &PARM05
973b42676b Ian *0150 bathyFile = 'bathy_1x1_1105m_testpool',
8f749e3979 Ian *0151 &
