/[MITgcm]/MITgcm_contrib/shelfice_remeshing/DIG/input/data
ViewVC logotype

Contents of /MITgcm_contrib/shelfice_remeshing/DIG/input/data

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (show annotations) (download)
Fri Apr 1 10:19:38 2016 UTC (9 years, 3 months ago) by dgoldberg
Branch: MAIN
CVS Tags: HEAD
Added rough code to dig ice shelf to make continuous ocean

1 #:::====================
2 # | Model parameters |
3 # ====================
4 #
5 # Continuous equation parameters
6 &PARM01
7 # momstepping=.FALSE.,
8 # tempstepping=.FALSE.,
9 # saltstepping=.FALSE.,
10 Tref = 50*-1.9,
11 Sref = 50*34.4,
12 viscAr=1.E-3,
13 #viscAh=600.0,
14 viscAhGrid=0.2,
15 viscA4Grid=0.02,
16 no_slip_sides=.FALSE.,
17 no_slip_bottom=.FALSE.,
18 diffKhT= 100.0,
19 #diffKhT= 600.0,
20 diffKrT=5.E-5,
21 diffKhS= 100.0,
22 #diffKhS= 600.0,
23 diffKrS=5.E-5,
24 bottomDragQuadratic=2.5E-3,
25 staggerTimestep = .true.,
26 tempAdvScheme=30,
27 saltAdvScheme=30,
28 vectorInvariantMomentum = .true.,
29 ### momImplVertAdv=.true.,
30 tempImplVertAdv=.true.,
31 saltImplVertAdv=.true.,
32 eosType='JMD95Z',
33 #tAlpha = 2.E-4,
34 #sBeta = 7.4E-4,
35 HeatCapacity_cp = 3974.0,
36 rhoConst=1000.,
37 rhoNil=1000.,
38 gravity=9.81,
39 gBaro=9.81,
40 implicitDiffusion = .true.,
41 implicitViscosity = .true.,
42 rigidLid=.FALSE.,
43 implicitFreeSurface=.TRUE.,
44 useRealFreshWaterFlux = .true.,
45 exactConserv=.TRUE.,
46 # start nlfs with rstar
47 #select_rStar=2,
48 nonlinFreeSurf=4,
49 hFacInf=0.2,
50 hFacSup=2.0,
51 # end nlfs with rstar
52 hFacMin=0.2,
53 #integr_GeoPot = 1,
54 #ivdc_kappa = 0.005.,
55 #implicitDiffusion=.true.,
56 #implicitViscosity=.true.,
57 #useCDscheme=.true.,
58 #nonHydrostatic=.true.,
59 #allowFreezing = .true.,
60 readBinaryPrec=64,
61 #writeBinaryPrec=64,
62 debuglevel = 1,
63 selectCoriMap = 0,
64 f0 = 0.0,
65 useSingleCpuIO=.TRUE.,
66 globalFiles=.TRUE.,
67 useCoriolis=.TRUE.,
68 &
69
70 # Elliptic solver parameters
71 &PARM02
72 cg2dMaxIters=300,
73 cg2dTargetResidual=1.E-11,
74 #cg3dMaxIters=40,
75 #cg3dTargetResidual=1.E-13,
76 &
77
78 #Time stepping parameters
79 &PARM03
80 niter0=1728,
81 nTimeSteps=1036800.,
82 # nTimeSteps=8640,
83 # startTime=51840.,
84 # startTime = 0.,
85 # startTime=62208000.,
86 # startTime=15552000.,
87 # endTime=31104000.,
88 deltaT=300.0,
89 # deltaT=1200.0,
90 forcing_In_AB = .false.,
91 #rCD = 0.9896,
92 abEps=0.1,
93 cAdjFreq = -1.,
94 pChkptFreq=15552000,
95 chkptFreq=15552000.0,
96 # pChkptFreq=600.0,
97 # chkptFreq=600.0,
98 monitorFreq=15552000.,
99 dumpFreq=15552000,
100 &
101
102 # Gridding parameters
103 &PARM04
104 usingCartesianGrid=.FALSE.,
105 usingSphericalPolarGrid=.TRUE.,
106 #delR = 10., 10., 10., 10., 10., 10., 10., 10.01, 10.03, 10.11,
107 #10.32, 10.8, 11.76, 13.42, 16.04, 19.82, 24.85, 31.1, 38.42, 46.5,
108 #55., 63.5, 71.58, 78.9, 85.15, 90.18, 93.96, 96.58, 98.25, 99.25,
109 #100.01, 101.33, 104.56, 111.33, 122.83, 139.09, 158.94, 180.83, 203.55, 226.5,
110 #249.5, 272.5, 295.5, 318.5, 341.5, 364.5, 387.5, 410.5, 433.5, 456.5,
111 delR=210*10.,
112 delX=3*.125,
113 delY=400*.0078125
114 # delYFile='dlat.bin',
115 # delY=30*10e3,
116 xgOrigin = -105.5,
117 ygOrigin = -75.4457,
118 &
119
120 # Input datasets
121 &PARM05
122 bathyFile='bathysteep.box',
123 hydrogThetaFile='theta.init',
124 hydrogSaltFile='salt.init',
125 pSurfInitFile='etainit.round.bin'
126
127 &

  ViewVC Help
Powered by ViewVC 1.1.22