| 1 |
# ==================== |
| 2 |
# | Model parameters | |
| 3 |
# ==================== |
| 4 |
# |
| 5 |
# Continuous equation parameters |
| 6 |
&PARM01 |
| 7 |
tRef = 42*20., |
| 8 |
sRef = 42*37., |
| 9 |
bottomDragQuadratic = 0.002, |
| 10 |
viscAr=1.5E-4, |
| 11 |
# viscAr=1.5E-6, |
| 12 |
# This is the minimum viscosity used total=viscah+viscleith |
| 13 |
# Often U*deltaX/2 for grid noise (Griffies and Hallberg,2000) |
| 14 |
# viscAh=1.1E4, |
| 15 |
### viscAh=1.E3, |
| 16 |
# viscA4=3.E2, |
| 17 |
# viscA4=1.E4, |
| 18 |
# This is the standard Leith Viscosity, prop to grad(curl v_h) |
| 19 |
## viscC2leith=1., |
| 20 |
# viscC4leith=.2, |
| 21 |
# This is an additional Leith Viscosity, prop to grad(div.v_h) |
| 22 |
# It traps instabilities that have no signal in grad(curl v_h) |
| 23 |
## viscC2leithD=1., |
| 24 |
# viscC4leithD=.2, |
| 25 |
# These are the minimum viscosities used with Smagorinsky |
| 26 |
# total=viscastrain+viscsmag, etc. |
| 27 |
# viscAstrain=3.E1, |
| 28 |
# viscAtension=3.E1, |
| 29 |
# This is harmonic Smagorinsky Coefficient |
| 30 |
# 0.15=value from 3-d turbulence theory; |
| 31 |
# 3-4=stable oceanic values (Griffies and Hallberg, 2000) |
| 32 |
# viscC2Smag=4., |
| 33 |
# This is the max Viscosity used, be it Smag or Leith, and it scales |
| 34 |
# with grid size and timestep |
| 35 |
# Often 0.5 for stability (Griffies and Hallberg,2000) |
| 36 |
# viscAhGridMax=0.5, |
| 37 |
# viscA4GridMax=.5, |
| 38 |
# viscA4max = 1.E11, |
| 39 |
## useFullLeith=.TRUE., |
| 40 |
# viscA4=1.4E10, |
| 41 |
viscA4=1.2E9, |
| 42 |
# ********** DIFFUSION ************* |
| 43 |
#- diffKhT=.5E3, |
| 44 |
#- diffKhS=.5E3, |
| 45 |
diffK4T=1.5E9, |
| 46 |
diffK4S=1.5E9, |
| 47 |
diffKrNrT=3*3.E-5,12*3.E-5,1*1.E-5,1*7.E-6,1*2.E-6,1*8.E-7,1*4.E-7,22*1.E-7, |
| 48 |
diffKrNrS=3*3.E-5,12*3.E-5,1*1.E-5,1*7.E-6,1*2.E-6,1*8.E-7,1*4.E-7,22*1.E-7, |
| 49 |
staggerTimeStep = .TRUE. |
| 50 |
tempAdvScheme = 33, |
| 51 |
saltAdvScheme = 33, |
| 52 |
implicitDiffusion=.TRUE., |
| 53 |
implicitViscosity=.TRUE., |
| 54 |
vectorInvariantMomentum=.TRUE., |
| 55 |
nonLinFreeSurf = 0, |
| 56 |
#exactConserv = .TRUE., |
| 57 |
rhonil=1035., |
| 58 |
gravity=9.81, |
| 59 |
# eosType = 'MDJWF', |
| 60 |
eosType = 'JMD95Z', |
| 61 |
ivdc_kappa=1., |
| 62 |
useOldFreezing=.TRUE., |
| 63 |
# useRealFreshWaterFlux=.TRUE., |
| 64 |
useNHMTerms=.TRUE., |
| 65 |
# turn on looped cells |
| 66 |
hFacMin=.1, |
| 67 |
hFacMindr=30., |
| 68 |
# set precision of data files |
| 69 |
readBinaryPrec=32, |
| 70 |
useSingleCpuIO=.TRUE., |
| 71 |
#debugMode = .TRUE., |
| 72 |
& |
| 73 |
|
| 74 |
# Elliptic solver parameters |
| 75 |
&PARM02 |
| 76 |
cg2dMaxIters=500, |
| 77 |
cg2dTargetResidual=1.E-13, |
| 78 |
& |
| 79 |
|
| 80 |
# Time stepping parameters |
| 81 |
&PARM03 |
| 82 |
nIter0 = 0, |
| 83 |
nTimeSteps = 720, |
| 84 |
#- 366.days ; 32.days (1 yr of integration will yield a reasonable flow field) |
| 85 |
#endTime = 31622400., |
| 86 |
#endTime = 2764800., |
| 87 |
deltaTmom = 600.0, |
| 88 |
deltaTtracer= 600.0, |
| 89 |
deltaTClock = 600.0, |
| 90 |
# if you are using a version later than checkpoint45d on the main branch |
| 91 |
# you can uncomment the following line and increase the time step |
| 92 |
# deltaTtracer and deltaTClock to 172800.0 as well to speed up the |
| 93 |
# asynchronous time stepping |
| 94 |
# deltaTfreesurf = 172800.0, |
| 95 |
abEps = 0.1, |
| 96 |
#pChkptFreq= 2628000., |
| 97 |
dumpFreq= 864000., |
| 98 |
taveFreq= 2628000., |
| 99 |
# taveFreq= 31536000., |
| 100 |
monitorFreq= 86400., |
| 101 |
monitorFreq= 3600., |
| 102 |
# 1 months restoring timescale for temperature |
| 103 |
# tauThetaClimRelax = 2592000., |
| 104 |
# 5 days restoring timescale for salinity |
| 105 |
# tauSaltClimRelax = 432000., |
| 106 |
# periodicExternalForcing=.TRUE., |
| 107 |
# externForcingPeriod=86400., |
| 108 |
# externForcingCycle=31104000., |
| 109 |
nTimeSteps = 36, |
| 110 |
monitorFreq= 1., |
| 111 |
& |
| 112 |
|
| 113 |
# Gridding parameters |
| 114 |
&PARM04 |
| 115 |
usingCartesianGrid=.FALSE., |
| 116 |
usingSphericalPolarGrid=.TRUE., |
| 117 |
delR = 10., 10., 10., 20., 20., |
| 118 |
20., 30., 30., 30., 40., |
| 119 |
40., 40., 40., 40., 40., |
| 120 |
40., 40., 50., 50., 50., |
| 121 |
60., 60., 70., 70., 80., |
| 122 |
80., 90., 90.,100.,120., |
| 123 |
130.,150.,150.,200.,200., |
| 124 |
200.,200.,200.,200.,300., |
| 125 |
300.,300., |
| 126 |
ygOrigin = 33.0625000000000000000000000000000, |
| 127 |
xgOrigin = -8.35416666666666666666666666666667, |
| 128 |
dySpacing=.041666666666666666666667, |
| 129 |
dxSpacing=.041666666666666666666667, |
| 130 |
& |
| 131 |
|
| 132 |
# Input datasets |
| 133 |
&PARM05 |
| 134 |
bathyFile= 'bathy_MITgcm.bin', |
| 135 |
hydrogThetaFile='init_temp.bin', |
| 136 |
hydrogSaltFile= 'init_sal.bin', |
| 137 |
# zonalWindFile= 'wind_stress_x.bin', |
| 138 |
# meridWindFile= 'wind_stress_y.bin', |
| 139 |
# thetaClimFile= 'lev_sst.bin', |
| 140 |
# saltClimFile= 'lev_sss.bin', |
| 141 |
# surfQFile= 'ncep_qnet.bin', |
| 142 |
# EmPmRFile= 'ncep_emp.bin', |
| 143 |
& |
Parent Directory
| 
Revision Log
| 
Revision Graph