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

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

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

revision 1.20 by adcroft, Thu Jun 18 13:56:06 1998 UTC revision 1.32 by cnh, Wed Sep 9 15:19:07 1998 UTC
# Line 26  C                dumping model state. Line 26  C                dumping model state.
26        PARAMETER ( precFloat32 = 0 )        PARAMETER ( precFloat32 = 0 )
27        INTEGER precFloat64        INTEGER precFloat64
28        PARAMETER ( precFloat64 = 1 )        PARAMETER ( precFloat64 = 1 )
29    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    
41  C     Checkpoint data  C     Checkpoint data
42        INTEGER maxNoChkptLev        INTEGER maxNoChkptLev
# Line 44  C     thetaClimFile   - File containing Line 55  C     thetaClimFile   - File containing
55  C                       in relaxation term -lambda(theta-theta*)  C                       in relaxation term -lambda(theta-theta*)
56  C     saltClimFile    - File containing salt climataology used  C     saltClimFile    - File containing salt climataology used
57  C                       in relaxation term -lambda(salt-salt*)  C                       in relaxation term -lambda(salt-salt*)
58    C     buoyancyRelation - Flag used to indicate which relation to use to
59    C                        get buoyancy.
60        COMMON /PARM_C/ checkPtSuff,        COMMON /PARM_C/ checkPtSuff,
61       &                bathyFile, hydrogThetaFile, hydrogSaltFile,       &                bathyFile, hydrogThetaFile, hydrogSaltFile,
62       &                zonalWindFile, meridWindFile, thetaClimFile,       &                zonalWindFile, meridWindFile, thetaClimFile,
63       &                saltClimFile       &                saltClimFile, buoyancyRelation
64        CHARACTER*(5) checkPtSuff(maxNoChkptLev)        CHARACTER*(5) checkPtSuff(maxNoChkptLev)
65        CHARACTER*(MAX_LEN_FNAM) bathyFile        CHARACTER*(MAX_LEN_FNAM) bathyFile
66        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile
# Line 56  C                       in relaxation te Line 69  C                       in relaxation te
69        CHARACTER*(MAX_LEN_FNAM) meridWindFile        CHARACTER*(MAX_LEN_FNAM) meridWindFile
70        CHARACTER*(MAX_LEN_FNAM) thetaClimFile        CHARACTER*(MAX_LEN_FNAM) thetaClimFile
71        CHARACTER*(MAX_LEN_FNAM) saltClimFile        CHARACTER*(MAX_LEN_FNAM) saltClimFile
72          CHARACTER*(MAX_LEN_FNAM) buoyancyRelation
73    
74  C--   COMMON /PARM_I/ Integer valued parameters used by the model.  C--   COMMON /PARM_I/ Integer valued parameters used by the model.
75  C     cg2dMaxIters        - Maximum number of iterations in the  C     cg2dMaxIters        - Maximum number of iterations in the
# Line 127  C                        climatology is Line 141  C                        climatology is
141  C     doSaltClimRelax  - Set true if relaxation to salinity  C     doSaltClimRelax  - Set true if relaxation to salinity
142  C                        climatology is required.  C                        climatology is required.
143  C     periodicExternalForcing - Set true if forcing is time-dependant  C     periodicExternalForcing - Set true if forcing is time-dependant
144    C     usingPCoords     - Set to indicate that we are working in pressure
145    C                        coords.
146    C     usingZCoords     - Set to indicate that we are working in height
147    C                        coords.
148        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,
149       & momViscosity, momAdvection, momForcing, useCoriolis, momPressureForcing,       & momViscosity, momAdvection, momForcing, useCoriolis, momPressureForcing,
150       & tempDiffusion, tempAdvection, tempForcing,       & tempDiffusion, tempAdvection, tempForcing,
# Line 136  C     periodicExternalForcing - Set true Line 154  C     periodicExternalForcing - Set true
154       & metricTerms, usingSphericalPolarMTerms,       & metricTerms, usingSphericalPolarMTerms,
155       & useConstantF, useBetaPlaneF, useSphereF,       & useConstantF, useBetaPlaneF, useSphereF,
156       & implicitDiffusion, doThetaClimRelax, doSaltClimRelax,       & implicitDiffusion, doThetaClimRelax, doSaltClimRelax,
157       & periodicExternalForcing       & periodicExternalForcing, usingPCoords, usingZCoords
158        LOGICAL usingCartesianGrid        LOGICAL usingCartesianGrid
159        LOGICAL usingSphericalPolarGrid        LOGICAL usingSphericalPolarGrid
160        LOGICAL usingSphericalPolarMTerms        LOGICAL usingSphericalPolarMTerms
# Line 164  C     periodicExternalForcing - Set true Line 182  C     periodicExternalForcing - Set true
182        LOGICAL doThetaClimRelax        LOGICAL doThetaClimRelax
183        LOGICAL doSaltClimRelax        LOGICAL doSaltClimRelax
184        LOGICAL periodicExternalForcing        LOGICAL periodicExternalForcing
185          LOGICAL usingPCoords
186          LOGICAL usingZCoords
187    
188  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.
189  C     cg2dTargetResidual  C     cg2dTargetResidual
# Line 182  C           number of iterations for con Line 202  C           number of iterations for con
202  C           drop form 192 -> 134! Need to investigate this further!  C           drop form 192 -> 134! Need to investigate this further!
203  C           For now I have introduced a parameter cg2dpcOffDFac which  C           For now I have introduced a parameter cg2dpcOffDFac which
204  C           defaults to 0.51 but can be set at runtime.  C           defaults to 0.51 but can be set at runtime.
205  C     delZ      - Vertical grid spacing ( m ) - delZ is the distance  C     delP      - Vertical grid spacing ( Pa ).
206  C                 between "w" surfaces.  C     delZ      - Vertical grid spacing ( m  ).
207    C     delR      - Vertical grid spacing ( units of r ).
208  C     delX      - Separation between cell faces (m) or (deg), depending  C     delX      - Separation between cell faces (m) or (deg), depending
209  C     delY        on input flags.  C     delY        on input flags.
210  C     gravity   - Accel. due to gravity ( m/s^2 )  C     gravity   - Accel. due to gravity ( m/s^2 )
211    C     recip_gravity and its inverse
212  C     gBaro     - Accel. due to gravity used in barotropic equation ( m/s^2 )  C     gBaro     - Accel. due to gravity used in barotropic equation ( m/s^2 )
213  C     ronil     - Reference density  C     ronil     - Reference density
214    C     rhoConst  - Vertically constant reference density
215  C     startTime - Start time for model ( s )  C     startTime - Start time for model ( s )
216  C     phiMin    - Latitude of southern most cell face.  C     phiMin    - Latitude of southern most cell face.
217  C     thetaMin  - Longitude of western most cell face (this  C     thetaMin  - Longitude of western most cell face (this
218  C                 is an "inert" parameter but it is included  C                 is an "inert" parameter but it is included
219  C                 to make geographical references simple.)  C                 to make geographical references simple.)
220  C     rSphere   - Radius of sphere for a spherical polar grid ( m ).  C     rSphere   - Radius of sphere for a spherical polar grid ( m ).
221  C     rRSphere  - Reciprocal radius of sphere for a spherical polar grid ( m ).  C     recip_RSphere  - Reciprocal radius of sphere ( m ).
222  C     f0        - Reference coriolis parameter ( 1/s )  C     f0        - Reference coriolis parameter ( 1/s )
223  C                 ( Southern edge f for beta plane )  C                 ( Southern edge f for beta plane )
224  C     beta      - df/dy ( s^-1.m^-1 )  C     beta      - df/dy ( s^-1.m^-1 )
# Line 204  C     viscAh    - Eddy viscosity coeff. Line 227  C     viscAh    - Eddy viscosity coeff.
227  C                 momentum laterally ( m^2/s )  C                 momentum laterally ( m^2/s )
228  C     viscAz    - Eddy viscosity coeff. for mixing of  C     viscAz    - Eddy viscosity coeff. for mixing of
229  C                 momentum vertically ( m^2/s )  C                 momentum vertically ( m^2/s )
230    C     viscAp    - Eddy viscosity coeff. for mixing of
231    C                 momentum vertically ( Pa^2/s )
232    C     viscAr    - Eddy viscosity coeff. for mixing of
233    C                 momentum vertically ( units of r^2/s )
234  C     viscA4    - Biharmonic viscosity coeff. for mixing of  C     viscA4    - Biharmonic viscosity coeff. for mixing of
235  C                 momentum laterally ( m^4/s )  C                 momentum laterally ( m^4/s )
236  C     diffKhT   - Laplacian diffusion coeff. for mixing of  C     diffKhT   - Laplacian diffusion coeff. for mixing of
237  C                 heat laterally ( m^2/s )  C                 heat laterally ( m^2/s )
238  C     diffKzT   - Laplacian diffusion coeff. for mixing of  C     diffKzT   - Laplacian diffusion coeff. for mixing of
239  C                 heat vertically ( m^2/s )  C                 heat vertically ( m^2/s )
240    C     diffKpT   - Laplacian diffusion coeff. for mixing of
241    C                 heat vertically ( Pa^2/s )
242    C     diffKrT   - Laplacian diffusion coeff. for mixing of
243    C                 heat vertically ( units of r^2/s )
244  C     diffK4T   - Biharmonic diffusion coeff. for mixing of  C     diffK4T   - Biharmonic diffusion coeff. for mixing of
245  C                 heat laterally ( m^4/s )  C                 heat laterally ( m^4/s )
246  C     diffKhS  -  Laplacian diffusion coeff. for mixing of  C     diffKhS  -  Laplacian diffusion coeff. for mixing of
247  C                 salt laterally ( m^2/s )  C                 salt laterally ( m^2/s )
248  C     diffKzS   - Laplacian diffusion coeff. for mixing of  C     diffKzS   - Laplacian diffusion coeff. for mixing of
249  C                 salt vertically ( m^2/s )  C                 salt vertically ( m^2/s )
250    C     diffKpS   - Laplacian diffusion coeff. for mixing of
251    C                 salt vertically ( Pa^2/s )
252    C     diffKrS   - Laplacian diffusion coeff. for mixing of
253    C                 salt vertically ( units of r^2/s )
254  C     diffK4S   - Biharmonic diffusion coeff. for mixing of  C     diffK4S   - Biharmonic diffusion coeff. for mixing of
255  C                 salt laterally ( m^4/s )  C                 salt laterally ( m^4/s )
256  C     deltaT    - Default timestep ( s )  C     deltaT    - Default timestep ( s )
# Line 230  C     deltaTtracer - Timestep for tracer Line 265  C     deltaTtracer - Timestep for tracer
265  C     freesurfFac  - Parameter to turn implicit free surface term on or off  C     freesurfFac  - Parameter to turn implicit free surface term on or off
266  C                    freesurfac = 1. uses implicit free surface  C                    freesurfac = 1. uses implicit free surface
267  C                    freesurfac = 0. uses rigid lid  C                    freesurfac = 0. uses rigid lid
268  C     tauCD     - CD scheme coupling timescale ( 1/s )  C     hFacMin      - Minimum fraction size of a cell (affects hFacC etc...)
269  C     rCD       - CD scheme normalised coupling parameter ( 0-1 )  C     hFacMinDz    - Minimum dimesional size of a cell (affects hFacC etc..., m)
270  C     GMmaxslope  - max. slope allowed in GM/Redi tensor  C     hFacMinDp    - Minimum dimesional size of a cell (affects hFacC etc..., Pa)
271  C     GMlength  - Length to use in Visbeck et al. formula for K (m)  C     hFacMinDr    - Minimum dimesional size of a cell (affects hFacC etc..., units of r)
272  C     GMalpha   - alpha to use in Visbeck et al. formula for K  C     tauCD        - CD scheme coupling timescale ( 1/s )
273  C     GMdepth   - Depth over which to integrate Richardson # (Visbeck et al.)  C     rCD          - CD scheme normalised coupling parameter ( 0-1 )
274  C     GMkbackground - background value of GM/Redi coefficient  C     GMmaxslope    - max. slope allowed in GM/Redi tensor
275  C     startTime - Starting time for this integration ( s ).  C     GMlength      - Length to use in Visbeck et al. formula for K (m)
276  C     endTime   - Ending time for this integration ( s ).  C     GMalpha       - alpha to use in Visbeck et al. formula for K
277  C     chkPtFreq  - Frequency of rolling check pointing ( s ).  C     GMdepth       - Depth over which to integrate Richardson # (Visbeck et al.)
278  C     pChkPtFreq - Frequency of permanent check pointing ( s ).  C     GMkbackground - background value of GM/Redi diffusion coefficient ( m^2/s )
279  C     dumpFreq  - Frequency with which model state is written to  C     GMmaxval      - max. value of KapGM allowed in GM/Redi scheme ( m^2/s )
280  C                 post-processing files ( s ).  C     startTime     - Starting time for this integration ( s ).
281  C     taveFreq  - Frequency with which time-averaged model state is written to  C     endTime       - Ending time for this integration ( s ).
282  C                 post-processing files ( s ).  C     chkPtFreq     - Frequency of rolling check pointing ( s ).
283  C     afFacMom  - Advection of momentum term scaling parameter  C     pChkPtFreq    - Frequency of permanent check pointing ( s ).
284  C     vfFacMom  - Momentum viscosity scaling parameter  C     dumpFreq      - Frequency with which model state is written to
285  C     pfFacMom  - Momentum pressure forcing parameter  C                     post-processing files ( s ).
286  C     cfFacMom  - Coriolis term scaling parameter  C     afFacMom      - Advection of momentum term tracer parameter
287  C     foFacMom  - Momentum forcing scaling parameter  C     vfFacMom      - Momentum viscosity tracer parameter
288  C     mtFacMom  - Metric terms scaling parameter  C     pfFacMom      - Momentum pressure forcing tracer parameter
289  C     cAdjFreq  - Frequency of convective adjustment  C     cfFacMom      - Coriolis term tracer parameter
290    C     foFacMom      - Momentum forcing tracer parameter
291    C     mtFacMom      - Metric terms tracer parameter
292    C     cAdjFreq      - Frequency of convective adjustment
293    C
294    C     taveFreq      - Frequency with which time-averaged model state is written to
295    C                     post-processing files ( s ).
296  C     tauThetaClimRelax - Relaxation to climatology time scale ( s ).  C     tauThetaClimRelax - Relaxation to climatology time scale ( s ).
297  C     lambdaThetaClimRelax - Inverse time scale for relaxation ( 1/s ).  C     lambdaThetaClimRelax - Inverse time scale for relaxation ( 1/s ).
298  C     tauSaltClimRelax - Relaxation to climatology time scale ( s ).  C     tauSaltClimRelax - Relaxation to climatology time scale ( s ).
# Line 260  C     externForcingPeriod - Is the perio Line 301  C     externForcingPeriod - Is the perio
301  C     externForcingCycle - Is the repeat time of the forcing (eg. 1 year)  C     externForcingCycle - Is the repeat time of the forcing (eg. 1 year)
302  C                          (note: externForcingCycle must be an integer  C                          (note: externForcingCycle must be an integer
303  C                           number times externForcingPeriod)  C                           number times externForcingPeriod)
304        COMMON /PARM_R/ cg2dTargetResidual, cg2dpcOffDFac, delZ, delX, delY,  C     horiVertRatio      - Ratio on units in vertical to units in horizontal.
305    C     recip_horiVertRatio  ( 1 if horiz in m and vertical in m ).
306    C                          ( g*rho if horiz in m and vertical in Pa ).
307          COMMON /PARM_R/ cg2dTargetResidual, cg2dpcOffDFac,
308         & delP, delZ, delR, delX, delY,
309       & deltaT,deltaTmom, deltaTtracer, deltaTClock,abeps, startTime, phiMin,       & deltaT,deltaTmom, deltaTtracer, deltaTClock,abeps, startTime, phiMin,
310       & thetaMin, rSphere, rRSphere, f0, fCori, beta, viscAh, viscAz, viscA4,       & thetaMin, rSphere, recip_RSphere, f0, fCori, beta,
311       & diffKhT, diffKzT, diffK4T, diffKhS, diffKzS, diffK4S, delT,       & viscAh,  viscAz,  viscA4,  viscAr,
312       & tauCD, rCD, freeSurfFac,       & diffKhT, diffKzT, diffK4T, diffKrT,
313       & GMmaxslope,GMlength,GMalpha,GMdepth,GMkbackground,       & diffKhS, diffKzS, diffK4S, diffKrS,
314       & gravity, gBaro, rhonil, tRef, sRef,       & delT, tauCD, rCD, freeSurfFac, hFacMin, hFacMinDz,
315         & GMmaxslope,GMlength,GMalpha,GMdepth,GMkbackground,GMmaxval,
316         & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,
317         & recip_rhoConst, rhoConst, tRef, sRef,
318       & endTime, chkPtFreq, pchkPtFreq, dumpFreq, taveFreq,       & endTime, chkPtFreq, pchkPtFreq, dumpFreq, taveFreq,
319       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
320       & cAdjFreq, omega, tauThetaClimRelax, lambdaThetaClimRelax,       & cAdjFreq, omega, tauThetaClimRelax, lambdaThetaClimRelax,
321       & tauSaltClimRelax, lambdaSaltClimRelax,       & tauSaltClimRelax, lambdaSaltClimRelax,
322       & externForcingCycle, externForcingPeriod       & externForcingCycle, externForcingPeriod,
323         & viscAp, diffKpT, diffKpS, hFacMinDr, hFacMinDp,
324         & theta_S, specVol_S, horiVertRatio, recip_horiVertRatio
325    
326        _RL cg2dTargetResidual        _RL cg2dTargetResidual
327        _RL cg2dpcOffDFac        _RL cg2dpcOffDFac
328        _RL delZ(Nz)        _RL delZ(Nr)
329          _RL delP(Nr)
330          _RL delR(Nr)
331        _RL delX(Nx)        _RL delX(Nx)
332        _RL delY(Ny)        _RL delY(Ny)
333        _RL deltaT        _RL deltaT
# Line 285  C                           number times Line 338  C                           number times
338        _RL phiMin        _RL phiMin
339        _RL thetaMin        _RL thetaMin
340        _RL rSphere        _RL rSphere
341        _RL rRSphere        _RL recip_RSphere
342        _RL f0        _RL f0
343        _RL freeSurfFac        _RL freeSurfFac
344          _RL hFacMin
345          _RL hFacMinDz
346          _RL hFacMinDp
347          _RL hFacMinDr
348        _RL beta        _RL beta
349        _RL viscAh        _RL viscAh
350        _RL viscAz        _RL viscAz
351          _RL viscAp
352          _RL viscAr
353        _RL viscA4        _RL viscA4
354        _RL diffKhT        _RL diffKhT
355          _RL diffKrT
356        _RL diffKzT        _RL diffKzT
357          _RL diffKpT
358        _RL diffK4T        _RL diffK4T
359        _RL diffKhS        _RL diffKhS
360          _RL diffKrS
361        _RL diffKzS        _RL diffKzS
362          _RL diffKpS
363        _RL diffK4S        _RL diffK4S
364        _RL delt        _RL delt
365        _RL tauCD        _RL tauCD
# Line 306  C                           number times Line 369  C                           number times
369        _RL GMalpha        _RL GMalpha
370        _RL GMdepth        _RL GMdepth
371        _RL GMkbackground        _RL GMkbackground
372          _RL GMmaxval
373        _RL gravity        _RL gravity
374          _RL recip_gravity
375        _RL gBaro        _RL gBaro
376        _RL rhonil        _RL rhonil
377        _RL tRef(Nz)        _RL recip_rhonil
378        _RL sRef(Nz)        _RL rhoConst
379          _RL recip_rhoConst
380          _RL specVol_S(Nr)
381          _RL tRef(Nr)
382          _RL theta_S(Nr)
383          _RL sRef(Nr)
384        _RS Fcori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS Fcori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
385        _RL startTime        _RL startTime
386        _RL endTime        _RL endTime
# Line 332  C                           number times Line 402  C                           number times
402        _RL lambdaSaltClimRelax        _RL lambdaSaltClimRelax
403        _RL externForcingCycle        _RL externForcingCycle
404        _RL externForcingPeriod        _RL externForcingPeriod
405          _RL horiVertRatio
406          _RL recip_horiVertRatio
407    
408        COMMON /PARM_A/ HeatCapacity_Cp,        COMMON /PARM_A/ HeatCapacity_Cp,
409       &                Lamba_theta       &                Lamba_theta
# Line 340  C                           number times Line 412  C                           number times
412    
413  C Equation of State (polynomial coeffients)  C Equation of State (polynomial coeffients)
414        COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS        COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS
415        _RL eosC(9,Nz+1),eosSig0(Nz+1),eosRefT(Nz+1),eosRefS(Nz+1)        _RL eosC(9,Nr+1),eosSig0(Nr+1),eosRefT(Nr+1),eosRefS(Nr+1)
416  C Linear equation of state  C Linear equation of state
417  C     tAlpha    - Linear EOS thermal expansion coefficient ( 1/degree ).  C     tAlpha    - Linear EOS thermal expansion coefficient ( 1/degree ).
418  C     sBeta     - Linear EOS haline contraction coefficient.  C     sBeta     - Linear EOS haline contraction coefficient.
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.32
  ViewVC Help
Powered by ViewVC 1.1.22