MITgcm_contrib/edge_problem/data.visca4grid

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef=15*20.,
 sRef=15*35.,
 viscAh =0,
 viscA4 =0,
 viscA4grid = 4e-3
 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,
 deltaT =1200.,
 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	# ====================
2	# \| Model parameters \|
3	# ====================
4	#
5	# Continuous equation parameters
6	&PARM01
7	tRef=15*20.,
8	sRef=15*35.,
9	viscAh =0,
10	viscA4 =0,
11	viscA4grid = 4e-3
12	viscAr =1.E-3,
13	diffKhT=0.,
14	diffK4T=0.,
15	diffKrT=3.E-5,
16	diffKhS=0.,
17	diffK4S=0.,
18	diffKrS=3.E-5,
19	ivdc_kappa=10.,
20	implicitDiffusion=.TRUE.,
21	gravity=9.81,
22	rhonil=1035.,
23	eosType='JMD95Z',
24	#staggerTimeStep=.TRUE.,
25	vectorInvariantMomentum=.TRUE.,
26	implicitFreeSurface=.TRUE.,
27	exactConserv=.TRUE.,
28	select_rStar=2,
29	nonlinFreeSurf=4,
30	hFacInf=0.2,
31	hFacSup=2.0,
32	useRealFreshWaterFlux=.TRUE.,
33	allowFreezing=.TRUE.,
34	hFacMin=.1,
35	hFacMinDr=20.,
36	readBinaryPrec=64,
37	globalfiles=.true.
38	&
39
40	# Elliptic solver parameters
41	&PARM02
42	cg2dMaxIters=200,
43	#cg2dTargetResidual=1.E-9,
44	cg2dTargetResWunit=1.E-14,
45	&
46
47	# Time stepping parameters
48	&PARM03
49	nIter0=72000,
50	nTimeSteps=3,
51	deltaT =1200.,
52	abEps = 0.1,
53	forcing_In_AB=.FALSE.,
54	pChkptFreq =1728000.,
55	#chkptFreq =311040000.,
56	taveFreq =1244160000.,
57	dumpFreq =622080000.,
58	monitorFreq =31104000.,
59	periodicExternalForcing=.TRUE.,
60	externForcingPeriod=2592000.,
61	externForcingCycle=31104000.,
62	# 2 months restoring timescale for temperature
63	tauThetaClimRelax = 5184000.,
64	# 2yrs restoring timescale for salinity
65	tauSaltClimRelax = 62208000.,
66	monitorFreq =1.,
67	&
68
69	# Gridding parameters
70	&PARM04
71	usingCartesianGrid=.FALSE.,
72	usingSphericalPolarGrid=.FALSE.,
73	usingCurvilinearGrid=.TRUE.,
74	delR= 50., 70., 100., 140., 190.,
75	240., 290., 340., 390., 440.,
76	490., 540., 590., 640., 690.,
77	&
78
79	# Input datasets
80	&PARM05
81	bathyFile ='bathy_Hmin50.bin',
82	hydrogThetaFile='lev_T_cs_15k.bin',
83	hydrogSaltFile ='lev_S_cs_15k.bin',
84	zonalWindFile ='trenberth_taux.bin',
85	meridWindFile ='trenberth_tauy.bin',
86	thetaClimFile ='lev_surfT_cs_12m.bin',
87	saltClimFile ='lev_surfS_cs_12m.bin',
88	surfQFile ='shiQnet_cs32.bin',
89	EmPmRFile ='shiEmPR_cs32.bin',
90	&