Warning, /verification/1D_ocean_ice_column/input/data is written in an unsupported language. File is not indexed.
view on githubraw file Latest commit e37fc501 on 2012-07-04 18:09:14 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
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,
84248c5b54 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.,
bcd6fab2ff Ian *0076 # debugLevel=1,
8f749e3979 Ian *0077 &
0078
0079 # Elliptic solver parameters
0080 #
0081 # cg2dMaxIters - Maximum number of 2d solver iterations
0082 # cg2dTargetResidual - Solver target residual
0083 #
0084 &PARM02
0085 cg2dMaxIters=1000,
0086 cg2dTargetResidual=1.D-13,
0087 &
0088
0089 # Time stepping parameters
0090 #
0091 # startTime - Integration starting time (s)
0092 # endTime - Integration ending time (s)
0093 # tauCD - CD scheme coupling timescale (s)
0094 # deltaTMom - Timestep for momemtum equations (s)
0095 # deltaTtracer - Tracer timestep (s)
0096 # deltaTClock - Timestep used as model "clock" (s)
0097 # abEps - Adams-Bashforth stabilising factor
0098 # pChkPtFreq - Frequency of permanent check pointing (s)
0099 # chkPtFreq - Frequency of rolling check pointing (s)
0100 # dumpFreq - Frequency at which model state is stored (s)
0101 # tauThetaClimRelax - Relaxation to climatology time scale (s)
0102 # tauSaltClimRelax - Relaxation to climatology time scale (s)
0103 #
0104 &PARM03
0105 startTime=0.0,
378593a5ba Ian *0106 nTimeSteps= 10,
0107 # nTimeSteps= 11000,
8f749e3979 Ian *0108 deltaTtracer=3600.0,
0109 deltaTClock =3600.0,
0110 cAdjFreq=0.,
0111 abEps=0.1,
e37fc501c5 Jean*0112 tracForcingOutAB=1,
8f749e3979 Ian *0113 pChkptFreq=0.,
0114 chkptFreq= 0.,
0115 dumpFreq = 0.,
e37fc501c5 Jean*0116 taveFreq = 0.,
0117 # monitorFreq=86400.,
b5c7abc5d3 Ian *0118 # adjDumpFreq=86400.,
0119 # adjMonitorFreq = 86400.,
e37fc501c5 Jean*0120 monitorFreq=1.,
0121 adjDumpFreq=1.,
0122 adjMonitorFreq = 1.,
8f749e3979 Ian *0123 &
0124
0125 # Gridding parameters
0126 #
0127 # usingSphericalPolarGrid - On/Off flag for spherical polar coordinates
0128 # delX - Zonal grid spacing (degrees)
0129 # delY - Meridional grid spacing (degrees)
0130 # delZ - Vertical grid spacing (m)
0131 # ygOrigin - Southern boundary latitude (degrees)
0132 #
0133 &PARM04
bdf2bf0470 Ian *0134 usingCartesianGrid=.TRUE.,
0135 dXspacing=5000.,
0136 dYspacing=5000.,
0137 # usingSphericalPolarGrid=.TRUE.,
0138 # delX=1*2.E0,
0139 # delY=1*2.E0,
8f749e3979 Ian *0140 delZ = 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.01,
63584c545e Jean*0141 10.03, 10.11, 10.32, 10.80, 11.76, 13.42, 16.04 , 19.82, 24.85,
0142 31.10, 38.42, 46.50, 55.00, 63.50, 71.58,
bdf2bf0470 Ian *0143 # ygOrigin=65.,
0144 # xgOrigin=280.,
0145 # rSphere = 6371.D3,
8f749e3979 Ian *0146 &
0147
0148 &PARM05
b5c7abc5d3 Ian *0149 bathyFile = 'bathy_1x1_1105m_testpool',
8f749e3979 Ian *0150 &
0151
