| 1 |
jmc |
1.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