GFS Hybrid Eddy-Diffusivity Mass-Flux (HEDMF) Scheme Module

This subroutine contains all of logic for the Hybrid EDMF PBL scheme except for the calculation of the updraft properties and mass flux. More...

Detailed Description

Argument Table

local_name	standard_name	long_name	units	type	dimensions	kind	intent	optional
im	horizontal_loop_extent	horizontal loop extent	count	integer	()		in	False
km	vertical_dimension	vertical layer dimension	count	integer	()		in	False
ntrac	number_of_vertical_diffusion_tracers	number of tracers to diffuse vertically	count	integer	()		in	False
ntcw	index_for_liquid_cloud_condensate	cloud condensate index in tracer array	index	integer	()		in	False
dv	tendency_of_y_wind_due_to_model_physics	updated tendency of the y wind	m s-2	real	(horizontal_dimension, vertical_dimension)	kind_phys	inout	False
du	tendency_of_x_wind_due_to_model_physics	updated tendency of the x wind	m s-2	real	(horizontal_dimension, vertical_dimension)	kind_phys	inout	False
tau	tendency_of_air_temperature_due_to_model_physics	updated tendency of the temperature	K s-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	inout	False
rtg	tendency_of_vertically_diffused_tracer_concentration	updated tendency of the tracers due to vertical diffusion in PBL scheme	kg kg-1 s-1	real	(horizontal_dimension, vertical_dimension, number_of_vertical_diffusion_tracers)	kind_phys	inout	False
u1	x_wind	x component of layer wind	m s-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	in	False
v1	y_wind	y component of layer wind	m s-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	in	False
t1	air_temperature	layer mean air temperature	K	real	(horizontal_dimension, vertical_dimension)	kind_phys	in	False
q1	vertically_diffused_tracer_concentration	tracer concentration diffused by PBL scheme	kg kg-1	real	(horizontal_dimension, vertical_dimension, number_of_vertical_diffusion_tracers)	kind_phys	in	False
swh	tendency_of_air_temperature_due_to_shortwave_heating_on_radiation_time_step	total sky shortwave heating rate	K s-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	in	False
hlw	tendency_of_air_temperature_due_to_longwave_heating_on_radiation_time_step	total sky longwave heating rate	K s-1	real	(horizontal_dimension, vertical_dimension)	kind_phys	in	False
xmu	zenith_angle_temporal_adjustment_factor_for_shortwave_fluxes	zenith angle temporal adjustment factor for shortwave	none	real	(horizontal_dimension)	kind_phys	in	False
psk	dimensionless_exner_function_at_lowest_model_interface	dimensionless Exner function at the surface interface	none	real	(horizontal_dimension)	kind_phys	in	False
rbsoil	bulk_richardson_number_at_lowest_model_level	bulk Richardson number at the surface	none	real	(horizontal_dimension)	kind_phys	in	False
zorl	surface_roughness_length	surface roughness length in cm	cm	real	(horizontal_dimension)	kind_phys	in	False
u10m	x_wind_at_10m	x component of wind at 10 m	m s-1	real	(horizontal_dimension)	kind_phys	in	False
v10m	y_wind_at_10m	y component of wind at 10 m	m s-1	real	(horizontal_dimension)	kind_phys	in	False
fm	Monin_Obukhov_similarity_function_for_momentum	Monin-Obukhov similarity function for momentum	none	real	(horizontal_dimension)	kind_phys	in	False
fh	Monin_Obukhov_similarity_function_for_heat	Monin-Obukhov similarity function for heat	none	real	(horizontal_dimension)	kind_phys	in	False
tsea	surface_skin_temperature	surface skin temperature	K	real	(horizontal_dimension)	kind_phys	in	False
heat	kinematic_surface_upward_sensible_heat_flux_reduced_by_surface_roughness	kinematic surface upward sensible heat flux	K m s-1	real	(horizontal_dimension)	kind_phys	in	False
evap	kinematic_surface_upward_latent_heat_flux_reduced_by_surface_roughness	kinematic surface upward latent heat flux	kg kg-1 m s-1	real	(horizontal_dimension)	kind_phys	in	False
stress	surface_wind_stress	surface wind stress	m2 s-2	real	(horizontal_dimension)	kind_phys	in	False
spd1	wind_speed_at_lowest_model_layer	wind speed at lowest model level	m s-1	real	(horizontal_dimension)	kind_phys	in	False
kpbl	vertical_index_at_top_of_atmosphere_boundary_layer	PBL top model level index	index	integer	(horizontal_dimension)		out	False
prsi	air_pressure_at_interface	air pressure at model layer interfaces	Pa	real	(horizontal_dimension, vertical_dimension_plus_one)	kind_phys	in	False
del	air_pressure_difference_between_midlayers	pres(k) - pres(k+1)	Pa	real	(horizontal_dimension, vertical_dimension)	kind_phys	in	False
prsl	air_pressure	mean layer pressure	Pa	real	(horizontal_dimension, vertical_dimension)	kind_phys	in	False
prslk	dimensionless_exner_function_at_model_layers	Exner function at layers	none	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
phil	geopotential	geopotential at model layer centers	m2 s-2	real	(horizontal_dimension, vertical_dimension)	kind_phys	in	False
delt	time_step_for_physics	time step for physics	s	real	()	kind_phys	in	False
dspheat	flag_TKE_dissipation_heating	flag for using TKE dissipation heating	flag	logical	()		in	False
dusfc	instantaneous_surface_x_momentum_flux	x momentum flux	Pa	real	(horizontal_dimension)	kind_phys	out	False
dvsfc	instantaneous_surface_y_momentum_flux	y momentum flux	Pa	real	(horizontal_dimension)	kind_phys	out	False
dtsfc	instantaneous_surface_upward_sensible_heat_flux	surface upward sensible heat flux	W m-2	real	(horizontal_dimension)	kind_phys	out	False
dqsfc	instantaneous_surface_upward_latent_heat_flux	surface upward latent heat flux	W m-2	real	(horizontal_dimension)	kind_phys	out	False
hpbl	atmosphere_boundary_layer_thickness	PBL thickness	m	real	(horizontal_dimension)	kind_phys	out	False
hgamt	countergradient_mixing_term_for_temperature	countergradient mixing term for temperature	K	real	(horizontal_dimension)	kind_phys	inout	False
hgamq	countergradient_mixing_term_for_water_vapor	countergradient mixing term for water vapor	kg kg-1	real	(horizontal_dimension)	kind_phys	inout	False
dkt	atmosphere_heat_diffusivity	diffusivity for heat	m2 s-1	real	(horizontal_dimension, vertical_dimension_minus_one)	kind_phys	out	False
kinver	index_of_highest_temperature_inversion	index of highest temperature inversion	index	integer	(horizontal_dimension)		in	False
xkzm_m	atmosphere_momentum_diffusivity_background	background value of momentum diffusivity	m2 s-1	real	()	kind_phys	in	False
xkzm_h	atmosphere_heat_diffusivity_background	background value of heat diffusivity	m2 s-1	real	()	kind_phys	in	False
xkzm_s	diffusivity_background_sigma_level	sigma level threshold for background diffusivity	none	real	()	kind_phys	in	False
lprnt	flag_print	flag for printing diagnostics to output	flag	logical	()		in	False
ipr	horizontal_index_of_printed_column	horizontal index of printed column	index	integer	()		in	False
xkzminv	atmosphere_heat_diffusivity_background_maximum	maximum background value of heat diffusivity	m2 s-1	real	()	kind_phys	in	False
moninq_fac	atmosphere_diffusivity_coefficient_factor	multiplicative constant for atmospheric diffusivities	none	real	()	kind_phys	in	False
lssav	flag_diagnostics	logical flag for storing diagnostics	flag	logical	()		in	False
ldiag3d	flag_diagnostics_3D	flag for 3d diagnostic fields	flag	logical	()		in	False
qdiag3d	flag_tracer_diagnostics_3D	flag for 3d tracer diagnostic fields	flag	logical	()		in	False
ntoz	index_for_ozone	tracer index for ozone mixing ratio	index	integer	()		in	False
du3dt_PBL	cumulative_change_in_x_wind_due_to_PBL	cumulative change in x wind due to PBL	m s-1	real	(horizontal_loop_extent, vertical_dimension)	kind_phys	inout	False
dv3dt_PBL	cumulative_change_in_y_wind_due_to_PBL	cumulative change in y wind due to PBL	m s-1	real	(horizontal_loop_extent, vertical_dimension)	kind_phys	inout	False
dt3dt_PBL	cumulative_change_in_temperature_due_to_PBL	cumulative change in temperature due to PBL	K	real	(horizontal_loop_extent, vertical_dimension)	kind_phys	inout	False
dq3dt_PBL	cumulative_change_in_water_vapor_specific_humidity_due_to_PBL	cumulative change in water vapor specific humidity due to PBL	kg kg-1	real	(horizontal_loop_extent, vertical_dimension)	kind_phys	inout	False
do3dt_PBL	cumulative_change_in_ozone_mixing_ratio_due_to_PBL	cumulative change in ozone mixing ratio due to PBL	kg kg-1	real	(horizontal_loop_extent, vertical_dimension)	kind_phys	inout	False
flag_for_pbl_generic_tend	flag_for_generic_planetary_boundary_layer_tendency	true if GFS_PBL_generic should calculate tendencies	flag	logical	()		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

GFS Hybrid EDMF General Algorithm

1. Compute preliminary variables from input arguments.
2. Calculate the first estimate of the PBL height ("Predictor step").
3. Calculate Monin-Obukhov similarity parameters.
4. Update thermal properties of surface parcel and recompute PBL height ("Corrector step").
5. Determine whether stratocumulus layers exist and compute quantities needed for enhanced diffusion.
6. Calculate the inverse Prandtl number.
7. Compute diffusion coefficients below the PBL top.
8. Compute diffusion coefficients above the PBL top.
9. If the PBL is convective, call the mass flux scheme to replace the countergradient terms.
10. Compute enhanced diffusion coefficients related to stratocumulus-topped PBLs.
11. Solve for the temperature and moisture tendencies due to vertical mixing.
12. Calculate heating due to TKE dissipation and add to the tendency for temperature.
13. Solve for the horizontal momentum tendencies and add them to output tendency terms.

GFS Hybrid HEDMF Detailed Algorithm

Functions/Subroutines
subroutine	hedmf::hedmf_run (im, km, ntrac, ntcw, dv, du, tau, rtg, u1, v1, t1, q1, swh, hlw, xmu, psk, rbsoil, zorl, u10m, v10m, fm, fh, tsea, heat, evap, stress, spd1, kpbl, prsi, del, prsl, prslk, phii, phil, delt, dspheat, dusfc, dvsfc, dtsfc, dqsfc, hpbl, hgamt, hgamq, dkt, kinver, xkzm_m, xkzm_h, xkzm_s, lprnt, ipr, xkzminv, moninq_fac, lssav, ldiag3d, qdiag3d, ntoz, du3dt_PBL, dv3dt_PBL, dt3dt_PBL, dq3dt_PBL, do3dt_PBL, flag_for_pbl_generic_tend, errmsg, errflg)
subroutine	tridi1 (l, n, cl, cm, cu, r1, au, a1)
	Routine to solve the tridiagonal system to calculate temperature and moisture at \( t + \Delta t \); part of two-part process to calculate time tendencies due to vertical diffusion. More...
subroutine	mfpbl (im, ix, km, ntrac, delt, cnvflg, zl, zm, thvx, q1, t1, u1, v1, hpbl, kpbl, sflx, ustar, wstar, xmf, tcko, qcko, ucko, vcko)
	This subroutine is used for calculating the mass flux and updraft properties. More...