parameters.pxi

#Adapated from PyCLES: https://github.com/pressel/pycles
cdef double pi = 3.14159265359
cdef double g = 9.80665
cdef double Rd = 287.1
cdef double Rv = 461.5
cdef double eps_v = 0.62210184182   # Rd / Rv
cdef double eps_vi = 1.60745384883  # Rv / Rd
cdef double cpd = 1004.0
cdef double cpv = 1859.0
cdef double cl = 4218.0
cdef double ci = 2106.0
cdef double kappa = 0.285956175299
cdef double Tf = 273.15
cdef double Tt = 273.16
cdef double T_tilde = 298.15
cdef double p_tilde = 100000.0
cdef double pv_star_t = 611.7
cdef double sd_tilde = 6864.8
cdef double sv_tilde = 10513.6
cdef double omega = 7.29211514671e-05
cdef double ql_threshold = 1e-08
cdef double vkb = 0.4
cdef double Pr0 = 0.74
cdef double beta_m = 4.8
cdef double beta_h = 6.5 # beta_h = beta_m / Pr0
cdef double gamma_m = 15.0
cdef double gamma_h = 9.0
# constants defined in Stevens et al 2005 (that are different from scampy)
# needed for DYCOMS case setup
cdef double dycoms_cp = 1015.
cdef double dycoms_L = 2.47 * 1e6
cdef double dycoms_Rd = 287.
# CLIMA microphysics parameters
cdef double C_drag = 0.55
cdef double rho_cloud_liq = 1e3
cdef double MP_n_0 = 16 * 1e6
cdef double tau_cond_evap = 10
cdef double q_liq_threshold = 5e-4 #5e-4
cdef double tau_acnv = 1e3
cdef double E_col = 0.8
cdef double nu_air = 1.6e-5
cdef double K_therm = 2.4e-2
cdef double D_vapor = 2.26e-5
cdef double a_vent = 1.5
cdef double b_vent = 0.53
# cutoff microphysics parameter
cdef double max_supersaturation = 0.02