Convective Heating Rate

gitagiriharta

New member
Hello all,

I'm curious if there's a way to extract the convective heating rate from WRF output data. Any help on this would be greatly appreciated. Thank you!
 

Attachments

  • IMG_8409.jpeg
    IMG_8409.jpeg
    184 KB · Views: 18
RTHCUTEN is the temperature tendency due to cumulus. Other than this, there is no output of convective heating rate in standard WRF.
 
Official WRF codes don't output this variable. You need to modify codes if you want to output this variable.
 
RTHCUTEN is the temperature tendency due to cumulus. Other than this, there is no output of convective heating rate in standard WRF.
I'm confused about the RTHCUTEN and RTHMPTEN, i wonder which represent the CONVECTIVE HEATING RATE and which represent the CONDENSATION HEATING RATE.
Besides, I'm confused with the meaning of RTHSHTEN.
Appreciate very much if you can reply!
 
RTHCUTEN is represent the CONVECTIVE HEATING
RTHMPTEN is CONDENSATION HEATING RATE
RTHSHTEN is SHALLOW CONVECTIVE HEATING
 
Hi Ming!
I didnot find "R"THCUTEN in my output files. However, there are "A"THCUTEN, "A"THMPTEN and "A"THSHTEN .
Is there a spell error?
I check the value of the 3 variables. Surprisingly, the ATHCUTEN and ATHSHTEN are all 0 during the entire simulation. Only ATHMPTEN is not zero. Can you tell me the reason?
Thank you very much! I upload my namelist file below:

&time_control
run_days = 6,
run_hours = 1,
run_minutes = 0,
run_seconds = 0,
start_year = 2024,
start_month = 01,
start_day =21,
start_hour = 0,
start_minute = 00,
start_second = 00,
history_interval = 60,
frames_per_outfile = 1,
restart = .true.,
restart_interval = 2880,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
/

&domains
time_step = 10,
time_step_fract_num = 0,
time_step_fract_den = 1,





max_dom = 1,
s_we = 1,
e_we = 385,
s_sn = 1,
e_sn = 385,
s_vert = 1,
e_vert = 50,
dx = 3000,
dy = 3000,
ztop = 28000.,
eta_levels = 1.000, 0.995, 0.990, 0.983, 0.975,
0.965, 0.953, 0.938, 0.920, 0.898,
0.876, 0.854, 0.832, 0.810, 0.788,
0.766, 0.744, 0.722, 0.700, 0.678,
0.656, 0.634, 0.612, 0.590, 0.570,
0.550, 0.530, 0.510, 0.490, 0.470,
0.450, 0.430, 0.410, 0.390, 0.370,
0.350, 0.330, 0.310, 0.290, 0.270,
0.250, 0.225, 0.195, 0.160, 0.120,
0.070, 0.035, 0.015, 0.005, 0.000,

/

&physics
sf_ocean_physics=1,
oml_hml0=20,
oml_relaxation_time=1296000,
oml_gamma=0.14,
acc_phy_tend = 1,
mp_physics = 6,
ra_lw_physics = 4,
ra_sw_physics = 4,
radt = 3,
sf_sfclay_physics = 1,
sf_surface_physics = 1,
bl_pbl_physics = 1,
bldt = 0,
cu_physics = 0,
cudt = 0,
bucket_mm =10.,
bucket_J =1.e8,
/

&fdda
/

&dynamics
hybrid_opt = 0,
rk_ord = 3,
diff_opt = 1, 1, 1,
km_opt = 4, 4, 4,
damp_opt = 3,
dampcoef = .2,
zdamp = 13000.,
w_damping=1,
khdif = 0,
kvdif = 0,
smdiv = 0.1,
emdiv = 0.01,
epssm = 0.4,
time_step_sound = 6,
h_mom_adv_order = 5,
v_mom_adv_order = 3,
h_sca_adv_order = 5,
v_sca_adv_order = 3,
mix_full_fields = .true.,
non_hydrostatic = .true.,
pert_coriolis = .true.,
/

&bdy_control
periodic_x = .true.,
symmetric_xs = .false.,
symmetric_xe = .false.,
open_xs = .false.,
open_xe = .false.,
periodic_y = .true.,
symmetric_ys = .false.,
symmetric_ye = .false.,
open_ys = .false.,
open_ye = .false.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
&ideal
ideal_case = 3
/
 
Back
Top