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

00z variables = 0 in everywhere in datastreams in restarts

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,
/
 
Hi,
I believe this may be due to the way the I/O is listed for these two variables in the WRF/Registry/Registry.EM_COMMON file. All the others include 'r' in the I/O column. For e.g., for RAINNC:

Code:
state  real  RAINNC   ij  misc    1   -   rh01du   "RAINNC"  "ACCUMULATED TOTAL GRID SCALE PRECIPITATION"              "mm"

Above, the 8th column (right after the "-" is the I/O. The 'r' in that string indicates it should be written to the restart file. But the 'r' is missing from Z and REFL_10CM I/O. I'm not sure the specific reason the code was written this way. If you want to try modifying it, you can add an 'r' to the I/O for those variables in the registry file. However, you will need to recompile the code, and because it's a registry change, you'll need to clean the code (./clean -a) and reconfigure before recompiling to incorporate those changes. Doing this won't remove any files you've generated. After it's built, you can try running the initial run again so that those variables will be included in the wrfrst* files, and then try running the restart again. If, however, you think the values of Z and REFL_10CM seem correct, even though they aren't included at the initial restart time, you can continue without doing this.
 
Top