[已解决]WRF运行结果后处理

短腿儿章鱼

请教各位大佬，我运行WRF模式一切顺利，没有报错，但是结果用ARWpost进行后处理就总有问题，之前模拟出来的文件后处理处理速度比较慢，每个时次都会在小屏幕上报出来，生成的ctl文件也可以用，但是最近模拟的三层嵌套，42小时，替换边界层方案的结果进行后处理速度特别快，嗖的一下就完事了，好像就处理了第一个时次，后边的直接跳过，但是也没有报错，甚至还提示成功了，但是再用grads提取数据时就只能提取第一个时次的，后边的全部都是替换值，模拟了2次了都是这样的，我也不知道wrf运行过程中哪块出问题了，请大家帮我看一下
这是namelist.input
time_control            
run_days = 0，
run_hours = 42，
run_minutes = 0，
run_seconds = 0，
start_year = 2017，
2017、2017 ，start_month = 7、7、7，
start_day =
7、7、7 ，start_hour = 0、0、0，
start_minute = 00、00、00，
start_second = 00，00，00，
end_year =
2017，2017，2017，
end_month =
7，7，7，
end_day =
8，8，8 ，end_hour = 18，18，18 ，end_minute = 00，00，00，end_second = 00、00、00，
interval_seconds = 21600，
input_from_file = .true。，。true。，。true。，
history_interval =
180、180、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 = 76，
time_step_fract_num = 0，
time_step_fract_den = 1
，
max_dom = 3，e_we =
100，88，76，
e_sn = 134，124，112 ，e_vert = 29， 29，29，
p_top_requested = 5000，
num_metgrid_levels = 32，
num_metgrid_soil_levels = 4，
dx = 25300，8433.333，2811.111，
dy = 25300，8433.333，2811.111，
grid_id = 1、2、3，
parent_id =
1、1、2，
i_parent_start =
1、36、32，
j_parent_start = 1、47、44 ，parent_grid_ratio = 1、3、3 ，parent_time_step_ratio = 1、3、3，
反馈= 1，
smooth_option = 0，
/

＆physics                  
mp_physics =
3、3、3，
ra_lw_physics =
1、1、1，
ra_sw_physics = 1、1、1 ，radt = 30、30、30 ，sf_sfclay_physics = 2、2、2
sf_surface_physics = 2，2，2，
bl_pbl_physics =
2，2，2，
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.，
wet_opt = 0，
zdamp = 5000.，5000.，5000.，
wetcoef =
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，
指定= .true。，。false。，。false。，
嵌套= .false。，。true。，。true。，
/

＆grib2                    
/ /

namelist_quilt           
nio_tasks_per_group = 0，
nio_groups = 1，
/

短腿儿章鱼

问题已解决，积分步长必须是数据输出时间的整数倍

短腿儿章鱼

这是模拟结果查询[root@localhost run]# tail rsl.out.0000
Timing for main: time 2017-07-08_17:59:57 on domain   3:    0.62131 elapsed seconds
Timing for main: time 2017-07-08_18:00:06 on domain   3:    0.61713 elapsed seconds
Timing for Writing wrfout_d03_2017-07-08_18:00:06 for domain        3:    0.06490 elapsed seconds
Timing for main: time 2017-07-08_18:00:14 on domain   3:    0.69141 elapsed seconds
Timing for main: time 2017-07-08_18:00:14 on domain   2:    2.95974 elapsed seconds
Timing for Writing wrfout_d02_2017-07-08_18:00:14 for domain        2:    0.08096 elapsed seconds
Timing for main: time 2017-07-08_18:00:40 on domain   2:    0.87850 elapsed seconds
Timing for main: time 2017-07-08_18:00:40 on domain   1:    8.05232 elapsed seconds
Timing for Writing wrfout_d01_2017-07-08_18:00:40 for domain        1:    0.13058 elapsed seconds
d01 2017-07-08_18:00:40 wrf: SUCCESS COMPLETE WRF
[root@localhost run]#

短腿儿章鱼

这个是namelist.ARWpost
&datetime
start_date = '2017-07-07_00:00:00',
end_date   = '2017-07-08_18:00:00',
interval_seconds = 3600,
tacc = 0,
debug_level = 0,
/

&io
input_root_name = '/opt/software/wrf-bins/wrf-test/MYJ/wrfout_d03_*'
output_root_name = '/opt/software/wrf-bins/wrf-test/1111'
plot = 'all_list'
fields = 'height,pressure,tk,tc,ws10,wspd,wd10,wdir,u10m,v10m'
mercator_defs = .true.
/
split_output = .true.
frames_per_outfile = 2



plot = 'all_list'
! Below is a list of all available diagnostics
fields = 'height,geopt,theta,tc,tk,td,td2,rh,rh2,umet,vmet,pressure,u10m,v10m,wdir,wspd,wd10,ws10,slp,mcape,mcin,lcl,lfc,cape,cin,dbz,max_dbz,clfr'


&interp
interp_method = 0,
interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
/
extrapolate = .true.

interp_method = 0,     ! 0 is model levels, -1 is nice height levels, 1 is user specified pressure/height

interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
interp_levels = 0.25, 0.50, 0.75, 1.00, 2.00, 3.00, 4.00, 5.00, 6.00, 7.00, 8.00, 9.00, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 17.0, 18.0, 19.0, 20.0,

短腿儿章鱼

这是运行arwpost的结果
[root@localhost ARWpost]# ./ARWpost.exe

!!!!!!!!!!!!!!!!
  ARWpost v3.1
!!!!!!!!!!!!!!!!

FOUND the following input files:
/opt/software/wrf-bins/wrf-test/MYJ/wrfout_d03_2017-07-07_00:00:00

START PROCESSING DATA

Processing  time --- 2017-07-07_00:00:00
   Found the right date - continue

Processing  time --- 2017-07-07_01:00:00
   Found  2017-07-08_18:00:06  before  2017-07-07_01:00:00

DONE Processing Data

CREATING .ctl file

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!  Successful completion of ARWpost  !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

[root@localhost ARWpost]#

短腿儿章鱼

坐等大佬前来提出解决方案

小其其格

Found  2017-07-08_18:00:06  before  2017-07-07_01:00:00

你转出来的这个文件，导致你的wrfout文件间隔不是3600秒了，所以后处理进行不下去
正确应该是wrfout_d03_2017-07-08_18:00:00

短腿儿章鱼

小其其格 发表于 2020-9-30 16:08
Found  2017-07-08_18:00:06  before  2017-07-07_01:00:00

你转出来的这个文件，导致你的wrfout文件间 ...

那请问我应该修改哪个文件呢，这和我WRF运行结果有关系吗，我两层嵌套就不会出现这个问题，但三层就老是这样，是我namelist里有啥地方没设置对吗？