stormchasegenie
Member
All,
I am using auxiliary datastreams to more efficiently output downscaled variables, and I am running into a strange issue (WRF v4.6.1).
For some variables, and only on restarts, the 00z value is 0 everywhere, but only in the 00z time slice. Other times in the file look as they should. Here is my datastream specification:
+:h:14:T2,RAINNC,U10,V10,Q2
+:h:22:OLR,SWUPT
-:h:22:XLAT,XLONG,T2,RAINC,RAINNC,U10,V10,C4F,C4H,C3F,C3H,C2F,C2H,C1F,C1H,Q2
+:h:22:SWDOWN,SWDDIF,SWDDNI,SWDDIR,ETRAN,EDIR,GRDFLX,LH,FIRA,HFX,TSK,FSA,PSFC,SNOW,PBLH
+:h:13:U,V,W,T,Z,QVAPOR,QICE,CLDFRA,QCLOUD,QRAIN,P,P_HYD,REFL_10CM,TSLB,SMOIS
I have underlined and bolded the problematic variables. All other variables look fantastic!
Initially, I tried changing write_hist_at_0h_rst = .false. to .true., but the problem remains.
See below my most recent namelist.input.
Best, and thank you,
-Stefan
&time_control
run_days = 1,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2025, 2025, 2025, 2025
start_month = 3, 4, 4, 4,
start_day = 31, 1, 1, 1,
start_hour = 0, 0, 0, 0,
end_year = 2000, 2000, 2000, 2000
end_month = 8, 8, 8, 8,
end_day = 1, 1, 1, 1,
end_hour = 1, 0, 0, 0,
interval_seconds = 21600
input_from_file = .true.,.true.,.true.,.true.,
history_interval = 360, 360, 360, 360,
frames_per_outfile = 1, 1, 1, 1,
restart = .true.,
restart_interval = 28800,
write_hist_at_0h_rst = .true.
adjust_output_times = .true.,
override_restart_timers = .true.
io_form_history = 0
io_form_restart = 102
io_form_input = 2
io_form_boundary = 2
io_form_auxinput4 = 2
auxinput4_inname = "wrflowinp_d<domain>",
auxinput4_interval_m = 360, 360, 360, 360,
iofields_filename = "vario.txt",
ignore_iofields_warning = .true.,
auxhist14_outname = "wrf5min_d<domain>_<date>",
auxhist14_interval_m = 5
frames_per_auxhist14 = 288
io_form_auxhist14 = 2
auxhist13_outname = "wrf3d_d<domain>_<date>",
auxhist13_interval_m = 360
frames_per_auxhist13 = 4
io_form_auxhist13 = 2
auxhist22_outname = "wrf2d_d<domain>_<date>",
auxhist22_interval = 60, 60, 60, 60
frames_per_auxhist22 = 24, 24, 24, 24,
io_form_auxhist22 = 2
debug_level = 0,
/
&diags
z_lev_diags = 1,
num_z_levels = 1,
z_levels = -100
/
&domains
time_step = 120,
use_adaptive_time_step = .true.
step_to_output_time = .true.
target_cfl = 1.2
max_step_increase_pct = 5, 51
starting_time_step = -1
max_time_step = -1
min_time_step = 1
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
parent_id = 1,
parent_grid_ratio = 1,
i_parent_start = 1,
j_parent_start = 1,
e_we = 1130,
e_sn = 800,
e_vert = 51, 30, 30,
p_top_requested = 5000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
sfcp_to_sfcp = .true.,
dx = 6000
dy = 6000
grid_id = 1, 2, 3, 4,
parent_time_step_ratio = 1, 5, 3, 3,
feedback = 0,
smooth_option = 0
smooth_cg_topo = .true.,
/
&physics
mp_physics = 51, 51, 51, 51,
cu_physics = 0, 6, 0, 0,
ra_lw_physics = 4, 4, 4, 4,
ra_sw_physics = 4, 4, 4, 4,
ghg_input = 1,
bl_pbl_physics = 1, 1, 1, 1,
sf_sfclay_physics = 1, 1, 1, 1,
sf_surface_physics = 4, 4, 4, 4,
radt = 5, 9, 3, 3,
bldt = 0, 0, 0, 0,
cudt = 0, 0, 0, 0,
icloud = 1,
num_land_cat = 28,
sf_urban_physics = 0, 0, 0, 0,
sst_update = 1,
sst_skin = 1,
sf_lake_physics = 1, 1, 1, 1,
use_lakedepth = 1, 1, 1, 1,
/
&fdda
grid_fdda = 2
gfdda_inname = "wrffdda_d<domain>"
gfdda_interval_m = 360,
gfdda_end_h = 4800000,
io_form_gfdda = 2
fgdt = 0,
if_no_pbl_nudging_uv = 1
if_no_pbl_nudging_t = 1
if_no_pbl_nudging_q = 1
if_no_pbl_nudging_ph = 1
if_zfac_uv = 0
if_zfac_t = 0
if_zfac_q = 0
guv = 0.00003
gt = 0.00003
gq = 0.0
gph = 0.0
xwavenum = 4
ywavenum = 3
/
&dynamics
hybrid_opt = 0,
w_damping = 1,
diff_opt = 2, 1, 1, 1,
km_opt = 4, 4, 4, 4,
diff_6th_opt = 0, 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 3,
zdamp = 10000., 10000., 10000.,
dampcoef = 0.2, 0.2, 0.2, 0.2,
khdif = 0, 0, 0, 0.,
kvdif = 0, 0, 0, 0.,
non_hydrostatic = .true., .true., .true., .true.,
moist_adv_opt = 1, 1, 1, 1,
scalar_adv_opt = 1, 1, 1, 1,
gwd_opt = 0,
/
&noah_mp
opt_run = 5,
dveg = 5,
/
&bdy_control
spec_bdy_width = 20,
specified = .true.
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
I am using auxiliary datastreams to more efficiently output downscaled variables, and I am running into a strange issue (WRF v4.6.1).
For some variables, and only on restarts, the 00z value is 0 everywhere, but only in the 00z time slice. Other times in the file look as they should. Here is my datastream specification:
+:h:14:T2,RAINNC,U10,V10,Q2
+:h:22:OLR,SWUPT
-:h:22:XLAT,XLONG,T2,RAINC,RAINNC,U10,V10,C4F,C4H,C3F,C3H,C2F,C2H,C1F,C1H,Q2
+:h:22:SWDOWN,SWDDIF,SWDDNI,SWDDIR,ETRAN,EDIR,GRDFLX,LH,FIRA,HFX,TSK,FSA,PSFC,SNOW,PBLH
+:h:13:U,V,W,T,Z,QVAPOR,QICE,CLDFRA,QCLOUD,QRAIN,P,P_HYD,REFL_10CM,TSLB,SMOIS
I have underlined and bolded the problematic variables. All other variables look fantastic!
Initially, I tried changing write_hist_at_0h_rst = .false. to .true., but the problem remains.
See below my most recent namelist.input.
Best, and thank you,
-Stefan
&time_control
run_days = 1,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2025, 2025, 2025, 2025
start_month = 3, 4, 4, 4,
start_day = 31, 1, 1, 1,
start_hour = 0, 0, 0, 0,
end_year = 2000, 2000, 2000, 2000
end_month = 8, 8, 8, 8,
end_day = 1, 1, 1, 1,
end_hour = 1, 0, 0, 0,
interval_seconds = 21600
input_from_file = .true.,.true.,.true.,.true.,
history_interval = 360, 360, 360, 360,
frames_per_outfile = 1, 1, 1, 1,
restart = .true.,
restart_interval = 28800,
write_hist_at_0h_rst = .true.
adjust_output_times = .true.,
override_restart_timers = .true.
io_form_history = 0
io_form_restart = 102
io_form_input = 2
io_form_boundary = 2
io_form_auxinput4 = 2
auxinput4_inname = "wrflowinp_d<domain>",
auxinput4_interval_m = 360, 360, 360, 360,
iofields_filename = "vario.txt",
ignore_iofields_warning = .true.,
auxhist14_outname = "wrf5min_d<domain>_<date>",
auxhist14_interval_m = 5
frames_per_auxhist14 = 288
io_form_auxhist14 = 2
auxhist13_outname = "wrf3d_d<domain>_<date>",
auxhist13_interval_m = 360
frames_per_auxhist13 = 4
io_form_auxhist13 = 2
auxhist22_outname = "wrf2d_d<domain>_<date>",
auxhist22_interval = 60, 60, 60, 60
frames_per_auxhist22 = 24, 24, 24, 24,
io_form_auxhist22 = 2
debug_level = 0,
/
&diags
z_lev_diags = 1,
num_z_levels = 1,
z_levels = -100
/
&domains
time_step = 120,
use_adaptive_time_step = .true.
step_to_output_time = .true.
target_cfl = 1.2
max_step_increase_pct = 5, 51
starting_time_step = -1
max_time_step = -1
min_time_step = 1
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
parent_id = 1,
parent_grid_ratio = 1,
i_parent_start = 1,
j_parent_start = 1,
e_we = 1130,
e_sn = 800,
e_vert = 51, 30, 30,
p_top_requested = 5000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
sfcp_to_sfcp = .true.,
dx = 6000
dy = 6000
grid_id = 1, 2, 3, 4,
parent_time_step_ratio = 1, 5, 3, 3,
feedback = 0,
smooth_option = 0
smooth_cg_topo = .true.,
/
&physics
mp_physics = 51, 51, 51, 51,
cu_physics = 0, 6, 0, 0,
ra_lw_physics = 4, 4, 4, 4,
ra_sw_physics = 4, 4, 4, 4,
ghg_input = 1,
bl_pbl_physics = 1, 1, 1, 1,
sf_sfclay_physics = 1, 1, 1, 1,
sf_surface_physics = 4, 4, 4, 4,
radt = 5, 9, 3, 3,
bldt = 0, 0, 0, 0,
cudt = 0, 0, 0, 0,
icloud = 1,
num_land_cat = 28,
sf_urban_physics = 0, 0, 0, 0,
sst_update = 1,
sst_skin = 1,
sf_lake_physics = 1, 1, 1, 1,
use_lakedepth = 1, 1, 1, 1,
/
&fdda
grid_fdda = 2
gfdda_inname = "wrffdda_d<domain>"
gfdda_interval_m = 360,
gfdda_end_h = 4800000,
io_form_gfdda = 2
fgdt = 0,
if_no_pbl_nudging_uv = 1
if_no_pbl_nudging_t = 1
if_no_pbl_nudging_q = 1
if_no_pbl_nudging_ph = 1
if_zfac_uv = 0
if_zfac_t = 0
if_zfac_q = 0
guv = 0.00003
gt = 0.00003
gq = 0.0
gph = 0.0
xwavenum = 4
ywavenum = 3
/
&dynamics
hybrid_opt = 0,
w_damping = 1,
diff_opt = 2, 1, 1, 1,
km_opt = 4, 4, 4, 4,
diff_6th_opt = 0, 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 3,
zdamp = 10000., 10000., 10000.,
dampcoef = 0.2, 0.2, 0.2, 0.2,
khdif = 0, 0, 0, 0.,
kvdif = 0, 0, 0, 0.,
non_hydrostatic = .true., .true., .true., .true.,
moist_adv_opt = 1, 1, 1, 1,
scalar_adv_opt = 1, 1, 1, 1,
gwd_opt = 0,
/
&noah_mp
opt_run = 5,
dveg = 5,
/
&bdy_control
spec_bdy_width = 20,
specified = .true.
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/