TKE output field is all zeros, although the TKE_PBL variable seems to be outputting normal values.

[Key]

Hi all,

I'm encountering an issue when running WRF from version 4.0 onwards: my TKE output field is all zeros, although the TKE_PBL variable seems to be outputting normal values.

I’ve checked the following:

• The model was compiled correctly and runs without errors.
• Other turbulence-related variables seem fine.
• I've tried different PBL schemes (e.g., MYJ, MYNN), but the issue persists.
• History output includes TKE (TKE) and TKE_PBL (TKE_PBL), but only TKE remains zero in the output files.

From what I understand, TKE typically represents the total turbulent kinetic energy, while TKE_PBL refers to the PBL-related component. However, if TKE is always zero while TKE_PBL is non-zero, it suggests either a change in how these fields are computed or written in newer versions of WRF.

Has anyone else encountered this issue with WRF v4.0+?
Is there a configuration setting or namelist flag that needs to be adjusted to get proper TKE output?

Any help or insight would be greatly appreciated!

Thanks in advance,


The following is my namelist.input,:
&physics
mp_physics = 2,2,
ra_lw_physics = 1,1,
ra_sw_physics = 1,1,
sf_surface_physics = 1,1
cu_physics = 5,5,1
num_land_cat = 21,
num_soil_layers = 4,
surface_input_source = 1,
cudt = 5,5,
icloud = 1,
sst_update = 1,
sf_sfclay_physics = 2, 2,
bl_pbl_physics = 5, 5,
bl_mynn_tkebudget = 1, 1,
bl_mynn_tkeadvect =.true., .true.
bl_mynn_edmf_tke = 1,1

/

&fdda
/

&dynamics
hybrid_opt = 0,
etac = 0.15,
w_damping = 1,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0, 0,
diff_6th_factor = 0.12, 0.12,
base_temp = 290.,
damp_opt = 0,
zdamp = 5000., 5000.,
dampcoef = 0.2, 0.2,
khdif = 0, 0,
kvdif = 0, 0,
non_hydrostatic = .true., .true.,
moist_adv_opt = 2, 1,
scalar_adv_opt = 2, 1,
h_sca_adv_order = 5,
tke_adv_opt = 1,1
time_step_sound = 0,

 

[Key]

all right
I've found that qke is twice TKE ,and I've decided to use this variable.