Dear all,
When I used the "iofields_filename= 'iofields_d01.txt', 'iofields_d02.txt'" in namelist.input to specify some variables output in a WRF restart running, and
the "iofields_d01.txt" is:
"# fields not required in wrfout files
-:h:0:HFX_FORCE,LH_FORCE,TSK_FORCE,HFX_FORCE_TEND,LH_FORCE_TEND,TSK_FORCE_TEND,THIS_IS_AN_IDEAL_RUN
"
the "iofields_d02.txt" is:
# wrfhrly fields
+:h:4:XLAT,XLONG,HAILCAST_DIAM_MAX
# fields not required in wrfout files
-:h:0:HFX_FORCE,LH_FORCE,TSK_FORCE,HFX_FORCE_TEND,LH_FORCE_TEND,TSK_FORCE_TEND,THIS_IS_AN_IDEAL_RUN
####--------------------------------
and the output stream 4 setting in namelist.input is :
auxhist4_outname = "wrfhrly_hail_d<domain>_<date>"
io_form_auxhist4 = 2
auxhist4_interval = 60, 60,
frames_per_auxhist4 = 24, 24,
###-----------------------------------
But the outputs of HAILCAST_DIAM_MAX from wrfhrly_hail_d02_<date> are different so much from the default outputs "wrfout_d02_<date>"
Seem that the values of HAILCAST_DIAM_MAX in the specifying files wrfhrly_hail_d02_<date> are the max values since the files started time, for example, the output file wrfhrly_hail_d02_2020-01-06_00:00 contains 24 values for 24 hours (2020-01-06_01:00:00~2020-01-07_00:00:00), and these values here represent the max values since 2020-01-06_01:00:00, but not the max values during the one hour it was produced.
While the default wrfout file, named wrfout_d02_2020-01-06_00:00:00, it produced the 24 values that correctly represented the max value during the one hour it was produced.
So I am wondering why it could happen ?
and if I wanted to output HAILCAST_DIAM_MAX only and improve the output efficiency in this restarted WRF simulation, what can I do ?
Some additional information：
The WRF version here is 4.3, and the related namelist.input are shown here:
&time_control
run_days = 6,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2020, 2020,
start_month = 1, 1
start_day = 1, 1,
start_hour = 0, 0,
start_minute = 0, 0,
start_second = 0, 0,
end_year = 2020, 2020,
end_month = 1, 1
end_day = 7, 7
end_hour = 0, 0,
end_minute = 0, 0,
end_second = 0, 0,
interval_seconds = 3600,
input_from_file = .True., .True.,
history_interval = 60, 60,
history_outname = "wrfout_d<domain>_<date>"
frames_per_outfile = 24, 24,
restart = .True.,
restart_interval = 7200,!!!
override_restart_timers = .true.,
write_hist_at_0h_rst = .False.,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
debug_level = 0,
!!!output_diagnostics = 1,
auxinput4_inname = "wrflowinp_d<domain>"
auxinput4_interval = 60, 60,
io_form_auxinput4 = 2,
!!!auxhist3_outname = "wrfxtrm_d<domain>_<date>"
!!!io_form_auxhist3 = 2
!!!auxhist3_interval = 1440, 1440,
!!!frames_per_auxhist3 = 1, 1,
auxhist4_outname = "wrfhrly_hail_d<domain>_<date>"
io_form_auxhist4 = 2
auxhist4_interval = 60, 60,
frames_per_auxhist4 = 24, 24,
!!!auxhist3_outname = "wrfdly_d<domain>_<date>"
!!!io_form_auxhist3 = 2
!!!auxhist3_interval = 1440, 1440,
!!!frames_per_auxhist3 = 1, 1,
iofields_filename = 'iofields_d01.txt', 'iofields_d02.txt'
use_netcdf_classic = .True.
/
&diagnostics
max_window = 302
max_rain_5m = 1,1
max_rain_10m = 1,1
max_rain_20m = 1,1
max_rain_30m = 1,1
max_rain_60m = 1,1
max_wind_5m = 1,1
max_wind_10m = 1,1
max_wind_20m = 1,1
max_wind_30m = 1,1
max_wind_60m = 1,1
/
&domains
time_step = 60,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
...
....
&physics
hailcast_opt = 0,1,
haildt = 0,0,
mp_physics = 8, 8,
ra_lw_physics = 4, 4,
ra_sw_physics = 4, 4,
o3input = 2
aer_opt = 1
radt = 10, 10,
cam_abs_freq_s = 10800
levsiz = 59
paerlev = 29
cam_abs_dim1 = 4
cam_abs_dim2 = 45
sf_sfclay_physics = 1, 1,
sf_surface_physics = 4, 4,
bl_pbl_physics = 7, 7,
bldt = 0, 0,
cu_physics = 6, 0,
cudt = 0, 0,
isfflx = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 1, 1,
sst_update = 1,
tmn_update = 1,
lagday = 150,
sst_skin = 1,
usemonalb = .true.,
rdmaxalb = .true.,
slope_rad = 1,
topo_shading = 1,
shadlen = 25000.,
bucket_mm = 1000,
prec_acc_dt = 60, 60,
/
&noah_mp
dveg = 2
opt_crs = 1
opt_sfc = 1
opt_btr = 1
opt_run = 3
opt_frz = 1
opt_inf = 1
opt_rad = 3
opt_alb = 2
opt_snf = 1
opt_tbot = 2
opt_stc = 1
opt_gla = 1
opt_rsf = 1
opt_soil = 1
opt_pedo = 1
opt_crop = 0
/
&dynamics
rk_ord = 3,
w_damping = 1,
diff_opt = 1, 1,
km_opt = 4, 4,
diff_6th_opt = 0, 0,
diff_6th_factor = 0.12,
base_temp = 290.
damp_opt = 1,
zdamp = 5000.0, 5000.0,
dampcoef = 0.01, 0.01,
khdif = 0, 0,
kvdif = 0, 0,
non_hydrostatic = .True., .True.,
moist_adv_opt = 1, 1,
scalar_adv_opt = 1, 1,
gwd_opt = 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .True., .False.,
nested = .False., .True.,
/
&namelist_quilt
poll_servers= .False.,
nio_tasks_per_group = 0,
nio_groups = 0,
/
When I used the "iofields_filename= 'iofields_d01.txt', 'iofields_d02.txt'" in namelist.input to specify some variables output in a WRF restart running, and
the "iofields_d01.txt" is:
"# fields not required in wrfout files
-:h:0:HFX_FORCE,LH_FORCE,TSK_FORCE,HFX_FORCE_TEND,LH_FORCE_TEND,TSK_FORCE_TEND,THIS_IS_AN_IDEAL_RUN
"
the "iofields_d02.txt" is:
# wrfhrly fields
+:h:4:XLAT,XLONG,HAILCAST_DIAM_MAX
# fields not required in wrfout files
-:h:0:HFX_FORCE,LH_FORCE,TSK_FORCE,HFX_FORCE_TEND,LH_FORCE_TEND,TSK_FORCE_TEND,THIS_IS_AN_IDEAL_RUN
####--------------------------------
and the output stream 4 setting in namelist.input is :
auxhist4_outname = "wrfhrly_hail_d<domain>_<date>"
io_form_auxhist4 = 2
auxhist4_interval = 60, 60,
frames_per_auxhist4 = 24, 24,
###-----------------------------------
But the outputs of HAILCAST_DIAM_MAX from wrfhrly_hail_d02_<date> are different so much from the default outputs "wrfout_d02_<date>"
Seem that the values of HAILCAST_DIAM_MAX in the specifying files wrfhrly_hail_d02_<date> are the max values since the files started time, for example, the output file wrfhrly_hail_d02_2020-01-06_00:00 contains 24 values for 24 hours (2020-01-06_01:00:00~2020-01-07_00:00:00), and these values here represent the max values since 2020-01-06_01:00:00, but not the max values during the one hour it was produced.
While the default wrfout file, named wrfout_d02_2020-01-06_00:00:00, it produced the 24 values that correctly represented the max value during the one hour it was produced.
So I am wondering why it could happen ?
and if I wanted to output HAILCAST_DIAM_MAX only and improve the output efficiency in this restarted WRF simulation, what can I do ?
Some additional information：
The WRF version here is 4.3, and the related namelist.input are shown here:
&time_control
run_days = 6,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2020, 2020,
start_month = 1, 1
start_day = 1, 1,
start_hour = 0, 0,
start_minute = 0, 0,
start_second = 0, 0,
end_year = 2020, 2020,
end_month = 1, 1
end_day = 7, 7
end_hour = 0, 0,
end_minute = 0, 0,
end_second = 0, 0,
interval_seconds = 3600,
input_from_file = .True., .True.,
history_interval = 60, 60,
history_outname = "wrfout_d<domain>_<date>"
frames_per_outfile = 24, 24,
restart = .True.,
restart_interval = 7200,!!!
override_restart_timers = .true.,
write_hist_at_0h_rst = .False.,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
debug_level = 0,
!!!output_diagnostics = 1,
auxinput4_inname = "wrflowinp_d<domain>"
auxinput4_interval = 60, 60,
io_form_auxinput4 = 2,
!!!auxhist3_outname = "wrfxtrm_d<domain>_<date>"
!!!io_form_auxhist3 = 2
!!!auxhist3_interval = 1440, 1440,
!!!frames_per_auxhist3 = 1, 1,
auxhist4_outname = "wrfhrly_hail_d<domain>_<date>"
io_form_auxhist4 = 2
auxhist4_interval = 60, 60,
frames_per_auxhist4 = 24, 24,
!!!auxhist3_outname = "wrfdly_d<domain>_<date>"
!!!io_form_auxhist3 = 2
!!!auxhist3_interval = 1440, 1440,
!!!frames_per_auxhist3 = 1, 1,
iofields_filename = 'iofields_d01.txt', 'iofields_d02.txt'
use_netcdf_classic = .True.
/
&diagnostics
max_window = 302
max_rain_5m = 1,1
max_rain_10m = 1,1
max_rain_20m = 1,1
max_rain_30m = 1,1
max_rain_60m = 1,1
max_wind_5m = 1,1
max_wind_10m = 1,1
max_wind_20m = 1,1
max_wind_30m = 1,1
max_wind_60m = 1,1
/
&domains
time_step = 60,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
...
....
&physics
hailcast_opt = 0,1,
haildt = 0,0,
mp_physics = 8, 8,
ra_lw_physics = 4, 4,
ra_sw_physics = 4, 4,
o3input = 2
aer_opt = 1
radt = 10, 10,
cam_abs_freq_s = 10800
levsiz = 59
paerlev = 29
cam_abs_dim1 = 4
cam_abs_dim2 = 45
sf_sfclay_physics = 1, 1,
sf_surface_physics = 4, 4,
bl_pbl_physics = 7, 7,
bldt = 0, 0,
cu_physics = 6, 0,
cudt = 0, 0,
isfflx = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 1, 1,
sst_update = 1,
tmn_update = 1,
lagday = 150,
sst_skin = 1,
usemonalb = .true.,
rdmaxalb = .true.,
slope_rad = 1,
topo_shading = 1,
shadlen = 25000.,
bucket_mm = 1000,
prec_acc_dt = 60, 60,
/
&noah_mp
dveg = 2
opt_crs = 1
opt_sfc = 1
opt_btr = 1
opt_run = 3
opt_frz = 1
opt_inf = 1
opt_rad = 3
opt_alb = 2
opt_snf = 1
opt_tbot = 2
opt_stc = 1
opt_gla = 1
opt_rsf = 1
opt_soil = 1
opt_pedo = 1
opt_crop = 0
/
&dynamics
rk_ord = 3,
w_damping = 1,
diff_opt = 1, 1,
km_opt = 4, 4,
diff_6th_opt = 0, 0,
diff_6th_factor = 0.12,
base_temp = 290.
damp_opt = 1,
zdamp = 5000.0, 5000.0,
dampcoef = 0.01, 0.01,
khdif = 0, 0,
kvdif = 0, 0,
non_hydrostatic = .True., .True.,
moist_adv_opt = 1, 1,
scalar_adv_opt = 1, 1,
gwd_opt = 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .True., .False.,
nested = .False., .True.,
/
&namelist_quilt
poll_servers= .False.,
nio_tasks_per_group = 0,
nio_groups = 0,
/
