/[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.38 - (hide annotations) (download)
Wed May 5 18:32:34 1999 UTC (25 years ago) by adcroft
Branch: MAIN
CVS Tags: checkpoint21, checkpoint22, checkpoint23
Changes since 1.37: +7 -2 lines
File MIME type: text/plain
Implemented new I/O package (mdsio.F). This package does parallel
I/O in much the same way as dfile.F used to except it uses "direct
access" rather than (f77) unformatted sequential access.

Problems with dfile.F package included:
  o unnecessary memory use (each process had two global sized buffers)
  o inability to read the files it had written without post-processing
  o "tiled" files were tiled by process/thread rather than actual tiles
  o created huge numbers of files with no alternatives

Features of the mdsio.F package:
  o direct-access binary writes
  o no excessive memory use
  o ability to read/write from multiple record files
  o "tiled" files are based on "WRAPPER" tiles so that the number
    and content of files is independent of the number of threads
    and/or processes
  o option to create single "global" files rather than "tiled" files
  o ability to read both "global" and "tiled" files
    [Caveat: the tiling of files must match the model tiles]
  o checkpoints now use a single file per model section
    ie.  one file for the hydrostatic model core, one file
    for the non-hydrostatic extensions and one file for the C-D
    extensions
  o the mid-level I/O routines now is broken into more source files
    read_write_fld.F supplies basic I/O routines with the same interface
                     as the original I/O package
    read_write_rec.F supplies I/O routines which allow multiple records
    write_state.F    writes the model state
    checkpoint.F     supplies the read/write checkpoint routines

All the example input data has had to be modified to be direct-access.
Otherwise only routines that used I/O have been affected and not
all of those have been due to the continuity of arguments in
the read_write_fld.F routines.

What needs to be done?  We have to create a suite of conversion
utilities for users with old-style data. Also supply the option
for using old-style I/O, not just for die-hards but for reading
data too extensive to be converted. And more...

1 adcroft 1.38 C $Header: /u/gcmpack/models/MITgcmUV/model/inc/PARAMS.h,v 1.37 1999/03/22 15:54:03 adcroft Exp $
2 cnh 1.1 C
3     C /==========================================================\
4     C | PARAMS.h |
5     C | o Header file defining model "parameters". |
6     C |==========================================================|
7 cnh 1.33 C | The values from the model standard input file are |
8 cnh 1.1 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 adcroft 1.36 PARAMETER ( precFloat32 = 32 )
27 cnh 1.7 INTEGER precFloat64
28 adcroft 1.36 PARAMETER ( precFloat64 = 64 )
29 cnh 1.28 C UNSET_xxx - Used to indicate variables that have not been given a value
30     Real*8 UNSET_FLOAT8
31     PARAMETER ( UNSET_FLOAT8 = 1.234567D5 )
32     Real*4 UNSET_FLOAT4
33     PARAMETER ( UNSET_FLOAT4 = 1.234567E5 )
34     _RL UNSET_RL
35     PARAMETER ( UNSET_RL = 1.234567D5 )
36     _RS UNSET_RS
37     PARAMETER ( UNSET_RS = 1.234567E5 )
38     INTEGER UNSET_I
39     PARAMETER ( UNSET_I = 123456789 )
40 cnh 1.7
41     C Checkpoint data
42     INTEGER maxNoChkptLev
43     PARAMETER ( maxNoChkptLev = 2 )
44    
45 cnh 1.1 C-- COMMON /PARM_C/ Character valued parameters used by the model.
46 cnh 1.7 C checkPtSuff - List of checkpoint file suffices
47 cnh 1.16 C bathyFile - File containing bathymetry. If not defined bathymetry
48     C is taken from inline function.
49     C hydrogThetaFile - File containing initial hydrographic data for potential
50     C temperature.
51     C hydrogSaltFile - File containing initial hydrographic data for salinity.
52     C zonalWindFile - File containing zonal wind data
53     C meridWindFile - File containing meridional wind data
54 cnh 1.18 C thetaClimFile - File containing theta climataology used
55     C in relaxation term -lambda(theta-theta*)
56     C saltClimFile - File containing salt climataology used
57     C in relaxation term -lambda(salt-salt*)
58 adcroft 1.37 C surfQfile - File containing surface heat flux
59     C EmPmRfile - File containing surface fresh water flux
60 cnh 1.26 C buoyancyRelation - Flag used to indicate which relation to use to
61     C get buoyancy.
62 cnh 1.16 COMMON /PARM_C/ checkPtSuff,
63     & bathyFile, hydrogThetaFile, hydrogSaltFile,
64 cnh 1.18 & zonalWindFile, meridWindFile, thetaClimFile,
65 adcroft 1.37 & saltClimFile, buoyancyRelation,
66     & EmPmRfile,surfQfile
67 cnh 1.7 CHARACTER*(5) checkPtSuff(maxNoChkptLev)
68 cnh 1.16 CHARACTER*(MAX_LEN_FNAM) bathyFile
69     CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile
70     CHARACTER*(MAX_LEN_FNAM) hydrogSaltFile
71     CHARACTER*(MAX_LEN_FNAM) zonalWindFile
72     CHARACTER*(MAX_LEN_FNAM) meridWindFile
73 cnh 1.18 CHARACTER*(MAX_LEN_FNAM) thetaClimFile
74     CHARACTER*(MAX_LEN_FNAM) saltClimFile
75 adcroft 1.37 CHARACTER*(MAX_LEN_FNAM) surfQfile
76     CHARACTER*(MAX_LEN_FNAM) EmPmRfile
77 cnh 1.26 CHARACTER*(MAX_LEN_FNAM) buoyancyRelation
78 cnh 1.7
79 cnh 1.1 C-- COMMON /PARM_I/ Integer valued parameters used by the model.
80     C cg2dMaxIters - Maximum number of iterations in the
81     C two-dimensional con. grad solver.
82     C cg2dChkResFreq - Frequency with which to check residual
83     C in con. grad solver.
84 cnh 1.33 C cg3dMaxIters - Maximum number of iterations in the
85     C three-dimensional con. grad solver.
86     C cg3dChkResFreq - Frequency with which to check residual
87     C in con. grad solver.
88 cnh 1.1 C nIter0 - Start time-step number of for this run
89     C nTimeSteps - Number of timesteps to execute
90     C numStepsPerPickup - For offline setup. Frequency of pickup
91     C of flow fields.
92 cnh 1.7 C writeStatePrec - Precision used for writing model state.
93     C writeBinaryPrec - Precision used for writing binary files
94     C readBinaryPrec - Precision used for reading binary files
95     C nCheckLev - Holds current checkpoint level
96 cnh 1.34 C nShap - "Order" of shapiro filter to apply to
97     C model prognositic fields.
98     C nShap == 1 => del2
99     C nShap == 2 => del4
100     C nShap == 3 => del6
101     C etc...
102    
103 cnh 1.1 COMMON /PARM_I/
104     & cg2dMaxIters,
105     & cg2dChkResFreq,
106 cnh 1.33 & cg3dMaxIters,
107     & cg3dChkResFreq,
108 cnh 1.1 & nIter0, nTimeSteps,
109 cnh 1.7 & numStepsPerPickup,
110     & writeStatePrec, nCheckLev,
111 cnh 1.34 & writeBinaryPrec, readBinaryPrec,
112     & nShap
113 cnh 1.1 INTEGER cg2dMaxIters
114     INTEGER cg2dChkResFreq
115 cnh 1.33 INTEGER cg3dMaxIters
116     INTEGER cg3dChkResFreq
117 cnh 1.1 INTEGER nIter0
118     INTEGER nTimeSteps
119     INTEGER numStepsPerPickup
120 cnh 1.7 INTEGER writeStatePrec
121     INTEGER writeBinaryPrec
122     INTEGER readBinaryPrec
123     INTEGER nCheckLev
124 cnh 1.34 INTEGER nShap
125 cnh 1.1
126     C-- COMMON /PARM_L/ Logical valued parameters used by the model.
127     C usingCartesianGrid - If TRUE grid generation will be in a cartesian
128     C coordinate frame.
129     C usingSphericalPolarGrid - If TRUE grid generation will be in a
130     C spherical polar frame.
131 adcroft 1.36 C no_slip_sides - Impose "no-slip" at lateral boundaries.
132     C no_slip_bottom- Impose "no-slip" at bottom boundary.
133 cnh 1.1 C momViscosity - Flag which turns momentum friction terms on and off.
134     C momAdvection - Flag which turns advection of momentum on and off.
135     C momForcing - Flag which turns external forcing of momentum on
136     C and off.
137 cnh 1.9 C momPressureForcing - Flag which turns pressure term in momentum equation
138     C on and off.
139 cnh 1.15 C metricTerms - Flag which turns metric terms on or off.
140     C usingSphericalPolarMTerms - If TRUE use spherical polar metric terms.
141 cnh 1.1 C useCoriolis - Flag which turns the coriolis terms on and off.
142     C tempDiffusion - Flag which turns diffusion of temperature on
143     C and off.
144     C tempAdvection - Flag which turns advection of temperature on
145     C and off.
146     C tempForcing - Flag which turns external forcing of temperature on
147     C and off.
148     C saltDiffusion - Flag which turns diffusion of salinit on
149     C and off.
150     C saltAdvection - Flag which turns advection of salinit on
151     C and off.
152     C saltForcing - Flag which turns external forcing of salinit on
153     C and off.
154 cnh 1.8 C implicitFreeSurface - Set to true to use implcit free surface
155     C rigidLid - Set to true to use rigid lid
156 cnh 1.10 C momStepping - Turns momentum equation time-stepping off
157 cnh 1.15 C tempStepping - Turns temperature equation time-stepping off
158 adcroft 1.17 C saltStepping - Turns salinity equation time-stepping off
159 cnh 1.15 C useConstantF - Coriolis parameter set to f0
160     C useBetaPlaneF - Coriolis parameter set to f0 + beta.y
161     C useSphereF - Coriolis parameter set to 2.omega.sin(phi)
162 adcroft 1.14 C implicitDiffusion - Turns implicit vertical diffusion on
163 cnh 1.18 C doThetaClimRelax - Set true if relaxation to temperature
164     C climatology is required.
165     C doSaltClimRelax - Set true if relaxation to salinity
166     C climatology is required.
167 adcroft 1.19 C periodicExternalForcing - Set true if forcing is time-dependant
168 cnh 1.29 C usingPCoords - Set to indicate that we are working in pressure
169     C coords.
170     C usingZCoords - Set to indicate that we are working in height
171     C coords.
172 adcroft 1.35 C openBoundaries - Using open-boundaries
173 adcroft 1.37 C nonHydrostatic - Using non-hydrostatic terms
174 adcroft 1.38 C globalFiles - Selects between "global" and "tiled" files
175     C allowFreezing - Allows water to freeze and form ice
176 cnh 1.1 COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,
177 adcroft 1.36 & no_slip_sides,no_slip_bottom,
178 cnh 1.33 & momViscosity, momAdvection, momForcing, useCoriolis,
179     & momPressureForcing,tempDiffusion, tempAdvection, tempForcing,
180 cnh 1.8 & saltDiffusion, saltAdvection, saltForcing,
181 adcroft 1.17 & implicitFreeSurface, rigidLid,
182     & momStepping, tempStepping, saltStepping,
183 cnh 1.15 & metricTerms, usingSphericalPolarMTerms,
184     & useConstantF, useBetaPlaneF, useSphereF,
185 adcroft 1.19 & implicitDiffusion, doThetaClimRelax, doSaltClimRelax,
186 adcroft 1.35 & periodicExternalForcing, usingPCoords, usingZCoords,
187 adcroft 1.38 & openBoundaries, nonHydrostatic, globalFiles,
188     & allowFreezing
189 cnh 1.1 LOGICAL usingCartesianGrid
190     LOGICAL usingSphericalPolarGrid
191 cnh 1.15 LOGICAL usingSphericalPolarMTerms
192 adcroft 1.36 LOGICAL no_slip_sides
193     LOGICAL no_slip_bottom
194 cnh 1.1 LOGICAL momViscosity
195     LOGICAL momAdvection
196     LOGICAL momForcing
197 cnh 1.9 LOGICAL momPressureForcing
198 cnh 1.1 LOGICAL useCoriolis
199     LOGICAL tempDiffusion
200     LOGICAL tempAdvection
201     LOGICAL tempForcing
202     LOGICAL saltDiffusion
203     LOGICAL saltAdvection
204     LOGICAL saltForcing
205 cnh 1.8 LOGICAL implicitFreeSurface
206     LOGICAL rigidLid
207 cnh 1.10 LOGICAL momStepping
208     LOGICAL tempStepping
209 adcroft 1.17 LOGICAL saltStepping
210 cnh 1.15 LOGICAL metricTerms
211     LOGICAL useConstantF
212     LOGICAL useBetaPlaneF
213     LOGICAL useSphereF
214 adcroft 1.14 LOGICAL implicitDiffusion
215 cnh 1.18 LOGICAL doThetaClimRelax
216     LOGICAL doSaltClimRelax
217 adcroft 1.19 LOGICAL periodicExternalForcing
218 cnh 1.29 LOGICAL usingPCoords
219     LOGICAL usingZCoords
220 adcroft 1.35 LOGICAL openBoundaries
221 adcroft 1.37 LOGICAL nonHydrostatic
222 adcroft 1.38 LOGICAL globalFiles
223     LOGICAL allowFreezing
224 cnh 1.1
225     C-- COMMON /PARM_R/ "Real" valued parameters used by the model.
226     C cg2dTargetResidual
227     C - Target residual for cg2d solver.
228 cnh 1.33 C cg3dTargetResidual
229     C - Target residual for cg3d solver.
230 cnh 1.7 C cg2dpcOffDFac - Averaging weight for preconditioner off-diagonal.
231     C Note. 20th May 1998
232     C I made a weird discovery! In the model paper we argue
233     C for the form of the preconditioner used here ( see
234     C A Finite-volume, Incompressible Navier-Stokes Model
235     C ...., Marshall et. al ). The algebra gives a simple
236     C 0.5 factor for the averaging of ac and aCw to get a
237     C symmettric pre-conditioner. By using a factor of 0.51
238     C i.e. scaling the off-diagonal terms in the
239     C preconditioner down slightly I managed to get the
240     C number of iterations for convergence in a test case to
241     C drop form 192 -> 134! Need to investigate this further!
242     C For now I have introduced a parameter cg2dpcOffDFac which
243     C defaults to 0.51 but can be set at runtime.
244 cnh 1.24 C delP - Vertical grid spacing ( Pa ).
245     C delZ - Vertical grid spacing ( m ).
246 cnh 1.23 C delR - Vertical grid spacing ( units of r ).
247 cnh 1.1 C delX - Separation between cell faces (m) or (deg), depending
248     C delY on input flags.
249     C gravity - Accel. due to gravity ( m/s^2 )
250 cnh 1.32 C recip_gravity and its inverse
251 cnh 1.8 C gBaro - Accel. due to gravity used in barotropic equation ( m/s^2 )
252 cnh 1.1 C ronil - Reference density
253 cnh 1.24 C rhoConst - Vertically constant reference density
254 cnh 1.3 C startTime - Start time for model ( s )
255 cnh 1.1 C phiMin - Latitude of southern most cell face.
256     C thetaMin - Longitude of western most cell face (this
257     C is an "inert" parameter but it is included
258     C to make geographical references simple.)
259 cnh 1.3 C rSphere - Radius of sphere for a spherical polar grid ( m ).
260 cnh 1.25 C recip_RSphere - Reciprocal radius of sphere ( m ).
261 cnh 1.1 C f0 - Reference coriolis parameter ( 1/s )
262     C ( Southern edge f for beta plane )
263     C beta - df/dy ( s^-1.m^-1 )
264 cnh 1.15 C omega - Angular velocity ( rad/s )
265 cnh 1.1 C viscAh - Eddy viscosity coeff. for mixing of
266     C momentum laterally ( m^2/s )
267     C viscAz - Eddy viscosity coeff. for mixing of
268     C momentum vertically ( m^2/s )
269 cnh 1.27 C viscAp - Eddy viscosity coeff. for mixing of
270     C momentum vertically ( Pa^2/s )
271 cnh 1.24 C viscAr - Eddy viscosity coeff. for mixing of
272     C momentum vertically ( units of r^2/s )
273 cnh 1.1 C viscA4 - Biharmonic viscosity coeff. for mixing of
274     C momentum laterally ( m^4/s )
275     C diffKhT - Laplacian diffusion coeff. for mixing of
276     C heat laterally ( m^2/s )
277     C diffKzT - Laplacian diffusion coeff. for mixing of
278     C heat vertically ( m^2/s )
279 cnh 1.27 C diffKpT - Laplacian diffusion coeff. for mixing of
280     C heat vertically ( Pa^2/s )
281 cnh 1.24 C diffKrT - Laplacian diffusion coeff. for mixing of
282     C heat vertically ( units of r^2/s )
283 cnh 1.1 C diffK4T - Biharmonic diffusion coeff. for mixing of
284     C heat laterally ( m^4/s )
285     C diffKhS - Laplacian diffusion coeff. for mixing of
286     C salt laterally ( m^2/s )
287     C diffKzS - Laplacian diffusion coeff. for mixing of
288     C salt vertically ( m^2/s )
289 cnh 1.27 C diffKpS - Laplacian diffusion coeff. for mixing of
290     C salt vertically ( Pa^2/s )
291 cnh 1.24 C diffKrS - Laplacian diffusion coeff. for mixing of
292     C salt vertically ( units of r^2/s )
293 cnh 1.1 C diffK4S - Biharmonic diffusion coeff. for mixing of
294     C salt laterally ( m^4/s )
295 cnh 1.3 C deltaT - Default timestep ( s )
296 cnh 1.7 C deltaTClock - Timestep used as model "clock". This determines the
297     C IO frequencies and is used in tagging output. It can
298     C be totally different to the dynamical time. Typically
299     C it will be the deep-water timestep for accelerated runs.
300     C Frequency of checkpointing and dumping of the model state
301     C are referenced to this clock. ( s )
302 cnh 1.3 C deltaTMom - Timestep for momemtum equations ( s )
303     C deltaTtracer - Timestep for tracer equations ( s )
304 cnh 1.8 C freesurfFac - Parameter to turn implicit free surface term on or off
305     C freesurfac = 1. uses implicit free surface
306     C freesurfac = 0. uses rigid lid
307 cnh 1.24 C hFacMin - Minimum fraction size of a cell (affects hFacC etc...)
308     C hFacMinDz - Minimum dimesional size of a cell (affects hFacC etc..., m)
309 cnh 1.28 C hFacMinDp - Minimum dimesional size of a cell (affects hFacC etc..., Pa)
310 cnh 1.24 C hFacMinDr - Minimum dimesional size of a cell (affects hFacC etc..., units of r)
311     C tauCD - CD scheme coupling timescale ( 1/s )
312     C rCD - CD scheme normalised coupling parameter ( 0-1 )
313 cnh 1.27 C GMmaxslope - max. slope allowed in GM/Redi tensor
314     C GMlength - Length to use in Visbeck et al. formula for K (m)
315     C GMalpha - alpha to use in Visbeck et al. formula for K
316     C GMdepth - Depth over which to integrate Richardson # (Visbeck et al.)
317 cnh 1.24 C GMkbackground - background value of GM/Redi diffusion coefficient ( m^2/s )
318     C GMmaxval - max. value of KapGM allowed in GM/Redi scheme ( m^2/s )
319     C startTime - Starting time for this integration ( s ).
320     C endTime - Ending time for this integration ( s ).
321     C chkPtFreq - Frequency of rolling check pointing ( s ).
322     C pChkPtFreq - Frequency of permanent check pointing ( s ).
323     C dumpFreq - Frequency with which model state is written to
324     C post-processing files ( s ).
325     C afFacMom - Advection of momentum term tracer parameter
326     C vfFacMom - Momentum viscosity tracer parameter
327     C pfFacMom - Momentum pressure forcing tracer parameter
328     C cfFacMom - Coriolis term tracer parameter
329     C foFacMom - Momentum forcing tracer parameter
330     C mtFacMom - Metric terms tracer parameter
331     C cAdjFreq - Frequency of convective adjustment
332     C
333 cnh 1.27 C taveFreq - Frequency with which time-averaged model state is written to
334     C post-processing files ( s ).
335 cnh 1.18 C tauThetaClimRelax - Relaxation to climatology time scale ( s ).
336     C lambdaThetaClimRelax - Inverse time scale for relaxation ( 1/s ).
337     C tauSaltClimRelax - Relaxation to climatology time scale ( s ).
338     C lambdaSaltClimRelax - Inverse time scale for relaxation ( 1/s ).
339 adcroft 1.19 C externForcingPeriod - Is the period of which forcing varies (eg. 1 month)
340     C externForcingCycle - Is the repeat time of the forcing (eg. 1 year)
341     C (note: externForcingCycle must be an integer
342     C number times externForcingPeriod)
343 cnh 1.30 C horiVertRatio - Ratio on units in vertical to units in horizontal.
344 cnh 1.31 C recip_horiVertRatio ( 1 if horiz in m and vertical in m ).
345 cnh 1.30 C ( g*rho if horiz in m and vertical in Pa ).
346 cnh 1.34 C latFFTFiltLo - Low latitude for FFT filtering of latitude
347     C circles ( see filter*.F )
348 cnh 1.24 COMMON /PARM_R/ cg2dTargetResidual, cg2dpcOffDFac,
349 cnh 1.33 & cg3dTargetResidual,
350 cnh 1.24 & delP, delZ, delR, delX, delY,
351 cnh 1.33 & deltaT,deltaTmom, deltaTtracer, deltaTClock,abeps, startTime,
352     & phiMin, thetaMin, rSphere, recip_RSphere, f0, fCori, beta,
353 cnh 1.24 & viscAh, viscAz, viscA4, viscAr,
354     & diffKhT, diffKzT, diffK4T, diffKrT,
355 cnh 1.25 & diffKhS, diffKzS, diffK4S, diffKrS,
356 cnh 1.27 & delT, tauCD, rCD, freeSurfFac, hFacMin, hFacMinDz,
357 adcroft 1.21 & GMmaxslope,GMlength,GMalpha,GMdepth,GMkbackground,GMmaxval,
358 cnh 1.32 & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,
359 cnh 1.27 & recip_rhoConst, rhoConst, tRef, sRef,
360 adcroft 1.20 & endTime, chkPtFreq, pchkPtFreq, dumpFreq, taveFreq,
361 cnh 1.15 & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
362 cnh 1.18 & cAdjFreq, omega, tauThetaClimRelax, lambdaThetaClimRelax,
363 adcroft 1.19 & tauSaltClimRelax, lambdaSaltClimRelax,
364 cnh 1.27 & externForcingCycle, externForcingPeriod,
365 cnh 1.28 & viscAp, diffKpT, diffKpS, hFacMinDr, hFacMinDp,
366 cnh 1.34 & theta_S, specVol_S, horiVertRatio, recip_horiVertRatio,
367     & latFFTFiltLo
368 cnh 1.27
369 adcroft 1.6 _RL cg2dTargetResidual
370 cnh 1.33 _RL cg3dTargetResidual
371 cnh 1.7 _RL cg2dpcOffDFac
372 cnh 1.25 _RL delZ(Nr)
373     _RL delP(Nr)
374     _RL delR(Nr)
375 adcroft 1.6 _RL delX(Nx)
376     _RL delY(Ny)
377     _RL deltaT
378 cnh 1.7 _RL deltaTClock
379 adcroft 1.6 _RL deltaTmom
380     _RL deltaTtracer
381     _RL abeps
382     _RL phiMin
383     _RL thetaMin
384     _RL rSphere
385 cnh 1.25 _RL recip_RSphere
386 adcroft 1.6 _RL f0
387 cnh 1.8 _RL freeSurfFac
388 adcroft 1.22 _RL hFacMin
389     _RL hFacMinDz
390 cnh 1.28 _RL hFacMinDp
391 cnh 1.27 _RL hFacMinDr
392 adcroft 1.6 _RL beta
393     _RL viscAh
394     _RL viscAz
395 cnh 1.27 _RL viscAp
396 cnh 1.25 _RL viscAr
397 adcroft 1.6 _RL viscA4
398     _RL diffKhT
399 cnh 1.25 _RL diffKrT
400 adcroft 1.6 _RL diffKzT
401 cnh 1.27 _RL diffKpT
402 adcroft 1.6 _RL diffK4T
403     _RL diffKhS
404 cnh 1.25 _RL diffKrS
405 adcroft 1.6 _RL diffKzS
406 cnh 1.27 _RL diffKpS
407 adcroft 1.6 _RL diffK4S
408     _RL delt
409     _RL tauCD
410     _RL rCD
411     _RL GMmaxslope
412     _RL GMlength
413     _RL GMalpha
414     _RL GMdepth
415     _RL GMkbackground
416 adcroft 1.21 _RL GMmaxval
417 adcroft 1.6 _RL gravity
418 cnh 1.32 _RL recip_gravity
419 cnh 1.8 _RL gBaro
420 adcroft 1.6 _RL rhonil
421 cnh 1.25 _RL recip_rhonil
422     _RL rhoConst
423 cnh 1.26 _RL recip_rhoConst
424     _RL specVol_S(Nr)
425 cnh 1.25 _RL tRef(Nr)
426 cnh 1.26 _RL theta_S(Nr)
427 cnh 1.25 _RL sRef(Nr)
428 cnh 1.11 _RS Fcori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
429 adcroft 1.6 _RL startTime
430     _RL endTime
431     _RL chkPtFreq
432 cnh 1.7 _RL pChkPtFreq
433 adcroft 1.6 _RL dumpFreq
434 adcroft 1.20 _RL taveFreq
435 cnh 1.9 _RL afFacMom
436     _RL vfFacMom
437     _RL pfFacMom
438     _RL cfFacMom
439     _RL foFacMom
440 cnh 1.15 _RL mTFacMom
441 cnh 1.9 _RL cAdjFreq
442 cnh 1.15 _RL omega
443 cnh 1.18 _RL tauThetaClimRelax
444     _RL lambdaThetaClimRelax
445     _RL tauSaltClimRelax
446     _RL lambdaSaltClimRelax
447 adcroft 1.19 _RL externForcingCycle
448     _RL externForcingPeriod
449 cnh 1.30 _RL horiVertRatio
450 cnh 1.31 _RL recip_horiVertRatio
451 cnh 1.34 _RL latFFTFiltLo
452 adcroft 1.6
453 adcroft 1.36 COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp,
454 cnh 1.1 & Lamba_theta
455 adcroft 1.6 _RL HeatCapacity_Cp
456     _RL Lamba_theta
457 adcroft 1.36 _RL recip_Cp
458 adcroft 1.5
459     C Equation of State (polynomial coeffients)
460     COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS
461 cnh 1.25 _RL eosC(9,Nr+1),eosSig0(Nr+1),eosRefT(Nr+1),eosRefS(Nr+1)
462 adcroft 1.5 C Linear equation of state
463     C tAlpha - Linear EOS thermal expansion coefficient ( 1/degree ).
464     C sBeta - Linear EOS haline contraction coefficient.
465 adcroft 1.12 COMMON /PARM_EOS_LIN/ tAlpha,sBeta,eosType
466 adcroft 1.6 _RL tAlpha
467     _RL sBeta
468 adcroft 1.12 character*(6) eosType
469 cnh 1.7
470 adcroft 1.35 C These are input arrays (of integers) that contain the *absolute*
471     C computational index of an open-boundary (OB) point.
472     C A zero (0) element means there is no corresponding OB in that column/row.
473     C The computational coordinate refers to "tracer" cells.
474     C For a northern/southern OB, the OB V point is to the south/north.
475     C For an eastern/western OB, the OB U point is to the west/east.
476     C eg.
477     C OB_Jnorth(3)=34 means that:
478     C T( 3 ,34) is a an OB point
479     C U(3:4,34) is a an OB point
480     C V( 4 ,34) is a an OB point
481     C while
482     C OB_Jsouth(3)=1 means that:
483     C T( 3 ,1) is a an OB point
484     C U(3:4,1) is a an OB point
485     C V( 4 ,2) is a an OB point
486     C
487     C For convenience, negative values for Jnorth/Ieast refer to
488     C points relative to the Northern/Eastern edges of the model
489     C eg. OB_Jnorth(3)=-1 means that the point (3,Ny-1) is a northern O-B.
490     C
491     COMMON /PARM_OB/
492     & OB_Jnorth,OB_Jsouth,OB_Ieast,OB_Iwest
493     INTEGER OB_Jnorth(Nx)
494     INTEGER OB_Jsouth(Nx)
495     INTEGER OB_Ieast(Ny)
496     INTEGER OB_Iwest(Ny)

  ViewVC Help
Powered by ViewVC 1.1.22