/[MITgcm]/MITgcm/model/inc/PARAMS.h
ViewVC logotype

Diff of /MITgcm/model/inc/PARAMS.h

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.213 by mlosch, Mon Feb 25 20:57:24 2008 UTC revision 1.214 by jmc, Sun Mar 30 21:44:09 2008 UTC
# Line 173  C     tempVertAdvScheme   :: Temp. Vert. Line 173  C     tempVertAdvScheme   :: Temp. Vert.
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
# Line 189  C     debugLevel          :: debug level Line 190  C     debugLevel          :: debug level
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
# Line 208  C     debugLevel          :: debug level Line 209  C     debugLevel          :: debug level
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    
# Line 274  C     use3dCoriolis :: Turns the 3-D cor Line 276  C     use3dCoriolis :: Turns the 3-D cor
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)
# Line 332  C                        & Last iteratio Line 333  C                        & Last iteratio
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,
# Line 354  C                        according to Eu Line 355  C                        according to Eu
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,
# Line 419  C                        according to Eu Line 420  C                        according to Eu
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
# Line 529  C                 ( Southern edge f for Line 529  C                 ( Southern edge f for
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
# Line 551  C     viscA4D   :: Biharmonic viscosity Line 551  C     viscA4D   :: Biharmonic viscosity
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
# Line 639  C     mtFacMom      :: Metric terms trac Line 639  C     mtFacMom      :: Metric terms trac
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 ).
# Line 654  C                          (note: extern Line 654  C                          (note: extern
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)
# Line 687  C     psiEuler      :: Euler angle, rota Line 687  C     psiEuler      :: Euler angle, rota
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,
# Line 748  C     psiEuler      :: Euler angle, rota Line 748  C     psiEuler      :: Euler angle, rota
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
# Line 771  C     psiEuler      :: Euler angle, rota Line 768  C     psiEuler      :: Euler angle, rota
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
# Line 861  C     atm_kappa :: kappa = R/Cp (R: cons Line 859  C     atm_kappa :: kappa = R/Cp (R: cons
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/
Legend:
Removed from v.1.213  
changed lines
  Added in v.1.214
  ViewVC Help
Powered by ViewVC 1.1.22