Hi, I am trying to retrieve the vertical and horizontal eddy viscosity from WRF v3.8.1. FAQs in WRF user page (Q16, http://www2.mmm.ucar.edu/wrf/users/FAQ_files/FAQ_wrf_physics.html) suggest that ‘the XKMV and XKHV arrays are used by other PBL options (not MYJ)’. Then, I modify the Registry.EM_COMMON adding “h” on the line of xkmv and xkmh. The corresponding fields are found in the output after recompiling. However, the values of xkmv are zero at any position and any time while the values of xkmh seem to be valid. Does any one have an explanation for that or another way to obtain eddy viscosity? Part of the namelist.input is listed below:
Thanks in advance.
Code:
&physics
mp_physics = 6, 6, 10
ra_lw_physics = 3, 3, 1,
ra_sw_physics = 3, 3, 1,
radt = 30, 30, 3 ,
sf_sfclay_physics = 7, 7, 7,
sf_surface_physics = 7, 7, 7,
bl_pbl_physics = 7, 7, 7,
bldt = 0, 0, 0,
cu_physics = 0, 0, 0,
cudt = 0, 0, 0,
isfflx = 1
ifsnow = 0
icloud = 1
surface_input_source = 1
num_soil_layers = 2
sf_urban_physics = 0, 0, 0,
NUM_LAND_CAT = 20
/
&dynamics
w_damping = 1,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.25, 0.25, 0.25,
base_temp = 290.
damp_opt = 3,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.05, 0.05,
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 2, 1, 1,
scalar_adv_opt = 1, 0, 0,
/
Thanks in advance.