CCPP SciDoc v7.0.0  v7.0.0
Common Community Physics Package Developed at DTC
 
Loading...
Searching...
No Matches
GFS Scale-aware TKE-based Moist Eddy-Diffusivity Mass-flux (TKE-EDMF) Scheme Module

This subroutine contains all of the logic for the scale-aware TKE-based moist eddy-diffusion mass-flux (TKE-EDMF) scheme. More...

subroutine mfpblt_mod::mfpblt (im, ix, km, kmpbl, ntcw, ntrac1, delt, cnvflg, zl, zm, q1, t1, u1, v1, plyr, pix, thlx, thvx, gdx, hpbl, kpbl, vpert, buo, xmf, tcko, qcko, ucko, vcko, xlamue)
 This subroutine computes mass flux and updraft parcel properties for thermals driven by surface heating.
 
subroutine mfscu_mod::mfscu (im, ix, km, kmscu, ntcw, ntrac1, delt, cnvflg, zl, zm, q1, t1, u1, v1, plyr, pix, thlx, thvx, thlvx, gdx, thetae, radj, krad, mrad, radmin, buo, xmfd, tcdo, qcdo, ucdo, vcdo, xlamde)
 This subroutine computes mass flux and downdraft parcel properties for stratocumulus-top-driven turbulence.
 

Detailed Description

Argument Table

satmedmfvdif_run argument table
local_namestandard_namelong_nameunitstypedimensionskindintent
imhorizontal_loop_extenthorizontal loop extentcountinteger()in
kmvertical_layer_dimensionvertical layer dimensioncountinteger()in
ntracnumber_of_vertical_diffusion_tracersnumber of tracers to diffuse verticallycountinteger()in
ntcwindex_for_liquid_cloud_condensate_vertical_diffusion_tracertracer index for cloud condensate (or liquid water)indexinteger()in
ntiwindex_for_ice_cloud_condensate_vertical_diffusion_tracertracer index for ice water in the vertically diffused tracer arrayindexinteger()in
ntkeindex_for_turbulent_kinetic_energy_vertical_diffusion_tracerindex for turbulent kinetic energy in the vertically diffused tracer arrayindexinteger()in
gravgravitational_accelerationgravitational accelerationm s-2real()kind_physin
rdgas_constant_of_dry_airideal gas constant for dry airJ kg-1 K-1real()kind_physin
cpspecific_heat_of_dry_air_at_constant_pressurespecific heat of dry air at constant pressureJ kg-1 K-1real()kind_physin
rvgas_constant_water_vaporideal gas constant for water vaporJ kg-1 K-1real()kind_physin
hvaplatent_heat_of_vaporization_of_water_at_0clatent heat of evaporation/sublimationJ kg-1real()kind_physin
hfuslatent_heat_of_fusion_of_water_at_0clatent heat of fusionJ kg-1real()kind_physin
fvratio_of_vapor_to_dry_air_gas_constants_minus_one(rv/rd) - 1 (rv = ideal gas constant for water vapor)nonereal()kind_physin
epsratio_of_dry_air_to_water_vapor_gas_constantsrd/rvnonereal()kind_physin
epsm1ratio_of_dry_air_to_water_vapor_gas_constants_minus_one(rd/rv) - 1nonereal()kind_physin
dvprocess_split_cumulative_tendency_of_y_windupdated tendency of the y windm s-2real(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)kind_physinout
duprocess_split_cumulative_tendency_of_x_windupdated tendency of the x windm s-2real(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)kind_physinout
tdtprocess_split_cumulative_tendency_of_air_temperatureupdated tendency of the temperatureK s-1real(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)kind_physinout
rtgtendency_of_vertically_diffused_tracer_concentrationupdated tendency of the tracers due to vertical diffusion in PBL schemekg kg-1 s-1real(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension, ccpp_constant_one:number_of_vertical_diffusion_tracers)kind_physinout
u1x_windx component of layer windm s-1real(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)kind_physin
v1y_windy component of layer windm s-1real(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)kind_physin
t1air_temperaturelayer mean air temperatureKreal(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)kind_physin
q1vertically_diffused_tracer_concentrationtracer concentration diffused by PBL schemekg kg-1real(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension, ccpp_constant_one:number_of_vertical_diffusion_tracers)kind_physin
swhtendency_of_air_temperature_due_to_shortwave_heating_on_radiation_timesteptotal sky shortwave heating rateK s-1real(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)kind_physin
hlwtendency_of_air_temperature_due_to_longwave_heating_on_radiation_timesteptotal sky longwave heating rateK s-1real(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)kind_physin
xmuzenith_angle_temporal_adjustment_factor_for_shortwave_fluxeszenith angle temporal adjustment factor for shortwavenonereal(ccpp_constant_one:horizontal_loop_extent)kind_physin
gareacell_areaarea of the grid cellm2real(ccpp_constant_one:horizontal_loop_extent)kind_physin
psksurface_dimensionless_exner_functiondimensionless Exner function at the surface interfacenonereal(ccpp_constant_one:horizontal_loop_extent)kind_physin
rbsoilbulk_richardson_number_at_lowest_model_levelbulk Richardson number at the surfacenonereal(ccpp_constant_one:horizontal_loop_extent)kind_physin
zorlsurface_roughness_lengthsurface roughness length in cmcmreal(ccpp_constant_one:horizontal_loop_extent)kind_physin
u10mx_wind_at_10mx component of wind at 10 mm s-1real(ccpp_constant_one:horizontal_loop_extent)kind_physin
v10my_wind_at_10my component of wind at 10 mm s-1real(ccpp_constant_one:horizontal_loop_extent)kind_physin
fmmonin_obukhov_similarity_function_for_momentumMonin-Obukhov similarity function for momentumnonereal(ccpp_constant_one:horizontal_loop_extent)kind_physin
fhmonin_obukhov_similarity_function_for_heatMonin-Obukhov similarity function for heatnonereal(ccpp_constant_one:horizontal_loop_extent)kind_physin
tseasurface_skin_temperaturesurface skin temperatureKreal(ccpp_constant_one:horizontal_loop_extent)kind_physin
heatkinematic_surface_upward_sensible_heat_flux_reduced_by_surface_roughness_and_vegetationkinematic surface upward sensible heat flux reduced by surface roughness and vegetationK m s-1real(ccpp_constant_one:horizontal_loop_extent)kind_physin
evapsurface_upward_specific_humidity_fluxkinematic surface upward latent heat fluxkg kg-1 m s-1real(ccpp_constant_one:horizontal_loop_extent)kind_physin
stresssurface_wind_stresssurface wind stressm2 s-2real(ccpp_constant_one:horizontal_loop_extent)kind_physin
spd1wind_speed_at_lowest_model_layerwind speed at lowest model levelm s-1real(ccpp_constant_one:horizontal_loop_extent)kind_physin
kpblvertical_index_at_top_of_atmosphere_boundary_layerPBL top model level indexindexinteger(ccpp_constant_one:horizontal_loop_extent)out
prsiair_pressure_at_interfaceair pressure at model layer interfacesPareal(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_interface_dimension)kind_physin
delair_pressure_difference_between_midlayerspres(k) - pres(k+1)Pareal(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)kind_physin
prslair_pressuremean layer pressurePareal(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)kind_physin
prslkdimensionless_exner_functionExner function at layersnonereal(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)kind_physin
phiigeopotential_at_interfacegeopotential at model layer interfacesm2 s-2real(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_interface_dimension)kind_physin
philgeopotentialgeopotential at model layer centersm2 s-2real(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension)kind_physin
delttimestep_for_physicstime step for physicssreal()kind_physin
dspheatflag_tke_dissipation_heatingflag for using TKE dissipation heatingflaglogical()in
dusfcinstantaneous_surface_x_momentum_fluxx momentum fluxPareal(ccpp_constant_one:horizontal_loop_extent)kind_physout
dvsfcinstantaneous_surface_y_momentum_fluxy momentum fluxPareal(ccpp_constant_one:horizontal_loop_extent)kind_physout
dtsfcinstantaneous_surface_upward_sensible_heat_fluxsurface upward sensible heat fluxW m-2real(ccpp_constant_one:horizontal_loop_extent)kind_physout
dqsfcinstantaneous_surface_upward_latent_heat_fluxsurface upward latent heat fluxW m-2real(ccpp_constant_one:horizontal_loop_extent)kind_physout
hpblatmosphere_boundary_layer_thicknessPBL thicknessmreal(ccpp_constant_one:horizontal_loop_extent)kind_physout
kinverindex_of_highest_temperature_inversionindex of highest temperature inversionindexinteger(ccpp_constant_one:horizontal_loop_extent)in
xkzm_matmosphere_momentum_diffusivity_due_to_backgroundbackground value of momentum diffusivitym2 s-1real()kind_physin
xkzm_hatmosphere_heat_diffusivity_due_to_backgroundbackground value of heat diffusivitym2 s-1real()kind_physin
xkzm_ssigma_pressure_threshold_at_upper_extent_of_background_diffusionsigma level threshold for background diffusivitynonereal()kind_physin
dtendcumulative_change_of_state_variablesdiagnostic tendencies for state variablesmixedreal(ccpp_constant_one:horizontal_loop_extent, ccpp_constant_one:vertical_layer_dimension, ccpp_constant_one:cumulative_change_of_state_variables_outer_index_max)kind_physinout
dtidxcumulative_change_of_state_variables_outer_indexindex of state-variable and process in last dimension of diagnostic tendencies array AKA cumulative_change_indexindexinteger(ccpp_constant_one:number_of_tracers_plus_one_hundred, ccpp_constant_one:number_of_cumulative_change_processes)in
index_of_temperatureindex_of_temperature_in_cumulative_change_indexindex of temperature in first dimension of array cumulative change indexindexinteger()in
index_of_x_windindex_of_x_wind_in_cumulative_change_indexindex of x-wind in first dimension of array cumulative change indexindexinteger()in
index_of_y_windindex_of_y_wind_in_cumulative_change_indexindex of x-wind in first dimension of array cumulative change indexindexinteger()in
index_of_process_pblindex_of_subgrid_scale_vertical_mixing_process_in_cumulative_change_indexindex of subgrid scale vertical mixing process in second dimension of array cumulative change indexindexinteger()in
ntqvindex_of_specific_humidity_in_tracer_concentration_arraytracer index for water vapor (specific humidity)indexinteger()in
ntozindex_of_ozone_mixing_ratio_in_tracer_concentration_arraytracer index for ozone mixing ratioindexinteger()in
gen_tendflag_for_generic_tendency_due_to_planetary_boundary_layertrue if GFS_PBL_generic should calculate tendenciesflaglogical()in
ldiag3dflag_for_diagnostics_3dflag for 3d diagnostic fieldsflaglogical()in
errmsgccpp_error_messageerror message for error handling in CCPPnonecharacter()len=*out
errflgccpp_error_codeerror code for error handling in CCPP1integer()out

GFS satmedmfvdif General Algorithm

satmedmfvdif_run() computes subgrid vertical turbulence mixing using the scale-aware TKE-based moist eddy-diffusion mass-flux (EDMF) parameterization of Han and Bretherton (2019) [75] .

  1. The local turbulent mixing is represented by an eddy-diffusivity scheme which is a function of a prognostic TKE.
  2. For the convective boundary layer, nonlocal transport by large eddies (mfpblt.f), is represented using a mass flux approach (Siebesma et al.(2007) [180] ).
  3. A mass-flux approach is also used to represent the stratocumulus-top-induced turbulence (mfscu.f).

GFS satmedmfvdif Detailed Algorithm