173 |
C saltAdvScheme :: Salt. Horiz.advection scheme selector |
C saltAdvScheme :: Salt. Horiz.advection scheme selector |
174 |
C saltVertAdvScheme :: Salt. Vert. Advection scheme selector |
C saltVertAdvScheme :: Salt. Vert. Advection scheme selector |
175 |
C selectKEscheme :: Kinetic Energy scheme selector (Vector Inv.) |
C selectKEscheme :: Kinetic Energy scheme selector (Vector Inv.) |
176 |
|
C selectVortScheme :: Scheme selector for Vorticity term (Vector Inv.) |
177 |
C monitorSelect :: select group of variables to monitor |
C monitorSelect :: select group of variables to monitor |
178 |
C =1 : dynvars ; =2 : + vort ; =3 : + surface |
C =1 : dynvars ; =2 : + vort ; =3 : + surface |
179 |
C debugLevel :: debug level selector: higher -> more writing |
C debugLevel :: debug level selector: higher -> more writing |
190 |
& momForcingOutAB, tracForcingOutAB, |
& momForcingOutAB, tracForcingOutAB, |
191 |
& tempAdvScheme, tempVertAdvScheme, |
& tempAdvScheme, tempVertAdvScheme, |
192 |
& saltAdvScheme, saltVertAdvScheme, |
& saltAdvScheme, saltVertAdvScheme, |
193 |
& selectKEscheme, |
& selectKEscheme, selectVortScheme, |
194 |
& monitorSelect, debugLevel |
& monitorSelect, debugLevel |
195 |
INTEGER cg2dMaxIters |
INTEGER cg2dMaxIters |
196 |
INTEGER cg2dChkResFreq |
INTEGER cg2dChkResFreq |
209 |
INTEGER tempAdvScheme, tempVertAdvScheme |
INTEGER tempAdvScheme, tempVertAdvScheme |
210 |
INTEGER saltAdvScheme, saltVertAdvScheme |
INTEGER saltAdvScheme, saltVertAdvScheme |
211 |
INTEGER selectKEscheme |
INTEGER selectKEscheme |
212 |
|
INTEGER selectVortScheme |
213 |
INTEGER monitorSelect |
INTEGER monitorSelect |
214 |
INTEGER debugLevel |
INTEGER debugLevel |
215 |
|
|
276 |
C useCDscheme :: use CD-scheme to calculate Coriolis terms. |
C useCDscheme :: use CD-scheme to calculate Coriolis terms. |
277 |
C useJamartWetPoints :: Use wet-point method for Coriolis (Jamart and Ozer, 1986) |
C useJamartWetPoints :: Use wet-point method for Coriolis (Jamart and Ozer, 1986) |
278 |
C useJamartMomAdv :: Use wet-point method for V.I. non-linear term |
C useJamartMomAdv :: Use wet-point method for V.I. non-linear term |
|
C SadournyCoriolis :: use the enstrophy conserving scheme by Sadourny |
|
279 |
C upwindVorticity :: bias interpolation of vorticity in the Coriolis term |
C upwindVorticity :: bias interpolation of vorticity in the Coriolis term |
280 |
C highOrderVorticity :: use 3rd/4th order interp. of vorticity (V.I., advection) |
C highOrderVorticity :: use 3rd/4th order interp. of vorticity (V.I., advection) |
281 |
C upwindShear :: use 1rst order upwind interp. (V.I., vertical advection) |
C upwindShear :: use 1rst order upwind interp. (V.I., vertical advection) |
333 |
C balanceEmPmR :: substract global mean of EmPmR at every time step |
C balanceEmPmR :: substract global mean of EmPmR at every time step |
334 |
C balanceQnet :: substract global mean of Qnet at every time step |
C balanceQnet :: substract global mean of Qnet at every time step |
335 |
C balancePrintMean:: print substracted global means to STDOUT |
C balancePrintMean:: print substracted global means to STDOUT |
336 |
C rotateGrid :: rotate grid coordinates to geographical coordinates |
C rotateGrid :: rotate grid coordinates to geographical coordinates |
337 |
C according to Euler angles phiEuler, thetaEuler, psiEuler |
C according to Euler angles phiEuler, thetaEuler, psiEuler |
338 |
|
|
339 |
COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid, |
COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid, |
355 |
& useConstantF, useBetaPlaneF, useSphereF, use3dCoriolis, |
& useConstantF, useBetaPlaneF, useSphereF, use3dCoriolis, |
356 |
& useCDscheme, |
& useCDscheme, |
357 |
& useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv, |
& useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv, |
358 |
& SadournyCoriolis, upwindVorticity, highOrderVorticity, |
& upwindVorticity, highOrderVorticity, |
359 |
& useAbsVorticity, upwindShear, |
& useAbsVorticity, upwindShear, |
360 |
& implicitDiffusion, implicitViscosity, |
& implicitDiffusion, implicitViscosity, |
361 |
& tempImplVertAdv, saltImplVertAdv, momImplVertAdv, |
& tempImplVertAdv, saltImplVertAdv, momImplVertAdv, |
420 |
LOGICAL useEnergyConservingCoriolis |
LOGICAL useEnergyConservingCoriolis |
421 |
LOGICAL useJamartWetPoints |
LOGICAL useJamartWetPoints |
422 |
LOGICAL useJamartMomAdv |
LOGICAL useJamartMomAdv |
|
LOGICAL SadournyCoriolis |
|
423 |
LOGICAL upwindVorticity |
LOGICAL upwindVorticity |
424 |
LOGICAL highOrderVorticity |
LOGICAL highOrderVorticity |
425 |
LOGICAL useAbsVorticity |
LOGICAL useAbsVorticity |
529 |
C beta :: df/dy ( s^-1.m^-1 ) |
C beta :: df/dy ( s^-1.m^-1 ) |
530 |
C omega :: Angular velocity ( rad/s ) |
C omega :: Angular velocity ( rad/s ) |
531 |
C rotationPeriod :: Rotation period (s) (= 2.pi/omega) |
C rotationPeriod :: Rotation period (s) (= 2.pi/omega) |
532 |
|
C viscAr :: Eddy viscosity coeff. for mixing of |
533 |
|
C momentum vertically ( units of r^2/s ) |
534 |
C viscAh :: Eddy viscosity coeff. for mixing of |
C viscAh :: Eddy viscosity coeff. for mixing of |
535 |
C momentum laterally ( m^2/s ) |
C momentum laterally ( m^2/s ) |
536 |
C viscAhW :: Eddy viscosity coeff. for mixing of vertical |
C viscAhW :: Eddy viscosity coeff. for mixing of vertical |
537 |
C momentum laterally, no effect for hydrostatic |
C momentum laterally, no effect for hydrostatic |
538 |
C model, defaults to viscAh if unset ( m^2/s ) |
C model, defaults to viscAh if unset ( m^2/s ) |
539 |
C Not used if variable horiz. viscosity is used. |
C Not used if variable horiz. viscosity is used. |
|
C viscAr :: Eddy viscosity coeff. for mixing of |
|
|
C momentum vertically ( units of r^2/s ) |
|
540 |
C viscA4 :: Biharmonic viscosity coeff. for mixing of |
C viscA4 :: Biharmonic viscosity coeff. for mixing of |
541 |
C momentum laterally ( m^4/s ) |
C momentum laterally ( m^4/s ) |
542 |
C viscA4W :: Biharmonic viscosity coeff. for mixing of vertical |
C viscA4W :: Biharmonic viscosity coeff. for mixing of vertical |
551 |
C (act on Divergence part) ( m^4/s ) |
C (act on Divergence part) ( m^4/s ) |
552 |
C viscA4Z :: Biharmonic viscosity coeff. for mixing of momentum laterally |
C viscA4Z :: Biharmonic viscosity coeff. for mixing of momentum laterally |
553 |
C (act on Vorticity part) ( m^4/s ) |
C (act on Vorticity part) ( m^4/s ) |
554 |
C viscC2leith :: Leith non-dimensional viscosity factor (grad(vort)) |
C viscC2leith :: Leith non-dimensional viscosity factor (grad(vort)) |
555 |
C viscC2leithD :: Modified Leith non-dimensional visc. factor (grad(div)) |
C viscC2leithD :: Modified Leith non-dimensional visc. factor (grad(div)) |
556 |
C viscC2smag :: Smagorinsky non-dimensional viscosity factor (harmonic) |
C viscC4leith :: Leith non-dimensional viscosity factor (grad(vort)) |
557 |
C viscC4smag :: Smagorinsky non-dimensional viscosity factor (biharmonic) |
C viscC4leithD :: Modified Leith non-dimensional viscosity factor (grad(div)) |
558 |
C viscAhMax :: Maximum eddy viscosity coeff. for mixing of |
C viscC2smag :: Smagorinsky non-dimensional viscosity factor (harmonic) |
559 |
C momentum laterally ( m^2/s ) |
C viscC4smag :: Smagorinsky non-dimensional viscosity factor (biharmonic) |
560 |
C viscAhReMax :: Maximum gridscale Reynolds number for eddy viscosity |
C viscAhMax :: Maximum eddy viscosity coeff. for mixing of |
561 |
C coeff. for mixing of momentum laterally (non-dim) |
C momentum laterally ( m^2/s ) |
562 |
|
C viscAhReMax :: Maximum gridscale Reynolds number for eddy viscosity |
563 |
|
C coeff. for mixing of momentum laterally (non-dim) |
564 |
|
C viscAhGrid :: non-dimensional grid-size dependent viscosity |
565 |
C viscAhGridMax:: maximum and minimum harmonic viscosity coefficients ... |
C viscAhGridMax:: maximum and minimum harmonic viscosity coefficients ... |
566 |
C viscAhGridMin:: in terms of non-dimensional grid-size dependent visc. |
C viscAhGridMin:: in terms of non-dimensional grid-size dependent visc. |
567 |
C viscA4Max :: Maximum biharmonic viscosity coeff. for mixing of |
C viscA4Max :: Maximum biharmonic viscosity coeff. for mixing of |
568 |
C momentum laterally ( m^4/s ) |
C momentum laterally ( m^4/s ) |
569 |
C viscA4ReMax :: Maximum Gridscale Reynolds number for |
C viscA4ReMax :: Maximum Gridscale Reynolds number for |
570 |
C biharmonic viscosity coeff. momentum laterally (non-dim) |
C biharmonic viscosity coeff. momentum laterally (non-dim) |
571 |
C viscAhGrid:: non-dimensional grid-size dependent viscosity |
C viscA4Grid :: non-dimensional grid-size dependent bi-harmonic viscosity |
|
C viscA4Grid:: non-dimensional grid-size dependent bi-harmonic viscosity |
|
572 |
C viscA4GridMax:: maximum and minimum biharmonic viscosity coefficients ... |
C viscA4GridMax:: maximum and minimum biharmonic viscosity coefficients ... |
573 |
C viscA4GridMin:: in terms of non-dimensional grid-size dependent viscosity |
C viscA4GridMin:: in terms of non-dimensional grid-size dependent viscosity |
|
C viscC4leith :: Leith non-dimensional viscosity factor (grad(vort)) |
|
|
C viscC4leithD :: Modified Leith non-dimensional viscosity factor (grad(div)) |
|
574 |
C diffKhT :: Laplacian diffusion coeff. for mixing of |
C diffKhT :: Laplacian diffusion coeff. for mixing of |
575 |
C heat laterally ( m^2/s ) |
C heat laterally ( m^2/s ) |
576 |
C diffKrNrT :: vertical profile of Laplacian diffusion coeff. |
C diffKrNrT :: vertical profile of Laplacian diffusion coeff. |
577 |
C for mixing of heat vertically ( units of r^2/s ) |
C for mixing of heat vertically ( units of r^2/s ) |
578 |
C diffK4T :: Biharmonic diffusion coeff. for mixing of |
C diffK4T :: Biharmonic diffusion coeff. for mixing of |
579 |
C heat laterally ( m^4/s ) |
C heat laterally ( m^4/s ) |
580 |
C diffKhS :: Laplacian diffusion coeff. for mixing of |
C diffKhS :: Laplacian diffusion coeff. for mixing of |
581 |
C salt laterally ( m^2/s ) |
C salt laterally ( m^2/s ) |
582 |
C diffKrNrS :: vertical profile of Laplacian diffusion coeff. |
C diffKrNrS :: vertical profile of Laplacian diffusion coeff. |
583 |
C for mixing of salt vertically ( units of r^2/s ), |
C for mixing of salt vertically ( units of r^2/s ), |
584 |
C diffK4S :: Biharmonic diffusion coeff. for mixing of |
C diffK4S :: Biharmonic diffusion coeff. for mixing of |
585 |
C salt laterally ( m^4/s ) |
C salt laterally ( m^4/s ) |
586 |
C diffKrBL79surf :: T/S surface diffusivity (m^2/s) Bryan and Lewis, 1979 |
C diffKrBL79surf :: T/S surface diffusivity (m^2/s) Bryan and Lewis, 1979 |
639 |
C cosPower :: Power of cosine of latitude to multiply viscosity |
C cosPower :: Power of cosine of latitude to multiply viscosity |
640 |
C cAdjFreq :: Frequency of convective adjustment |
C cAdjFreq :: Frequency of convective adjustment |
641 |
C |
C |
642 |
C taveFreq :: Frequency with which time-averaged model state |
C taveFreq :: Frequency with which time-averaged model state |
643 |
C is written to post-processing files ( s ). |
C is written to post-processing files ( s ). |
644 |
C tave_lastIter :: (for state variable only) fraction of the last time |
C tave_lastIter :: (for state variable only) fraction of the last time |
645 |
C step (of each taveFreq period) put in the time average. |
C step (of each taveFreq period) put in the time average. |
646 |
C (fraction for 1rst iter = 1 - tave_lastIter) |
C (fraction for 1rst iter = 1 - tave_lastIter) |
647 |
C tauThetaClimRelax :: Relaxation to climatology time scale ( s ). |
C tauThetaClimRelax :: Relaxation to climatology time scale ( s ). |
648 |
C tauSaltClimRelax :: Relaxation to climatology time scale ( s ). |
C tauSaltClimRelax :: Relaxation to climatology time scale ( s ). |
654 |
C number times externForcingPeriod) |
C number times externForcingPeriod) |
655 |
C convertFW2Salt :: salinity, used to convert Fresh-Water Flux to Salt Flux |
C convertFW2Salt :: salinity, used to convert Fresh-Water Flux to Salt Flux |
656 |
C (use model surface (local) value if set to -1) |
C (use model surface (local) value if set to -1) |
657 |
C temp_EvPrRn :: temperature of Rain & Evap. |
C temp_EvPrRn :: temperature of Rain & Evap. |
658 |
C salt_EvPrRn :: salinity of Rain & Evap. |
C salt_EvPrRn :: salinity of Rain & Evap. |
659 |
C (notes: a) tracer content of Rain/Evap only used if both |
C (notes: a) tracer content of Rain/Evap only used if both |
660 |
C NonLin_FrSurf & useRealFreshWater are set. |
C NonLin_FrSurf & useRealFreshWater are set. |
661 |
C b) use model surface (local) value if set to UNSET_RL) |
C b) use model surface (local) value if set to UNSET_RL) |
662 |
C hMixCrit :: criteria for mixed-layer diagnostic |
C hMixCrit :: criteria for mixed-layer diagnostic |
663 |
C ivdc_kappa :: implicit vertical diffusivity for convection [m^2/s] |
C ivdc_kappa :: implicit vertical diffusivity for convection [m^2/s] |
664 |
C Ro_SeaLevel :: standard position of Sea-Level in "R" coordinate, used as |
C Ro_SeaLevel :: standard position of Sea-Level in "R" coordinate, used as |
665 |
C starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel) |
C starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel) |
666 |
C sideDragFactor :: side-drag scaling factor (used only if no_slip_sides) |
C sideDragFactor :: side-drag scaling factor (used only if no_slip_sides) |
667 |
C (default=2: full drag ; =1: gives half-slip BC) |
C (default=2: full drag ; =1: gives half-slip BC) |
668 |
C bottomDragLinear :: Linear bottom-drag coefficient (units of [r]/s) |
C bottomDragLinear :: Linear bottom-drag coefficient (units of [r]/s) |
669 |
C bottomDragQuadratic :: Quadratic bottom-drag coefficient (units of [r]/m) |
C bottomDragQuadratic :: Quadratic bottom-drag coefficient (units of [r]/m) |
687 |
& viscC2leith, viscC2leithD, |
& viscC2leith, viscC2leithD, |
688 |
& viscC2smag, viscC4smag, |
& viscC2smag, viscC4smag, |
689 |
& viscAhD, viscAhZ, viscA4D, viscA4Z, |
& viscAhD, viscAhZ, viscA4D, viscA4Z, |
690 |
& viscA4, viscA4W, |
& viscA4, viscA4W, viscA4Max, |
691 |
& viscA4Max, viscA4Grid, viscA4GridMax, viscA4GridMin, |
& viscA4Grid, viscA4GridMax, viscA4GridMin, |
692 |
& viscAhRemax, viscA4Remax, |
& viscAhReMax, viscA4ReMax, |
693 |
& viscC4leith, viscC4leithD, viscAr, |
& viscC4leith, viscC4leithD, viscAr, |
694 |
& diffKhT, diffK4T, diffKrNrT, |
& diffKhT, diffK4T, diffKrNrT, |
695 |
& diffKhS, diffK4S, diffKrNrS, |
& diffKhS, diffK4S, diffKrNrS, |
748 |
_RL hFacInf |
_RL hFacInf |
749 |
_RL hFacSup |
_RL hFacSup |
750 |
_RL beta |
_RL beta |
751 |
|
_RL viscAr |
752 |
_RL viscAh |
_RL viscAh |
753 |
_RL viscAhW |
_RL viscAhW |
754 |
_RL viscAhD |
_RL viscAhD |
755 |
_RL viscAhZ |
_RL viscAhZ |
756 |
_RL viscAhMax |
_RL viscAhMax |
757 |
_RL viscAhReMax |
_RL viscAhReMax |
758 |
_RL viscAhGrid |
_RL viscAhGrid, viscAhGridMax, viscAhGridMin |
|
_RL viscAhGridMax |
|
|
_RL viscAhGridMin |
|
759 |
_RL viscC2leith |
_RL viscC2leith |
760 |
_RL viscC2leithD |
_RL viscC2leithD |
761 |
_RL viscC2smag |
_RL viscC2smag |
|
_RL viscC4smag |
|
|
_RL viscAr |
|
762 |
_RL viscA4 |
_RL viscA4 |
763 |
_RL viscA4W |
_RL viscA4W |
764 |
_RL viscA4D |
_RL viscA4D |
768 |
_RL viscA4Grid, viscA4GridMax, viscA4GridMin |
_RL viscA4Grid, viscA4GridMax, viscA4GridMin |
769 |
_RL viscC4leith |
_RL viscC4leith |
770 |
_RL viscC4leithD |
_RL viscC4leithD |
771 |
|
_RL viscC4smag |
772 |
_RL diffKhT |
_RL diffKhT |
773 |
_RL diffKrNrT(Nr) |
_RL diffKrNrT(Nr) |
774 |
_RL diffK4T |
_RL diffK4T |
859 |
C atm_Rq :: water vapour specific volume anomaly relative to dry air |
C atm_Rq :: water vapour specific volume anomaly relative to dry air |
860 |
C (e.g. typical value = (29/18 -1) 10^-3 with q [g/kg]) |
C (e.g. typical value = (29/18 -1) 10^-3 with q [g/kg]) |
861 |
C integr_GeoPot :: option to select the way we integrate the geopotential |
C integr_GeoPot :: option to select the way we integrate the geopotential |
862 |
C (still a subject of discussions ...) |
C (still a subject of discussions ...) |
863 |
C selectFindRoSurf :: select the way surf. ref. pressure (=Ro_surf) is |
C selectFindRoSurf :: select the way surf. ref. pressure (=Ro_surf) is |
864 |
C derived from the orography. Implemented: 0,1 (see INI_P_GROUND) |
C derived from the orography. Implemented: 0,1 (see INI_P_GROUND) |
865 |
COMMON /PARM_ATM/ |
COMMON /PARM_ATM/ |