micro_mg_utils Module Reference

Data Types
interface	mghydrometeorprops
interface	rising_factorial
interface	size_dist_param_basic
interface	size_dist_param_ice
interface	size_dist_param_liq
interface	var_coef

Functions/Subroutines
subroutine, public	micro_mg_utils_init (kind, rair, rh2o, cpair, tmelt_in, latvap, latice, dcs)
	Initialize module variables. More...
type(mghydrometeorprops) function	newmghydrometeorprops (rho, eff_dim, lambda_bounds, min_mean_mass)
	Constructor for a constituent property object. More...
pure real(r8) function	rising_factorial_r8 (x, n)
pure real(r8) function	rising_factorial_integer (x, n)
elemental real(r8) function	calc_ab (t, qv, xxl)
elemental subroutine	size_dist_param_liq_line (props, qcic, ncic, rho, pgam, lamc)
	get cloud droplet size distribution parameters More...
subroutine	size_dist_param_liq_vect (props, qcic, ncic, rho, pgam, lamc, mgncol)
	This subroutine gets cloud droplet size distribution parameters. More...
elemental subroutine	size_dist_param_basic_line (props, qic, nic, lam, n0)
	Basic routine for getting size distribution parameters. More...
subroutine	size_dist_param_basic_vect (props, qic, nic, lam, mgncol, n0)
	This subroutine calculates. More...
elemental subroutine	size_dist_param_ice_line (props, qic, nic, lam, n0)
	ice routine for getting size distribution parameters. More...
subroutine	size_dist_param_ice_vect (props, qic, nic, lam, mgncol, n0)
	This subroutine. More...
real(r8) elemental function, public	avg_diameter (q, n, rho_air, rho_sub)
	Finds the average diameter of particles given their density, and mass/number concentrations in the air. Assumes that diameter follows an exponential distribution. More...
elemental real(r8) function	var_coef_r8 (relvar, a)
	Finds a coefficient for process rates based on the relative variance of cloud water. More...
elemental real(r8) function	var_coef_integer (relvar, a)
	Finds a coefficient for process rates based on the relative variance of cloud water. More...
subroutine, public	ice_deposition_sublimation (t, qv, qi, ni, icldm, rho, dv, qvl, qvi, berg, vap_dep, ice_sublim, mgncol)
	Initial ice deposition and sublimation loop. Run before the main loop This subroutine written by Peter Caldwell. More...
subroutine, public	kk2000_liq_autoconversion (microp_uniform, qcic, ncic, rho, relvar, prc, nprc, nprc1, mgncol)
	autoconversion of cloud liquid water to rain formula from Khrouditnov and Kogan (2000), modified for sub-grid distribution of qc minimum qc of 1 x 10^-8 prevents floating point error More...
subroutine, public	sb2001v2_liq_autoconversion (pgam, qc, nc, qr, rho, relvar, au, nprc, nprc1, mgncol)
	This subroutine. More...
subroutine, public	liu_liq_autoconversion (pgam, qc, nc, qr, rho, relvar, au, nprc, nprc1, mgncol)
	Anning Cheng 10/5/2017 add Liu et al. autoconversion. More...
subroutine, public	sb2001v2_accre_cld_water_rain (qc, nc, qr, rho, relvar, pra, npra, mgncol)
subroutine, public	ice_autoconversion (t, qiic, lami, n0i, dcs, ac_time, prci, nprci, mgncol)
	Autoconversion of cloud ice to snow similar to Ferrier (1994) More...
subroutine, public	gmao_ice_autoconversion (t, qiic, niic, lami, n0i, dcs, ac_time, prci, nprci, mgncol)
	GMAO ice autoconversion. More...
subroutine, public	immersion_freezing (microp_uniform, t, pgam, lamc, qcic, ncic, relvar, mnuccc, nnuccc, mgncol)
	immersion freezing (Bigg, 1953) More...
subroutine, public	contact_freezing (microp_uniform, t, p, rndst, nacon, pgam, lamc, qcic, ncic, relvar, mnucct, nnucct, mgncol, mdust)
	contact freezing (-40<T<-3 C) (Young, 1974) with hooks into simulated dust dust size and number in multiple bins are read in from companion routine More...
subroutine, public	snow_self_aggregation (t, rho, asn, rhosn, qsic, nsic, nsagg, mgncol)
	snow self-aggregation from passarelli, 1978, used by reisner, 1998 More...
subroutine, public	accrete_cloud_water_snow (t, rho, asn, uns, mu, qcic, ncic, qsic, pgam, lamc, lams, n0s, psacws, npsacws, mgncol)
	accretion of cloud droplets onto snow/graupel More...
subroutine, public	secondary_ice_production (t, psacws, msacwi, nsacwi, mgncol)
	add secondary ice production due to accretion of droplets by snow More...
subroutine, public	accrete_rain_snow (t, rho, umr, ums, unr, uns, qric, qsic, lamr, n0r, lams, n0s, pracs, npracs, mgncol)
	accretion of rain water by snow More...
subroutine, public	heterogeneous_rain_freezing (t, qric, nric, lamr, mnuccr, nnuccr, mgncol)
	heterogeneous freezing of rain drops More...
subroutine, public	accrete_cloud_water_rain (microp_uniform, qric, qcic, ncic, relvar, accre_enhan, pra, npra, mgncol)
	accretion of cloud liquid water by rain formula from Khrouditnov and Kogan (2000) More...
subroutine, public	self_collection_rain (rho, qric, nric, nragg, mgncol)
	Self-collection of rain drops from Beheng(1994) More...
subroutine, public	accrete_cloud_ice_snow (t, rho, asn, qiic, niic, qsic, lams, n0s, prai, nprai, mgncol)
	Accretion of cloud ice by snow. More...
subroutine, public	evaporate_sublimate_precip (t, rho, dv, mu, sc, q, qvl, qvi, lcldm, precip_frac, arn, asn, qcic, qiic, qric, qsic, lamr, n0r, lams, n0s, pre, prds, am_evp_st, mgncol)
	calculate evaporation/sublimation of rain and snow More...
subroutine	evaporate_sublimate_precip_graupel (t, rho, dv, mu, sc, q, qvl, qvi, lcldm, precip_frac, arn, asn, agn, bg, qcic, qiic, qric, qsic, qgic, lamr, n0r, lams, n0s, lamg, n0g, pre, prds, prdg, am_evp_st, mgncol)
	evaporation/sublimation of rain, snow and graupel More...
subroutine, public	bergeron_process_snow (t, rho, dv, mu, sc, qvl, qvi, asn, qcic, qsic, lams, n0s, bergs, mgncol)
	bergeron process - evaporation of droplets and deposition onto snow More...
subroutine	graupel_collecting_snow (qsic, qric, umr, ums, rho, lamr, n0r, lams, n0s, psacr, mgncol)
	Collection of snow by rain to form graupel. More...
subroutine	graupel_collecting_cld_water (qgic, qcic, ncic, rho, n0g, lamg, bg, agn, psacwg, npsacwg, mgncol)
	Collection of cloud water by graupel. More...
subroutine	graupel_riming_liquid_snow (psacws, qsic, qcic, nsic, rho, rhosn, rhog, asn, lams, n0s, dtime, pgsacw, nscng, mgncol)
	Conversion of rimed cloud water onto snow to graupel/hail. More...
subroutine	graupel_collecting_rain (qric, qgic, umg, umr, ung, unr, rho, n0r, lamr, n0g, lamg, pracg, npracg, mgncol)
	CHANGE IN Q,N COLLECTION RAIN BY GRAUPEL. More...
subroutine	graupel_rain_riming_snow (pracs, npracs, psacr, qsic, qric, nric, nsic, n0s, lams, n0r, lamr, dtime, pgracs, ngracs, mgncol)
	Rain riming snow to graupel. More...
subroutine	graupel_rime_splintering (t, qcic, qric, qgic, psacwg, pracg, qmultg, nmultg, qmultrg, nmultrg, mgncol)
	Rime splintering. More...
pure real(r8) function	no_limiter ()
pure logical function	limiter_is_on (lim)
real(r8) function	gamma_incomp (muice, x)

Variables
type(mghydrometeorprops), public	mg_liq_props
type(mghydrometeorprops), public	mg_ice_props
type(mghydrometeorprops), public	mg_rain_props
type(mghydrometeorprops), public	mg_snow_props
type(mghydrometeorprops), public	mg_graupel_props
real(r8), parameter, public	pi = 3.14159265358979323846_r8
	Pi to 20 digits; more than enough to reach the limit of double precision. More...
real(r8), parameter, public	omsm = 1._r8 - 1.e-6_r8
	"One minus small number": number near unity for round-off issues. More...
real(r8), parameter, public	qsmall = 1.e-18_r8
	Smallest mixing ratio considered in microphysics. More...
real(r8), parameter, public	mincld = 0.000001_r8
	minimum allowed cloud fraction More...
real(r8), parameter, public	rhosn = 250._r8
	bulk density snow More...
real(r8), parameter, public	rhoi = 500._r8
	bulk density ice More...
real(r8), parameter, public	rhow = 1000._r8
	bulk density liquid More...
real(r8), parameter, public	rhows = 917._r8
	bulk density water solid More...
real(r8), parameter, public	rhog = 500._r8
real(r8), parameter, public	rhoh = 400._r8
real(r8), parameter, public	ac = 3.e7_r8
real(r8), parameter, public	bc = 2._r8
real(r8), parameter, public	as = 11.72_r8
real(r8), parameter, public	bs = 0.41_r8
real(r8), parameter, public	ai = 700._r8
real(r8), parameter, public	bi = 1._r8
real(r8), parameter, public	aj = ac*((rhoi/rhows)**(bc/3._r8))*rhows/rhow
real(r8), parameter, public	bj = bc
real(r8), parameter, public	ar = 841.99667_r8
real(r8), parameter, public	br = 0.8_r8
real(r8), parameter, public	ag = 19.3_r8
real(r8), parameter, public	bg = 0.37_r8
real(r8), parameter, public	ah = 114.5_r8
real(r8), parameter, public	bh = 0.5_r8
real(r8), parameter, public	mi0 = 4._r8/3._r8*pi*rhoi*(1.e-6_r8)**3
	mass of new crystal due to aerosol freezing and growth (kg) Make this consistent with the lower bound, to support UTLS and stratospheric ice, and the smaller ice size limit. More...
real(r8), parameter, public	mg0 = 1.6e-10_r8
real(r8), parameter, public	mmult = 4._r8/3._r8*pi*rhoi*(5.e-6_r8)**3
integer(i8), parameter	limiter_off = int(Z'7FF1111111111111', i8)
real(r8), parameter	icsmall = 1.e-8_r8
real(r8), parameter	dsph = 3._r8
real(r8), dimension(2), parameter	lam_bnd_rain = 1._r8/[500.e-6_r8, 20.e-6_r8]
real(r8), dimension(2), parameter	lam_bnd_snow = 1._r8/[2000.e-6_r8, 10.e-6_r8]
real(r8), parameter	min_mean_mass_liq = 1.e-20_r8
real(r8), parameter	min_mean_mass_ice = 1.e-20_r8
real(r8), parameter	f1s = 0.86_r8
real(r8), parameter	f2s = 0.28_r8
real(r8), parameter	f1r = 0.78_r8
real(r8), parameter	f2r = 0.308_r8
real(r8), parameter	eii = 0.2_r8
real(r8), parameter, public	ecid = 0.7_r8
real(r8), parameter, public	ecr = 1.0_r8
real(r8), parameter	bimm = 100._r8
real(r8), parameter	aimm = 0.66_r8
real(r8), parameter	droplet_mass_25um = 4._r8/3._r8*pi*rhow*(25.e-6_r8)**3
real(r8), parameter	droplet_mass_40um = 4._r8/3._r8*pi*rhow*(40.e-6_r8)**3
real(r8), parameter	droplet_mass_40umi = 1._r8/droplet_mass_40um
real(r8)	rv
real(r8)	cpp
real(r8)	tmelt
real(r8)	ra
real(r8)	xxlv
real(r8)	xlf
real(r8)	xxls
real(r8)	gamma_bs_plus3
real(r8)	gamma_half_br_plus5
real(r8)	gamma_half_bs_plus5
real(r8)	gamma_2bs_plus2
real(r8), parameter	zero = 0._r8
real(r8), parameter	one = 1._r8
real(r8), parameter	two = 2._r8
real(r8), parameter	three = 3._r8
real(r8), parameter	four = 4._r8
real(r8), parameter	five = 5._r8
real(r8), parameter	six = 6._r8
real(r8), parameter	pio6 = pi/six
real(r8), parameter	pio3 = pi/three
real(r8), parameter	half = 0.5_r8
real(r8), parameter	oneo3 = one/three
real(r8), parameter	twopi = pi + pi