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 |
C useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux |
177 |
C as a real Fresh Water (=> changes the seal level) |
C as a real Fresh Water (=> changes the Sea Level) |
178 |
C if false, converts P+R-E to virtual salt flux |
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, |
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 |
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 ) |
367 |
C are referenced to this clock. ( s ) |
C are referenced to this clock. ( s ) |
368 |
C deltaTMom :: Timestep for momemtum equations ( s ) |
C deltaTMom :: Timestep for momemtum equations ( s ) |
369 |
C deltaTtracer :: Timestep for tracer equations ( s ) |
C deltaTtracer :: Timestep for tracer equations ( s ) |
370 |
|
C deltaTfreesurf :: Timestep for free-surface equation ( s ) |
371 |
C freesurfFac :: Parameter to turn implicit free surface term on or off |
C freesurfFac :: Parameter to turn implicit free surface term on or off |
372 |
C freesurfac = 1. uses implicit free surface |
C freesurfac = 1. uses implicit free surface |
373 |
C freesurfac = 0. uses rigid lid |
C freesurfac = 0. uses rigid lid |
427 |
C --"-"-- Quadratic ( linear: 1/s, quadratic: 1/m ) |
C --"-"-- Quadratic ( linear: 1/s, quadratic: 1/m ) |
428 |
COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit, |
COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit, |
429 |
& cg2dpcOffDFac, cg3dTargetResidual, |
& cg2dpcOffDFac, cg3dTargetResidual, |
430 |
& delP, delZ, delR, delX, delY, |
& delP, delZ, delR, delRc, delX, delY, |
431 |
& deltaT,deltaTmom, deltaTtracer, deltaTClock,abeps, startTime, |
& deltaT, deltaTmom, deltaTtracer, deltaTfreesurf, deltaTClock, |
432 |
|
& abeps, startTime, |
433 |
& phiMin, thetaMin, rSphere, recip_RSphere, f0, beta, |
& phiMin, thetaMin, rSphere, recip_RSphere, f0, beta, |
434 |
& fCori, fCoriG, |
& fCori, fCoriG, |
435 |
& viscAh, viscAz, viscA4, viscAr, viscAstrain, viscAtension, |
& viscAh, viscAz, viscA4, viscAr, viscAstrain, viscAtension, |
460 |
_RL delZ(Nr) |
_RL delZ(Nr) |
461 |
_RL delP(Nr) |
_RL delP(Nr) |
462 |
_RL delR(Nr) |
_RL delR(Nr) |
463 |
|
_RL delRc(Nr+1) |
464 |
_RL delX(Nx) |
_RL delX(Nx) |
465 |
_RL delY(Ny) |
_RL delY(Ny) |
466 |
_RL deltaT |
_RL deltaT |
467 |
_RL deltaTClock |
_RL deltaTClock |
468 |
_RL deltaTmom |
_RL deltaTmom |
469 |
_RL deltaTtracer |
_RL deltaTtracer |
470 |
|
_RL deltaTfreesurf |
471 |
_RL abeps |
_RL abeps |
472 |
_RL phiMin |
_RL phiMin |
473 |
_RL thetaMin |
_RL thetaMin |