/[MITgcm]/MITgcm/verification/MLAdjust/input.1.leith/data
ViewVC logotype

Annotation of /MITgcm/verification/MLAdjust/input.1.leith/data

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

Revision 1.6 - (hide annotations) (download)
Wed Oct 12 20:58:27 2005 UTC (18 years, 8 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint58l_post, checkpoint58e_post, checkpoint57v_post, checkpoint58u_post, checkpoint60, checkpoint61, checkpoint58r_post, checkpoint57y_post, checkpoint58g_post, checkpoint57x_post, checkpoint58n_post, checkpoint58x_post, checkpoint58t_post, checkpoint58h_post, checkpoint58w_post, checkpoint58j_post, checkpoint57y_pre, checkpoint58q_post, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint59j, checkpoint59, checkpoint58, checkpoint58f_post, checkpoint58d_post, checkpoint58c_post, checkpoint57w_post, checkpoint58a_post, checkpoint58i_post, checkpoint58o_post, checkpoint57z_post, checkpoint58y_post, checkpoint58k_post, checkpoint58v_post, checkpoint58s_post, checkpoint61f, checkpoint61g, checkpoint61d, checkpoint61e, checkpoint61b, checkpoint61c, checkpoint58p_post, checkpoint61a, checkpoint61n, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint58b_post, checkpoint58m_post
Changes since 1.5: +0 -1 lines 
return to default tRef profile (now uniform, =20 deg.C)
1 baylor 1.1 # ====================
2     # | Model parameters |
3     # ====================
4     #
5     # Continuous equation parameters
6     &PARM01
7 baylor 1.3 viscAr=5.E-5,
8     # This is the background viscosity used total=viscah+viscleith
9     # viscAh=1.E2,
10     # viscA4=1.E2,
11     # This is the background viscosity used on vorticity points
12     # viscAhZ=1.E2,
13     # viscA4Z=1.E2,
14     # This is the background viscosity used on divergence points
15     # viscAhD=1.E2,
16     # viscA4D=1.E2,
17     # viscAhd and viscahz are set to viscah if they are unset initially
18     #
19     # This is a maximum gridscale Reynolds number limiter
20     # ViscAhRemax=10,
21     # ViscA4Remax=10,
22     # This is a minimum autoscaling viscosity
23     # ViscAhGridMin=.01,
24     # ViscA4GridMin=.01,
25 baylor 1.1 # This is the standard Leith Viscosity, prop to grad(curl v_h)
26 baylor 1.4 # viscC2leith=1.85,
27     viscC4leith=1.85,
28 baylor 1.1 # This is an additional Leith Viscosity, prop to grad(div.v_h)
29     # It traps instabilities that have no signal in grad(curl v_h)
30 baylor 1.4 # viscC2leithD=1.85,
31     viscC4leithD=1.85,
32 baylor 1.1 # This is harmonic Smagorinsky Coefficient
33 baylor 1.3 # 0.2-0.9=value from 3-d turbulence theory; Smagorinsky (1993)
34     # 2.2-4=stable oceanic values (Griffies and Hallberg, 2000)
35 baylor 1.4 # 1.85 smag is approximately the same as 1.85 leith+leithd
36 baylor 1.3 # viscC2Smag=1.0,
37     viscC4Smag=1.0,
38 baylor 1.1 # This is the max Viscosity used, be it Smag or Leith, and it scales
39     # with grid size and timestep
40 baylor 1.3 # viscAhGridMax=1.0,
41     viscA4GridMax=1.0,
42     # This uses the full grad(vort) and grad(div)
43 baylor 1.1 useFullLeith=.TRUE.,
44 baylor 1.3 # This uses the Strain-Tension form for calculating viscous terms
45     # useStrainTensionVisc=.TRUE.,
46 baylor 1.1 no_slip_sides=.FALSE.,
47     no_slip_bottom=.FALSE.,
48     ivdc_kappa=10.,
49     implicitDiffusion=.TRUE.,
50     implicitViscosity=.TRUE.,
51     f0=7.29e-5,
52     beta=0.E-11,
53     eosType='LINEAR',
54     # Using Temperature as only Active Tracer
55     tAlpha=2.E-4,
56     # Using Salinity as Spice
57     sBeta =0.E-4,
58     saltStepping=.TRUE.,
59     gravity=9.81,
60     rhonil=1035.,
61     implicitFreeSurface=.TRUE.,
62     usejamartwetpoints=.TRUE.,
63     rigidLid=.FALSE.,
64     staggerTimeStep=.TRUE.,
65     readBinaryPrec=64,
66     writeBinaryPrec=64,
67     tempAdvScheme=33,
68     saltAdvScheme=33,
69     nonHydrostatic=.FALSE.,
70     vectorInvariantMomentum=.TRUE.,
71     /
72    
73     # Elliptic solver parameters
74     &PARM02
75     cg2dMaxIters=1000,
76     cg2dTargetResidual=1.E-11,
77     cg3dMaxIters=40,
78     cg3dTargetResidual=1.E-9,
79     /
80    
81     # Time stepping parameters
82     &PARM03
83     niter0=36,
84 baylor 1.3 #endTime=86400.,
85     #- for testing, only run for 12 iterations:
86     nTimeSteps=12,
87 baylor 1.1 deltaT=1200.,
88     abEps=0.01,
89 baylor 1.3 #- for testing, reduce amount of output by commenting out output frequencies:
90     # pchkptFreq=43200.,
91     #taveFreq=7200.,
92     #dumpFreq=7200.,
93     #diagFreq=7200.,
94 baylor 1.1 monitorFreq=1200.,
95     cadjfreq=0.,
96     /
97    
98     # Gridding parameters
99     &PARM04
100     usingCartesianGrid=.TRUE.,
101     usingSphericalPolarGrid=.FALSE.,
102     delR= 5d0, 5d0, 5d0, 5d0, 5d0,
103     5d0, 5d0, 5d0, 5d0, 5d0,
104     5d0, 5d0, 5d0, 5d0, 5d0,
105     5d0, 5d0, 5d0, 5d0, 5d0,
106     5d0, 5d0, 5d0, 5d0, 5d0,
107     5d0, 5d0, 5d0, 5d0, 5d0,
108     5d0, 5d0, 5d0, 5d0, 5d0,
109     5d0, 5d0, 5d0, 5d0, 5d0,
110     dxspacing=1d3,
111     dyspacing=1d3,
112     /
113    
114     # Input datasets
115     &PARM05
116     bathyFile='topo_sl.bin',
117     hydrogThetaFile='thetaInitial.bin',
118     hydrogSaltFile='spiceInitial.bin',
119     /
  ViewVC Help
Powered by ViewVC 1.1.22