global_ocean.cs32x15/input.thsice/data

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef=15*20.,
 sRef=15*35.,
 viscAh =3.E5,
 viscAr =1.E-3,
 diffKrT=3.E-5,
 diffKrS=3.E-5,
 ivdc_kappa=10.,
 implicitDiffusion=.TRUE.,
 implicitViscosity=.TRUE.,
 selectImplicitDrag=2,
 gravity=9.81,
 rhoConst=1035.,
 rhoConstFresh=1000.,
 eosType='JMD95Z',
 staggerTimeStep=.TRUE.,
 vectorInvariantMomentum=.TRUE.,
 implicitFreeSurface=.TRUE.,
 exactConserv=.TRUE.,
 select_rStar=2,
 nonlinFreeSurf=4,
 hFacInf=0.2,
 hFacSup=2.0,
 useRealFreshWaterFlux=.TRUE.,
 hFacMin=.1,
 hFacMinDr=20.,
 readBinaryPrec=64,
#writeBinaryPrec=64,
 plotLevel=0,
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=200,
#cg2dTargetResidual=1.E-9,
 cg2dTargetResWunit=1.E-14,
 &

# Time stepping parameters
 &PARM03
 nIter0=36000,
 nTimeSteps=20,
 deltaTMom   =1200.,
 deltaTtracer=86400.,
 deltaTFreeSurf=86400.,
 deltaTClock =86400.,
 abEps = 0.1,
#forcing_In_AB=.FALSE.,
 tracForcingOutAB=1,
 pChkptFreq  =31104000.,
 taveFreq    =31104000.,
 dumpFreq    =31104000.,
 monitorFreq =2592000.,
#tave_lastIter=0.,
 periodicExternalForcing=.TRUE.,
 externForcingPeriod=2592000.,
 externForcingCycle=31104000.,
# 2 months restoring timescale for temperature
 tauThetaClimRelax = 5184000.,
# 6 months restoring timescale for salinity
 tauSaltClimRelax = 15552000.,
 latBandClimRelax=50.,
 monitorFreq =1.,
 &

# Gridding parameters
 &PARM04
 usingCurvilinearGrid=.TRUE.,
 horizGridFile='grid_cs32',
 radius_fromHorizGrid=6370.E3,
 delR= 50., 70., 100., 140., 190.,
       240., 290., 340., 390., 440.,
       490., 540., 590., 640., 690.,
 &

# Input datasets
 &PARM05
 bathyFile      ='bathy_Hmin50.bin',
 hydrogThetaFile='lev_T_cs_15k.bin',
 hydrogSaltFile ='lev_S_cs_15k.bin',
 zonalWindFile  ='trenberth_taux.bin',
 meridWindFile  ='trenberth_tauy.bin',
 thetaClimFile  ='lev_surfT_cs_12m.bin',
 saltClimFile   ='lev_surfS_cs_12m.bin',
 &
1	jmc	1.1	# ====================
2			# \| Model parameters \|
3			# ====================
4			#
5			# Continuous equation parameters
6			&PARM01
7			tRef=15*20.,
8			sRef=15*35.,
9			viscAh =3.E5,
10			viscAr =1.E-3,
11			diffKrT=3.E-5,
12			diffKrS=3.E-5,
13			ivdc_kappa=10.,
14			implicitDiffusion=.TRUE.,
15	jmc	1.11	implicitViscosity=.TRUE.,
16			selectImplicitDrag=2,
17	jmc	1.1	gravity=9.81,
18	jmc	1.11	rhoConst=1035.,
19	jmc	1.1	rhoConstFresh=1000.,
20			eosType='JMD95Z',
21	jmc	1.4	staggerTimeStep=.TRUE.,
22	jmc	1.1	vectorInvariantMomentum=.TRUE.,
23			implicitFreeSurface=.TRUE.,
24			exactConserv=.TRUE.,
25	jmc	1.4	select_rStar=2,
26			nonlinFreeSurf=4,
27	jmc	1.1	hFacInf=0.2,
28			hFacSup=2.0,
29	jmc	1.4	useRealFreshWaterFlux=.TRUE.,
30	jmc	1.1	hFacMin=.1,
31			hFacMinDr=20.,
32			readBinaryPrec=64,
33	jmc	1.4	#writeBinaryPrec=64,
34	jmc	1.12	plotLevel=0,
35	jmc	1.1	&
36
37			# Elliptic solver parameters
38			&PARM02
39			cg2dMaxIters=200,
40			#cg2dTargetResidual=1.E-9,
41			cg2dTargetResWunit=1.E-14,
42			&
43
44			# Time stepping parameters
45			&PARM03
46	jmc	1.4	nIter0=36000,
47	jmc	1.1	nTimeSteps=20,
48	jmc	1.11	deltaTMom =1200.,
49	jmc	1.1	deltaTtracer=86400.,
50	jmc	1.11	deltaTFreeSurf=86400.,
51	jmc	1.1	deltaTClock =86400.,
52			abEps = 0.1,
53	jmc	1.8	#forcing_In_AB=.FALSE.,
54			tracForcingOutAB=1,
55	jmc	1.4	pChkptFreq =31104000.,
56			taveFreq =31104000.,
57			dumpFreq =31104000.,
58	jmc	1.1	monitorFreq =2592000.,
59			#tave_lastIter=0.,
60			periodicExternalForcing=.TRUE.,
61			externForcingPeriod=2592000.,
62			externForcingCycle=31104000.,
63			# 2 months restoring timescale for temperature
64			tauThetaClimRelax = 5184000.,
65			# 6 months restoring timescale for salinity
66			tauSaltClimRelax = 15552000.,
67			latBandClimRelax=50.,
68			monitorFreq =1.,
69			&
70
71			# Gridding parameters
72			&PARM04
73			usingCurvilinearGrid=.TRUE.,
74	jmc	1.9	horizGridFile='grid_cs32',
75	jmc	1.11	radius_fromHorizGrid=6370.E3,
76	jmc	1.1	delR= 50., 70., 100., 140., 190.,
77			240., 290., 340., 390., 440.,
78			490., 540., 590., 640., 690.,
79			&
80
81			# Input datasets
82			&PARM05
83			bathyFile ='bathy_Hmin50.bin',
84			hydrogThetaFile='lev_T_cs_15k.bin',
85			hydrogSaltFile ='lev_S_cs_15k.bin',
86			zonalWindFile ='trenberth_taux.bin',
87			meridWindFile ='trenberth_tauy.bin',
88			thetaClimFile ='lev_surfT_cs_12m.bin',
89			saltClimFile ='lev_surfS_cs_12m.bin',
90			&