subroutine, public mp_thompson::mp_thompson_run	(	integer, intent(in)	ncol,
		integer, intent(in)	nlev,
		real(kind_phys), intent(in)	con_g,
		real(kind_phys), intent(in)	con_rd,
		real(kind_phys), dimension(1:ncol,1:nlev), intent(inout)	spechum,
		real(kind_phys), dimension(1:ncol,1:nlev), intent(inout)	qc,
		real(kind_phys), dimension(1:ncol,1:nlev), intent(inout)	qr,
		real(kind_phys), dimension(1:ncol,1:nlev), intent(inout)	qi,
		real(kind_phys), dimension(1:ncol,1:nlev), intent(inout)	qs,
		real(kind_phys), dimension(1:ncol,1:nlev), intent(inout)	qg,
		real(kind_phys), dimension(1:ncol,1:nlev), intent(inout)	ni,
		real(kind_phys), dimension(1:ncol,1:nlev), intent(inout)	nr,
		logical, intent(in)	is_aerosol_aware,
		real(kind_phys), dimension(:,:), intent(inout), optional	nc,
		real(kind_phys), dimension(:,:), intent(inout), optional	nwfa,
		real(kind_phys), dimension(:,:), intent(inout), optional	nifa,
		real(kind_phys), dimension(:), intent(in), optional	nwfa2d,
		real(kind_phys), dimension(:), intent(in), optional	nifa2d,
		real(kind_phys), dimension(1:ncol,1:nlev), intent(inout)	tgrs,
		real(kind_phys), dimension(1:ncol,1:nlev), intent(in)	prsl,
		real(kind_phys), dimension(1:ncol,1:nlev+1), intent(in)	phii,
		real(kind_phys), dimension(1:ncol,1:nlev), intent(in)	omega,
		real(kind_phys), intent(in)	dtp,
		real(kind_phys), dimension(1:ncol), intent(out)	prcp,
		real(kind_phys), dimension(1:ncol), intent(out)	rain,
		real(kind_phys), dimension(1:ncol), intent(out)	graupel,
		real(kind_phys), dimension(1:ncol), intent(out)	ice,
		real(kind_phys), dimension(1:ncol), intent(out)	snow,
		real(kind_phys), dimension(1:ncol), intent(out)	sr,
		real(kind_phys), dimension(1:ncol,1:nlev), intent(out)	refl_10cm,
		logical, intent(in)	reset,
		logical, intent(in), optional	do_radar_ref,
		real(kind_phys), dimension(1:ncol,1:nlev), intent(out), optional	re_cloud,
		real(kind_phys), dimension(1:ncol,1:nlev), intent(out), optional	re_ice,
		real(kind_phys), dimension(1:ncol,1:nlev), intent(out), optional	re_snow,
		integer, intent(in)	mpicomm,
		integer, intent(in)	mpirank,
		integer, intent(in)	mpiroot,
		character(len=*), intent(out)	errmsg,
		integer, intent(out)	errflg
	)

Argument Table

local_name	standard_name	long_name	units	type	dimensions	kind	intent	optional
ncol	horizontal_loop_extent	horizontal loop extent	count	integer	()		in	False
nlev	vertical_dimension	number of vertical levels	count	integer	()		in	False
con_g	gravitational_acceleration	gravitational acceleration	m s-2	real	()	kind_phys	in	False
con_rd	gas_constant_dry_air	ideal gas constant for dry air	J kg-1 K-1	real	()	kind_phys	in	False
spechum	water_vapor_specific_humidity_updated_by_physics	water vapor specific humidity	kg kg-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	inout	False
qc	cloud_condensed_water_mixing_ratio_updated_by_physics	cloud water mixing ratio wrt dry+vapor (no condensates)	kg kg-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	inout	False
qr	rain_water_mixing_ratio_updated_by_physics	rain water mixing ratio wrt dry+vapor (no condensates)	kg kg-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	inout	False
qi	ice_water_mixing_ratio_updated_by_physics	ice water mixing ratio wrt dry+vapor (no condensates)	kg kg-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	inout	False
qs	snow_water_mixing_ratio_updated_by_physics	snow water mixing ratio wrt dry+vapor (no condensates)	kg kg-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	inout	False
qg	graupel_mixing_ratio_updated_by_physics	graupel mixing ratio wrt dry+vapor (no condensates)	kg kg-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	inout	False
ni	ice_number_concentration_updated_by_physics	ice number concentration	kg-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	inout	False
nr	rain_number_concentration_updated_by_physics	rain number concentration	kg-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	inout	False
is_aerosol_aware	flag_for_aerosol_physics	flag for aerosol-aware physics	flag	logical	()		in	False
nc	cloud_droplet_number_concentration_updated_by_physics	cloud droplet number concentration	kg-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	inout	True
nwfa	water_friendly_aerosol_number_concentration_updated_by_physics	number concentration of water-friendly aerosols	kg-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	inout	True
nifa	ice_friendly_aerosol_number_concentration_updated_by_physics	number concentration of ice-friendly aerosols	kg-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	inout	True
nwfa2d	tendency_of_water_friendly_aerosols_at_surface	instantaneous fake water-friendly surface aerosol source	kg-1 s-1	real	(horizontal_dimension)	kind_phys	in	True
nifa2d	tendency_of_ice_friendly_aerosols_at_surface	instantaneous fake ice-friendly surface aerosol source	kg-1 s-1	real	(horizontal_dimension)	kind_phys	in	True
tgrs	air_temperature_updated_by_physics	model layer mean temperature	K	real	(horizontal_dimension, vertical_dimension)	kind_phys	inout	False
prsl	air_pressure	mean layer pressure	Pa	real	(horizontal_dimension, vertical_dimension)	kind_phys	in	False
phii	geopotential_at_interface	geopotential at model layer interfaces	m2 s-2	real	(horizontal_dimension, vertical_dimension_plus_one)	kind_phys	in	False
omega	omega	layer mean vertical velocity	Pa s-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	in	False
dtp	time_step_for_physics	physics timestep	s	real	()	kind_phys	in	False
prcp	lwe_thickness_of_explicit_precipitation_amount	explicit precipitation (rain, ice, snow, graupel) on physics timestep	m	real	(horizontal_dimension)	kind_phys	inout	False
rain	lwe_thickness_of_explicit_rain_amount	explicit rain fall on physics timestep	m	real	(horizontal_dimension)	kind_phys	inout	False
graupel	lwe_thickness_of_graupel_amount	graupel fall on physics timestep	m	real	(horizontal_dimension)	kind_phys	inout	False
ice	lwe_thickness_of_ice_amount	ice fall on physics timestep	m	real	(horizontal_dimension)	kind_phys	inout	False
snow	lwe_thickness_of_snow_amount	snow fall on physics timestep	m	real	(horizontal_dimension)	kind_phys	inout	False
sr	ratio_of_snowfall_to_rainfall	ratio of snowfall to large-scale rainfall	frac	real	(horizontal_dimension)	kind_phys	out	False
refl_10cm	radar_reflectivity_10cm	instantaneous refl_10cm	dBZ	real	(horizontal_dimension, vertical_dimension)	kind_phys	out	False
reset	flag_for_resetting_radar_reflectivity_calculation	flag for resetting radar reflectivity calculation	flag	logical	()		in	False
do_radar_ref	flag_for_radar_reflectivity	flag for radar reflectivity	flag	logical	()		in	False
re_cloud	effective_radius_of_stratiform_cloud_liquid_water_particle_in_um	eff. radius of cloud liquid water particle in micrometer (meter here)	m	real	(horizontal_dimension, vertical_dimension)	kind_phys	out	True
re_ice	effective_radius_of_stratiform_cloud_ice_particle_in_um	eff. radius of cloud ice water particle in micrometer (meter here)	m	real	(horizontal_dimension, vertical_dimension)	kind_phys	out	True
re_snow	effective_radius_of_stratiform_cloud_snow_particle_in_um	effective radius of cloud snow particle in micrometer (meter here)	m	real	(horizontal_dimension, vertical_dimension)	kind_phys	out	True
mpicomm	mpi_comm	MPI communicator	index	integer	()		in	False
mpirank	mpi_rank	current MPI-rank	index	integer	()		in	False
mpiroot	mpi_root	master MPI-rank	index	integer	()		in	False
errmsg	ccpp_error_message	error message for error handling in CCPP	none	character	()	len=*	out	False
errflg	ccpp_error_flag	error flag for error handling in CCPP	flag	integer	()		out	False

Thompson MP General Algorithm

• Convert specific humidity/moist mixing ratios to dry mixing ratios
• Convert number concentrations from moist to dry
• Density of air in kg m-3
• Convert omega in Pa s-1 to vertical velocity w in m s-1
• Layer width in m from geopotential in m2 s-2
• Call mp_gt_driver() with or without aerosols
• Convert dry mixing ratios to specific humidity/moist mixing ratios
• Convert number concentrations from dry to moist
• Convert rainfall deltas from mm to m (on physics timestep); add to inout variables

References is_initialized.