global_ocean.cs32x15/inp_thsice/data

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef=15*20.,
 sRef=15*35.,
 viscAh =3.E5,
 viscAr =1.E-3,
 diffKhT=0., 
 diffK4T=0., 
 diffKrT=3.E-5,
 diffKhS=0.,  
 diffK4S=0., 
 diffKrS=3.E-5,
 ivdc_kappa=10.,
 implicitDiffusion=.TRUE.,
 gravity=9.81,
 rhonil=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.,
#allowFreezing=.TRUE.,
 hFacMin=.1,
 hFacMinDr=20.,
 readBinaryPrec=64,
 &

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

# Time stepping parameters
 &PARM03
 nIter0=72000,
 nTimeSteps=20,
 deltaTmom   =1200.,
 deltaTtracer=86400.,
 deltaTfreesurf=86400.,
 deltaTClock =86400.,
 abEps = 0.1,
 forcing_In_AB=.FALSE.,
 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
 usingCartesianGrid=.FALSE.,
 usingSphericalPolarGrid=.FALSE.,
 usingCurvilinearGrid=.TRUE.,
 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',
 surfQFile      ='shiQnet_cs32.bin',
 EmPmRFile      ='shiEmPR_cs32.bin',
 &
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	diffKhT=0.,
12	diffK4T=0.,
13	diffKrT=3.E-5,
14	diffKhS=0.,
15	diffK4S=0.,
16	diffKrS=3.E-5,
17	ivdc_kappa=10.,
18	implicitDiffusion=.TRUE.,
19	gravity=9.81,
20	rhonil=1035.,
21	rhoConstFresh=1000.,
22	eosType='JMD95Z',
23	#staggerTimeStep=.TRUE.,
24	vectorInvariantMomentum=.TRUE.,
25	implicitFreeSurface=.TRUE.,
26	exactConserv=.TRUE.,
27	#select_rStar=2,
28	#nonlinFreeSurf=4,
29	hFacInf=0.2,
30	hFacSup=2.0,
31	#useRealFreshWaterFlux=.TRUE.,
32	#allowFreezing=.TRUE.,
33	hFacMin=.1,
34	hFacMinDr=20.,
35	readBinaryPrec=64,
36	&
37
38	# Elliptic solver parameters
39	&PARM02
40	cg2dMaxIters=200,
41	#cg2dTargetResidual=1.E-9,
42	cg2dTargetResWunit=1.E-14,
43	&
44
45	# Time stepping parameters
46	&PARM03
47	nIter0=72000,
48	nTimeSteps=20,
49	deltaTmom =1200.,
50	deltaTtracer=86400.,
51	deltaTfreesurf=86400.,
52	deltaTClock =86400.,
53	abEps = 0.1,
54	forcing_In_AB=.FALSE.,
55	pChkptFreq =31104000.,
56	taveFreq =31104000.,
57	dumpFreq =31104000.,
58	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	usingCartesianGrid=.FALSE.,
74	usingSphericalPolarGrid=.FALSE.,
75	usingCurvilinearGrid=.TRUE.,
76	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	surfQFile ='shiQnet_cs32.bin',
91	EmPmRFile ='shiEmPR_cs32.bin',
92	&