MITgcm_contrib/edge_problem/data.base

# ====================
# | 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.,
 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,
 globalfiles=.true.
 &

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

# Time stepping parameters
 &PARM03
 nIter0=72000,
 nTimeSteps=3,
 deltaTmom   =1200.,
 deltaTtracer=86400.,
 deltaTfreesurf=86400.,
 deltaTClock =86400.,
 abEps = 0.1,
 forcing_In_AB=.FALSE.,
 pChkptFreq  =1728000., 
#chkptFreq   =311040000., 
 taveFreq    =1244160000., 
 dumpFreq    =622080000.,
 monitorFreq =31104000.,
 periodicExternalForcing=.TRUE.,
 externForcingPeriod=2592000.,
 externForcingCycle=31104000.,
# 2 months restoring timescale for temperature
 tauThetaClimRelax = 5184000.,
# 2yrs restoring timescale for salinity
 tauSaltClimRelax = 62208000.,
 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	dimitri	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			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			eosType='JMD95Z',
22			#staggerTimeStep=.TRUE.,
23			vectorInvariantMomentum=.TRUE.,
24			implicitFreeSurface=.TRUE.,
25			exactConserv=.TRUE.,
26			select_rStar=2,
27			nonlinFreeSurf=4,
28			hFacInf=0.2,
29			hFacSup=2.0,
30			useRealFreshWaterFlux=.TRUE.,
31			allowFreezing=.TRUE.,
32			hFacMin=.1,
33			hFacMinDr=20.,
34			readBinaryPrec=64,
35			globalfiles=.true.
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=3,
49			deltaTmom =1200.,
50			deltaTtracer=86400.,
51			deltaTfreesurf=86400.,
52			deltaTClock =86400.,
53			abEps = 0.1,
54			forcing_In_AB=.FALSE.,
55			pChkptFreq =1728000.,
56			#chkptFreq =311040000.,
57			taveFreq =1244160000.,
58			dumpFreq =622080000.,
59			monitorFreq =31104000.,
60			periodicExternalForcing=.TRUE.,
61			externForcingPeriod=2592000.,
62			externForcingCycle=31104000.,
63			# 2 months restoring timescale for temperature
64			tauThetaClimRelax = 5184000.,
65			# 2yrs restoring timescale for salinity
66			tauSaltClimRelax = 62208000.,
67			monitorFreq =1.,
68			&
69
70			# Gridding parameters
71			&PARM04
72			usingCartesianGrid=.FALSE.,
73			usingSphericalPolarGrid=.FALSE.,
74			usingCurvilinearGrid=.TRUE.,
75			delR= 50., 70., 100., 140., 190.,
76			240., 290., 340., 390., 440.,
77			490., 540., 590., 640., 690.,
78			&
79
80			# Input datasets
81			&PARM05
82			bathyFile ='bathy_Hmin50.bin',
83			hydrogThetaFile='lev_T_cs_15k.bin',
84			hydrogSaltFile ='lev_S_cs_15k.bin',
85			zonalWindFile ='trenberth_taux.bin',
86			meridWindFile ='trenberth_tauy.bin',
87			thetaClimFile ='lev_surfT_cs_12m.bin',
88			saltClimFile ='lev_surfS_cs_12m.bin',
89			surfQFile ='shiQnet_cs32.bin',
90			EmPmRFile ='shiEmPR_cs32.bin',
91			&