/[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.18 by cnh, Mon Jun 15 05:13:55 1998 UTC revision 1.27 by cnh, Mon Aug 24 02:25:01 1998 UTC
# Line 44  C     thetaClimFile   - File containing Line 44  C     thetaClimFile   - File containing
44  C                       in relaxation term -lambda(theta-theta*)  C                       in relaxation term -lambda(theta-theta*)
45  C     saltClimFile    - File containing salt climataology used  C     saltClimFile    - File containing salt climataology used
46  C                       in relaxation term -lambda(salt-salt*)  C                       in relaxation term -lambda(salt-salt*)
47    C     buoyancyRelation - Flag used to indicate which relation to use to
48    C                        get buoyancy.
49        COMMON /PARM_C/ checkPtSuff,        COMMON /PARM_C/ checkPtSuff,
50       &                bathyFile, hydrogThetaFile, hydrogSaltFile,       &                bathyFile, hydrogThetaFile, hydrogSaltFile,
51       &                zonalWindFile, meridWindFile, thetaClimFile,       &                zonalWindFile, meridWindFile, thetaClimFile,
52       &                saltClimFile       &                saltClimFile, buoyancyRelation
53        CHARACTER*(5) checkPtSuff(maxNoChkptLev)        CHARACTER*(5) checkPtSuff(maxNoChkptLev)
54        CHARACTER*(MAX_LEN_FNAM) bathyFile        CHARACTER*(MAX_LEN_FNAM) bathyFile
55        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile
# Line 56  C                       in relaxation te Line 58  C                       in relaxation te
58        CHARACTER*(MAX_LEN_FNAM) meridWindFile        CHARACTER*(MAX_LEN_FNAM) meridWindFile
59        CHARACTER*(MAX_LEN_FNAM) thetaClimFile        CHARACTER*(MAX_LEN_FNAM) thetaClimFile
60        CHARACTER*(MAX_LEN_FNAM) saltClimFile        CHARACTER*(MAX_LEN_FNAM) saltClimFile
61          CHARACTER*(MAX_LEN_FNAM) buoyancyRelation
62    
63  C--   COMMON /PARM_I/ Integer valued parameters used by the model.  C--   COMMON /PARM_I/ Integer valued parameters used by the model.
64  C     cg2dMaxIters        - Maximum number of iterations in the  C     cg2dMaxIters        - Maximum number of iterations in the
# Line 126  C     doThetaClimRelax - Set true if rel Line 129  C     doThetaClimRelax - Set true if rel
129  C                        climatology is required.  C                        climatology is required.
130  C     doSaltClimRelax  - Set true if relaxation to salinity  C     doSaltClimRelax  - Set true if relaxation to salinity
131  C                        climatology is required.  C                        climatology is required.
132    C     periodicExternalForcing - Set true if forcing is time-dependant
133        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,
134       & momViscosity, momAdvection, momForcing, useCoriolis, momPressureForcing,       & momViscosity, momAdvection, momForcing, useCoriolis, momPressureForcing,
135       & tempDiffusion, tempAdvection, tempForcing,       & tempDiffusion, tempAdvection, tempForcing,
# Line 134  C                        climatology is Line 138  C                        climatology is
138       & momStepping, tempStepping, saltStepping,       & momStepping, tempStepping, saltStepping,
139       & metricTerms, usingSphericalPolarMTerms,       & metricTerms, usingSphericalPolarMTerms,
140       & useConstantF, useBetaPlaneF, useSphereF,       & useConstantF, useBetaPlaneF, useSphereF,
141       & implicitDiffusion, doThetaClimRelax, doSaltClimRelax       & implicitDiffusion, doThetaClimRelax, doSaltClimRelax,
142         & periodicExternalForcing
143        LOGICAL usingCartesianGrid        LOGICAL usingCartesianGrid
144        LOGICAL usingSphericalPolarGrid        LOGICAL usingSphericalPolarGrid
145        LOGICAL usingSphericalPolarMTerms        LOGICAL usingSphericalPolarMTerms
# Line 161  C                        climatology is Line 166  C                        climatology is
166        LOGICAL implicitDiffusion        LOGICAL implicitDiffusion
167        LOGICAL doThetaClimRelax        LOGICAL doThetaClimRelax
168        LOGICAL doSaltClimRelax        LOGICAL doSaltClimRelax
169          LOGICAL periodicExternalForcing
170    
171  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.
172  C     cg2dTargetResidual  C     cg2dTargetResidual
# Line 179  C           number of iterations for con Line 185  C           number of iterations for con
185  C           drop form 192 -> 134! Need to investigate this further!  C           drop form 192 -> 134! Need to investigate this further!
186  C           For now I have introduced a parameter cg2dpcOffDFac which  C           For now I have introduced a parameter cg2dpcOffDFac which
187  C           defaults to 0.51 but can be set at runtime.  C           defaults to 0.51 but can be set at runtime.
188  C     delZ      - Vertical grid spacing ( m ) - delZ is the distance  C     delP      - Vertical grid spacing ( Pa ).
189  C                 between "w" surfaces.  C     delZ      - Vertical grid spacing ( m  ).
190    C     delR      - Vertical grid spacing ( units of r ).
191  C     delX      - Separation between cell faces (m) or (deg), depending  C     delX      - Separation between cell faces (m) or (deg), depending
192  C     delY        on input flags.  C     delY        on input flags.
193  C     gravity   - Accel. due to gravity ( m/s^2 )  C     gravity   - Accel. due to gravity ( m/s^2 )
194  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 )
195  C     ronil     - Reference density  C     ronil     - Reference density
196    C     rhoConst  - Vertically constant reference density
197  C     startTime - Start time for model ( s )  C     startTime - Start time for model ( s )
198  C     phiMin    - Latitude of southern most cell face.  C     phiMin    - Latitude of southern most cell face.
199  C     thetaMin  - Longitude of western most cell face (this  C     thetaMin  - Longitude of western most cell face (this
200  C                 is an "inert" parameter but it is included  C                 is an "inert" parameter but it is included
201  C                 to make geographical references simple.)  C                 to make geographical references simple.)
202  C     rSphere   - Radius of sphere for a spherical polar grid ( m ).  C     rSphere   - Radius of sphere for a spherical polar grid ( m ).
203  C     rRSphere  - Reciprocal radius of sphere for a spherical polar grid ( m ).  C     recip_RSphere  - Reciprocal radius of sphere ( m ).
204  C     f0        - Reference coriolis parameter ( 1/s )  C     f0        - Reference coriolis parameter ( 1/s )
205  C                 ( Southern edge f for beta plane )  C                 ( Southern edge f for beta plane )
206  C     beta      - df/dy ( s^-1.m^-1 )  C     beta      - df/dy ( s^-1.m^-1 )
# Line 201  C     viscAh    - Eddy viscosity coeff. Line 209  C     viscAh    - Eddy viscosity coeff.
209  C                 momentum laterally ( m^2/s )  C                 momentum laterally ( m^2/s )
210  C     viscAz    - Eddy viscosity coeff. for mixing of  C     viscAz    - Eddy viscosity coeff. for mixing of
211  C                 momentum vertically ( m^2/s )  C                 momentum vertically ( m^2/s )
212    C     viscAp    - Eddy viscosity coeff. for mixing of
213    C                 momentum vertically ( Pa^2/s )
214    C     viscAr    - Eddy viscosity coeff. for mixing of
215    C                 momentum vertically ( units of r^2/s )
216  C     viscA4    - Biharmonic viscosity coeff. for mixing of  C     viscA4    - Biharmonic viscosity coeff. for mixing of
217  C                 momentum laterally ( m^4/s )  C                 momentum laterally ( m^4/s )
218  C     diffKhT   - Laplacian diffusion coeff. for mixing of  C     diffKhT   - Laplacian diffusion coeff. for mixing of
219  C                 heat laterally ( m^2/s )  C                 heat laterally ( m^2/s )
220  C     diffKzT   - Laplacian diffusion coeff. for mixing of  C     diffKzT   - Laplacian diffusion coeff. for mixing of
221  C                 heat vertically ( m^2/s )  C                 heat vertically ( m^2/s )
222    C     diffKpT   - Laplacian diffusion coeff. for mixing of
223    C                 heat vertically ( Pa^2/s )
224    C     diffKrT   - Laplacian diffusion coeff. for mixing of
225    C                 heat vertically ( units of r^2/s )
226  C     diffK4T   - Biharmonic diffusion coeff. for mixing of  C     diffK4T   - Biharmonic diffusion coeff. for mixing of
227  C                 heat laterally ( m^4/s )  C                 heat laterally ( m^4/s )
228  C     diffKhS  -  Laplacian diffusion coeff. for mixing of  C     diffKhS  -  Laplacian diffusion coeff. for mixing of
229  C                 salt laterally ( m^2/s )  C                 salt laterally ( m^2/s )
230  C     diffKzS   - Laplacian diffusion coeff. for mixing of  C     diffKzS   - Laplacian diffusion coeff. for mixing of
231  C                 salt vertically ( m^2/s )  C                 salt vertically ( m^2/s )
232    C     diffKpS   - Laplacian diffusion coeff. for mixing of
233    C                 salt vertically ( Pa^2/s )
234    C     diffKrS   - Laplacian diffusion coeff. for mixing of
235    C                 salt vertically ( units of r^2/s )
236  C     diffK4S   - Biharmonic diffusion coeff. for mixing of  C     diffK4S   - Biharmonic diffusion coeff. for mixing of
237  C                 salt laterally ( m^4/s )  C                 salt laterally ( m^4/s )
238  C     deltaT    - Default timestep ( s )  C     deltaT    - Default timestep ( s )
# Line 227  C     deltaTtracer - Timestep for tracer Line 247  C     deltaTtracer - Timestep for tracer
247  C     freesurfFac  - Parameter to turn implicit free surface term on or off  C     freesurfFac  - Parameter to turn implicit free surface term on or off
248  C                    freesurfac = 1. uses implicit free surface  C                    freesurfac = 1. uses implicit free surface
249  C                    freesurfac = 0. uses rigid lid  C                    freesurfac = 0. uses rigid lid
250  C     tauCD     - CD scheme coupling timescale ( 1/s )  C     hFacMin      - Minimum fraction size of a cell (affects hFacC etc...)
251  C     rCD       - CD scheme normalised coupling parameter ( 0-1 )  C     hFacMinDz    - Minimum dimesional size of a cell (affects hFacC etc..., m)
252  C     GMmaxslope  - max. slope allowed in GM/Redi tensor  C     hFacMinDr    - Minimum dimesional size of a cell (affects hFacC etc..., units of r)
253  C     GMlength  - Length to use in Visbeck et al. formula for K (m)  C     tauCD        - CD scheme coupling timescale ( 1/s )
254  C     GMalpha   - alpha to use in Visbeck et al. formula for K  C     rCD          - CD scheme normalised coupling parameter ( 0-1 )
255  C     GMdepth   - Depth over which to integrate Richardson # (Visbeck et al.)  C     GMmaxslope    - max. slope allowed in GM/Redi tensor
256  C     GMkbackground - background value of GM/Redi coefficient  C     GMlength      - Length to use in Visbeck et al. formula for K (m)
257  C     startTime - Starting time for this integration ( s ).  C     GMalpha       - alpha to use in Visbeck et al. formula for K
258  C     endTime   - Ending time for this integration ( s ).  C     GMdepth       - Depth over which to integrate Richardson # (Visbeck et al.)
259  C     chkPtFreq  - Frequency of rolling check pointing ( s ).  C     GMkbackground - background value of GM/Redi diffusion coefficient ( m^2/s )
260  C     pChkPtFreq - Frequency of permanent check pointing ( s ).  C     GMmaxval      - max. value of KapGM allowed in GM/Redi scheme ( m^2/s )
261  C     dumpFreq  - Frequency with which model state is written to  C     startTime     - Starting time for this integration ( s ).
262  C                 post-processing files ( s ).  C     endTime       - Ending time for this integration ( s ).
263  C     afFacMom  - Advection of momentum term scaling parameter  C     chkPtFreq     - Frequency of rolling check pointing ( s ).
264  C     vfFacMom  - Momentum viscosity scaling parameter  C     pChkPtFreq    - Frequency of permanent check pointing ( s ).
265  C     pfFacMom  - Momentum pressure forcing parameter  C     dumpFreq      - Frequency with which model state is written to
266  C     cfFacMom  - Coriolis term scaling parameter  C                     post-processing files ( s ).
267  C     foFacMom  - Momentum forcing scaling parameter  C     afFacMom      - Advection of momentum term tracer parameter
268  C     mtFacMom  - Metric terms scaling parameter  C     vfFacMom      - Momentum viscosity tracer parameter
269  C     cAdjFreq  - Frequency of convective adjustment  C     pfFacMom      - Momentum pressure forcing tracer parameter
270    C     cfFacMom      - Coriolis term tracer parameter
271    C     foFacMom      - Momentum forcing tracer parameter
272    C     mtFacMom      - Metric terms tracer parameter
273    C     cAdjFreq      - Frequency of convective adjustment
274    C
275    C     taveFreq      - Frequency with which time-averaged model state is written to
276    C                     post-processing files ( s ).
277  C     tauThetaClimRelax - Relaxation to climatology time scale ( s ).  C     tauThetaClimRelax - Relaxation to climatology time scale ( s ).
278  C     lambdaThetaClimRelax - Inverse time scale for relaxation ( 1/s ).  C     lambdaThetaClimRelax - Inverse time scale for relaxation ( 1/s ).
279  C     tauSaltClimRelax - Relaxation to climatology time scale ( s ).  C     tauSaltClimRelax - Relaxation to climatology time scale ( s ).
280  C     lambdaSaltClimRelax - Inverse time scale for relaxation ( 1/s ).  C     lambdaSaltClimRelax - Inverse time scale for relaxation ( 1/s ).
281        COMMON /PARM_R/ cg2dTargetResidual, cg2dpcOffDFac, delZ, delX, delY,  C     externForcingPeriod - Is the period of which forcing varies (eg. 1 month)
282    C     externForcingCycle - Is the repeat time of the forcing (eg. 1 year)
283    C                          (note: externForcingCycle must be an integer
284    C                           number times externForcingPeriod)
285          COMMON /PARM_R/ cg2dTargetResidual, cg2dpcOffDFac,
286         & delP, delZ, delR, delX, delY,
287       & deltaT,deltaTmom, deltaTtracer, deltaTClock,abeps, startTime, phiMin,       & deltaT,deltaTmom, deltaTtracer, deltaTClock,abeps, startTime, phiMin,
288       & thetaMin, rSphere, rRSphere, f0, fCori, beta, viscAh, viscAz, viscA4,       & thetaMin, rSphere, recip_RSphere, f0, fCori, beta,
289       & diffKhT, diffKzT, diffK4T, diffKhS, diffKzS, diffK4S, delT,       & viscAh,  viscAz,  viscA4,  viscAr,
290       & tauCD, rCD, freeSurfFac,       & diffKhT, diffKzT, diffK4T, diffKrT,
291       & GMmaxslope,GMlength,GMalpha,GMdepth,GMkbackground,       & diffKhS, diffKzS, diffK4S, diffKrS,
292       & gravity, gBaro, rhonil, tRef, sRef,       & delT, tauCD, rCD, freeSurfFac, hFacMin, hFacMinDz,
293       & endTime, chkPtFreq, pchkPtFreq, dumpFreq,       & GMmaxslope,GMlength,GMalpha,GMdepth,GMkbackground,GMmaxval,
294         & gravity, gBaro, rhonil, recip_rhonil,
295         & recip_rhoConst, rhoConst, tRef, sRef,
296         & endTime, chkPtFreq, pchkPtFreq, dumpFreq, taveFreq,
297       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
298       & cAdjFreq, omega, tauThetaClimRelax, lambdaThetaClimRelax,       & cAdjFreq, omega, tauThetaClimRelax, lambdaThetaClimRelax,
299       & tauSaltClimRelax, lambdaSaltClimRelax       & tauSaltClimRelax, lambdaSaltClimRelax,
300         & externForcingCycle, externForcingPeriod,
301         & viscAp, diffKpT, diffKpS, hFacMinDr,
302         & theta_S, specVol_S
303    
304        _RL cg2dTargetResidual        _RL cg2dTargetResidual
305        _RL cg2dpcOffDFac        _RL cg2dpcOffDFac
306        _RL delZ(Nz)        _RL delZ(Nr)
307          _RL delP(Nr)
308          _RL delR(Nr)
309        _RL delX(Nx)        _RL delX(Nx)
310        _RL delY(Ny)        _RL delY(Ny)
311        _RL deltaT        _RL deltaT
# Line 275  C     lambdaSaltClimRelax - Inverse time Line 316  C     lambdaSaltClimRelax - Inverse time
316        _RL phiMin        _RL phiMin
317        _RL thetaMin        _RL thetaMin
318        _RL rSphere        _RL rSphere
319        _RL rRSphere        _RL recip_RSphere
320        _RL f0        _RL f0
321        _RL freeSurfFac        _RL freeSurfFac
322          _RL hFacMin
323          _RL hFacMinDz
324          _RL hFacMinDr
325        _RL beta        _RL beta
326        _RL viscAh        _RL viscAh
327        _RL viscAz        _RL viscAz
328          _RL viscAp
329          _RL viscAr
330        _RL viscA4        _RL viscA4
331        _RL diffKhT        _RL diffKhT
332          _RL diffKrT
333        _RL diffKzT        _RL diffKzT
334          _RL diffKpT
335        _RL diffK4T        _RL diffK4T
336        _RL diffKhS        _RL diffKhS
337          _RL diffKrS
338        _RL diffKzS        _RL diffKzS
339          _RL diffKpS
340        _RL diffK4S        _RL diffK4S
341        _RL delt        _RL delt
342        _RL tauCD        _RL tauCD
# Line 296  C     lambdaSaltClimRelax - Inverse time Line 346  C     lambdaSaltClimRelax - Inverse time
346        _RL GMalpha        _RL GMalpha
347        _RL GMdepth        _RL GMdepth
348        _RL GMkbackground        _RL GMkbackground
349          _RL GMmaxval
350        _RL gravity        _RL gravity
351        _RL gBaro        _RL gBaro
352        _RL rhonil        _RL rhonil
353        _RL tRef(Nz)        _RL recip_rhonil
354        _RL sRef(Nz)        _RL rhoConst
355          _RL recip_rhoConst
356          _RL specVol_S(Nr)
357          _RL tRef(Nr)
358          _RL theta_S(Nr)
359          _RL sRef(Nr)
360        _RS Fcori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS Fcori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
361        _RL startTime        _RL startTime
362        _RL endTime        _RL endTime
363        _RL chkPtFreq        _RL chkPtFreq
364        _RL pChkPtFreq        _RL pChkPtFreq
365        _RL dumpFreq        _RL dumpFreq
366          _RL taveFreq
367        _RL afFacMom        _RL afFacMom
368        _RL vfFacMom        _RL vfFacMom
369        _RL pfFacMom        _RL pfFacMom
# Line 319  C     lambdaSaltClimRelax - Inverse time Line 376  C     lambdaSaltClimRelax - Inverse time
376        _RL lambdaThetaClimRelax        _RL lambdaThetaClimRelax
377        _RL tauSaltClimRelax        _RL tauSaltClimRelax
378        _RL lambdaSaltClimRelax        _RL lambdaSaltClimRelax
379          _RL externForcingCycle
380          _RL externForcingPeriod
381    
382        COMMON /PARM_A/ HeatCapacity_Cp,        COMMON /PARM_A/ HeatCapacity_Cp,
383       &                Lamba_theta       &                Lamba_theta
# Line 327  C     lambdaSaltClimRelax - Inverse time Line 386  C     lambdaSaltClimRelax - Inverse time
386    
387  C Equation of State (polynomial coeffients)  C Equation of State (polynomial coeffients)
388        COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS        COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS
389        _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)
390  C Linear equation of state  C Linear equation of state
391  C     tAlpha    - Linear EOS thermal expansion coefficient ( 1/degree ).  C     tAlpha    - Linear EOS thermal expansion coefficient ( 1/degree ).
392  C     sBeta     - Linear EOS haline contraction coefficient.  C     sBeta     - Linear EOS haline contraction coefficient.
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.27
  ViewVC Help
Powered by ViewVC 1.1.22