module_mp_thompson.F90 File Reference

This file contains the entity of GSD Thompson MP scheme. More...

Modules
module	module_mp_thompson

Functions/Subroutines
subroutine	module_mp_thompson::thompson_finalize ()
subroutine	module_mp_thompson::qr_acr_qg
	Rain collecting graupel (and inverse). Explicit CE integration. More...
subroutine	module_mp_thompson::qr_acr_qs
	Rain collecting snow (and inverse). Explicit CE integration. More...
subroutine	module_mp_thompson::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. More...
subroutine	module_mp_thompson::qi_aut_qs
	Cloud ice converting to snow since portion greater than min snow size. Given cloud ice content (kg/m**3), number concentration (#/m**3) 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). More...
subroutine	module_mp_thompson::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. More...
subroutine	module_mp_thompson::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.". More...
real function	module_mp_thompson::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). More...
subroutine	module_mp_thompson::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. More...
subroutine	module_mp_thompson::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. More...
real function	module_mp_thompson::activ_ncloud (Tt, Ww, NCCN)
	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. More...
subroutine	module_mp_thompson::gcf (GAMMCF, A, X, GLN)
	Returns the incomplete gamma function q(a,x) evaluated by its continued fraction representation as gammcf. More...
subroutine	module_mp_thompson::gser (GAMSER, A, X, GLN)
	Returns the incomplete gamma function p(a,x) evaluated by its series representation as gamser. More...
real function	module_mp_thompson::gammln (XX)
	Returns the value ln(gamma(xx)) for xx > 0. More...
real function	module_mp_thompson::gammp (A, X)
real function	module_mp_thompson::wgamma (y)
real function	module_mp_thompson::rslf (P, T)
	THIS FUNCTION CALCULATES THE LIQUID SATURATION VAPOR MIXING RATIO AS A FUNCTION OF TEMPERATURE AND PRESSURE. More...
real function	module_mp_thompson::rsif (P, T)
	THIS FUNCTION CALCULATES THE ICE SATURATION VAPOR MIXING RATIO AS A FUNCTION OF TEMPERATURE AND PRESSURE. More...
real function	module_mp_thompson::icedemott (tempc, qv, qvs, qvsi, rho, nifa)
real function	module_mp_thompson::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. More...
real function	module_mp_thompson::delta_p (yy, y1, y2, aa, bb)
	Helper routine for Phillips et al (2008) ice nucleation. Trude. More...
subroutine	module_mp_thompson::calc_effectrad (t1d, p1d, qv1d, qc1d, nc1d, qi1d, ni1d, qs1d, re_qc1d, re_qi1d, re_qs1d, 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. More...
subroutine	module_mp_thompson::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. More...
subroutine	module_mp_thompson::readwrite_tables (mode, mpicomm, mpirank, mpiroot, ierr)
real *8 function	calculate_checksum ()
subroutine	module_mp_thompson::thompson_init (nwfa2d, nifa2d, nwfa, nifa, mpicomm, mpirank, mpiroot, threads, errmsg, errflg)
	This subroutine calculates simplified cloud species equations and create lookup tables in Thomspson scheme. More...
subroutine	module_mp_thompson::mp_gt_driver (WRF_CHEM )
	This is a wrapper routine designed to transfer values from 3D to 1D. More...
subroutine	module_mp_thompson::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)[172] [173]. More...

Variables
logical, parameter, private	module_mp_thompson::iiwarm = .false.
logical, private	module_mp_thompson::is_aerosol_aware = .false.
logical, parameter, private	module_mp_thompson::dustyice = .true.
logical, parameter, private	module_mp_thompson::homogice = .true.
integer, parameter, private	module_mp_thompson::ifdry = 0
real, parameter, private	module_mp_thompson::t_0 = 273.15
real, parameter, private	module_mp_thompson::pi = 3.1415926536
real, parameter, private	module_mp_thompson::rho_w = 1000.0
real, parameter, private	module_mp_thompson::rho_s = 100.0
real, parameter, private	module_mp_thompson::rho_g = 500.0
real, parameter, private	module_mp_thompson::rho_i = 890.0
real, parameter	module_mp_thompson::nt_c = 100.E6
real, parameter, private	module_mp_thompson::nt_c_max = 1999.E6
real, parameter	module_mp_thompson::nain0 = 1.5E6
real, parameter	module_mp_thompson::nain1 = 0.5E6
real, parameter	module_mp_thompson::naccn0 = 300.0E6
real, parameter	module_mp_thompson::naccn1 = 50.0E6
real, parameter, private	module_mp_thompson::mu_r = 0.0
real, parameter, private	module_mp_thompson::mu_g = 0.0
real, parameter, private	module_mp_thompson::mu_i = 0.0
real, private	module_mp_thompson::mu_c
real, parameter, private	module_mp_thompson::mu_s = 0.6357
real, parameter, private	module_mp_thompson::kap0 = 490.6
real, parameter, private	module_mp_thompson::kap1 = 17.46
real, parameter, private	module_mp_thompson::lam0 = 20.78
real, parameter, private	module_mp_thompson::lam1 = 3.29
real, parameter, private	module_mp_thompson::gonv_min = 1.E4
real, parameter, private	module_mp_thompson::gonv_max = 3.E6
real, parameter, private	module_mp_thompson::am_r = PI*rho_w/6.0
real, parameter, private	module_mp_thompson::bm_r = 3.0
real, parameter, private	module_mp_thompson::am_s = 0.069
real, parameter, private	module_mp_thompson::bm_s = 2.0
real, parameter, private	module_mp_thompson::am_g = PI*rho_g/6.0
real, parameter, private	module_mp_thompson::bm_g = 3.0
real, parameter, private	module_mp_thompson::am_i = PI*rho_i/6.0
real, parameter, private	module_mp_thompson::bm_i = 3.0
real, parameter, private	module_mp_thompson::av_r = 4854.0
real, parameter, private	module_mp_thompson::bv_r = 1.0
real, parameter, private	module_mp_thompson::fv_r = 195.0
real, parameter, private	module_mp_thompson::av_s = 40.0
real, parameter, private	module_mp_thompson::bv_s = 0.55
real, parameter, private	module_mp_thompson::fv_s = 100.0
real, parameter, private	module_mp_thompson::av_g = 442.0
real, parameter, private	module_mp_thompson::bv_g = 0.89
real, parameter, private	module_mp_thompson::av_i = 1847.5
real, parameter, private	module_mp_thompson::bv_i = 1.0
real, parameter, private	module_mp_thompson::av_c = 0.316946E8
real, parameter, private	module_mp_thompson::bv_c = 2.0
real, parameter, private	module_mp_thompson::c_cube = 0.5
real, parameter, private	module_mp_thompson::c_sqrd = 0.15
real, parameter, private	module_mp_thompson::ef_si = 0.05
real, parameter, private	module_mp_thompson::ef_rs = 0.95
real, parameter, private	module_mp_thompson::ef_rg = 0.75
real, parameter, private	module_mp_thompson::ef_ri = 0.95
real, parameter, private	module_mp_thompson::r1 = 1.E-12
real, parameter, private	module_mp_thompson::r2 = 1.E-6
real, parameter	module_mp_thompson::eps = 1.E-15
real, parameter, private	module_mp_thompson::tno = 5.0
real, parameter, private	module_mp_thompson::ato = 0.304
real, parameter, private	module_mp_thompson::rho_not = 101325.0/(287.05*298.0)
real, parameter, private	module_mp_thompson::sc = 0.632
real, private	module_mp_thompson::sc3
real, parameter, private	module_mp_thompson::hgfr = 235.16
	Homogeneous freezing temperature. More...
real, parameter, private	module_mp_thompson::rv = 461.5
real, parameter, private	module_mp_thompson::orv = 1./Rv
real, parameter, private	module_mp_thompson::r = 287.04
real, parameter, private	module_mp_thompson::cp = 1004.0
real, parameter, private	module_mp_thompson::r_uni = 8.314
	J (mol K)-1. More...
double precision, parameter, private	module_mp_thompson::k_b = 1.38065E-23
	Boltzmann constant [J/K]. More...
double precision, parameter, private	module_mp_thompson::m_w = 18.01528E-3
	molecular mass of water [kg/mol] More...
double precision, parameter, private	module_mp_thompson::m_a = 28.96E-3
	molecular mass of air [kg/mol] More...
double precision, parameter, private	module_mp_thompson::n_avo = 6.022E23
	Avogadro number [1/mol]. More...
double precision, parameter, private	module_mp_thompson::ma_w = M_w / N_avo
	mass of water molecule [kg] More...
real, parameter, private	module_mp_thompson::ar_volume = 4./3.*PI*(2.5e-6)**3
	assume radius of 0.025 micrometer, 2.5e-6 cm More...
real, parameter, private	module_mp_thompson::lsub = 2.834E6
real, parameter, private	module_mp_thompson::lvap0 = 2.5E6
real, parameter, private	module_mp_thompson::lfus = lsub - lvap0
real, parameter, private	module_mp_thompson::olfus = 1./lfus
real, parameter, private	module_mp_thompson::xm0i = 1.E-12
real, parameter, private	module_mp_thompson::d0c = 1.E-6
real, parameter, private	module_mp_thompson::d0r = 50.E-6
real, parameter, private	module_mp_thompson::d0s = 200.E-6
real, parameter, private	module_mp_thompson::d0g = 250.E-6
real, private	module_mp_thompson::d0i
real, private	module_mp_thompson::xm0s
real, private	module_mp_thompson::xm0g
integer, parameter, private	module_mp_thompson::nbins = 100
integer, parameter, private	module_mp_thompson::nbc = nbins
integer, parameter, private	module_mp_thompson::nbi = nbins
integer, parameter, private	module_mp_thompson::nbr = nbins
integer, parameter, private	module_mp_thompson::nbs = nbins
integer, parameter, private	module_mp_thompson::nbg = nbins
integer, parameter, private	module_mp_thompson::ntb_c = 37
integer, parameter, private	module_mp_thompson::ntb_i = 64
integer, parameter, private	module_mp_thompson::ntb_r = 37
integer, parameter, private	module_mp_thompson::ntb_s = 28
integer, parameter, private	module_mp_thompson::ntb_g = 28
integer, parameter, private	module_mp_thompson::ntb_g1 = 28
integer, parameter, private	module_mp_thompson::ntb_r1 = 37
integer, parameter, private	module_mp_thompson::ntb_i1 = 55
integer, parameter, private	module_mp_thompson::ntb_t = 9
integer, private	module_mp_thompson::nic1
integer, private	module_mp_thompson::nic2
integer, private	module_mp_thompson::nii2
integer, private	module_mp_thompson::nii3
integer, private	module_mp_thompson::nir2
integer, private	module_mp_thompson::nir3
integer, private	module_mp_thompson::nis2
integer, private	module_mp_thompson::nig2
integer, private	module_mp_thompson::nig3
integer, parameter, private	module_mp_thompson::ntb_arc = 7
integer, parameter, private	module_mp_thompson::ntb_arw = 9
integer, parameter, private	module_mp_thompson::ntb_art = 7
integer, parameter, private	module_mp_thompson::ntb_arr = 5
integer, parameter, private	module_mp_thompson::ntb_ark = 4
integer, parameter, private	module_mp_thompson::ntb_in = 55
integer, private	module_mp_thompson::niin2
double precision, dimension(nbins+1)	module_mp_thompson::xdx
double precision, dimension(nbc)	module_mp_thompson::dc
double precision, dimension(nbc)	module_mp_thompson::dtc
double precision, dimension(nbi)	module_mp_thompson::di
double precision, dimension(nbi)	module_mp_thompson::dti
double precision, dimension(nbr)	module_mp_thompson::dr
double precision, dimension(nbr)	module_mp_thompson::dtr
double precision, dimension(nbs)	module_mp_thompson::ds
double precision, dimension(nbs)	module_mp_thompson::dts
double precision, dimension(nbg)	module_mp_thompson::dg
double precision, dimension(nbg)	module_mp_thompson::dtg
double precision, dimension(nbc)	module_mp_thompson::t_nc
real, dimension(ntb_c), parameter, private	module_mp_thompson::r_c = (/1.e-6,2.e-6,3.e-6,4.e-6,5.e-6,6.e-6,7.e-6,8.e-6,9.e-6, 1.e-5,2.e-5,3.e-5,4.e-5,5.e-5,6.e-5,7.e-5,8.e-5,9.e-5, 1.e-4,2.e-4,3.e-4,4.e-4,5.e-4,6.e-4,7.e-4,8.e-4,9.e-4, 1.e-3,2.e-3,3.e-3,4.e-3,5.e-3,6.e-3,7.e-3,8.e-3,9.e-3, 1.e-2/)
	Lookup tables for cloud water content (kg/m**3). More...
real, dimension(ntb_i), parameter, private	module_mp_thompson::r_i = (/1.e-10,2.e-10,3.e-10,4.e-10, 5.e-10,6.e-10,7.e-10,8.e-10,9.e-10, 1.e-9,2.e-9,3.e-9,4.e-9,5.e-9,6.e-9,7.e-9,8.e-9,9.e-9, 1.e-8,2.e-8,3.e-8,4.e-8,5.e-8,6.e-8,7.e-8,8.e-8,9.e-8, 1.e-7,2.e-7,3.e-7,4.e-7,5.e-7,6.e-7,7.e-7,8.e-7,9.e-7, 1.e-6,2.e-6,3.e-6,4.e-6,5.e-6,6.e-6,7.e-6,8.e-6,9.e-6, 1.e-5,2.e-5,3.e-5,4.e-5,5.e-5,6.e-5,7.e-5,8.e-5,9.e-5, 1.e-4,2.e-4,3.e-4,4.e-4,5.e-4,6.e-4,7.e-4,8.e-4,9.e-4, 1.e-3/)
	Lookup tables for cloud ice content (kg/m**3). More...
real, dimension(ntb_r), parameter, private	module_mp_thompson::r_r = (/1.e-6,2.e-6,3.e-6,4.e-6,5.e-6,6.e-6,7.e-6,8.e-6,9.e-6, 1.e-5,2.e-5,3.e-5,4.e-5,5.e-5,6.e-5,7.e-5,8.e-5,9.e-5, 1.e-4,2.e-4,3.e-4,4.e-4,5.e-4,6.e-4,7.e-4,8.e-4,9.e-4, 1.e-3,2.e-3,3.e-3,4.e-3,5.e-3,6.e-3,7.e-3,8.e-3,9.e-3, 1.e-2/)
	Lookup tables for rain content (kg/m**3). More...
real, dimension(ntb_g), parameter, private	module_mp_thompson::r_g = (/1.e-5,2.e-5,3.e-5,4.e-5,5.e-5,6.e-5,7.e-5,8.e-5,9.e-5, 1.e-4,2.e-4,3.e-4,4.e-4,5.e-4,6.e-4,7.e-4,8.e-4,9.e-4, 1.e-3,2.e-3,3.e-3,4.e-3,5.e-3,6.e-3,7.e-3,8.e-3,9.e-3, 1.e-2/)
	Lookup tables for graupel content (kg/m**3). More...
real, dimension(ntb_s), parameter, private	module_mp_thompson::r_s = (/1.e-5,2.e-5,3.e-5,4.e-5,5.e-5,6.e-5,7.e-5,8.e-5,9.e-5, 1.e-4,2.e-4,3.e-4,4.e-4,5.e-4,6.e-4,7.e-4,8.e-4,9.e-4, 1.e-3,2.e-3,3.e-3,4.e-3,5.e-3,6.e-3,7.e-3,8.e-3,9.e-3, 1.e-2/)
	Lookup tables for snow content (kg/m**3). More...
real, dimension(ntb_r1), parameter, private	module_mp_thompson::n0r_exp = (/1.e6,2.e6,3.e6,4.e6,5.e6,6.e6,7.e6,8.e6,9.e6, 1.e7,2.e7,3.e7,4.e7,5.e7,6.e7,7.e7,8.e7,9.e7, 1.e8,2.e8,3.e8,4.e8,5.e8,6.e8,7.e8,8.e8,9.e8, 1.e9,2.e9,3.e9,4.e9,5.e9,6.e9,7.e9,8.e9,9.e9, 1.e10/)
	Lookup tables for rain y-intercept parameter (/m**4). More...
real, dimension(ntb_g1), parameter, private	module_mp_thompson::n0g_exp = (/1.e4,2.e4,3.e4,4.e4,5.e4,6.e4,7.e4,8.e4,9.e4, 1.e5,2.e5,3.e5,4.e5,5.e5,6.e5,7.e5,8.e5,9.e5, 1.e6,2.e6,3.e6,4.e6,5.e6,6.e6,7.e6,8.e6,9.e6, 1.e7/)
	Lookup tables for graupel y-intercept parameter (/m**4). More...
real, dimension(ntb_i1), parameter, private	module_mp_thompson::nt_i = (/1.0,2.0,3.0,4.0,5.0,6.0,7.0,8.0,9.0, 1.e1,2.e1,3.e1,4.e1,5.e1,6.e1,7.e1,8.e1,9.e1, 1.e2,2.e2,3.e2,4.e2,5.e2,6.e2,7.e2,8.e2,9.e2, 1.e3,2.e3,3.e3,4.e3,5.e3,6.e3,7.e3,8.e3,9.e3, 1.e4,2.e4,3.e4,4.e4,5.e4,6.e4,7.e4,8.e4,9.e4, 1.e5,2.e5,3.e5,4.e5,5.e5,6.e5,7.e5,8.e5,9.e5, 1.e6/)
	Lookup tables for ice number concentration (/m**3). More...
real, dimension(ntb_arc), parameter, private	module_mp_thompson::ta_na = (/10.0, 31.6, 100.0, 316.0, 1000.0, 3160.0, 10000.0/)
real, dimension(ntb_arw), parameter, private	module_mp_thompson::ta_ww = (/0.01, 0.0316, 0.1, 0.316, 1.0, 3.16, 10.0, 31.6, 100.0/)
real, dimension(ntb_art), parameter, private	module_mp_thompson::ta_tk = (/243.15, 253.15, 263.15, 273.15, 283.15, 293.15, 303.15/)
real, dimension(ntb_arr), parameter, private	module_mp_thompson::ta_ra = (/0.01, 0.02, 0.04, 0.08, 0.16/)
real, dimension(ntb_ark), parameter, private	module_mp_thompson::ta_ka = (/0.2, 0.4, 0.6, 0.8/)
real, dimension(ntb_in), parameter, private	module_mp_thompson::nt_in = (/1.0,2.0,3.0,4.0,5.0,6.0,7.0,8.0,9.0, 1.e1,2.e1,3.e1,4.e1,5.e1,6.e1,7.e1,8.e1,9.e1, 1.e2,2.e2,3.e2,4.e2,5.e2,6.e2,7.e2,8.e2,9.e2, 1.e3,2.e3,3.e3,4.e3,5.e3,6.e3,7.e3,8.e3,9.e3, 1.e4,2.e4,3.e4,4.e4,5.e4,6.e4,7.e4,8.e4,9.e4, 1.e5,2.e5,3.e5,4.e5,5.e5,6.e5,7.e5,8.e5,9.e5, 1.e6/)
	Lookup tables for IN concentration (/m**3) from 0.001 to 1000/Liter. More...
real, dimension(10), parameter, private	module_mp_thompson::sa = (/ 5.065339, -0.062659, -3.032362, 0.029469, -0.000285, 0.31255, 0.000204, 0.003199, 0.0, -0.015952/)
	For snow moments conversions (from Field et al. 2005) More...
real, dimension(10), parameter, private	module_mp_thompson::sb = (/ 0.476221, -0.015896, 0.165977, 0.007468, -0.000141, 0.060366, 0.000079, 0.000594, 0.0, -0.003577/)
real, dimension(ntb_t), parameter, private	module_mp_thompson::tc = (/-0.01, -5., -10., -15., -20., -25., -30., -35., -40./)
	Temperatures (5 C interval 0 to -40) used in lookup tables. More...
integer, parameter, private	module_mp_thompson::r8size = 8
integer, parameter, private	module_mp_thompson::r4size = 4
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tcg_racg
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tmr_racg
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tcr_gacr
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tmg_gacr
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tnr_racg
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tnr_gacr
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tcs_racs1
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tmr_racs1
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tcs_racs2
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tmr_racs2
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tcr_sacr1
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tms_sacr1
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tcr_sacr2
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tms_sacr2
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tnr_racs1
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tnr_racs2
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tnr_sacr1
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tnr_sacr2
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tpi_qcfz
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tni_qcfz
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tpi_qrfz
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tpg_qrfz
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tni_qrfz
real(kind=r8size), dimension(:,:,:,:), allocatable	module_mp_thompson::tnr_qrfz
real(kind=r8size), dimension(:,:), allocatable	module_mp_thompson::tps_iaus
real(kind=r8size), dimension(:,:), allocatable	module_mp_thompson::tni_iaus
real(kind=r8size), dimension(:,:), allocatable	module_mp_thompson::tpi_ide
real(kind=r8size), dimension(:,:), allocatable	module_mp_thompson::t_efrw
real(kind=r8size), dimension(:,:), allocatable	module_mp_thompson::t_efsw
real(kind=r8size), dimension(:,:,:), allocatable	module_mp_thompson::tnr_rev
real(kind=r8size), dimension(:,:,:), allocatable	module_mp_thompson::tpc_wev
real(kind=r8size), dimension(:,:,:), allocatable	module_mp_thompson::tnc_wev
real(kind=r4size), dimension(:,:,:,:,:), allocatable	module_mp_thompson::tnccn_act
real, dimension(5, 15), private	module_mp_thompson::cce
real, dimension(5, 15), private	module_mp_thompson::ccg
real, dimension(15), private	module_mp_thompson::ocg1
real, dimension(15), private	module_mp_thompson::ocg2
real, dimension(7), private	module_mp_thompson::cie
real, dimension(7), private	module_mp_thompson::cig
real, private	module_mp_thompson::oig1
real, private	module_mp_thompson::oig2
real, private	module_mp_thompson::obmi
real, dimension(13), private	module_mp_thompson::cre
real, dimension(13), private	module_mp_thompson::crg
real, private	module_mp_thompson::ore1
real, private	module_mp_thompson::org1
real, private	module_mp_thompson::org2
real, private	module_mp_thompson::org3
real, private	module_mp_thompson::obmr
real, dimension(18), private	module_mp_thompson::cse
real, dimension(18), private	module_mp_thompson::csg
real, private	module_mp_thompson::oams
real, private	module_mp_thompson::obms
real, private	module_mp_thompson::ocms
real, dimension(12), private	module_mp_thompson::cge
real, dimension(12), private	module_mp_thompson::cgg
real, private	module_mp_thompson::oge1
real, private	module_mp_thompson::ogg1
real, private	module_mp_thompson::ogg2
real, private	module_mp_thompson::ogg3
real, private	module_mp_thompson::oamg
real, private	module_mp_thompson::obmg
real, private	module_mp_thompson::ocmg
real	module_mp_thompson::t1_qr_qc
real	module_mp_thompson::t1_qr_qi
real	module_mp_thompson::t2_qr_qi
real	module_mp_thompson::t1_qg_qc
real	module_mp_thompson::t1_qs_qc
real	module_mp_thompson::t1_qs_qi
real	module_mp_thompson::t1_qr_ev
real	module_mp_thompson::t2_qr_ev
real	module_mp_thompson::t1_qs_sd
real	module_mp_thompson::t2_qs_sd
real	module_mp_thompson::t1_qg_sd
real	module_mp_thompson::t2_qg_sd
real	module_mp_thompson::t1_qs_me
real	module_mp_thompson::t2_qs_me
real	module_mp_thompson::t1_qg_me
real	module_mp_thompson::t2_qg_me
integer	module_mp_thompson::mpi_communicator
logical	module_mp_thompson::thompson_table_writer