Common Community Physics Package (CCPP) Scientific Documentation  Version 2.0
GFS moninedmf Main

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 rank type kind intent optional
ix horizontal_dimension horizontal dimension count 0 integer in F
im horizontal_loop_extent horizontal loop extent count 0 integer in F
km vertical_dimension vertical layer dimension count 0 integer in F
ntrac number_of_vertical_diffusion_tracers number of tracers to diffuse vertically count 0 integer in F
ntcw index_for_liquid_cloud_condensate cloud condensate index in tracer array index 0 integer in F
dv tendency_of_y_wind_due_to_model_physics updated tendency of the y wind m s-2 2 real kind_phys inout F
du tendency_of_x_wind_due_to_model_physics updated tendency of the x wind m s-2 2 real kind_phys inout F
tau tendency_of_air_temperature_due_to_model_physics updated tendency of the temperature K s-1 2 real kind_phys inout F
rtg tendency_of_vertically_diffused_tracer_concentration updated tendency of the tracers due to vertical diffusion in PBL scheme kg kg-1 s-1 3 real kind_phys inout F
u1 x_wind x component of layer wind m s-1 2 real kind_phys in F
v1 y_wind y component of layer wind m s-1 2 real kind_phys in F
t1 air_temperature layer mean air temperature K 2 real kind_phys in F
q1 vertically_diffused_tracer_concentration tracer concentration diffused by PBL scheme kg kg-1 3 real kind_phys in F
swh tendency_of_air_temperature_due_to_shortwave_heating_on_radiation_time_step total sky shortwave heating rate K s-1 2 real kind_phys in F
hlw tendency_of_air_temperature_due_to_longwave_heating_on_radiation_time_step total sky longwave heating rate K s-1 2 real kind_phys in F
xmu zenith_angle_temporal_adjustment_factor_for_shortwave_fluxes zenith angle temporal adjustment factor for shortwave none 1 real kind_phys in F
psk dimensionless_exner_function_at_lowest_model_interface dimensionless Exner function at the surface interface none 1 real kind_phys in F
rbsoil bulk_richardson_number_at_lowest_model_level bulk Richardson number at the surface none 1 real kind_phys in F
zorl surface_roughness_length surface roughness length in cm cm 1 real kind_phys in F
u10m x_wind_at_10m x component of wind at 10 m m s-1 1 real kind_phys in F
v10m y_wind_at_10m y component of wind at 10 m m s-1 1 real kind_phys in F
fm Monin-Obukhov_similarity_function_for_momentum Monin-Obukhov similarity function for momentum none 1 real kind_phys in F
fh Monin-Obukhov_similarity_function_for_heat Monin-Obukhov similarity function for heat none 1 real kind_phys in F
tsea surface_skin_temperature surface skin temperature K 1 real kind_phys in F
heat kinematic_surface_upward_sensible_heat_flux kinematic surface upward sensible heat flux K m s-1 1 real kind_phys in F
evap kinematic_surface_upward_latent_heat_flux kinematic surface upward latent heat flux kg kg-1 m s-1 1 real kind_phys in F
stress surface_wind_stress surface wind stress m2 s-2 1 real kind_phys in F
spd1 wind_speed_at_lowest_model_layer wind speed at lowest model level m s-1 1 real kind_phys in F
kpbl vertical_index_at_top_of_atmosphere_boundary_layer PBL top model level index index 1 integer out F
prsi air_pressure_at_interface air pressure at model layer interfaces Pa 2 real kind_phys in F
del air_pressure_difference_between_midlayers pres(k) - pres(k+1) Pa 2 real kind_phys in F
prsl air_pressure mean layer pressure Pa 2 real kind_phys in F
prslk dimensionless_exner_function_at_model_layers Exner function at layers none 2 real kind_phys in F
phii geopotential_at_interface geopotential at model layer interfaces m2 s-2 2 real kind_phys in F
phil geopotential geopotential at model layer centers m2 s-2 2 real kind_phys in F
delt time_step_for_physics time step for physics s 0 real kind_phys in F
dspheat flag_TKE_dissipation_heating flag for using TKE dissipation heating flag 0 logical in F
dusfc instantaneous_surface_x_momentum_flux x momentum flux Pa 1 real kind_phys out F
dvsfc instantaneous_surface_y_momentum_flux y momentum flux Pa 1 real kind_phys out F
dtsfc instantaneous_surface_upward_sensible_heat_flux surface upward sensible heat flux W m-2 1 real kind_phys out F
dqsfc instantaneous_surface_upward_latent_heat_flux surface upward latent heat flux W m-2 1 real kind_phys out F
hpbl atmosphere_boundary_layer_thickness PBL thickness m 1 real kind_phys out F
hgamt countergradient_mixing_term_for_temperature countergradient mixing term for temperature K 1 real kind_phys inout F
hgamq countergradient_mixing_term_for_water_vapor countergradient mixing term for water vapor kg kg-1 1 real kind_phys inout F
dkt atmosphere_heat_diffusivity diffusivity for heat m2 s-1 2 real kind_phys out F
kinver index_of_highest_temperature_inversion index of highest temperature inversion index 1 integer in F
xkzm_m atmosphere_momentum_diffusivity_background background value of momentum diffusivity m2 s-1 0 real kind_phys in F
xkzm_h atmosphere_heat_diffusivity_background background value of heat diffusivity m2 s-1 0 real kind_phys in F
xkzm_s diffusivity_background_sigma_level sigma level threshold for background diffusivity none 0 real kind_phys in F
lprnt flag_print flag for printing diagnostics to output flag 0 logical in F
ipr horizontal_index_of_printed_column horizontal index of printed column index 0 integer in F
xkzminv atmosphere_heat_diffusivity_background_maximum maximum background value of heat diffusivity m2 s-1 0 real kind_phys in F
moninq_fac atmosphere_diffusivity_coefficient_factor multiplicative constant for atmospheric diffusivities none 0 real kind_phys in F
errmsg ccpp_error_message error message for error handling in CCPP none 0 character len=* out F
errflg ccpp_error_flag error flag for error handling in CCPP flag 0 integer out F

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.

Detailed Algorithm

Functions/Subroutines

subroutine hedmf::hedmf_run (ix, 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, errmsg, errflg)
 
subroutine tridi2 (l, n, cl, cm, cu, r1, r2, au, a1, a2)
 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 tridin (l, n, nt, cl, cm, cu, r1, r2, au, a1, a2)
 Routine to solve the tridiagonal system to calculate u- and v-momentum 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...