mp_thompson_run()

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), intent(in)	con_eps,
		logical, intent(in)	convert_dry_rho,
		real(kind_phys), dimension(:,:), intent(inout)	spechum,
		real(kind_phys), dimension(:,:), intent(inout)	qc,
		real(kind_phys), dimension(:,:), intent(inout)	qr,
		real(kind_phys), dimension(:,:), intent(inout)	qi,
		real(kind_phys), dimension(:,:), intent(inout)	qs,
		real(kind_phys), dimension(:,:), intent(inout)	qg,
		real(kind_phys), dimension(:,:), intent(inout)	ni,
		real(kind_phys), dimension(:,:), 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(:,:), intent(inout)	tgrs,
		real(kind_phys), dimension(:,:), intent(in)	prsl,
		real(kind_phys), dimension(:,:), intent(in)	phii,
		real(kind_phys), dimension(:,:), intent(in)	omega,
		logical, intent(in)	sedi_semi,
		integer, intent(in)	decfl,
		real(kind_phys), intent(in)	dtp,
		real, intent(in)	dt_inner,
		logical, intent(in)	first_time_step,
		integer, intent(in)	istep,
		integer, intent(in)	nsteps,
		real(kind_phys), dimension(:), intent(inout)	prcp,
		real(kind_phys), dimension(:), intent(inout)	rain,
		real(kind_phys), dimension(:), intent(inout)	graupel,
		real(kind_phys), dimension(:), intent(inout)	ice,
		real(kind_phys), dimension(:), intent(inout)	snow,
		real(kind_phys), dimension(:), intent(out)	sr,
		real(kind_phys), dimension(:,:), intent(inout)	refl_10cm,
		logical, intent(in)	reset_dBZ,
		logical, intent(in)	do_radar_ref,
		integer, intent(in)	mpicomm,
		integer, intent(in)	mpirank,
		integer, intent(in)	mpiroot,
		integer, intent(in)	blkno,
		logical, intent(in)	ext_diag,
		real(kind_phys), dimension(:,:,:), intent(inout), target	diag3d,
		logical, intent(in)	reset_diag3d,
		real(kind_phys), dimension(:,:), intent(in)	spp_wts_mp,
		integer, intent(in)	spp_mp,
		integer, intent(in)	n_var_spp,
		real(kind_phys), dimension(:), intent(in)	spp_prt_list,
		character(len=3), dimension(:), intent(in)	spp_var_list,
		real(kind_phys), dimension(:), intent(in)	spp_stddev_cutoff,
		character(len=*), intent(out)	errmsg,
		integer, intent(out)	errflg
	)

Argument Table

local_name	standard_name	long_name	units	type	dimensions	kind	intent
ncol	horizontal_loop_extent	horizontal loop extent	count	integer	()		in
nlev	vertical_layer_dimension	number of vertical levels	count	integer	()		in
con_g	gravitational_acceleration	gravitational acceleration	m s-2	real	()	kind_phys	in
con_rd	gas_constant_of_dry_air	ideal gas constant for dry air	J kg-1 K-1	real	()	kind_phys	in
con_eps	ratio_of_dry_air_to_water_vapor_gas_constants	rd/rv	none	real	()	kind_phys	in
convert_dry_rho	flag_for_converting_hydrometeors_from_moist_to_dry_air	flag for converting hydrometeors from moist to dry air	flag	logical	()		in
spechum	specific_humidity_of_new_state	water vapor specific humidity	kg kg-1	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	inout
qc	cloud_liquid_water_mixing_ratio_of_new_state	cloud water mixing ratio wrt dry+vapor (no condensates)	kg kg-1	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	inout
qr	rain_mixing_ratio_of_new_state	rain water mixing ratio wrt dry+vapor (no condensates)	kg kg-1	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	inout
qi	cloud_ice_mixing_ratio_of_new_state	ice water mixing ratio wrt dry+vapor (no condensates)	kg kg-1	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	inout
qs	snow_mixing_ratio_of_new_state	snow water mixing ratio wrt dry+vapor (no condensates)	kg kg-1	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	inout
qg	graupel_mixing_ratio_of_new_state	graupel mixing ratio wrt dry+vapor (no condensates)	kg kg-1	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	inout
ni	mass_number_concentration_of_cloud_ice_water_crystals_in_air_of_new_state	ice number concentration	kg-1	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	inout
nr	mass_number_concentration_of_rain_of_new_state	rain number concentration	kg-1	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	inout
is_aerosol_aware	flag_for_aerosol_physics	flag for aerosol-aware physics	flag	logical	()		in
nc	mass_number_concentration_of_cloud_liquid_water_particles_in_air_of_new_state	cloud droplet number concentration	kg-1	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	inout
nwfa	mass_number_concentration_of_hygroscopic_aerosols_of_new_state	number concentration of water-friendly aerosols	kg-1	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	inout
nifa	mass_number_concentration_of_nonhygroscopic_ice_nucleating_aerosols_of_new_state	number concentration of ice-friendly aerosols	kg-1	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	inout
nwfa2d	tendency_of_hygroscopic_aerosols_at_surface_adjacent_layer	instantaneous fake water-friendly surface aerosol source	kg-1 s-1	real	(ccpp_constant_one:horizontal_loop_extent)	kind_phys	in
nifa2d	tendency_of_nonhygroscopic_ice_nucleating_aerosols_at_surface_adjacent_layer	instantaneous fake ice-friendly surface aerosol source	kg-1 s-1	real	(ccpp_constant_one:horizontal_loop_extent)	kind_phys	in
tgrs	air_temperature_of_new_state	model layer mean temperature	K	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	inout
prsl	air_pressure	mean layer pressure	Pa	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	in
phii	geopotential_at_interface	geopotential at model layer interfaces	m2 s-2	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_interface_dimension)	kind_phys	in
omega	lagrangian_tendency_of_air_pressure	layer mean vertical velocity	Pa s-1	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	in
sedi_semi	flag_for_semi_lagrangian_sedi_rain	flag for semi Lagrangian sedi of rain	flag	logical	()		in
decfl	deformed_cfl_factor	deformed CFL factor	count	integer	()		in
dtp	timestep_for_physics	physics timestep	s	real	()	kind_phys	in
dt_inner	time_step_for_inner_loop	time step for inner loop	s	real	()	kind_phys	in
first_time_step	flag_for_first_timestep	flag for first time step for time integration loop (cold/warmstart)	flag	logical	()		in
istep	ccpp_loop_counter	loop counter for subcycling loops in CCPP	index	integer	()		in
nsteps	ccpp_loop_extent	loop extent for subcycling loops in CCPP	count	integer	()		in
prcp	lwe_thickness_of_explicit_precipitation_amount	explicit precipitation (rain, ice, snow, graupel) on physics timestep	m	real	(ccpp_constant_one:horizontal_loop_extent)	kind_phys	inout
rain	lwe_thickness_of_explicit_rain_amount	explicit rain fall on physics timestep	m	real	(ccpp_constant_one:horizontal_loop_extent)	kind_phys	inout
graupel	lwe_thickness_of_graupel_amount	graupel fall on physics timestep	m	real	(ccpp_constant_one:horizontal_loop_extent)	kind_phys	inout
ice	lwe_thickness_of_ice_amount	ice fall on physics timestep	m	real	(ccpp_constant_one:horizontal_loop_extent)	kind_phys	inout
snow	lwe_thickness_of_snow_amount	snow fall on physics timestep	m	real	(ccpp_constant_one:horizontal_loop_extent)	kind_phys	inout
sr	ratio_of_snowfall_to_rainfall	ratio of snowfall to large-scale rainfall	frac	real	(ccpp_constant_one:horizontal_loop_extent)	kind_phys	out
refl_10cm	radar_reflectivity_10cm	instantaneous refl_10cm	dBZ	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	out
reset_dBZ	flag_for_resetting_radar_reflectivity_calculation	flag for resetting radar reflectivity calculation	flag	logical	()		in
do_radar_ref	flag_for_radar_reflectivity	flag for radar reflectivity	flag	logical	()		in
mpicomm	mpi_communicator	MPI communicator	index	integer	()		in
mpirank	mpi_rank	current MPI-rank	index	integer	()		in
mpiroot	mpi_root	master MPI-rank	index	integer	()		in
blkno	ccpp_block_number	number of block for explicit data blocking in CCPP	index	integer	()		in
ext_diag	flag_for_extended_diagnostic_output_from_thompson_microphysics	flag for extended diagnostic output from thompson microphysics	flag	logical	()		in
diag3d	extended_diagnostics_output_from_thompson_microphysics	set of 3d arrays for extended diagnostics output from thompson microphysics	none	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension, ccpp_constant_one:number_of_3d_diagnostic_output_arrays_from_thompson_microphysics)	kind_phys	inout
reset_diag3d	flag_reset_extended_diagnostics_output_arrays_from_thompson_microphysics	flag for resetting extended diagnostics output arrays from thompson microphysics	flag	logical	()		in
spp_wts_mp	spp_weights_for_microphysics_scheme	spp weights for microphysics scheme	1	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	in
spp_mp	control_for_microphysics_spp_perturbations	control for microphysics spp perturbations	count	integer	()		in
n_var_spp	number_of_perturbed_spp_schemes	number of perturbed spp schemes	count	integer	()		in
spp_prt_list	magnitude_of_spp_perturbations	magnitude of spp perturbations	1	real	(ccpp_constant_one:number_of_perturbed_spp_schemes)	kind_phys	in
spp_stddev_cutoff	magnitude_of_spp_standard_deviation_cutoff	magnitude of spp standard deviation cutoff	1	real	(ccpp_constant_one:number_of_perturbed_spp_schemes)	kind_phys	in
spp_var_list	perturbed_spp_schemes	perturbed spp schemes	none	character	(ccpp_constant_one:number_of_perturbed_spp_schemes)	len=3	in
errmsg	ccpp_error_message	error message for error handling in CCPP	none	character	()	len=*	out
errflg	ccpp_error_code	error code for error handling in CCPP	1	integer	()		out

Thompson MP General Algorithm

• Convert specific humidity to water vapor mixing ratio.
  • Also, hydrometeor variables are mass or number mixing ratio
  • either kg of species per kg of dry air, or per kg of (dry + vapor).
• 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 water vapor mixing ratio back to specific humidity
• Convert rainfall deltas from mm to m (on physics timestep); add to inout variables

References module_mp_thompson::mp_gt_driver().

Here is the call graph for this function: