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,
		logical, intent(in)	merra2_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,
		logical, intent(in), optional	aero_ind_fdb,
		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,
		integer, dimension(:), intent(in)	islmsk,
		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), optional	rain,
		real(kind_phys), dimension(:), intent(inout), optional	graupel,
		real(kind_phys), dimension(:), intent(inout), optional	ice,
		real(kind_phys), dimension(:), intent(inout), optional	snow,
		real(kind_phys), dimension(:), intent(out)	sr,
		real(kind_phys), dimension(:,:), intent(inout)	refl_10cm,
		logical, intent(in)	fullradar_diag,
		real(kind_phys), dimension(:), intent(inout)	max_hail_diam_sfc,
		logical, intent(in)	do_radar_ref,
		real(kind_phys), dimension(:,:,:), intent(in)	aerfld,
		type(mpi_comm), 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), optional, target	diag3d,
		logical, intent(in)	reset_diag3d,
		real(kind_phys), dimension(:,:), intent(in), optional	spp_wts_mp,
		integer, intent(in)	spp_mp,
		integer, intent(in)	n_var_spp,
		real(kind_phys), dimension(:), intent(in), optional	spp_prt_list,
		character(len=10), dimension(:), intent(in), optional	spp_var_list,
		real(kind_phys), dimension(:), intent(in)	spp_stddev_cutoff,
		logical, intent(in)	cplchm,
		real(kind=kind_phys), dimension(:,:), intent(inout), optional	pfi_lsan,
		real(kind=kind_phys), dimension(:,:), intent(inout), optional	pfl_lsan,
		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
merra2_aerosol_aware	do_merra2_aerosol_awareness	flag for merra2 aerosol-aware physics for example the thompson microphysics	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
aero_ind_fdb	do_smoke_aerosol_indirect_feedback	flag for wfa ifa emission indirect feedback	flag	logical	()		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
islmsk	sea_land_ice_mask	sea/land/ice mask (=0/1/2)	flag	integer	(ccpp_constant_one:horizontal_loop_extent)		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
max_hail_diam_sfc	max_hail_diameter_sfc	instantaneous maximum hail diameter at lowest model level	m	real	(ccpp_constant_one:horizontal_loop_extent)	kind_phys	inout
fullradar_diag	do_full_radar_reflectivity	flag for computing full radar reflectivity	flag	logical	()		in
do_radar_ref	flag_for_radar_reflectivity	flag for radar reflectivity	flag	logical	()		in
aerfld	mass_mixing_ratio_of_aerosol_from_gocart_or_merra2	mass mixing ratio of aerosol from gocart or merra2	kg kg-1	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension, ccpp_constant_one:number_of_aerosol_tracers_MG)	kind_phys	in
mpicomm	mpi_communicator	MPI communicator	index	MPI_Comm	()	MPI_Comm	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=10	in
cplchm	flag_for_chemistry_coupling	flag controlling cplchm collection (default off)	flag	logical	()		in
pfi_lsan	ice_flux_due_to_large_scale_precipitation	instantaneous 3D flux of ice from nonconvective precipitation	kg m-2 s-1	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	inout
pfl_lsan	liquid_flux_due_to_large_scale_precipitation	instantaneous 3D flux of liquid water from nonconvective precipitation	kg m-2 s-1	real	(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)	kind_phys	inout
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, with or without effective radii, ...

• Convert water vapor mixing ratio back to specific humidity

• Convert rainfall deltas from mm to m (on physics timestep); add to inout variables

Definition at line 322 of file mp_thompson.F90.

References module_mp_thompson::mp_gt_driver().

Here is the call graph for this function: