Skip to content
Snippets Groups Projects
Normal view Permalink
sfx_esm_param.f90 1.29 KiB
Newer Older
Evgeny Mortikov's avatar
major code update
Evgeny Mortikov committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 
module sfx_esm_param
      !> @brief ESM surface flux model parameters
      !> @details  all in SI units

      ! modules used
      ! --------------------------------------------------------------------------------
      use sfx_phys_const
      ! --------------------------------------------------------------------------------

      ! directives list
      ! --------------------------------------------------------------------------------
      implicit none
      ! --------------------------------------------------------------------------------


      !> von Karman constant [n/d]
      real, parameter :: kappa = 0.40
      !> inverse Prandtl (turbulent) number in neutral conditions [n/d]
      real, parameter :: Pr_t_0_inv = 1.15
      !> inverse Prandtl (turbulent) number in free convection [n/d]
      real, parameter :: Pr_t_inf_inv = 3.5

      !> stability function coeff. (unstable) [= g4 & g10 in deprecated code]
      real, parameter :: alpha_m = 16.0
      real, parameter :: alpha_h = 16.0
      real, parameter :: alpha_h_fix = 1.2

      !> stability function coeff. (stable)
      real, parameter :: beta_m = 4.7
      real, parameter :: beta_h = beta_m

      !> --- max Ri-bulk value in stable case ( < 1 / beta_m )
      real, parameter ::  Rib_max = 0.9 / beta_m

end module sfx_esm_param