/[MITgcm]/MITgcm_contrib/high_res_cube/input/data_hr
ViewVC logotype

Contents of /MITgcm_contrib/high_res_cube/input/data_hr

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


Revision 1.23 - (show annotations) (download)
Sat Feb 25 16:59:00 2006 UTC (19 years, 5 months ago) by dimitri
Branch: MAIN
Changes since 1.22: +7 -5 lines
updated instructions for cube29

1 # ====================
2 # | Model parameters |
3 # ====================
4 #
5 # Continuous equation parameters
6 &PARM01
7 tRef = 3*23.,3*22.,21.,2*20.,19.,2*18.,17.,2*16.,15.,14.,13.,
8 12.,11.,2*9.,8.,7.,2*6.,2*5.,3*4.,3*3.,4*2.,12*1.,
9 sRef = 50*34.5,
10 viscAr= 1e-3,
11 no_slip_sides=.FALSE.,
12 no_slip_bottom=.TRUE.,
13 diffKhT=0.,
14 diffK4T=0.,
15 diffKrT=1.51e-5,
16 diffKhS=0.,
17 diffK4S=0.,
18 diffKrS=1.51e-5,
19 tAlpha=2.E-4,
20 sBeta =7.4E-4,
21 gravity=9.81,
22 eosType='JMD95Z',
23 readBinaryPrec=32,
24 writeBinaryPrec=32,
25 useSingleCPUio=.TRUE.
26 hFacMinDr=50.,
27 hFacMin=0.3,
28 hFacInf=0.1,
29 hFacSup=5.,
30 # z* - begin
31 select_rStar=2,
32 nonlinFreeSurf=4,
33 # * - end
34 implicitDiffusion=.true.,
35 implicitViscosity=.true.,
36 allowFreezing=.false.,
37 viscAh=0.,
38 viscAhGrid=0.,
39 # Following change took place starting in checkpoint57v_post
40 # Redefine viscC?Leith*. Change the coefficients of Leith viscosity
41 # to be proportional to enstrophy dissipation length. Thus,
42 # viscC?leith?_new=pi*(viscC?leith?_old)^(1/3)
43 # This is now consistent with Smagorinsky and with the manual.
44 # viscC4Leith=1.0,
45 # viscC4Leithd=1.0,
46 # viscC4Leith=3.1416,
47 # viscC4Leithd=3.1416,
48 viscC4Leith=1.5,
49 viscC4Leithd=1.5,
50 viscA4GridMax=0.5,
51 ##########################################################
52 viscC2Leith=0.,
53 viscA4=0.,
54 useAreaViscLength=.TRUE.,
55 sideDragFactor=0.,
56 SadournyCoriolis=.FALSE.,
57 useAbsVorticity=.FALSE.,
58 highOrderVorticity = .TRUE.
59 bottomDragQuadratic = 0.002,
60 tempAdvScheme=33,
61 saltAdvScheme=33,
62 StaggerTimeStep=.true.,
63 multiDimAdvection=.true.,
64 vectorInvariantMomentum=.TRUE.,
65 rigidLid=.FALSE.,
66 implicitFreeSurface=.TRUE.,
67 exactConserv=.TRUE.,
68 debuglevel=-1
69 &
70
71 # Elliptic solver parameters
72 &PARM02
73 cg2dMaxIters=100,
74 cg2dTargetResidual=1.E-5,
75 &
76
77 # Time stepping parameters
78 &PARM03
79 niter0=216,
80 endTime=421459200,
81 abEps=0.1,
82 deltaT=1200,
83 cAdjFreq=0.,
84 pChkptFreq=31536000,
85 chkptFreq=2592000
86 tauSaltClimRelax = 3842380.,
87 monitorFreq=86400
88 forcing_In_AB = .FALSE.
89 calendarDumps = .TRUE.
90 dumpInitAndLast=.FALSE.
91 &
92
93 # Gridding parameters
94 &PARM04
95 usingCartesianGrid=.FALSE.,
96 usingSphericalPolarGrid=.FALSE.,
97 usingCurvilinearGrid=.TRUE.,
98 delR = 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.01,
99 10.03, 10.11, 10.32, 10.80, 11.76, 13.42, 16.04 , 19.82, 24.85,
100 31.10, 38.42, 46.50, 55.00, 63.50, 71.58, 78.90, 85.15, 90.18,
101 93.96, 96.58, 98.25, 99.25,100.01,101.33,104.56,111.33,122.83,
102 139.09,158.94,180.83,203.55,226.50,249.50,272.50,295.50,318.50,
103 341.50,364.50,387.50,410.50,433.50,456.50,
104 &
105
106 # Input datasets
107 &PARM05
108 bathyFile ='ETOPO2_510x6x510_ver5e_filled.bin',
109 hydrogThetaFile='LEVT01_JAN_510x6x510x50.bin',
110 hydrogSaltFile ='LEVS01_JAN_510x6x510x50.bin',
111 &

  ViewVC Help
Powered by ViewVC 1.1.22