module sfx_sheba_noit_param
    !< @brief SHEBA 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

    real, parameter :: Rib_max = 0.5

    !< stability function coeff. (unstable)
    real, parameter :: alpha_m = 16.0
    real, parameter :: alpha_h = 16.0
    real, parameter :: alpha_h_fix = 16.0

    !< stability function coeff. (stable)
    real, parameter :: a_m = 5.0
    real, parameter :: b_m = a_m / 6.5
    real, parameter :: a_h = 5.0
    real, parameter :: b_h = 5.0
    real, parameter :: c_h = 3.0
    
    real, parameter  ::  gamma = 2.91, zeta_a = 3.6 ! for stable psi

        	!< snow parameters
    real, parameter :: rho_s =900
    real, parameter :: d_s=0.0000886
    real, parameter :: h_salt_const=0.05
    real, parameter :: u_thsnow=0.35   
    real, parameter :: S_salt_const=0.000004   
    real, parameter :: rho_air=1.2
    real, parameter :: Csn=3.25


end module sfx_sheba_noit_param