This module computes the moisture tendencies of water vapor, cloud droplets, rain, cloud ice (pristine), snow, and graupel. Prior to WRFv2.2 this code was based on Reisner et al (1998), but few of those pieces remain. A complete description is now found in Thompson, G., P. R. Field, R. M. Rasmussen, and W. D. Hall, 2008: Explicit Forecasts of winter precipitation using an improved bulk microphysics scheme. Part II: Implementation of a new snow parameterization. Mon. Wea. Rev., 136, 5095-5115. More...

Functions/Subroutines
subroutine	thompson_final ()

subroutine	qr_acr_qg
	Rain collecting graupel (and inverse). Explicit CE integration.

subroutine	qr_acr_qs
	Rain collecting snow (and inverse). Explicit CE integration.

subroutine	freezeh2o (threads)
	This is a literal adaptation of Bigg (1954) probability of drops of a particular volume freezing. Given this probability, simply freeze the proportion of drops summing their masses.

subroutine	qi_aut_qs
	Cloud ice converting to snow since portion greater than min snow size. Given cloud ice content (kg/m3), number concentration (#/m3) and gamma shape parameter, mu_i, break the distrib into bins and figure out the mass/number of ice with sizes larger than D0s. Also, compute incomplete gamma function for the integration of ice depositional growth from diameter=0 to D0s. Amount of ice depositional growth is this portion of distrib while larger diameters contribute to snow growth (as in Harrington et al. 1995).

subroutine	table_efrw
	Variable collision efficiency for rain collecting cloud water using method of Beard and Grover, 1974 if a/A less than 0.25; otherwise uses polynomials to get close match of Pruppacher & Klett Fig 14-9.

subroutine	table_efsw
	Variable collision efficiency for snow collecting cloud water using method of Wang and Ji, 2000 except equate melted snow diameter to their "effective collision cross-section.".

real function	eff_aero (d, da, visc, rhoa, temp, species)
	Function to compute collision efficiency of collector species (rain, snow, graupel) of aerosols. Follows Wang et al, 2010, ACP, which follows Slinn (1983).

subroutine	table_dropevap
	Integrate rain size distribution from zero to D-star to compute the number of drops smaller than D-star that evaporate in a single timestep. Drops larger than D-star dont evaporate entirely so do not affect number concentration.

subroutine	table_ccnact (errmess, errflag)
	Fill the table of CCN activation data created from parcel model run by Trude Eidhammer with inputs of aerosol number concentration, vertical velocity, temperature, lognormal mean aerosol radius, and hygroscopicity, kappa. The data are read from external file and contain activated fraction of CCN for given conditions.

real function	activ_ncloud (tt, ww, nccn, lsm_in)
	Retrieve fraction of CCN that gets activated given the model temp, vertical velocity, and available CCN concentration. The lookup table (read from external file) has CCN concentration varying the quickest, then updraft, then temperature, then mean aerosol radius, and finally hygroscopicity, kappa.

subroutine	gcf (gammcf, a, x, gln)
	Returns the incomplete gamma function q(a,x) evaluated by its continued fraction representation as gammcf.

subroutine	gser (gamser, a, x, gln)
	Returns the incomplete gamma function p(a,x) evaluated by its series representation as gamser.

real function	gammln (xx)
	Returns the value ln(gamma(xx)) for xx > 0.

real function	gammp (a, x)

real function	wgamma (y)

real function	rslf (p, t)
	THIS FUNCTION CALCULATES THE LIQUID SATURATION VAPOR MIXING RATIO AS A FUNCTION OF TEMPERATURE AND PRESSURE.

real function	rsif (p, t)
	THIS FUNCTION CALCULATES THE ICE SATURATION VAPOR MIXING RATIO AS A FUNCTION OF TEMPERATURE AND PRESSURE.

real function	icedemott (tempc, qv, qvs, qvsi, rho, nifa)

real function	icekoop (temp, qv, qvs, naero, dt)
	Newer research since Koop et al (2001) suggests that the freezing rate should be lower than original paper, so J_rate is reduced by two orders of magnitude.

real function	delta_p (yy, y1, y2, aa, bb)
	Helper routine for Phillips et al (2008) ice nucleation.

subroutine	calc_effectrad (t1d, p1d, qv1d, qc1d, nc1d, qi1d, ni1d, qs1d, re_qc1d, re_qi1d, re_qs1d, lsml, kts, kte)
	Compute radiation effective radii of cloud water, ice, and snow. These are entirely consistent with microphysics assumptions, not constant or otherwise ad hoc as is internal to most radiation schemes. Since only the smallest snowflakes should impact radiation, compute from first portion of complicated Field number distribution, not the second part, which is the larger sizes.

subroutine	calc_refl10cm (qv1d, qc1d, qr1d, nr1d, qs1d, qg1d, t1d, p1d, dbz, rand1, kts, kte, ii, jj, melti, vt_dbz, first_time_step)
	Compute radar reflectivity assuming 10 cm wavelength radar and using Rayleigh approximation. Only complication is melted snow/graupel which we treat as water-coated ice spheres and use Uli Blahak's library of routines. The meltwater fraction is simply the amount of frozen species remaining from what initially existed at the melting level interface.

subroutine	semi_lagrange_sedim (km, dzl, wwl, rql, precip, pfsan, dt, r1)
	This routine is a semi-Lagrangian forward advection for hydrometeors with mass conservation and positive definite advection 2nd order interpolation with monotonic piecewise parabolic method is used. This routine is under assumption of decfl < 1 for semi_Lagrangian (Juang and Hong, 2010 [87]).

subroutine	graupel_psd_parameters (kts, kte, rand1, rg, ilamg, n0_g)
	Calculates graupel size distribution parameters.

real(wp) function	hail_mass_99th_percentile (kts, kte, qg, temperature, pressure, qv)
	Calculates graupel/hail maximum diameter.


subroutine	thompson_init (is_aerosol_aware_in, merra2_aerosol_aware_in, mpicomm, mpirank, mpiroot, threads, errmsg, errflg)
	This subroutine calculates simplified cloud species equations and create lookup tables in Thomspson scheme.


subroutine	mp_gt_driver (wrf_chem)
	This is a wrapper routine designed to transfer values from 3D to 1D.


subroutine	mp_thompson (wrf_chem)
	This subroutine computes the moisture tendencies of water vapor, cloud droplets, rain, cloud ice (pristine), snow, and graupel. Previously this code was based on Reisner et al (1998), but few of those pieces remain. A complete description is now found in Thompson et al. (2004, 2008)[191] [192].

Variables
real(wp)	t_0

real(wp)	pi

real(wp)	nt_c_l = 150.e6

real(wp)	nt_c_o = 50.e6

real(wp)	av_i

real(wp)	xnc_max = 1000.e3

real(wp)	ssati_min = 0.15

real(wp)	nt_i_max = 4999.e3_dp

real(wp)	rr_min = 1000.0

real(wp), parameter	nain0 = 1.5e6

real(wp), parameter	nain1 = 0.5e6

real(wp), parameter	naccn0 = 300.0e6

real(wp), parameter	naccn1 = 50.0e6

real(wp), parameter	av_s = 40.0

real(wp), parameter	bv_s = 0.55

real(wp), parameter	bv_i = 1.0

real(wp), parameter	eps = 1.E-15

real(wp)	rho_not

real(wp)	rv

real(wp)	r

real(wp)	roverrv

real(wp)	cp

real(wp)	r_uni

real(dp)	k_b

real(dp)	m_w

real(dp)	m_a

real(dp)	n_avo

real(wp)	lvap0

real(wp)	lfus

real(wp), parameter	d0s = 300.e-6

real(wp), parameter	re_qc_min = 2.50e-6

real(wp), parameter	re_qc_max = 50.0e-6

real(wp), parameter	re_qi_min = 2.50e-6

real(wp), parameter	re_qi_max = 125.0e-6

real(wp), parameter	re_qs_min = 5.00e-6

real(wp), parameter	re_qs_max = 999.0e-6

real(dp), dimension(nbins+1)	xdx

real(dp), dimension(nbc)	dc

real(dp), dimension(nbc)	dtc

real(dp), dimension(nbi)	di

real(dp), dimension(nbi)	dti

real(dp), dimension(nbr)	dr

real(dp), dimension(nbr)	dtr

real(dp), dimension(nbs)	ds

real(dp), dimension(nbs)	dts

real(dp), dimension(nbg)	dg

real(dp), dimension(nbg)	dtg

real(dp), dimension(nbc)	t_nc

character(len= *), parameter	thomp_table_file = 'thompson_tables_precomp_v2.sl'

character(len= *), parameter	qr_acr_qg_file = 'qr_acr_qgV2.dat'

character(len= *), parameter	qr_acr_qs_file = 'qr_acr_qsV2.dat'

character(len= *), parameter	freeze_h2o_file = 'freezeH2O.dat'

real(dp), dimension(:,:,:,:), allocatable	tcg_racg

real(dp), dimension(:,:,:,:), allocatable	tmr_racg

real(dp), dimension(:,:,:,:), allocatable	tcr_gacr

real(dp), dimension(:,:,:,:), allocatable	tmg_gacr

real(dp), dimension(:,:,:,:), allocatable	tnr_racg

real(dp), dimension(:,:,:,:), allocatable	tnr_gacr

real(dp), dimension(:,:,:,:), allocatable	tcs_racs1

real(dp), dimension(:,:,:,:), allocatable	tmr_racs1

real(dp), dimension(:,:,:,:), allocatable	tcs_racs2

real(dp), dimension(:,:,:,:), allocatable	tmr_racs2

real(dp), dimension(:,:,:,:), allocatable	tcr_sacr1

real(dp), dimension(:,:,:,:), allocatable	tms_sacr1

real(dp), dimension(:,:,:,:), allocatable	tcr_sacr2

real(dp), dimension(:,:,:,:), allocatable	tms_sacr2

real(dp), dimension(:,:,:,:), allocatable	tnr_racs1

real(dp), dimension(:,:,:,:), allocatable	tnr_racs2

real(dp), dimension(:,:,:,:), allocatable	tnr_sacr1

real(dp), dimension(:,:,:,:), allocatable	tnr_sacr2

real(dp), dimension(:,:,:,:), allocatable	tpi_qcfz

real(dp), dimension(:,:,:,:), allocatable	tni_qcfz

real(dp), dimension(:,:,:,:), allocatable	tpi_qrfz

real(dp), dimension(:,:,:,:), allocatable	tpg_qrfz

real(dp), dimension(:,:,:,:), allocatable	tni_qrfz

real(dp), dimension(:,:,:,:), allocatable	tnr_qrfz

real(dp), dimension(:,:), allocatable	tps_iaus

real(dp), dimension(:,:), allocatable	tni_iaus

real(dp), dimension(:,:), allocatable	tpi_ide

real(dp), dimension(:,:), allocatable	t_efrw

real(dp), dimension(:,:), allocatable	t_efsw

real(dp), dimension(:,:,:), allocatable	tnr_rev

real(dp), dimension(:,:,:), allocatable	tpc_wev

real(dp), dimension(:,:,:), allocatable	tnc_wev

real(sp), dimension(:,:,:,:,:), allocatable	tnccn_act

real(wp)	t1_qr_qc

real(wp)	t1_qr_qi

real(wp)	t2_qr_qi

real(wp)	t1_qg_qc

real(wp)	t1_qs_qc

real(wp)	t1_qs_qi

real(wp)	t1_qr_ev

real(wp)	t2_qr_ev

real(wp)	t1_qs_sd

real(wp)	t2_qs_sd

real(wp)	t1_qg_sd

real(wp)	t2_qg_sd

real(wp)	t1_qs_me

real(wp)	t2_qs_me

real(wp)	t1_qg_me

real(wp)	t2_qg_me

type(mpi_comm)	mpi_communicator

logical	thompson_table_writer

Detailed Description

This module computes the moisture tendencies of water vapor, cloud droplets, rain, cloud ice (pristine), snow, and graupel. Prior to WRFv2.2 this code was based on Reisner et al (1998), but few of those pieces remain. A complete description is now found in Thompson, G., P. R. Field, R. M. Rasmussen, and W. D. Hall, 2008: Explicit Forecasts of winter precipitation using an improved bulk microphysics scheme. Part II: Implementation of a new snow parameterization. Mon. Wea. Rev., 136, 5095-5115.

Prior to WRFv3.1, this code was single-moment rain prediction as described in the reference above, but in v3.1 and higher, the scheme is two-moment rain (predicted rain number concentration).

Beginning with WRFv3.6, this is also the "aerosol-aware" scheme as described in Thompson, G. and T. Eidhammer, 2014: A study of aerosol impacts on clouds and precipitation development in a large winter cyclone. J. Atmos. Sci., 71, 3636-3658. Setting WRF namelist option mp_physics=8 utilizes the older one-moment cloud water with constant droplet concentration set as Nt_c (found below) while mp_physics=28 uses double-moment cloud droplet number concentration, which is not permitted to exceed Nt_c_max below.

Most importantly, users may wish to modify the prescribed number of cloud droplets (Nt_c; see guidelines mentioned below). Otherwise, users may alter the rain and graupel size distribution parameters to use exponential (Marshal-Palmer) or generalized gamma shape. The snow field assumes a combination of two gamma functions (from Field et al. 2005) and would require significant modifications throughout the entire code to alter its shape as well as accretion rates. Users may also alter the constants used for density of rain, graupel, ice, and snow, but the latter is not constant when using Paul Field's snow distribution and moments methods. Other values users can modify include the constants for mass and/or velocity power law relations and assumed capacitances used in deposition/ sublimation/evaporation/melting. Remaining values should probably be left alone.

Author: Greg Thompson, NCAR-RAL, gthom.nosp@m.psn@.nosp@m.ucar..nosp@m.edu, 303-497-2805

Last modified: 24 Jan 2018 Aerosol additions to v3.5.1 code 9/2013 Cloud fraction additions 11/2014 part of pre-v3.7
Imported in CCPP by: Dom Heinzeller, NOAA/ESRL/GSD, dom.h.nosp@m.einz.nosp@m.eller.nosp@m.@noa.nosp@m.a.gov
Last modified: 6 Aug 2018 Update of initial import to WRFV4.0
Last modified: 13 Mar 2020 Add logic to turtn on/off the calculation of melting layer in radar reflectivity routine
Last modified: 2 Jun 2020 Add option to rollback to version 3.8.1 used in RAPv5/HRRRv4, include stochastic physics perturbations to the graupel intercept parameter, the cloud water shape parameter, and the number concentration of nucleated aerosols.
Last modified: 12 Feb 2021 G. Thompson updated to align more closely with his WRF version, including bug fixes and designed changes.

Function/Subroutine Documentation

◆ semi_lagrange_sedim()

subroutine module_mp_thompson::semi_lagrange_sedim	(	integer, intent(in)	km,
		real(wp), dimension(km), intent(in)	dzl,
		real(wp), dimension(km), intent(in)	wwl,
		real(wp), dimension(km), intent(inout)	rql,
		real(wp), intent(out)	precip,
		real(wp), dimension(km), intent(out)	pfsan,
		real(wp), intent(in)	dt,
		real(wp), intent(in)	r1
	)

This routine is a semi-Lagrangian forward advection for hydrometeors with mass conservation and positive definite advection 2nd order interpolation with monotonic piecewise parabolic method is used. This routine is under assumption of decfl < 1 for semi_Lagrangian (Juang and Hong, 2010 [87]).

Definition at line 6231 of file module_mp_thompson.F90.

Variable Documentation

◆ av_i

real(wp) module_mp_thompson::av_i

Definition at line 98 of file module_mp_thompson.F90.

◆ av_s

real(wp), parameter module_mp_thompson::av_s = 40.0

Definition at line 154 of file module_mp_thompson.F90.

◆ bv_i

real(wp), parameter module_mp_thompson::bv_i = 1.0

Definition at line 159 of file module_mp_thompson.F90.

◆ bv_s

real(wp), parameter module_mp_thompson::bv_s = 0.55

Definition at line 155 of file module_mp_thompson.F90.

◆ cp

real(wp) module_mp_thompson::cp

Definition at line 203 of file module_mp_thompson.F90.

◆ d0s

real(wp), parameter module_mp_thompson::d0s = 300.e-6

Definition at line 224 of file module_mp_thompson.F90.

◆ dc

real(dp), dimension(nbc) module_mp_thompson::dc

Definition at line 263 of file module_mp_thompson.F90.

◆ dg

real(dp), dimension(nbg) module_mp_thompson::dg

Definition at line 267 of file module_mp_thompson.F90.

◆ di

real(dp), dimension(nbi) module_mp_thompson::di

Definition at line 264 of file module_mp_thompson.F90.

◆ dr

real(dp), dimension(nbr) module_mp_thompson::dr

Definition at line 265 of file module_mp_thompson.F90.

◆ ds

real(dp), dimension(nbs) module_mp_thompson::ds

Definition at line 266 of file module_mp_thompson.F90.

◆ dtc

real(dp), dimension(nbc) module_mp_thompson::dtc

Definition at line 263 of file module_mp_thompson.F90.

◆ dtg

real(dp), dimension(nbg) module_mp_thompson::dtg

Definition at line 267 of file module_mp_thompson.F90.

◆ dti

real(dp), dimension(nbi) module_mp_thompson::dti

Definition at line 264 of file module_mp_thompson.F90.

◆ dtr

real(dp), dimension(nbr) module_mp_thompson::dtr

Definition at line 265 of file module_mp_thompson.F90.

◆ dts

real(dp), dimension(nbs) module_mp_thompson::dts

Definition at line 266 of file module_mp_thompson.F90.

◆ eps

real(wp), parameter module_mp_thompson::eps = 1.E-15

Definition at line 182 of file module_mp_thompson.F90.

◆ freeze_h2o_file

character(len=*), parameter module_mp_thompson::freeze_h2o_file = 'freezeH2O.dat'

Definition at line 385 of file module_mp_thompson.F90.

◆ k_b

real(dp) module_mp_thompson::k_b

Definition at line 206 of file module_mp_thompson.F90.

◆ lfus

real(wp) module_mp_thompson::lfus

Definition at line 216 of file module_mp_thompson.F90.

◆ lvap0

real(wp) module_mp_thompson::lvap0

Definition at line 215 of file module_mp_thompson.F90.

◆ m_a

real(dp) module_mp_thompson::m_a

Definition at line 208 of file module_mp_thompson.F90.

◆ m_w

real(dp) module_mp_thompson::m_w

Definition at line 207 of file module_mp_thompson.F90.

◆ mpi_communicator

type(mpi_comm) module_mp_thompson::mpi_communicator

Definition at line 427 of file module_mp_thompson.F90.

◆ n_avo

real(dp) module_mp_thompson::n_avo

Definition at line 209 of file module_mp_thompson.F90.

◆ naccn0

real(wp), parameter module_mp_thompson::naccn0 = 300.0e6

Definition at line 109 of file module_mp_thompson.F90.

◆ naccn1

real(wp), parameter module_mp_thompson::naccn1 = 50.0e6

Definition at line 110 of file module_mp_thompson.F90.

◆ nain0

real(wp), parameter module_mp_thompson::nain0 = 1.5e6

Definition at line 107 of file module_mp_thompson.F90.

◆ nain1

real(wp), parameter module_mp_thompson::nain1 = 0.5e6

Definition at line 108 of file module_mp_thompson.F90.

◆ nt_c_l

real(wp) module_mp_thompson::nt_c_l = 150.e6

Definition at line 96 of file module_mp_thompson.F90.

◆ nt_c_o

real(wp) module_mp_thompson::nt_c_o = 50.e6

Definition at line 97 of file module_mp_thompson.F90.

◆ nt_i_max

real(wp) module_mp_thompson::nt_i_max = 4999.e3_dp

Definition at line 101 of file module_mp_thompson.F90.

◆ pi

real(wp) module_mp_thompson::pi

Definition at line 77 of file module_mp_thompson.F90.

◆ qr_acr_qg_file

character(len=*), parameter module_mp_thompson::qr_acr_qg_file = 'qr_acr_qgV2.dat'

Definition at line 383 of file module_mp_thompson.F90.

◆ qr_acr_qs_file

character(len=*), parameter module_mp_thompson::qr_acr_qs_file = 'qr_acr_qsV2.dat'

Definition at line 384 of file module_mp_thompson.F90.

◆ r

real(wp) module_mp_thompson::r

Definition at line 201 of file module_mp_thompson.F90.

◆ r_uni

real(wp) module_mp_thompson::r_uni

Definition at line 204 of file module_mp_thompson.F90.

◆ re_qc_max

real(wp), parameter module_mp_thompson::re_qc_max = 50.0e-6

Definition at line 231 of file module_mp_thompson.F90.

◆ re_qc_min

real(wp), parameter module_mp_thompson::re_qc_min = 2.50e-6

Definition at line 230 of file module_mp_thompson.F90.

◆ re_qi_max

real(wp), parameter module_mp_thompson::re_qi_max = 125.0e-6

Definition at line 233 of file module_mp_thompson.F90.

◆ re_qi_min

real(wp), parameter module_mp_thompson::re_qi_min = 2.50e-6

Definition at line 232 of file module_mp_thompson.F90.

◆ re_qs_max

real(wp), parameter module_mp_thompson::re_qs_max = 999.0e-6

Definition at line 235 of file module_mp_thompson.F90.

◆ re_qs_min

real(wp), parameter module_mp_thompson::re_qs_min = 5.00e-6

Definition at line 234 of file module_mp_thompson.F90.

◆ rho_not

real(wp) module_mp_thompson::rho_not

Definition at line 189 of file module_mp_thompson.F90.

◆ roverrv

real(wp) module_mp_thompson::roverrv

Definition at line 202 of file module_mp_thompson.F90.

◆ rr_min

real(wp) module_mp_thompson::rr_min = 1000.0

Definition at line 102 of file module_mp_thompson.F90.

◆ rv

real(wp) module_mp_thompson::rv

Definition at line 199 of file module_mp_thompson.F90.

◆ ssati_min

real(wp) module_mp_thompson::ssati_min = 0.15

Definition at line 100 of file module_mp_thompson.F90.

◆ t1_qg_me

real(wp) module_mp_thompson::t1_qg_me

Definition at line 424 of file module_mp_thompson.F90.

◆ t1_qg_qc

real(wp) module_mp_thompson::t1_qg_qc

Definition at line 421 of file module_mp_thompson.F90.

◆ t1_qg_sd

real(wp) module_mp_thompson::t1_qg_sd

Definition at line 423 of file module_mp_thompson.F90.

◆ t1_qr_ev

real(wp) module_mp_thompson::t1_qr_ev

Definition at line 422 of file module_mp_thompson.F90.

◆ t1_qr_qc

real(wp) module_mp_thompson::t1_qr_qc

Definition at line 421 of file module_mp_thompson.F90.

◆ t1_qr_qi

real(wp) module_mp_thompson::t1_qr_qi

Definition at line 421 of file module_mp_thompson.F90.

◆ t1_qs_me

real(wp) module_mp_thompson::t1_qs_me

Definition at line 424 of file module_mp_thompson.F90.

◆ t1_qs_qc

real(wp) module_mp_thompson::t1_qs_qc

Definition at line 421 of file module_mp_thompson.F90.

◆ t1_qs_qi

real(wp) module_mp_thompson::t1_qs_qi

Definition at line 421 of file module_mp_thompson.F90.

◆ t1_qs_sd

real(wp) module_mp_thompson::t1_qs_sd

Definition at line 423 of file module_mp_thompson.F90.

◆ t2_qg_me

real(wp) module_mp_thompson::t2_qg_me

Definition at line 424 of file module_mp_thompson.F90.

◆ t2_qg_sd

real(wp) module_mp_thompson::t2_qg_sd

Definition at line 423 of file module_mp_thompson.F90.

◆ t2_qr_ev

real(wp) module_mp_thompson::t2_qr_ev

Definition at line 422 of file module_mp_thompson.F90.

◆ t2_qr_qi

real(wp) module_mp_thompson::t2_qr_qi

Definition at line 421 of file module_mp_thompson.F90.

◆ t2_qs_me

real(wp) module_mp_thompson::t2_qs_me

Definition at line 424 of file module_mp_thompson.F90.

◆ t2_qs_sd

real(wp) module_mp_thompson::t2_qs_sd

Definition at line 423 of file module_mp_thompson.F90.

◆ t_0

real(wp) module_mp_thompson::t_0

Definition at line 76 of file module_mp_thompson.F90.

◆ t_efrw

real (dp), dimension(:,:), allocatable module_mp_thompson::t_efrw

Definition at line 400 of file module_mp_thompson.F90.

◆ t_efsw

real (dp), dimension(:,:), allocatable module_mp_thompson::t_efsw

Definition at line 401 of file module_mp_thompson.F90.

◆ t_nc

real(dp), dimension(nbc) module_mp_thompson::t_nc

Definition at line 268 of file module_mp_thompson.F90.

◆ tcg_racg

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tcg_racg

Definition at line 387 of file module_mp_thompson.F90.

◆ tcr_gacr

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tcr_gacr

Definition at line 387 of file module_mp_thompson.F90.

◆ tcr_sacr1

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tcr_sacr1

Definition at line 390 of file module_mp_thompson.F90.

◆ tcr_sacr2

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tcr_sacr2

Definition at line 390 of file module_mp_thompson.F90.

◆ tcs_racs1

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tcs_racs1

Definition at line 390 of file module_mp_thompson.F90.

◆ tcs_racs2

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tcs_racs2

Definition at line 390 of file module_mp_thompson.F90.

◆ thomp_table_file

character(len=*), parameter module_mp_thompson::thomp_table_file = 'thompson_tables_precomp_v2.sl'

Definition at line 382 of file module_mp_thompson.F90.

◆ thompson_table_writer

logical module_mp_thompson::thompson_table_writer

Definition at line 430 of file module_mp_thompson.F90.

◆ tmg_gacr

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tmg_gacr

Definition at line 387 of file module_mp_thompson.F90.

◆ tmr_racg

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tmr_racg

Definition at line 387 of file module_mp_thompson.F90.

◆ tmr_racs1

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tmr_racs1

Definition at line 390 of file module_mp_thompson.F90.

◆ tmr_racs2

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tmr_racs2

Definition at line 390 of file module_mp_thompson.F90.

◆ tms_sacr1

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tms_sacr1

Definition at line 390 of file module_mp_thompson.F90.

◆ tms_sacr2

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tms_sacr2

Definition at line 390 of file module_mp_thompson.F90.

◆ tnc_wev

real (dp), dimension(:,:,:), allocatable module_mp_thompson::tnc_wev

Definition at line 403 of file module_mp_thompson.F90.

◆ tnccn_act

real (sp), dimension(:,:,:,:,:), allocatable module_mp_thompson::tnccn_act

Definition at line 405 of file module_mp_thompson.F90.

◆ tni_iaus

real (dp), dimension(:,:), allocatable module_mp_thompson::tni_iaus

Definition at line 398 of file module_mp_thompson.F90.

◆ tni_qcfz

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tni_qcfz

Definition at line 394 of file module_mp_thompson.F90.

◆ tni_qrfz

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tni_qrfz

Definition at line 396 of file module_mp_thompson.F90.

◆ tnr_gacr

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tnr_gacr

Definition at line 387 of file module_mp_thompson.F90.

◆ tnr_qrfz

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tnr_qrfz

Definition at line 396 of file module_mp_thompson.F90.

◆ tnr_racg

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tnr_racg

Definition at line 387 of file module_mp_thompson.F90.

◆ tnr_racs1

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tnr_racs1

Definition at line 390 of file module_mp_thompson.F90.

◆ tnr_racs2

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tnr_racs2

Definition at line 390 of file module_mp_thompson.F90.

◆ tnr_rev

real (dp), dimension(:,:,:), allocatable module_mp_thompson::tnr_rev

Definition at line 402 of file module_mp_thompson.F90.

◆ tnr_sacr1

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tnr_sacr1

Definition at line 390 of file module_mp_thompson.F90.

◆ tnr_sacr2

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tnr_sacr2

Definition at line 390 of file module_mp_thompson.F90.

◆ tpc_wev

real (dp), dimension(:,:,:), allocatable module_mp_thompson::tpc_wev

Definition at line 403 of file module_mp_thompson.F90.

◆ tpg_qrfz

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tpg_qrfz

Definition at line 396 of file module_mp_thompson.F90.

◆ tpi_ide

real (dp), dimension(:,:), allocatable module_mp_thompson::tpi_ide

Definition at line 398 of file module_mp_thompson.F90.

◆ tpi_qcfz

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tpi_qcfz

Definition at line 394 of file module_mp_thompson.F90.

◆ tpi_qrfz

real (dp), dimension(:,:,:,:), allocatable module_mp_thompson::tpi_qrfz

Definition at line 396 of file module_mp_thompson.F90.

◆ tps_iaus

real (dp), dimension(:,:), allocatable module_mp_thompson::tps_iaus

Definition at line 398 of file module_mp_thompson.F90.

◆ xdx

real(dp), dimension(nbins+1) module_mp_thompson::xdx

Definition at line 262 of file module_mp_thompson.F90.

◆ xnc_max

real(wp) module_mp_thompson::xnc_max = 1000.e3

Definition at line 99 of file module_mp_thompson.F90.

Functions/Subroutines

Variables

Detailed Description

Function/Subroutine Documentation

◆ semi_lagrange_sedim()

Variable Documentation

◆ av_i

◆ av_s

◆ bv_i

◆ bv_s

◆ cp

◆ d0s

◆ dc

◆ dg

◆ di

◆ dr

◆ ds

◆ dtc

◆ dtg

◆ dti

◆ dtr

◆ dts

◆ eps

◆ freeze_h2o_file

◆ k_b

◆ lfus

◆ lvap0

◆ m_a

◆ m_w

◆ mpi_communicator

◆ n_avo

◆ naccn0

◆ naccn1

◆ nain0

◆ nain1

◆ nt_c_l

◆ nt_c_o

◆ nt_i_max

◆ pi

◆ qr_acr_qg_file

◆ qr_acr_qs_file

◆ r

◆ r_uni

◆ re_qc_max

◆ re_qc_min

◆ re_qi_max

◆ re_qi_min

◆ re_qs_max

◆ re_qs_min

◆ rho_not

◆ roverrv

◆ rr_min

◆ rv

◆ ssati_min

◆ t1_qg_me

◆ t1_qg_qc

◆ t1_qg_sd

◆ t1_qr_ev

◆ t1_qr_qc

◆ t1_qr_qi

◆ t1_qs_me

◆ t1_qs_qc

◆ t1_qs_qi

◆ t1_qs_sd

◆ t2_qg_me

◆ t2_qg_sd

◆ t2_qr_ev

◆ t2_qr_qi

◆ t2_qs_me

◆ t2_qs_sd

◆ t_0

◆ t_efrw

◆ t_efsw

◆ t_nc

◆ tcg_racg

◆ tcr_gacr

◆ tcr_sacr1

◆ tcr_sacr2

◆ tcs_racs1

◆ tcs_racs2