global_ocean.cs32x15/input_ad.thsice/data

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef=15*20.,
 sRef=15*35.,
 viscAh =3.E5,
#- biharmonic Viscosity: 3.e15 is close to the stability limit with deltaTmom=20mn
#viscA4 =3.E15,
 viscAr =1.E-3,
 diffKhT=0.,
 diffK4T=0.,
#- diffKrT unused when compiled with ALLOW_3D_DIFFKR
#diffKrT=3.E-5,
 diffKhS=0.,
 diffK4S=0.,
 diffKrS=3.E-5,
 ivdc_kappa=10.,
 implicitDiffusion=.TRUE.,
 gravity=9.81,
 rhoConst=1035.,
 rhoConstFresh=1000.,
 eosType='JMD95Z',
 staggerTimeStep=.TRUE.,
 vectorInvariantMomentum=.TRUE.,
 implicitFreeSurface=.TRUE.,
 tempAdvScheme=30,
 saltAdvScheme=30,
 exactConserv=.TRUE.,
 select_rStar=1,
 nonlinFreeSurf=2,
 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=5,
 deltaTMom   =1200.,
 deltaTtracer=86400.,
 deltaTFreeSurf=86400.,
 deltaTClock =86400.,
 abEps = 0.1,
 alph_AB=0.6,
 beta_AB=0.0,
 forcing_In_AB=.FALSE.,
#momDissip_In_AB=.FALSE.,
 pChkptFreq  =311040000.,
 chkptFreq   = 31104000.,
#taveFreq    =311040000.,
#dumpFreq    = 31104000.,
#adjDumpFreq = 31104000.,
#monitorFreq = 31104000.,
#- forcing is set by EXF
# periodicExternalForcing=.TRUE.,
# externForcingPeriod=2592000.,
# externForcingCycle=31104000.,
# 2 months restoring timescale for temperature
# tauThetaClimRelax = 5184000.,
# 2yrs restoring timescale for salinity
# tauSaltClimRelax = 62208000.,
 monitorFreq   =1.,
 adjMonitorFreq=1.,
 dumpFreq    = 432000.,
 adjDumpFreq = 432000.,
 pickupStrictlyMatch=.FALSE.,
 &

# Gridding parameters
 &PARM04
 usingCurvilinearGrid=.TRUE.,
 horizGridFile='grid_cs32',
 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',
#- forcing is set by EXF
# 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	#- biharmonic Viscosity: 3.e15 is close to the stability limit with deltaTmom=20mn
11	#viscA4 =3.E15,
12	viscAr =1.E-3,
13	diffKhT=0.,
14	diffK4T=0.,
15	#- diffKrT unused when compiled with ALLOW_3D_DIFFKR
16	#diffKrT=3.E-5,
17	diffKhS=0.,
18	diffK4S=0.,
19	diffKrS=3.E-5,
20	ivdc_kappa=10.,
21	implicitDiffusion=.TRUE.,
22	gravity=9.81,
23	rhoConst=1035.,
24	rhoConstFresh=1000.,
25	eosType='JMD95Z',
26	staggerTimeStep=.TRUE.,
27	vectorInvariantMomentum=.TRUE.,
28	implicitFreeSurface=.TRUE.,
29	tempAdvScheme=30,
30	saltAdvScheme=30,
31	exactConserv=.TRUE.,
32	select_rStar=1,
33	nonlinFreeSurf=2,
34	hFacInf=0.2,
35	hFacSup=2.0,
36	useRealFreshWaterFlux=.TRUE.,
37	### allowFreezing=.TRUE.,
38	hFacMin=.1,
39	hFacMinDr=20.,
40	readBinaryPrec=64,
41	&
42
43	# Elliptic solver parameters
44	&PARM02
45	cg2dMaxIters=200,
46	#cg2dTargetResidual=1.E-9,
47	cg2dTargetResWunit=1.E-14,
48	&
49
50	# Time stepping parameters
51	&PARM03
52	nIter0=72000,
53	nTimeSteps=5,
54	deltaTMom =1200.,
55	deltaTtracer=86400.,
56	deltaTFreeSurf=86400.,
57	deltaTClock =86400.,
58	abEps = 0.1,
59	alph_AB=0.6,
60	beta_AB=0.0,
61	forcing_In_AB=.FALSE.,
62	#momDissip_In_AB=.FALSE.,
63	pChkptFreq =311040000.,
64	chkptFreq = 31104000.,
65	#taveFreq =311040000.,
66	#dumpFreq = 31104000.,
67	#adjDumpFreq = 31104000.,
68	#monitorFreq = 31104000.,
69	#- forcing is set by EXF
70	# periodicExternalForcing=.TRUE.,
71	# externForcingPeriod=2592000.,
72	# externForcingCycle=31104000.,
73	# 2 months restoring timescale for temperature
74	# tauThetaClimRelax = 5184000.,
75	# 2yrs restoring timescale for salinity
76	# tauSaltClimRelax = 62208000.,
77	monitorFreq =1.,
78	adjMonitorFreq=1.,
79	dumpFreq = 432000.,
80	adjDumpFreq = 432000.,
81	pickupStrictlyMatch=.FALSE.,
82	&
83
84	# Gridding parameters
85	&PARM04
86	usingCurvilinearGrid=.TRUE.,
87	horizGridFile='grid_cs32',
88	delR= 50., 70., 100., 140., 190.,
89	240., 290., 340., 390., 440.,
90	490., 540., 590., 640., 690.,
91	&
92
93	# Input datasets
94	&PARM05
95	bathyFile ='bathy_Hmin50.bin',
96	hydrogThetaFile='lev_T_cs_15k.bin',
97	hydrogSaltFile ='lev_S_cs_15k.bin',
98	#- forcing is set by EXF
99	# zonalWindFile ='trenberth_taux.bin',
100	# meridWindFile ='trenberth_tauy.bin',
101	# thetaClimFile ='lev_surfT_cs_12m.bin',
102	# saltClimFile ='lev_surfS_cs_12m.bin',
103	# surfQFile ='shiQnet_cs32.bin',
104	# EmPmRFile ='shiEmPR_cs32.bin',
105	&