WRF 3.7计算过程中，异常中断，请各位大神看看

Torey

用WRF 3.7.1进行气象数据计算处理过程中，会出现异常中断，如下图所示：

error

Namelist.input文件如下：
&time_control
run_days                            = 0,
run_hours                           = 72,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          =  2015, 2015, 2015,,
start_month                         =  12, 12, 12,,
start_day                           =  21, 21, 21,,
start_hour                          =  00, 00, 00,,
start_minute                        =  00, 00, 00,,
start_second                        =  00, 00, 00,,
end_year                            =  2015, 2015, 2015,,
end_month                           =  12, 12, 12,,
end_day                             =  24, 24, 24,,
end_hour                            =  0, 0, 0,,
end_second                          =  0, 0, 0,,
interval_seconds                    =  10800,
input_from_file                     =  .true., .true., .true.,,
history_interval                    =  60,  60,   60,
frames_per_outfile                  =  1,    1,    1,
restart                             = .false.,
restart_interval                    = 7200,
reset_simulation_start              = F,
io_form_input                       = 2,
io_form_history                     = 2,
io_form_restart                     = 2,
io_form_boundary                    = 2,
io_form_auxinput4                   = 2,
auxinput4_inname                    = "wrflowinp_d<domain>",
auxinput4_interval                  = 180,
io_form_auxinput4                   = 2,
write_hist_at_0h_rst                = .false.,
debug_level                         = 0,
/

&fdda
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

&domains
time_step                           = 150,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 3,
e_we                                =  485,547,547,
e_sn                                =  555,544,913,
e_vert                              =   45,      45,    45,
eta_levels                          = 1.000, 0.999, 0.998, 0.997, 0.994,
                                       0.992, 0.989, 0.985, 0.981, 0.976,
                                       0.971, 0.965, 0.958, 0.952, 0.944,
                                       0.936, 0.927, 0.918, 0.908, 0.898,
                                       0.888, 0.877, 0.865, 0.853, 0.841,
                                       0.828, 0.814, 0.800, 0.783, 0.763,
                                       0.733, 0.694, 0.646, 0.591, 0.530,
                                       0.464, 0.396, 0.327, 0.259, 0.195,
                                       0.138, 0.090, 0.052, 0.022, 0.000,
p_top_requested                     = 5000,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
dx                                  =   18000.,  6000.,  2000.,
dy                                  =   18000.,  6000.,  2000.,

grid_id                             =  1, 2, 3,,
parent_id                           =  0, 1, 2,,
i_parent_start                      = 1,131,170,,
j_parent_start                      = 1,190,100,,
parent_grid_ratio                   = 1,3,3,,
parent_time_step_ratio              = 1,3,3,,
feedback                            = 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                   = 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,
/

&dynamics
w_damping                           = 0,
diff_opt                            = 1,      1,      1,
km_opt                              = 4,      4,      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.,
/


尝试把time_step参数减小，更改到120/90/60都会异常中断，但修改到30后，就可以正常计算完成了
有没有哪位大神知道为什么time_step在30以上都会导致计算异常终止，有没有方法可以在time_step为150的情况下，保证计算可以顺利完成?
感谢大家了

andrewsoong

18km X 6 =多少？？？？不是150 啊

Torey

andrewsoong 发表于 2016-3-3 09:48
18km X 6 =多少？？？？不是150 啊

我试过用timestep 为 108、54来测试，同样会异常中断
在36的时候可以正常计算完成，不知道是不是time step和其他参数不匹配导致的，time step只由dx来决定吗？

andrewsoong

Torey 发表于 2016-3-3 10:33
我试过用timestep 为 108、54来测试，同样会异常中断
在36的时候可以正常计算完成，不知道是不是time st ...

一般来说是dx的公里数的6倍，小于6倍，如4倍也可以。

wong

我也很好奇，坐等大神来回答

Lighting

可能是由于内存不足导致的。看你采用的是三层嵌套，你可以尝试使用两层嵌套运行试一下。你可以减小积分步长试一下。如果是在大型机上提交作业，可以更改一下节点数和核数试试，具体的每个节点的核数要求可以询问管理员。

星雨

积分步长取dx的6倍是WRF手册提供的参考值，但根据我的经验，这个值通常过大，适当的取小一点。比如我dx为18km，积分步长我只取90s。还有，你知不知道你模拟区域的数据量会有多大，我最内层dx=2km、垂直42层、网格数约400×400时，只模拟18h，d03的结果就有几十个G，你500*1000的格点，还积分3天，数据量得多少。
不知道你刚开始用WRF还是模拟区域设置确实要这么大，本人感觉你这样设置太浪费资源。
关于错误本身，我看你的设置下，根本没积分几分钟，重设置参数，减小积分步长试试吧。

fhmny

楼主的DOMIN确实牛逼，这覆盖了有半个地球了。
这么大的区域，Boundary是否有效呢。