167 |
C and off. |
C and off. |
168 |
C tempForcing :: Flag which turns external forcing of temperature on |
C tempForcing :: Flag which turns external forcing of temperature on |
169 |
C and off. |
C and off. |
170 |
C saltDiffusion :: Flag which turns diffusion of salinit on |
C saltDiffusion :: Flag which turns diffusion of salinity on |
171 |
C and off. |
C and off. |
172 |
C saltAdvection :: Flag which turns advection of salinit on |
C saltAdvection :: Flag which turns advection of salinity on |
173 |
C and off. |
C and off. |
174 |
C saltForcing :: Flag which turns external forcing of salinit on |
C saltForcing :: Flag which turns external forcing of salinity on |
175 |
C and off. |
C and off. |
176 |
|
C useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux |
177 |
|
C as a real Fresh Water (=> changes the Sea Level) |
178 |
|
C if F, converts P+R-E to salt flux (no SL effect) |
179 |
C rigidLid :: Set to true to use rigid lid |
C rigidLid :: Set to true to use rigid lid |
180 |
C implicitFreeSurface :: Set to true to use implcit free surface |
C implicitFreeSurface :: Set to true to use implcit free surface |
181 |
C exactConserv :: Set to true to conserve exactly the total Volume |
C exactConserv :: Set to true to conserve exactly the total Volume |
188 |
C useConstantF :: Coriolis parameter set to f0 |
C useConstantF :: Coriolis parameter set to f0 |
189 |
C useBetaPlaneF :: Coriolis parameter set to f0 + beta.y |
C useBetaPlaneF :: Coriolis parameter set to f0 + beta.y |
190 |
C useSphereF :: Coriolis parameter set to 2.omega.sin(phi) |
C useSphereF :: Coriolis parameter set to 2.omega.sin(phi) |
191 |
|
C useJamartWetPoints :: Use wet-point method for Coriolis (Jamart and Ozer, 1986) |
192 |
C implicitDiffusion :: Turns implicit vertical diffusion on |
C implicitDiffusion :: Turns implicit vertical diffusion on |
193 |
C implicitViscosity :: Turns implicit vertical viscosity on |
C implicitViscosity :: Turns implicit vertical viscosity on |
194 |
|
C multiDimAdvection :: Flag that enable multi-dimension advection |
195 |
|
C forcing_In_AB :: if False, put forcing (Temp,Salt,Tracers) contribution |
196 |
|
C out off Adams-Bashforth time stepping. |
197 |
C doThetaClimRelax :: Set true if relaxation to temperature |
C doThetaClimRelax :: Set true if relaxation to temperature |
198 |
C climatology is required. |
C climatology is required. |
199 |
C doSaltClimRelax :: Set true if relaxation to salinity |
C doSaltClimRelax :: Set true if relaxation to salinity |
203 |
C coords. |
C coords. |
204 |
C usingZCoords :: Set to indicate that we are working in height |
C usingZCoords :: Set to indicate that we are working in height |
205 |
C coords. |
C coords. |
206 |
|
C setCenterDr :: set cell Center depth and put Interface at the middle |
207 |
C nonHydrostatic :: Using non-hydrostatic terms |
C nonHydrostatic :: Using non-hydrostatic terms |
208 |
C globalFiles :: Selects between "global" and "tiled" files |
C globalFiles :: Selects between "global" and "tiled" files |
209 |
C allowFreezing :: Allows water to freeze and form ice |
C allowFreezing :: Allows water to freeze and form ice |
210 |
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) |
|
211 |
COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid, |
COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid, |
212 |
& usingCurvilinearGrid, |
& usingCurvilinearGrid, |
213 |
& no_slip_sides,no_slip_bottom, |
& no_slip_sides,no_slip_bottom, |
216 |
& momPressureForcing, vectorInvariantMomentum, |
& momPressureForcing, vectorInvariantMomentum, |
217 |
& tempDiffusion, tempAdvection, tempForcing, |
& tempDiffusion, tempAdvection, tempForcing, |
218 |
& saltDiffusion, saltAdvection, saltForcing, |
& saltDiffusion, saltAdvection, saltForcing, |
219 |
|
& useRealFreshWaterFlux, |
220 |
& rigidLid, implicitFreeSurface, exactConserv, uniformLin_PhiSurf, |
& rigidLid, implicitFreeSurface, exactConserv, uniformLin_PhiSurf, |
221 |
& momStepping, tempStepping, saltStepping, tr1Stepping, |
& momStepping, tempStepping, saltStepping, tr1Stepping, |
222 |
& metricTerms, usingSphericalPolarMTerms, |
& metricTerms, usingSphericalPolarMTerms, |
223 |
& useConstantF, useBetaPlaneF, useSphereF, |
& useConstantF, useBetaPlaneF, useSphereF, |
224 |
|
& useEnergyConservingCoriolis, useJamartWetPoints, |
225 |
& implicitDiffusion, implicitViscosity, |
& implicitDiffusion, implicitViscosity, |
226 |
|
& multiDimAdvection, forcing_In_AB, |
227 |
& doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax, |
& doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax, |
228 |
& periodicExternalForcing, usingPCoords, usingZCoords, |
& periodicExternalForcing, |
229 |
|
& usingPCoords, usingZCoords, setCenterDr, |
230 |
& nonHydrostatic, globalFiles, |
& nonHydrostatic, globalFiles, |
231 |
& allowFreezing, groundAtK1, |
& allowFreezing, groundAtK1, |
232 |
& usePickupBeforeC35, debugMode, |
& usePickupBeforeC35, debugMode, |
233 |
& readPickupWithTracer, writePickupWithTracer, |
& readPickupWithTracer, writePickupWithTracer |
|
& multiDimAdvection, useEnergyConservingCoriolis, |
|
|
& useJamartWetPoints |
|
234 |
LOGICAL usingCartesianGrid |
LOGICAL usingCartesianGrid |
235 |
LOGICAL usingSphericalPolarGrid |
LOGICAL usingSphericalPolarGrid |
236 |
LOGICAL usingCurvilinearGrid |
LOGICAL usingCurvilinearGrid |
250 |
LOGICAL saltDiffusion |
LOGICAL saltDiffusion |
251 |
LOGICAL saltAdvection |
LOGICAL saltAdvection |
252 |
LOGICAL saltForcing |
LOGICAL saltForcing |
253 |
|
LOGICAL useRealFreshWaterFlux |
254 |
LOGICAL rigidLid |
LOGICAL rigidLid |
255 |
LOGICAL implicitFreeSurface |
LOGICAL implicitFreeSurface |
256 |
LOGICAL exactConserv |
LOGICAL exactConserv |
263 |
LOGICAL useConstantF |
LOGICAL useConstantF |
264 |
LOGICAL useBetaPlaneF |
LOGICAL useBetaPlaneF |
265 |
LOGICAL useSphereF |
LOGICAL useSphereF |
266 |
|
LOGICAL useEnergyConservingCoriolis |
267 |
|
LOGICAL useJamartWetPoints |
268 |
LOGICAL implicitDiffusion |
LOGICAL implicitDiffusion |
269 |
LOGICAL implicitViscosity |
LOGICAL implicitViscosity |
270 |
|
LOGICAL multiDimAdvection |
271 |
|
LOGICAL forcing_In_AB |
272 |
LOGICAL doThetaClimRelax |
LOGICAL doThetaClimRelax |
273 |
LOGICAL doSaltClimRelax |
LOGICAL doSaltClimRelax |
274 |
LOGICAL doTr1ClimRelax |
LOGICAL doTr1ClimRelax |
275 |
LOGICAL periodicExternalForcing |
LOGICAL periodicExternalForcing |
276 |
LOGICAL usingPCoords |
LOGICAL usingPCoords |
277 |
LOGICAL usingZCoords |
LOGICAL usingZCoords |
278 |
|
LOGICAL setCenterDr |
279 |
LOGICAL nonHydrostatic |
LOGICAL nonHydrostatic |
280 |
LOGICAL globalFiles |
LOGICAL globalFiles |
281 |
LOGICAL allowFreezing |
LOGICAL allowFreezing |
284 |
LOGICAL debugMode |
LOGICAL debugMode |
285 |
LOGICAL readPickupWithTracer |
LOGICAL readPickupWithTracer |
286 |
LOGICAL writePickupWithTracer |
LOGICAL writePickupWithTracer |
|
LOGICAL multiDimAdvection |
|
|
LOGICAL useEnergyConservingCoriolis |
|
|
LOGICAL useJamartWetPoints |
|
287 |
|
|
288 |
C-- COMMON /PARM_R/ "Real" valued parameters used by the model. |
C-- COMMON /PARM_R/ "Real" valued parameters used by the model. |
289 |
C gg2dTargetResidual |
C gg2dTargetResidual |
309 |
C delP :: Vertical grid spacing ( Pa ). |
C delP :: Vertical grid spacing ( Pa ). |
310 |
C delZ :: Vertical grid spacing ( m ). |
C delZ :: Vertical grid spacing ( m ). |
311 |
C delR :: Vertical grid spacing ( units of r ). |
C delR :: Vertical grid spacing ( units of r ). |
312 |
|
C delRc :: Vertical grid spacing between cell centers (r unit). |
313 |
C delX :: Separation between cell faces (m) or (deg), depending |
C delX :: Separation between cell faces (m) or (deg), depending |
314 |
C delY on input flags. |
C delY on input flags. |
315 |
C gravity :: Accel. due to gravity ( m/s^2 ) |
C gravity :: Accel. due to gravity ( m/s^2 ) |
409 |
C externForcingCycle :: Is the repeat time of the forcing (eg. 1 year) |
C externForcingCycle :: Is the repeat time of the forcing (eg. 1 year) |
410 |
C (note: externForcingCycle must be an integer |
C (note: externForcingCycle must be an integer |
411 |
C number times externForcingPeriod) |
C number times externForcingPeriod) |
412 |
|
C convertFW2Salt :: salinity, used to convert Fresh-Water Flux to Salt Flux |
413 |
|
C (use model surface (local) value if set to -1) |
414 |
|
C temp_EvPrRn :: temperature of Rain & Evap. |
415 |
|
C salt_EvPrRn :: salinity of Rain & Evap. |
416 |
|
C trac_EvPrRn :: tracer concentration in Rain & Evap. |
417 |
|
C (notes: a) tracer content of Rain/Evap only used if both |
418 |
|
C NonLin_FrSurf & useRealFreshWater are set. |
419 |
|
C b) use model surface (local) value if set to UNSET_RL) |
420 |
C horiVertRatio :: Ratio on units in vertical to units in horizontal. |
C horiVertRatio :: Ratio on units in vertical to units in horizontal. |
421 |
C recip_horiVertRatio ( 1 if horiz in m and vertical in m ). |
C recip_horiVertRatio ( 1 if horiz in m and vertical in m ). |
422 |
C ( g*rho if horiz in m and vertical in Pa ). |
C ( g*rho if horiz in m and vertical in Pa ). |
426 |
C --"-"-- Quadratic ( linear: 1/s, quadratic: 1/m ) |
C --"-"-- Quadratic ( linear: 1/s, quadratic: 1/m ) |
427 |
COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit, |
COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit, |
428 |
& cg2dpcOffDFac, cg3dTargetResidual, |
& cg2dpcOffDFac, cg3dTargetResidual, |
429 |
& delP, delZ, delR, delX, delY, |
& delP, delZ, delR, delRc, delX, delY, |
430 |
& deltaT,deltaTmom, deltaTtracer, deltaTClock,abeps, startTime, |
& deltaT,deltaTmom, deltaTtracer, deltaTClock,abeps, startTime, |
431 |
& phiMin, thetaMin, rSphere, recip_RSphere, f0, beta, |
& phiMin, thetaMin, rSphere, recip_RSphere, f0, beta, |
432 |
& fCori, fCoriG, |
& fCori, fCoriG, |
445 |
& tauSaltClimRelax, lambdaSaltClimRelax, |
& tauSaltClimRelax, lambdaSaltClimRelax, |
446 |
& tauTr1ClimRelax, lambdaTr1ClimRelax, |
& tauTr1ClimRelax, lambdaTr1ClimRelax, |
447 |
& externForcingCycle, externForcingPeriod, |
& externForcingCycle, externForcingPeriod, |
448 |
|
& convertFW2Salt, temp_EvPrRn, salt_EvPrRn, trac_EvPrRn, |
449 |
& viscAp, diffKpT, diffKpS, hFacMinDr, hFacMinDp, |
& viscAp, diffKpT, diffKpS, hFacMinDr, hFacMinDp, |
450 |
& theta_S, specVol_S, horiVertRatio, recip_horiVertRatio, |
& horiVertRatio, recip_horiVertRatio, |
451 |
& ivdc_kappa, Ro_SeaLevel, |
& ivdc_kappa, Ro_SeaLevel, |
452 |
& bottomDragLinear,bottomDragQuadratic |
& bottomDragLinear,bottomDragQuadratic |
453 |
|
|
458 |
_RL delZ(Nr) |
_RL delZ(Nr) |
459 |
_RL delP(Nr) |
_RL delP(Nr) |
460 |
_RL delR(Nr) |
_RL delR(Nr) |
461 |
|
_RL delRc(Nr+1) |
462 |
_RL delX(Nx) |
_RL delX(Nx) |
463 |
_RL delY(Ny) |
_RL delY(Ny) |
464 |
_RL deltaT |
_RL deltaT |
508 |
_RL recip_rhonil |
_RL recip_rhonil |
509 |
_RL rhoConst |
_RL rhoConst |
510 |
_RL recip_rhoConst |
_RL recip_rhoConst |
|
_RL specVol_S(Nr) |
|
511 |
_RL tRef(Nr) |
_RL tRef(Nr) |
|
_RL theta_S(Nr) |
|
512 |
_RL sRef(Nr) |
_RL sRef(Nr) |
513 |
_RS fCori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS fCori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
514 |
_RS fCoriG(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS fCoriG(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
537 |
_RL lambdaTr1ClimRelax |
_RL lambdaTr1ClimRelax |
538 |
_RL externForcingCycle |
_RL externForcingCycle |
539 |
_RL externForcingPeriod |
_RL externForcingPeriod |
540 |
|
_RL convertFW2Salt |
541 |
|
_RL temp_EvPrRn |
542 |
|
_RL salt_EvPrRn |
543 |
|
_RL trac_EvPrRn |
544 |
_RL horiVertRatio |
_RL horiVertRatio |
545 |
_RL recip_horiVertRatio |
_RL recip_horiVertRatio |
546 |
_RL ivdc_kappa |
_RL ivdc_kappa |
586 |
LOGICAL useSHAP_FILT |
LOGICAL useSHAP_FILT |
587 |
LOGICAL useZONAL_FILT |
LOGICAL useZONAL_FILT |
588 |
LOGICAL useFLT |
LOGICAL useFLT |
589 |
|
LOGICAL usePTRACERS |
590 |
COMMON /PARM_PACKAGES/ |
COMMON /PARM_PACKAGES/ |
591 |
& useKPP, useGMRedi, useOBCS, useAIM, useECCO, |
& useKPP, useGMRedi, useOBCS, useAIM, useECCO, |
592 |
& useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT |
& useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT, |
593 |
|
& usePTRACERS |
594 |
|
|