72 |
C pLoadFile :: File containing pressure loading |
C pLoadFile :: File containing pressure loading |
73 |
C buoyancyRelation :: Flag used to indicate which relation to use to |
C buoyancyRelation :: Flag used to indicate which relation to use to |
74 |
C get buoyancy. |
C get buoyancy. |
75 |
|
C eosType :: choose the equation of state: |
76 |
|
C LINEAR, POLY3, UNESCO, JMD95Z, JMD95P |
77 |
COMMON /PARM_C/ checkPtSuff, |
COMMON /PARM_C/ checkPtSuff, |
78 |
& bathyFile, topoFile, |
& bathyFile, topoFile, |
79 |
& hydrogThetaFile, hydrogSaltFile, |
& hydrogThetaFile, hydrogSaltFile, |
81 |
& saltClimFile, buoyancyRelation, |
& saltClimFile, buoyancyRelation, |
82 |
& EmPmRfile, surfQfile, surfQswfile, |
& EmPmRfile, surfQfile, surfQswfile, |
83 |
& uVelInitFile, vVelInitFile, pSurfInitFile, |
& uVelInitFile, vVelInitFile, pSurfInitFile, |
84 |
& dQdTfile, ploadFile |
& dQdTfile, ploadFile, |
85 |
|
& eosType |
86 |
CHARACTER*(5) checkPtSuff(maxNoChkptLev) |
CHARACTER*(5) checkPtSuff(maxNoChkptLev) |
87 |
CHARACTER*(MAX_LEN_FNAM) bathyFile, topoFile |
CHARACTER*(MAX_LEN_FNAM) bathyFile, topoFile |
88 |
CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile |
CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile |
100 |
CHARACTER*(MAX_LEN_FNAM) pSurfInitFile |
CHARACTER*(MAX_LEN_FNAM) pSurfInitFile |
101 |
CHARACTER*(MAX_LEN_FNAM) dQdTfile |
CHARACTER*(MAX_LEN_FNAM) dQdTfile |
102 |
CHARACTER*(MAX_LEN_FNAM) ploadFile |
CHARACTER*(MAX_LEN_FNAM) ploadFile |
103 |
|
CHARACTER*(6) eosType |
104 |
|
|
105 |
C-- COMMON /PARM_I/ Integer valued parameters used by the model. |
C-- COMMON /PARM_I/ Integer valued parameters used by the model. |
106 |
C cg2dMaxIters :: Maximum number of iterations in the |
C cg2dMaxIters :: Maximum number of iterations in the |
166 |
C on and off. |
C on and off. |
167 |
C metricTerms :: Flag which turns metric terms on or off. |
C metricTerms :: Flag which turns metric terms on or off. |
168 |
C usingSphericalPolarMTerms :: If TRUE use spherical polar metric terms. |
C usingSphericalPolarMTerms :: If TRUE use spherical polar metric terms. |
169 |
|
C useNHMTerms :: If TRUE use non-hydrostatic metric terms. |
170 |
C useCoriolis :: Flag which turns the coriolis terms on and off. |
C useCoriolis :: Flag which turns the coriolis terms on and off. |
171 |
C tempDiffusion :: Flag which turns diffusion of temperature on |
C tempDiffusion :: Flag which turns diffusion of temperature on |
172 |
C and off. |
C and off. |
212 |
C coords. |
C coords. |
213 |
C setCenterDr :: set cell Center depth and put Interface at the middle |
C setCenterDr :: set cell Center depth and put Interface at the middle |
214 |
C nonHydrostatic :: Using non-hydrostatic terms |
C nonHydrostatic :: Using non-hydrostatic terms |
215 |
|
C quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm |
216 |
C globalFiles :: Selects between "global" and "tiled" files |
C globalFiles :: Selects between "global" and "tiled" files |
217 |
C allowFreezing :: Allows water to freeze and form ice |
C allowFreezing :: Allows water to freeze and form ice |
218 |
C groundAtK1 :: put the surface(k=1) at the Lower Boundary (=ground) |
C groundAtK1 :: put the surface(k=1) at the Lower Boundary (=ground) |
227 |
& useRealFreshWaterFlux, |
& useRealFreshWaterFlux, |
228 |
& rigidLid, implicitFreeSurface, exactConserv, uniformLin_PhiSurf, |
& rigidLid, implicitFreeSurface, exactConserv, uniformLin_PhiSurf, |
229 |
& momStepping, tempStepping, saltStepping, tr1Stepping, |
& momStepping, tempStepping, saltStepping, tr1Stepping, |
230 |
& metricTerms, usingSphericalPolarMTerms, |
& metricTerms, usingSphericalPolarMTerms, useNHMTerms, |
231 |
& useConstantF, useBetaPlaneF, useSphereF, |
& useConstantF, useBetaPlaneF, useSphereF, |
232 |
& useEnergyConservingCoriolis, useJamartWetPoints, |
& useEnergyConservingCoriolis, useJamartWetPoints, |
233 |
& implicitDiffusion, implicitViscosity, |
& implicitDiffusion, implicitViscosity, |
235 |
& doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax, |
& doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax, |
236 |
& periodicExternalForcing, |
& periodicExternalForcing, |
237 |
& usingPCoords, usingZCoords, setCenterDr, |
& usingPCoords, usingZCoords, setCenterDr, |
238 |
& nonHydrostatic, globalFiles, |
& nonHydrostatic, quasiHydrostatic, globalFiles, |
239 |
& allowFreezing, groundAtK1, |
& allowFreezing, groundAtK1, |
240 |
& usePickupBeforeC35, debugMode, |
& usePickupBeforeC35, debugMode, |
241 |
& readPickupWithTracer, writePickupWithTracer |
& readPickupWithTracer, writePickupWithTracer |
243 |
LOGICAL usingSphericalPolarGrid |
LOGICAL usingSphericalPolarGrid |
244 |
LOGICAL usingCurvilinearGrid |
LOGICAL usingCurvilinearGrid |
245 |
LOGICAL usingSphericalPolarMTerms |
LOGICAL usingSphericalPolarMTerms |
246 |
|
LOGICAL useNHMTerms |
247 |
LOGICAL no_slip_sides |
LOGICAL no_slip_sides |
248 |
LOGICAL no_slip_bottom |
LOGICAL no_slip_bottom |
249 |
LOGICAL staggerTimeStep |
LOGICAL staggerTimeStep |
286 |
LOGICAL usingZCoords |
LOGICAL usingZCoords |
287 |
LOGICAL setCenterDr |
LOGICAL setCenterDr |
288 |
LOGICAL nonHydrostatic |
LOGICAL nonHydrostatic |
289 |
|
LOGICAL quasiHydrostatic |
290 |
LOGICAL globalFiles |
LOGICAL globalFiles |
291 |
LOGICAL allowFreezing |
LOGICAL allowFreezing |
292 |
LOGICAL groundAtK1 |
LOGICAL groundAtK1 |
325 |
C gravity :: Accel. due to gravity ( m/s^2 ) |
C gravity :: Accel. due to gravity ( m/s^2 ) |
326 |
C recip_gravity and its inverse |
C recip_gravity and its inverse |
327 |
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 ) |
328 |
C ronil :: Reference density |
C rhoNil :: Reference density for the linear equation of state |
329 |
C rhoConst :: Vertically constant reference density |
C rhoConst :: Vertically constant reference density |
330 |
|
C rhoConstFresh :: Constant reference density for fresh water (rain) |
331 |
C startTime :: Start time for model ( s ) |
C startTime :: Start time for model ( s ) |
332 |
C phiMin :: Latitude of southern most cell face. |
C phiMin :: Latitude of southern most cell face. |
333 |
C thetaMin :: Longitude of western most cell face (this |
C thetaMin :: Longitude of western most cell face (this |
411 |
C cosPower :: Power of cosine of latitude to multiply viscosity |
C cosPower :: Power of cosine of latitude to multiply viscosity |
412 |
C cAdjFreq :: Frequency of convective adjustment |
C cAdjFreq :: Frequency of convective adjustment |
413 |
C |
C |
414 |
C taveFreq :: Frequency with which time-averaged model state is written to |
C taveFreq :: Frequency with which time-averaged model state |
415 |
C post-processing files ( s ). |
C is written to post-processing files ( s ). |
416 |
|
C tave_lastIter :: (for state variable only) fraction of the last time |
417 |
|
C step (of each taveFreq period) put in the time average. |
418 |
|
C (fraction for 1rst iter = 1 - tave_lastIter) |
419 |
C tauThetaClimRelax :: Relaxation to climatology time scale ( s ). |
C tauThetaClimRelax :: Relaxation to climatology time scale ( s ). |
420 |
C lambdaThetaClimRelax :: Inverse time scale for relaxation ( 1/s ). |
C lambdaThetaClimRelax :: Inverse time scale for relaxation ( 1/s ). |
421 |
C tauSaltClimRelax :: Relaxation to climatology time scale ( s ). |
C tauSaltClimRelax :: Relaxation to climatology time scale ( s ). |
445 |
& deltaT, deltaTmom, deltaTtracer, deltaTfreesurf, deltaTClock, |
& deltaT, deltaTmom, deltaTtracer, deltaTfreesurf, deltaTClock, |
446 |
& abeps, startTime, |
& abeps, startTime, |
447 |
& phiMin, thetaMin, rSphere, recip_RSphere, f0, beta, |
& phiMin, thetaMin, rSphere, recip_RSphere, f0, beta, |
448 |
& fCori, fCoriG, |
& fCori, fCoriG, fCoriCos, |
449 |
& viscAh, viscAz, viscA4, viscAr, viscAstrain, viscAtension, |
& viscAh, viscAz, viscA4, viscAr, viscAstrain, viscAtension, |
450 |
& diffKhT, diffKzT, diffK4T, diffKrT, |
& diffKhT, diffKzT, diffK4T, diffKrT, |
451 |
& diffKhS, diffKzS, diffK4S, diffKrS, |
& diffKhS, diffKzS, diffK4S, diffKrS, |
452 |
& delT, tauCD, rCD, freeSurfFac, implicSurfPress, implicDiv2Dflow, |
& delT, tauCD, rCD, freeSurfFac, implicSurfPress, implicDiv2Dflow, |
453 |
& hFacMin, hFacMinDz, hFacInf, hFacSup, |
& hFacMin, hFacMinDz, hFacInf, hFacSup, |
454 |
& gravity, recip_Gravity, gBaro, rhonil, recip_rhonil, |
& gravity, recip_Gravity, gBaro, rhonil, recip_rhonil, |
455 |
& recip_rhoConst, rhoConst, tRef, sRef, |
& recip_rhoConst, rhoConst, |
456 |
|
& rhoConstFresh, convertEmP2rUnit, tRef, sRef, |
457 |
& endTime, chkPtFreq, pchkPtFreq, dumpFreq, |
& endTime, chkPtFreq, pchkPtFreq, dumpFreq, |
458 |
& diagFreq, taveFreq, monitorFreq, |
& diagFreq, taveFreq, tave_lastIter, monitorFreq, |
459 |
& afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom, |
& afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom, |
460 |
& cosPower, cAdjFreq, omega, |
& cosPower, cAdjFreq, omega, |
461 |
& tauThetaClimRelax, lambdaThetaClimRelax, |
& tauThetaClimRelax, lambdaThetaClimRelax, |
526 |
_RL recip_rhonil |
_RL recip_rhonil |
527 |
_RL rhoConst |
_RL rhoConst |
528 |
_RL recip_rhoConst |
_RL recip_rhoConst |
529 |
|
_RL rhoConstFresh |
530 |
|
_RL convertEmP2rUnit |
531 |
_RL tRef(Nr) |
_RL tRef(Nr) |
532 |
_RL sRef(Nr) |
_RL sRef(Nr) |
533 |
_RS fCori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS fCori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
534 |
_RS fCoriG(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS fCoriG(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
535 |
|
_RS fCoriCos(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
536 |
_RL startTime |
_RL startTime |
537 |
_RL endTime |
_RL endTime |
538 |
_RL chkPtFreq |
_RL chkPtFreq |
540 |
_RL dumpFreq |
_RL dumpFreq |
541 |
_RL diagFreq |
_RL diagFreq |
542 |
_RL taveFreq |
_RL taveFreq |
543 |
|
_RL tave_lastIter |
544 |
_RL monitorFreq |
_RL monitorFreq |
545 |
_RL afFacMom |
_RL afFacMom |
546 |
_RL vfFacMom |
_RL vfFacMom |
576 |
_RL Lamba_theta |
_RL Lamba_theta |
577 |
_RL recip_Cp |
_RL recip_Cp |
578 |
|
|
579 |
C Equation of State (polynomial coeffients) |
CmlC Equation of State (polynomial coeffients) |
580 |
COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS |
Cml COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS |
581 |
_RL eosC(9,Nr+1),eosSig0(Nr+1),eosRefT(Nr+1),eosRefS(Nr+1) |
Cml _RL eosC(9,Nr+1),eosSig0(Nr+1),eosRefT(Nr+1),eosRefS(Nr+1) |
582 |
C Linear equation of state |
CmlC Linear equation of state |
583 |
C tAlpha :: Linear EOS thermal expansion coefficient ( 1/degree ). |
CmlC tAlpha :: Linear EOS thermal expansion coefficient ( 1/degree ). |
584 |
C sBeta :: Linear EOS haline contraction coefficient. |
CmlC sBeta :: Linear EOS haline contraction coefficient. |
585 |
COMMON /PARM_EOS_LIN/ tAlpha,sBeta,eosType |
Cml COMMON /PARM_EOS_LIN/ tAlpha,sBeta,eosType |
586 |
_RL tAlpha |
Cml _RL tAlpha |
587 |
_RL sBeta |
Cml _RL sBeta |
|
character*(6) eosType |
|
588 |
|
|
589 |
C Atmospheric physical parameters (Ideal Gas EOS, ...) |
C Atmospheric physical parameters (Ideal Gas EOS, ...) |
590 |
C atm_po :: standard reference pressure |
C atm_po :: standard reference pressure |
608 |
LOGICAL useZONAL_FILT |
LOGICAL useZONAL_FILT |
609 |
LOGICAL useFLT |
LOGICAL useFLT |
610 |
LOGICAL usePTRACERS |
LOGICAL usePTRACERS |
611 |
|
LOGICAL useSEAICE |
612 |
|
cswdblk -- add --- |
613 |
|
LOGICAL useBulkf |
614 |
|
LOGICAL useTSEAICE |
615 |
|
cswdblk --- end add --- |
616 |
COMMON /PARM_PACKAGES/ |
COMMON /PARM_PACKAGES/ |
617 |
& useKPP, useGMRedi, useOBCS, useAIM, useECCO, |
& useKPP, useGMRedi, useOBCS, useAIM, useECCO, |
618 |
& useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT, |
& useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT, |
619 |
& usePTRACERS |
& usePTRACERS, useSEAICE, |
620 |
|
cswdblk -- add --- |
621 |
|
& useTSEAICE, useBulkf |
622 |
|
cswdblk --- end add --- |