Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

How to add advection and tendency variables to wrfout file in WRFV4.4

yliu

New member
Hi,
I am new to WRF, I have some questions about how to add advection and tendency variables to wrfout file.
I found these variables in Registry.EM_COMMON, and I have read the "How to Use the WRF Registry". I notice that the I/O of these variables is "h", but there is no these variables in my wrfout.nc file.
How can I output these variables to wrfout?

Best wishes,
Yan Liu

The fig below show the variables in the Registry.
registry.PNG
registry_1.PNG

and below is my namelist.input
&time_control
run_days = 0,
run_hours = 24,
run_minutes = 0,
run_seconds = 0,
start_year = 2016, 2016, 2016,
start_month = 06, 06, 06,
start_day = 01, 01, 01,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2016, 2016, 2016,
end_month = 06, 06, 06,
end_day = 02, 02, 02,
end_hour = 00, 00, 00,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 10800,
input_from_file = .true., .true., .true.,
history_interval = 60, 60, 60,
frames_per_outfile = 1000, 1000, 1000,
restart = .false.,
restart_interval = 5000,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
debug_level = 0,
/

&domains
time_step = 162,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 45, 31, 19,
e_sn = 60, 46, 34,
e_vert = 45, 45, 45,
p_top_requested = 5000,
num_metgrid_levels = 32,
num_metgrid_soil_levels = 4,
dx = 27000, 9000, 3000,
dy = 27000, 9000, 3000,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 18, 13,
j_parent_start = 1, 23, 18,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 0,
smooth_option = 0,
/

&physics
mp_physics = 3, 3, 3,
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 1, 1, 1,
radt = 27, 27, 27,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 1, 1, 0,
cudt = 5, 5, 5,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0, 0, 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
/

&fdda
/

&dynamics
w_damping = 0,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.,
damp_opt = 0,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2,
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 1, 1, 1,
scalar_adv_opt = 1, 1, 1,
/

&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false., .false.,
nested = .false., .true., .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
 
Hi, Yan,
What you list in Registry are the forcing for WRF run.I suppose you are looking for tendencies from WRF physics/dynamics. There are various variables that represent the tendencies of variables like U, V, Theta, etc. Below is just an example of that tendencies:
"RTHCUTEN": this is from cumulus scheme
"RTHSHTEN": this is from shallow convection
"RTHRATEN": this is from radaition
"RTHBLTEN": this is from boundary layer
"RTHFTEN" is the advective tendency

Please look at Registry.EM_COMMON and find the tendency terms you are looking for.
 
Hi, Yan,
What you list in Registry are the forcing for WRF run.I suppose you are looking for tendencies from WRF physics/dynamics. There are various variables that represent the tendencies of variables like U, V, Theta, etc. Below is just an example of that tendencies:
"RTHCUTEN": this is from cumulus scheme
"RTHSHTEN": this is from shallow convection
"RTHRATEN": this is from radaition
"RTHBLTEN": this is from boundary layer
"RTHFTEN" is the advective tendency

Please look at Registry.EM_COMMON and find the tendency terms you are looking for.
Thanks for helping me figure out that these variables are used to force WRF. I flowed your suggest and succeed to output tendencies as you listed. But I confused to get these forcing variables output from WRF em_real to as input for forcing WRF SCM. Can you give me some advices to get these forcing variables from WRF em_real.
 
Top