69 |
C surfQswfile :: File containing surface shortwave radiation |
C surfQswfile :: File containing surface shortwave radiation |
70 |
C dQdTfile :: File containing thermal relaxation coefficient |
C dQdTfile :: File containing thermal relaxation coefficient |
71 |
C EmPmRfile :: File containing surface fresh water flux |
C EmPmRfile :: File containing surface fresh water flux |
72 |
|
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 |
& 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 |
99 |
CHARACTER*(MAX_LEN_FNAM) vVelInitFile |
CHARACTER*(MAX_LEN_FNAM) vVelInitFile |
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 |
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 |
173 |
C and off. |
C and off. |
174 |
C tempForcing :: Flag which turns external forcing of temperature on |
C tempForcing :: Flag which turns external forcing of temperature on |
175 |
C and off. |
C and off. |
176 |
C saltDiffusion :: Flag which turns diffusion of salinit on |
C saltDiffusion :: Flag which turns diffusion of salinity on |
177 |
C and off. |
C and off. |
178 |
C saltAdvection :: Flag which turns advection of salinit on |
C saltAdvection :: Flag which turns advection of salinity on |
179 |
C and off. |
C and off. |
180 |
C saltForcing :: Flag which turns external forcing of salinit on |
C saltForcing :: Flag which turns external forcing of salinity on |
181 |
C and off. |
C and off. |
182 |
C useRealFreshWaterFlux :: if true (=Natural BCS), treats P+R-E flux |
C useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux |
183 |
C as a real Fresh Water (=> changes the seal level) |
C as a real Fresh Water (=> changes the Sea Level) |
184 |
C if false, converts P+R-E to virtual salt flux |
C if F, converts P+R-E to salt flux (no SL effect) |
185 |
C rigidLid :: Set to true to use rigid lid |
C rigidLid :: Set to true to use rigid lid |
186 |
C implicitFreeSurface :: Set to true to use implcit free surface |
C implicitFreeSurface :: Set to true to use implcit free surface |
187 |
C exactConserv :: Set to true to conserve exactly the total Volume |
C exactConserv :: Set to true to conserve exactly the total Volume |
194 |
C useConstantF :: Coriolis parameter set to f0 |
C useConstantF :: Coriolis parameter set to f0 |
195 |
C useBetaPlaneF :: Coriolis parameter set to f0 + beta.y |
C useBetaPlaneF :: Coriolis parameter set to f0 + beta.y |
196 |
C useSphereF :: Coriolis parameter set to 2.omega.sin(phi) |
C useSphereF :: Coriolis parameter set to 2.omega.sin(phi) |
197 |
|
C useJamartWetPoints :: Use wet-point method for Coriolis (Jamart and Ozer, 1986) |
198 |
C implicitDiffusion :: Turns implicit vertical diffusion on |
C implicitDiffusion :: Turns implicit vertical diffusion on |
199 |
C implicitViscosity :: Turns implicit vertical viscosity on |
C implicitViscosity :: Turns implicit vertical viscosity on |
200 |
|
C multiDimAdvection :: Flag that enable multi-dimension advection |
201 |
|
C forcing_In_AB :: if False, put forcing (Temp,Salt,Tracers) contribution |
202 |
|
C out off Adams-Bashforth time stepping. |
203 |
C doThetaClimRelax :: Set true if relaxation to temperature |
C doThetaClimRelax :: Set true if relaxation to temperature |
204 |
C climatology is required. |
C climatology is required. |
205 |
C doSaltClimRelax :: Set true if relaxation to salinity |
C doSaltClimRelax :: Set true if relaxation to salinity |
209 |
C coords. |
C coords. |
210 |
C usingZCoords :: Set to indicate that we are working in height |
C usingZCoords :: Set to indicate that we are working in height |
211 |
C coords. |
C coords. |
212 |
|
C setCenterDr :: set cell Center depth and put Interface at the middle |
213 |
C nonHydrostatic :: Using non-hydrostatic terms |
C nonHydrostatic :: Using non-hydrostatic terms |
214 |
C globalFiles :: Selects between "global" and "tiled" files |
C globalFiles :: Selects between "global" and "tiled" files |
215 |
C allowFreezing :: Allows water to freeze and form ice |
C allowFreezing :: Allows water to freeze and form ice |
216 |
C groundAtK1 :: put the surface(k=1) at the Lower Boundary (=ground) |
C groundAtK1 :: put the surface(k=1) at the Lower Boundary (=ground) |
|
C useJamartWetPoints :: Use wet-point method for Coriolis (Jamart and Ozer, 1986) |
|
217 |
COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid, |
COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid, |
218 |
& usingCurvilinearGrid, |
& usingCurvilinearGrid, |
219 |
& no_slip_sides,no_slip_bottom, |
& no_slip_sides,no_slip_bottom, |
227 |
& momStepping, tempStepping, saltStepping, tr1Stepping, |
& momStepping, tempStepping, saltStepping, tr1Stepping, |
228 |
& metricTerms, usingSphericalPolarMTerms, |
& metricTerms, usingSphericalPolarMTerms, |
229 |
& useConstantF, useBetaPlaneF, useSphereF, |
& useConstantF, useBetaPlaneF, useSphereF, |
230 |
|
& useEnergyConservingCoriolis, useJamartWetPoints, |
231 |
& implicitDiffusion, implicitViscosity, |
& implicitDiffusion, implicitViscosity, |
232 |
|
& multiDimAdvection, forcing_In_AB, |
233 |
& doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax, |
& doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax, |
234 |
& periodicExternalForcing, usingPCoords, usingZCoords, |
& periodicExternalForcing, |
235 |
|
& usingPCoords, usingZCoords, setCenterDr, |
236 |
& nonHydrostatic, globalFiles, |
& nonHydrostatic, globalFiles, |
237 |
& allowFreezing, groundAtK1, |
& allowFreezing, groundAtK1, |
238 |
& usePickupBeforeC35, debugMode, |
& usePickupBeforeC35, debugMode, |
239 |
& readPickupWithTracer, writePickupWithTracer, |
& readPickupWithTracer, writePickupWithTracer |
|
& multiDimAdvection, useEnergyConservingCoriolis, |
|
|
& useJamartWetPoints |
|
240 |
LOGICAL usingCartesianGrid |
LOGICAL usingCartesianGrid |
241 |
LOGICAL usingSphericalPolarGrid |
LOGICAL usingSphericalPolarGrid |
242 |
LOGICAL usingCurvilinearGrid |
LOGICAL usingCurvilinearGrid |
269 |
LOGICAL useConstantF |
LOGICAL useConstantF |
270 |
LOGICAL useBetaPlaneF |
LOGICAL useBetaPlaneF |
271 |
LOGICAL useSphereF |
LOGICAL useSphereF |
272 |
|
LOGICAL useEnergyConservingCoriolis |
273 |
|
LOGICAL useJamartWetPoints |
274 |
LOGICAL implicitDiffusion |
LOGICAL implicitDiffusion |
275 |
LOGICAL implicitViscosity |
LOGICAL implicitViscosity |
276 |
|
LOGICAL multiDimAdvection |
277 |
|
LOGICAL forcing_In_AB |
278 |
LOGICAL doThetaClimRelax |
LOGICAL doThetaClimRelax |
279 |
LOGICAL doSaltClimRelax |
LOGICAL doSaltClimRelax |
280 |
LOGICAL doTr1ClimRelax |
LOGICAL doTr1ClimRelax |
281 |
LOGICAL periodicExternalForcing |
LOGICAL periodicExternalForcing |
282 |
LOGICAL usingPCoords |
LOGICAL usingPCoords |
283 |
LOGICAL usingZCoords |
LOGICAL usingZCoords |
284 |
|
LOGICAL setCenterDr |
285 |
LOGICAL nonHydrostatic |
LOGICAL nonHydrostatic |
286 |
LOGICAL globalFiles |
LOGICAL globalFiles |
287 |
LOGICAL allowFreezing |
LOGICAL allowFreezing |
290 |
LOGICAL debugMode |
LOGICAL debugMode |
291 |
LOGICAL readPickupWithTracer |
LOGICAL readPickupWithTracer |
292 |
LOGICAL writePickupWithTracer |
LOGICAL writePickupWithTracer |
|
LOGICAL multiDimAdvection |
|
|
LOGICAL useEnergyConservingCoriolis |
|
|
LOGICAL useJamartWetPoints |
|
293 |
|
|
294 |
C-- COMMON /PARM_R/ "Real" valued parameters used by the model. |
C-- COMMON /PARM_R/ "Real" valued parameters used by the model. |
295 |
C gg2dTargetResidual |
C gg2dTargetResidual |
315 |
C delP :: Vertical grid spacing ( Pa ). |
C delP :: Vertical grid spacing ( Pa ). |
316 |
C delZ :: Vertical grid spacing ( m ). |
C delZ :: Vertical grid spacing ( m ). |
317 |
C delR :: Vertical grid spacing ( units of r ). |
C delR :: Vertical grid spacing ( units of r ). |
318 |
|
C delRc :: Vertical grid spacing between cell centers (r unit). |
319 |
C delX :: Separation between cell faces (m) or (deg), depending |
C delX :: Separation between cell faces (m) or (deg), depending |
320 |
C delY on input flags. |
C delY on input flags. |
321 |
C gravity :: Accel. due to gravity ( m/s^2 ) |
C gravity :: Accel. due to gravity ( m/s^2 ) |
373 |
C are referenced to this clock. ( s ) |
C are referenced to this clock. ( s ) |
374 |
C deltaTMom :: Timestep for momemtum equations ( s ) |
C deltaTMom :: Timestep for momemtum equations ( s ) |
375 |
C deltaTtracer :: Timestep for tracer equations ( s ) |
C deltaTtracer :: Timestep for tracer equations ( s ) |
376 |
|
C deltaTfreesurf :: Timestep for free-surface equation ( s ) |
377 |
C freesurfFac :: Parameter to turn implicit free surface term on or off |
C freesurfFac :: Parameter to turn implicit free surface term on or off |
378 |
C freesurfac = 1. uses implicit free surface |
C freesurfac = 1. uses implicit free surface |
379 |
C freesurfac = 0. uses rigid lid |
C freesurfac = 0. uses rigid lid |
433 |
C --"-"-- Quadratic ( linear: 1/s, quadratic: 1/m ) |
C --"-"-- Quadratic ( linear: 1/s, quadratic: 1/m ) |
434 |
COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit, |
COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit, |
435 |
& cg2dpcOffDFac, cg3dTargetResidual, |
& cg2dpcOffDFac, cg3dTargetResidual, |
436 |
& delP, delZ, delR, delX, delY, |
& delP, delZ, delR, delRc, delX, delY, |
437 |
& deltaT,deltaTmom, deltaTtracer, deltaTClock,abeps, startTime, |
& deltaT, deltaTmom, deltaTtracer, deltaTfreesurf, deltaTClock, |
438 |
|
& abeps, startTime, |
439 |
& phiMin, thetaMin, rSphere, recip_RSphere, f0, beta, |
& phiMin, thetaMin, rSphere, recip_RSphere, f0, beta, |
440 |
& fCori, fCoriG, |
& fCori, fCoriG, |
441 |
& viscAh, viscAz, viscA4, viscAr, viscAstrain, viscAtension, |
& viscAh, viscAz, viscA4, viscAr, viscAstrain, viscAtension, |
466 |
_RL delZ(Nr) |
_RL delZ(Nr) |
467 |
_RL delP(Nr) |
_RL delP(Nr) |
468 |
_RL delR(Nr) |
_RL delR(Nr) |
469 |
|
_RL delRc(Nr+1) |
470 |
_RL delX(Nx) |
_RL delX(Nx) |
471 |
_RL delY(Ny) |
_RL delY(Ny) |
472 |
_RL deltaT |
_RL deltaT |
473 |
_RL deltaTClock |
_RL deltaTClock |
474 |
_RL deltaTmom |
_RL deltaTmom |
475 |
_RL deltaTtracer |
_RL deltaTtracer |
476 |
|
_RL deltaTfreesurf |
477 |
_RL abeps |
_RL abeps |
478 |
_RL phiMin |
_RL phiMin |
479 |
_RL thetaMin |
_RL thetaMin |
563 |
_RL Lamba_theta |
_RL Lamba_theta |
564 |
_RL recip_Cp |
_RL recip_Cp |
565 |
|
|
566 |
C Equation of State (polynomial coeffients) |
CmlC Equation of State (polynomial coeffients) |
567 |
COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS |
Cml COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS |
568 |
_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) |
569 |
C Linear equation of state |
CmlC Linear equation of state |
570 |
C tAlpha :: Linear EOS thermal expansion coefficient ( 1/degree ). |
CmlC tAlpha :: Linear EOS thermal expansion coefficient ( 1/degree ). |
571 |
C sBeta :: Linear EOS haline contraction coefficient. |
CmlC sBeta :: Linear EOS haline contraction coefficient. |
572 |
COMMON /PARM_EOS_LIN/ tAlpha,sBeta,eosType |
Cml COMMON /PARM_EOS_LIN/ tAlpha,sBeta,eosType |
573 |
_RL tAlpha |
Cml _RL tAlpha |
574 |
_RL sBeta |
Cml _RL sBeta |
|
character*(6) eosType |
|
575 |
|
|
576 |
C Atmospheric physical parameters (Ideal Gas EOS, ...) |
C Atmospheric physical parameters (Ideal Gas EOS, ...) |
577 |
C atm_po :: standard reference pressure |
C atm_po :: standard reference pressure |