/[MITgcm]/MITgcm/model/inc/PARAMS.h
ViewVC logotype

Annotation of /MITgcm/model/inc/PARAMS.h

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


Revision 1.25 - (hide annotations) (download)
Sat Aug 22 17:51:06 1998 UTC (25 years, 8 months ago) by cnh
Branch: MAIN
Changes since 1.24: +17 -14 lines
File MIME type: text/plain
Isomorphism consistency changes

1 cnh 1.25 C $Header: /u/gcmpack/models/MITgcmUV/model/inc/PARAMS.h,v 1.24 1998/08/21 12:29:55 cnh Exp $
2 cnh 1.1 C
3     C /==========================================================\
4     C | PARAMS.h |
5     C | o Header file defining model "parameters". |
6     C |==========================================================|
7     C | The values from the model's standard input file are |
8     C | stored into the variables held here. Notes describing |
9     C | the parameters can also be found here. |
10     C \==========================================================/
11 cnh 1.13
12     C Macros for special grid options
13     #include "PARAMS_MACROS.h"
14 cnh 1.1
15     C-- Contants
16     C Useful physical values
17     Real*8 PI
18     PARAMETER ( PI = 3.14159265358979323844D0 )
19     Real*8 deg2rad
20     PARAMETER ( deg2rad = 2.D0*PI/360.D0 )
21    
22 cnh 1.7 C Symbolic values
23     C precXXXX - Used to indicate what precision to use for
24     C dumping model state.
25     INTEGER precFloat32
26     PARAMETER ( precFloat32 = 0 )
27     INTEGER precFloat64
28     PARAMETER ( precFloat64 = 1 )
29    
30     C Checkpoint data
31     INTEGER maxNoChkptLev
32     PARAMETER ( maxNoChkptLev = 2 )
33    
34 cnh 1.1 C-- COMMON /PARM_C/ Character valued parameters used by the model.
35 cnh 1.7 C checkPtSuff - List of checkpoint file suffices
36 cnh 1.16 C bathyFile - File containing bathymetry. If not defined bathymetry
37     C is taken from inline function.
38     C hydrogThetaFile - File containing initial hydrographic data for potential
39     C temperature.
40     C hydrogSaltFile - File containing initial hydrographic data for salinity.
41     C zonalWindFile - File containing zonal wind data
42     C meridWindFile - File containing meridional wind data
43 cnh 1.18 C thetaClimFile - File containing theta climataology used
44     C in relaxation term -lambda(theta-theta*)
45     C saltClimFile - File containing salt climataology used
46     C in relaxation term -lambda(salt-salt*)
47 cnh 1.16 COMMON /PARM_C/ checkPtSuff,
48     & bathyFile, hydrogThetaFile, hydrogSaltFile,
49 cnh 1.18 & zonalWindFile, meridWindFile, thetaClimFile,
50     & saltClimFile
51 cnh 1.7 CHARACTER*(5) checkPtSuff(maxNoChkptLev)
52 cnh 1.16 CHARACTER*(MAX_LEN_FNAM) bathyFile
53     CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile
54     CHARACTER*(MAX_LEN_FNAM) hydrogSaltFile
55     CHARACTER*(MAX_LEN_FNAM) zonalWindFile
56     CHARACTER*(MAX_LEN_FNAM) meridWindFile
57 cnh 1.18 CHARACTER*(MAX_LEN_FNAM) thetaClimFile
58     CHARACTER*(MAX_LEN_FNAM) saltClimFile
59 cnh 1.7
60 cnh 1.1 C-- COMMON /PARM_I/ Integer valued parameters used by the model.
61     C cg2dMaxIters - Maximum number of iterations in the
62     C two-dimensional con. grad solver.
63     C cg2dChkResFreq - Frequency with which to check residual
64     C in con. grad solver.
65     C nIter0 - Start time-step number of for this run
66     C nTimeSteps - Number of timesteps to execute
67     C numStepsPerPickup - For offline setup. Frequency of pickup
68     C of flow fields.
69 cnh 1.7 C writeStatePrec - Precision used for writing model state.
70     C writeBinaryPrec - Precision used for writing binary files
71     C readBinaryPrec - Precision used for reading binary files
72     C nCheckLev - Holds current checkpoint level
73 cnh 1.1 COMMON /PARM_I/
74     & cg2dMaxIters,
75     & cg2dChkResFreq,
76     & nIter0, nTimeSteps,
77 cnh 1.7 & numStepsPerPickup,
78     & writeStatePrec, nCheckLev,
79     & writeBinaryPrec, readBinaryPrec
80 cnh 1.1 INTEGER cg2dMaxIters
81     INTEGER cg2dChkResFreq
82     INTEGER nIter0
83     INTEGER nTimeSteps
84     INTEGER numStepsPerPickup
85 cnh 1.7 INTEGER writeStatePrec
86     INTEGER writeBinaryPrec
87     INTEGER readBinaryPrec
88     INTEGER nCheckLev
89 cnh 1.1
90     C-- COMMON /PARM_L/ Logical valued parameters used by the model.
91     C usingCartesianGrid - If TRUE grid generation will be in a cartesian
92     C coordinate frame.
93     C usingSphericalPolarGrid - If TRUE grid generation will be in a
94     C spherical polar frame.
95     C momViscosity - Flag which turns momentum friction terms on and off.
96     C momAdvection - Flag which turns advection of momentum on and off.
97     C momForcing - Flag which turns external forcing of momentum on
98     C and off.
99 cnh 1.9 C momPressureForcing - Flag which turns pressure term in momentum equation
100     C on and off.
101 cnh 1.15 C metricTerms - Flag which turns metric terms on or off.
102     C usingSphericalPolarMTerms - If TRUE use spherical polar metric terms.
103 cnh 1.1 C useCoriolis - Flag which turns the coriolis terms on and off.
104     C tempDiffusion - Flag which turns diffusion of temperature on
105     C and off.
106     C tempAdvection - Flag which turns advection of temperature on
107     C and off.
108     C tempForcing - Flag which turns external forcing of temperature on
109     C and off.
110     C saltDiffusion - Flag which turns diffusion of salinit on
111     C and off.
112     C saltAdvection - Flag which turns advection of salinit on
113     C and off.
114     C saltForcing - Flag which turns external forcing of salinit on
115     C and off.
116 cnh 1.8 C implicitFreeSurface - Set to true to use implcit free surface
117     C rigidLid - Set to true to use rigid lid
118 cnh 1.10 C momStepping - Turns momentum equation time-stepping off
119 cnh 1.15 C tempStepping - Turns temperature equation time-stepping off
120 adcroft 1.17 C saltStepping - Turns salinity equation time-stepping off
121 cnh 1.15 C useConstantF - Coriolis parameter set to f0
122     C useBetaPlaneF - Coriolis parameter set to f0 + beta.y
123     C useSphereF - Coriolis parameter set to 2.omega.sin(phi)
124 adcroft 1.14 C implicitDiffusion - Turns implicit vertical diffusion on
125 cnh 1.18 C doThetaClimRelax - Set true if relaxation to temperature
126     C climatology is required.
127     C doSaltClimRelax - Set true if relaxation to salinity
128     C climatology is required.
129 adcroft 1.19 C periodicExternalForcing - Set true if forcing is time-dependant
130 cnh 1.1 COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,
131 cnh 1.9 & momViscosity, momAdvection, momForcing, useCoriolis, momPressureForcing,
132 cnh 1.1 & tempDiffusion, tempAdvection, tempForcing,
133 cnh 1.8 & saltDiffusion, saltAdvection, saltForcing,
134 adcroft 1.17 & implicitFreeSurface, rigidLid,
135     & momStepping, tempStepping, saltStepping,
136 cnh 1.15 & metricTerms, usingSphericalPolarMTerms,
137     & useConstantF, useBetaPlaneF, useSphereF,
138 adcroft 1.19 & implicitDiffusion, doThetaClimRelax, doSaltClimRelax,
139     & periodicExternalForcing
140 cnh 1.1 LOGICAL usingCartesianGrid
141     LOGICAL usingSphericalPolarGrid
142 cnh 1.15 LOGICAL usingSphericalPolarMTerms
143 cnh 1.1 LOGICAL momViscosity
144     LOGICAL momAdvection
145     LOGICAL momForcing
146 cnh 1.9 LOGICAL momPressureForcing
147 cnh 1.1 LOGICAL useCoriolis
148     LOGICAL tempDiffusion
149     LOGICAL tempAdvection
150     LOGICAL tempForcing
151     LOGICAL saltDiffusion
152     LOGICAL saltAdvection
153     LOGICAL saltForcing
154 cnh 1.8 LOGICAL implicitFreeSurface
155     LOGICAL rigidLid
156 cnh 1.10 LOGICAL momStepping
157     LOGICAL tempStepping
158 adcroft 1.17 LOGICAL saltStepping
159 cnh 1.15 LOGICAL metricTerms
160     LOGICAL useConstantF
161     LOGICAL useBetaPlaneF
162     LOGICAL useSphereF
163 adcroft 1.14 LOGICAL implicitDiffusion
164 cnh 1.18 LOGICAL doThetaClimRelax
165     LOGICAL doSaltClimRelax
166 adcroft 1.19 LOGICAL periodicExternalForcing
167 cnh 1.1
168     C-- COMMON /PARM_R/ "Real" valued parameters used by the model.
169     C cg2dTargetResidual
170     C - Target residual for cg2d solver.
171 cnh 1.7 C cg2dpcOffDFac - Averaging weight for preconditioner off-diagonal.
172     C Note. 20th May 1998
173     C I made a weird discovery! In the model paper we argue
174     C for the form of the preconditioner used here ( see
175     C A Finite-volume, Incompressible Navier-Stokes Model
176     C ...., Marshall et. al ). The algebra gives a simple
177     C 0.5 factor for the averaging of ac and aCw to get a
178     C symmettric pre-conditioner. By using a factor of 0.51
179     C i.e. scaling the off-diagonal terms in the
180     C preconditioner down slightly I managed to get the
181     C number of iterations for convergence in a test case to
182     C drop form 192 -> 134! Need to investigate this further!
183     C For now I have introduced a parameter cg2dpcOffDFac which
184     C defaults to 0.51 but can be set at runtime.
185 cnh 1.24 C delP - Vertical grid spacing ( Pa ).
186     C delZ - Vertical grid spacing ( m ).
187 cnh 1.23 C delR - Vertical grid spacing ( units of r ).
188 cnh 1.1 C delX - Separation between cell faces (m) or (deg), depending
189     C delY on input flags.
190     C gravity - Accel. due to gravity ( m/s^2 )
191 cnh 1.8 C gBaro - Accel. due to gravity used in barotropic equation ( m/s^2 )
192 cnh 1.1 C ronil - Reference density
193 cnh 1.24 C rhoConst - Vertically constant reference density
194 cnh 1.3 C startTime - Start time for model ( s )
195 cnh 1.1 C phiMin - Latitude of southern most cell face.
196     C thetaMin - Longitude of western most cell face (this
197     C is an "inert" parameter but it is included
198     C to make geographical references simple.)
199 cnh 1.3 C rSphere - Radius of sphere for a spherical polar grid ( m ).
200 cnh 1.25 C recip_RSphere - Reciprocal radius of sphere ( m ).
201 cnh 1.1 C f0 - Reference coriolis parameter ( 1/s )
202     C ( Southern edge f for beta plane )
203     C beta - df/dy ( s^-1.m^-1 )
204 cnh 1.15 C omega - Angular velocity ( rad/s )
205 cnh 1.1 C viscAh - Eddy viscosity coeff. for mixing of
206     C momentum laterally ( m^2/s )
207     C viscAz - Eddy viscosity coeff. for mixing of
208     C momentum vertically ( m^2/s )
209 cnh 1.24 C viscAr - Eddy viscosity coeff. for mixing of
210     C momentum vertically ( units of r^2/s )
211 cnh 1.1 C viscA4 - Biharmonic viscosity coeff. for mixing of
212     C momentum laterally ( m^4/s )
213     C diffKhT - Laplacian diffusion coeff. for mixing of
214     C heat laterally ( m^2/s )
215     C diffKzT - Laplacian diffusion coeff. for mixing of
216     C heat vertically ( m^2/s )
217 cnh 1.24 C diffKrT - Laplacian diffusion coeff. for mixing of
218     C heat vertically ( units of r^2/s )
219 cnh 1.1 C diffK4T - Biharmonic diffusion coeff. for mixing of
220     C heat laterally ( m^4/s )
221     C diffKhS - Laplacian diffusion coeff. for mixing of
222     C salt laterally ( m^2/s )
223     C diffKzS - Laplacian diffusion coeff. for mixing of
224     C salt vertically ( m^2/s )
225 cnh 1.24 C diffKrS - Laplacian diffusion coeff. for mixing of
226     C salt vertically ( units of r^2/s )
227 cnh 1.1 C diffK4S - Biharmonic diffusion coeff. for mixing of
228     C salt laterally ( m^4/s )
229 cnh 1.3 C deltaT - Default timestep ( s )
230 cnh 1.7 C deltaTClock - Timestep used as model "clock". This determines the
231     C IO frequencies and is used in tagging output. It can
232     C be totally different to the dynamical time. Typically
233     C it will be the deep-water timestep for accelerated runs.
234     C Frequency of checkpointing and dumping of the model state
235     C are referenced to this clock. ( s )
236 cnh 1.3 C deltaTMom - Timestep for momemtum equations ( s )
237     C deltaTtracer - Timestep for tracer equations ( s )
238 cnh 1.8 C freesurfFac - Parameter to turn implicit free surface term on or off
239     C freesurfac = 1. uses implicit free surface
240     C freesurfac = 0. uses rigid lid
241 cnh 1.24 C hFacMin - Minimum fraction size of a cell (affects hFacC etc...)
242     C hFacMinDz - Minimum dimesional size of a cell (affects hFacC etc..., m)
243     C hFacMinDr - Minimum dimesional size of a cell (affects hFacC etc..., units of r)
244     C tauCD - CD scheme coupling timescale ( 1/s )
245     C rCD - CD scheme normalised coupling parameter ( 0-1 )
246     C GMmaxslope - max. slope allowed in GM/Redi tensor
247     C GMlength - Length to use in Visbeck et al. formula for K (m)
248     C GMalpha - alpha to use in Visbeck et al. formula for K
249     C GMdepth - Depth over which to integrate Richardson # (Visbeck et al.)
250     C GMkbackground - background value of GM/Redi diffusion coefficient ( m^2/s )
251     C GMmaxval - max. value of KapGM allowed in GM/Redi scheme ( m^2/s )
252     C startTime - Starting time for this integration ( s ).
253     C endTime - Ending time for this integration ( s ).
254     C chkPtFreq - Frequency of rolling check pointing ( s ).
255     C pChkPtFreq - Frequency of permanent check pointing ( s ).
256     C dumpFreq - Frequency with which model state is written to
257     C post-processing files ( s ).
258     C afFacMom - Advection of momentum term tracer parameter
259     C vfFacMom - Momentum viscosity tracer parameter
260     C pfFacMom - Momentum pressure forcing tracer parameter
261     C cfFacMom - Coriolis term tracer parameter
262     C foFacMom - Momentum forcing tracer parameter
263     C mtFacMom - Metric terms tracer parameter
264     C cAdjFreq - Frequency of convective adjustment
265     C
266 adcroft 1.20 C taveFreq - Frequency with which time-averaged model state is written to
267     C post-processing files ( s ).
268 cnh 1.18 C tauThetaClimRelax - Relaxation to climatology time scale ( s ).
269     C lambdaThetaClimRelax - Inverse time scale for relaxation ( 1/s ).
270     C tauSaltClimRelax - Relaxation to climatology time scale ( s ).
271     C lambdaSaltClimRelax - Inverse time scale for relaxation ( 1/s ).
272 adcroft 1.19 C externForcingPeriod - Is the period of which forcing varies (eg. 1 month)
273     C externForcingCycle - Is the repeat time of the forcing (eg. 1 year)
274     C (note: externForcingCycle must be an integer
275     C number times externForcingPeriod)
276 cnh 1.24 COMMON /PARM_R/ cg2dTargetResidual, cg2dpcOffDFac,
277     & delP, delZ, delR, delX, delY,
278 cnh 1.7 & deltaT,deltaTmom, deltaTtracer, deltaTClock,abeps, startTime, phiMin,
279 cnh 1.25 & thetaMin, rSphere, recip_RSphere, f0, fCori, beta,
280 cnh 1.24 & viscAh, viscAz, viscA4, viscAr,
281     & diffKhT, diffKzT, diffK4T, diffKrT,
282 cnh 1.25 & diffKhS, diffKzS, diffK4S, diffKrS,
283 cnh 1.24 & delT, tauCD, rCD, freeSurfFac, hFacMin, hFacMinDz,
284 adcroft 1.21 & GMmaxslope,GMlength,GMalpha,GMdepth,GMkbackground,GMmaxval,
285 cnh 1.25 & gravity, gBaro, rhonil, recip_rhonil, rhoConst, tRef, sRef,
286 adcroft 1.20 & endTime, chkPtFreq, pchkPtFreq, dumpFreq, taveFreq,
287 cnh 1.15 & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
288 cnh 1.18 & cAdjFreq, omega, tauThetaClimRelax, lambdaThetaClimRelax,
289 adcroft 1.19 & tauSaltClimRelax, lambdaSaltClimRelax,
290     & externForcingCycle, externForcingPeriod
291 adcroft 1.6 _RL cg2dTargetResidual
292 cnh 1.7 _RL cg2dpcOffDFac
293 cnh 1.25 _RL delZ(Nr)
294     _RL delP(Nr)
295     _RL delR(Nr)
296 adcroft 1.6 _RL delX(Nx)
297     _RL delY(Ny)
298     _RL deltaT
299 cnh 1.7 _RL deltaTClock
300 adcroft 1.6 _RL deltaTmom
301     _RL deltaTtracer
302     _RL abeps
303     _RL phiMin
304     _RL thetaMin
305     _RL rSphere
306 cnh 1.25 _RL recip_RSphere
307 adcroft 1.6 _RL f0
308 cnh 1.8 _RL freeSurfFac
309 adcroft 1.22 _RL hFacMin
310     _RL hFacMinDz
311 adcroft 1.6 _RL beta
312     _RL viscAh
313     _RL viscAz
314 cnh 1.25 _RL viscAr
315 adcroft 1.6 _RL viscA4
316     _RL diffKhT
317 cnh 1.25 _RL diffKrT
318 adcroft 1.6 _RL diffKzT
319     _RL diffK4T
320     _RL diffKhS
321 cnh 1.25 _RL diffKrS
322 adcroft 1.6 _RL diffKzS
323     _RL diffK4S
324     _RL delt
325     _RL tauCD
326     _RL rCD
327     _RL GMmaxslope
328     _RL GMlength
329     _RL GMalpha
330     _RL GMdepth
331     _RL GMkbackground
332 adcroft 1.21 _RL GMmaxval
333 adcroft 1.6 _RL gravity
334 cnh 1.8 _RL gBaro
335 adcroft 1.6 _RL rhonil
336 cnh 1.25 _RL recip_rhonil
337     _RL rhoConst
338     _RL tRef(Nr)
339     _RL sRef(Nr)
340 cnh 1.11 _RS Fcori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
341 adcroft 1.6 _RL startTime
342     _RL endTime
343     _RL chkPtFreq
344 cnh 1.7 _RL pChkPtFreq
345 adcroft 1.6 _RL dumpFreq
346 adcroft 1.20 _RL taveFreq
347 cnh 1.9 _RL afFacMom
348     _RL vfFacMom
349     _RL pfFacMom
350     _RL cfFacMom
351     _RL foFacMom
352 cnh 1.15 _RL mTFacMom
353 cnh 1.9 _RL cAdjFreq
354 cnh 1.15 _RL omega
355 cnh 1.18 _RL tauThetaClimRelax
356     _RL lambdaThetaClimRelax
357     _RL tauSaltClimRelax
358     _RL lambdaSaltClimRelax
359 adcroft 1.19 _RL externForcingCycle
360     _RL externForcingPeriod
361 adcroft 1.6
362 cnh 1.1 COMMON /PARM_A/ HeatCapacity_Cp,
363     & Lamba_theta
364 adcroft 1.6 _RL HeatCapacity_Cp
365     _RL Lamba_theta
366 adcroft 1.5
367     C Equation of State (polynomial coeffients)
368     COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS
369 cnh 1.25 _RL eosC(9,Nr+1),eosSig0(Nr+1),eosRefT(Nr+1),eosRefS(Nr+1)
370 adcroft 1.5 C Linear equation of state
371     C tAlpha - Linear EOS thermal expansion coefficient ( 1/degree ).
372     C sBeta - Linear EOS haline contraction coefficient.
373 adcroft 1.12 COMMON /PARM_EOS_LIN/ tAlpha,sBeta,eosType
374 adcroft 1.6 _RL tAlpha
375     _RL sBeta
376 adcroft 1.12 character*(6) eosType
377 cnh 1.7

  ViewVC Help
Powered by ViewVC 1.1.22