| 1 |
#:::==================== |
| 2 |
# | Model parameters | |
| 3 |
# ==================== |
| 4 |
# |
| 5 |
# Continuous equation parameters |
| 6 |
&PARM01 |
| 7 |
# momstepping=.FALSE., |
| 8 |
# tempstepping=.FALSE., |
| 9 |
# saltstepping=.FALSE., |
| 10 |
Tref = 90*-1.9, |
| 11 |
Sref = 90*34.4, |
| 12 |
viscAr=1.E-3, |
| 13 |
#viscAh=600.0, |
| 14 |
viscAhGrid=0.2, |
| 15 |
viscA4Grid=0.02, |
| 16 |
no_slip_sides=.FALSE., |
| 17 |
no_slip_bottom=.FALSE., |
| 18 |
diffKhT= 100.0, |
| 19 |
#diffKhT= 600.0, |
| 20 |
diffKrT=5.E-5, |
| 21 |
diffKhS= 100.0, |
| 22 |
#diffKhS= 600.0, |
| 23 |
diffKrS=5.E-5, |
| 24 |
bottomDragQuadratic=2.5E-3, |
| 25 |
staggerTimestep = .true., |
| 26 |
tempAdvScheme=30, |
| 27 |
saltAdvScheme=30, |
| 28 |
vectorInvariantMomentum = .true., |
| 29 |
### momImplVertAdv=.true., |
| 30 |
tempImplVertAdv=.true., |
| 31 |
saltImplVertAdv=.true., |
| 32 |
eosType='JMD95Z', |
| 33 |
#tAlpha = 2.E-4, |
| 34 |
#sBeta = 7.4E-4, |
| 35 |
HeatCapacity_cp = 3974.0, |
| 36 |
rhoConst=1000., |
| 37 |
rhoNil=1000., |
| 38 |
gravity=9.81, |
| 39 |
gBaro=9.81, |
| 40 |
implicitDiffusion = .true., |
| 41 |
implicitViscosity = .true., |
| 42 |
rigidLid=.FALSE., |
| 43 |
implicitFreeSurface=.TRUE., |
| 44 |
useRealFreshWaterFlux = .true., |
| 45 |
exactConserv=.TRUE., |
| 46 |
# start nlfs with rstar |
| 47 |
#select_rStar=2, |
| 48 |
nonlinFreeSurf=4, |
| 49 |
hFacInf=0.2, |
| 50 |
hFacSup=2.0, |
| 51 |
# end nlfs with rstar |
| 52 |
hFacMin=0.2, |
| 53 |
#integr_GeoPot = 1, |
| 54 |
#ivdc_kappa = 0.005., |
| 55 |
#implicitDiffusion=.true., |
| 56 |
#implicitViscosity=.true., |
| 57 |
#useCDscheme=.true., |
| 58 |
#nonHydrostatic=.true., |
| 59 |
#allowFreezing = .true., |
| 60 |
readBinaryPrec=64, |
| 61 |
#writeBinaryPrec=64, |
| 62 |
debuglevel = 1, |
| 63 |
selectCoriMap = 0, |
| 64 |
f0 = 0.0, |
| 65 |
useSingleCpuIO=.TRUE., |
| 66 |
globalFiles=.TRUE., |
| 67 |
useCoriolis=.TRUE., |
| 68 |
& |
| 69 |
|
| 70 |
# Elliptic solver parameters |
| 71 |
&PARM02 |
| 72 |
cg2dMaxIters=300, |
| 73 |
cg2dTargetResidual=1.E-11, |
| 74 |
# cg2dMinColumnEps = 1. |
| 75 |
# pReleaseVisc = 1. |
| 76 |
#cg3dMaxIters=40, |
| 77 |
#cg3dTargetResidual=1.E-13, |
| 78 |
& |
| 79 |
|
| 80 |
#Time stepping parameters |
| 81 |
&PARM03 |
| 82 |
niter0=8640, |
| 83 |
# nTimeSteps=25920000., |
| 84 |
nTimeSteps=20, |
| 85 |
# startTime=0., |
| 86 |
# startTime = 0., |
| 87 |
#startTime=62208000., |
| 88 |
#startTime=15552000., |
| 89 |
#endTime=31104000., |
| 90 |
deltaT=300.0, |
| 91 |
### deltaT=1200.0, |
| 92 |
forcing_In_AB = .false., |
| 93 |
#rCD = 0.9896, |
| 94 |
abEps=0.1, |
| 95 |
cAdjFreq = -1., |
| 96 |
pChkptFreq=3000, |
| 97 |
chkptFreq=3000.0, |
| 98 |
# pChkptFreq=600.0, |
| 99 |
# chkptFreq=600.0, |
| 100 |
dumpFreq=3000, |
| 101 |
# monitorFreq=15552000., |
| 102 |
# monitorSelect=0, |
| 103 |
# monitorFreq=31104000., |
| 104 |
& |
| 105 |
|
| 106 |
# Gridding parameters |
| 107 |
&PARM04 |
| 108 |
usingCartesianGrid=.FALSE., |
| 109 |
usingSphericalPolarGrid=.TRUE., |
| 110 |
#delR = 10., 10., 10., 10., 10., 10., 10., 10.01, 10.03, 10.11, |
| 111 |
#10.32, 10.8, 11.76, 13.42, 16.04, 19.82, 24.85, 31.1, 38.42, 46.5, |
| 112 |
#55., 63.5, 71.58, 78.9, 85.15, 90.18, 93.96, 96.58, 98.25, 99.25, |
| 113 |
#100.01, 101.33, 104.56, 111.33, 122.83, 139.09, 158.94, 180.83, 203.55, 226.5, |
| 114 |
#249.5, 272.5, 295.5, 318.5, 341.5, 364.5, 387.5, 410.5, 433.5, 456.5, |
| 115 |
delR=90*10., |
| 116 |
delX=3*.125, |
| 117 |
delY=200*.0078125, |
| 118 |
# delYFile='dlat.bin', |
| 119 |
# delY=30*10e3, |
| 120 |
xgOrigin = -105.5, |
| 121 |
ygOrigin = -75.4457, |
| 122 |
& |
| 123 |
|
| 124 |
# Input datasets |
| 125 |
&PARM05 |
| 126 |
bathyFile='bathy.box', |
| 127 |
hydrogThetaFile='theta.init', |
| 128 |
hydrogSaltFile='salt.init', |
| 129 |
pSurfInitFile='etainit.round.bin', |
| 130 |
& |