WRF-Chem WRF run在第一个history_interval时停止

quwj

WRF-Chem WRF run在第一个history_interval时停止

刚接触WRF-Chem，第一次run的时候出现问题，总是在该write第一个时次模拟结果（wrfout文件）时终止

测试了一下，发现如果模拟时间段设短一点（比如6小时）可以正常运行结束

请各位江湖高手帮忙看看，有可能什么地方出现问题？非常感谢！


下面是namelist.input文件
&time_control
run_days                            = 20,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2018, 2010, 2000,
start_month                         = 12,   07,   01,
start_day                           = 28,   14,   24,
start_hour                          = 00,   00,   12,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2019, 2010, 2000,
end_month                           = 01,   07,   01,
end_day                             = 17,   16,   25,
end_hour                            = 00,   00,   12,
end_minute                          = 00,   00,   00,
end_second                          = 00,   00,   00,
interval_seconds                    = 21600,
input_from_file                     = .true.,.true.,.false.,
history_interval                    = 60,  60,   60,
frames_per_outfile                  = 1,       !1000, 1000, 1000,
restart                             = .false.,
restart_interval                    =    0,
io_form_history                     = 2,
io_form_restart                     = 2,
io_form_input                       = 2,
io_form_boundary                    = 2,
auxinput6_inname                    = 'wrfbiochemi_d01',
auxinput7_inname                    = 'wrffirechemi_d<domain>',
auxinput8_inname                    = 'wrfchemi_gocart_bg_d<domain>',
auxinput12_inname                   = 'wrf_chem_input',
auxinput13_inname                   = 'wrfchemv_d<domain>',
auxinput5_interval_m                = 2592000,         !86400, 86400, 60,
auxinput7_interval_m                = 2592000,         !86400, 86400, 60,
auxinput8_interval_m                = 2592000,         !86400, 86400, 60,
auxinput13_interval_m               = 2592000,         !86400, 86400, 60,
io_form_auxinput2                   = 2,
io_form_auxinput5                   = 2,
io_form_auxinput6                   = 0,
io_form_auxinput7                   = 2,
io_form_auxinput8                   = 2,
io_form_auxinput12                  = 2,
io_form_auxinput13                  = 0,
debug_level                         = 1000,    !300,
auxinput1_inname                    = "met_em.d<domain>.<date>",
/

&dfi_control
/

&domains
time_step                           = 600,     !600,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
s_we                                =   1,    1,     1,
e_we                                = 41,    41,    94,
s_sn                                =   1,    1,     1,
e_sn                                = 41,    41,    91,
e_vert                              = 32,    31,    28,
num_metgrid_levels                  = 32,
num_metgrid_soil_levels             = 4,
dx                                  = 111000, 4000,  3333.33,
dy                                  = 111000, 4000,  3333.33,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,     6,    30,
j_parent_start                      = 1,     65,    30,
parent_grid_ratio                   = 1,     5,     3,
parent_time_step_ratio              = 1,     5,     3,
p_top_requested                     = 5000,
feedback                            = 1,
smooth_option                       = 0
p_top_requested                     = 5000
zap_close_levels                    = 50
interp_type                         = 1
t_extrap_type                       = 2
force_sfc_in_vinterp                = 0
use_levels_below_ground             = .true.
use_surface                         = .true.
lagrange_order                      = 1
/
sfcp_to_sfcp                        = .true.

&physics
mp_physics                          = 4,     2,     2,
progn                               = 0,
ra_lw_physics                       = 1,     1,     1,
ra_sw_physics                       = 2,     1,     1,
radt                                = 30,    30,    30,
sf_sfclay_physics                   = 1,     1,     1,
sf_surface_physics                  = 2,     2,     1,
bl_pbl_physics                      = 1,     1,     1,
bldt                                = 0,     0,     0,
cu_physics                          = 5,     5,     0,
cu_diag                             = 1,     0,     0,
cudt                                = 0,     0,     0,
ishallow                            = 0,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
sf_urban_physics                    = 0,     0,     0,
mp_zero_out                         = 2,
mp_zero_out_thresh                  = 1.e-12
maxiens                             = 1,
maxens                              = 3,
maxens2                             = 3,
maxens3                             = 16,
ensdim                              = 144,
cu_rad_feedback                     = .true.,
/

&fdda
/

&dynamics
rk_ord                              = 3,
w_damping                           = 1,
diff_opt                            = 1,
km_opt                              = 4,
diff_6th_opt                        = 0,
diff_6th_factor                     = 0.12,
base_temp                           = 290.
damp_opt                            = 0,
zdamp                               = 5000.,  5000.,  5000.,
dampcoef                            = 0.01,   0.01,   0.01
khdif                               = 0,      0,      0,
kvdif                               = 0,      0,      0,
non_hydrostatic                     = .true., .true., .true.,
moist_adv_opt                       = 2,      0,      0,
scalar_adv_opt                      = 2,      0,      0,
chem_adv_opt                        = 2,      0,      0,
tke_adv_opt                         = 2,      0,      0,
time_step_sound                     = 4,      4,      4,
h_mom_adv_order                     = 5,      5,      5,
v_mom_adv_order                     = 3,      3,      3,
h_sca_adv_order                     = 5,      5,      5,
v_sca_adv_order                     = 3,      3,      3,
/

&bdy_control
spec_bdy_width                      = 5,
spec_zone                           = 1,
relax_zone                          = 4,
specified                           = .true., .false.,.false.,
nested                              = .false.,.false.,.false.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

&chem
kemit                               = 1,
chem_opt                            = 301,     !301,
bioemdt                             = 30,
photdt                              = 30,
chemdt                              = 10,
io_style_emissions                  = 2,
emiss_opt                           = 5,          0,
emiss_opt_vol                       = 0,          0,
emiss_ash_hgt                       = 20000.,
chem_in_opt                         = 0,          0,
phot_opt                            = 1,       !1,          0,
gas_drydep_opt                      = 1,          1,
aer_drydep_opt                      = 1,          1,
bio_emiss_opt                       = 1,          0,
ne_area                             = 1,
depo_fact                           = 0.25,
dust_opt                            = 1,
dmsemis_opt                         = 1,
seas_opt                            = 1,
gas_bc_opt                          = 1,          0,
gas_ic_opt                          = 1,          0,
aer_bc_opt                          = 1,          0,
aer_ic_opt                          = 1,          0,
gaschem_onoff                       = 1,          0,
aerchem_onoff                       = 1,          0,
wetscav_onoff                       = 0,          0,
cldchem_onoff                       = 0,          0,
vertmix_onoff                       = 1,          0,
chem_conv_tr                        = 1,          0,
conv_tr_wetscav                     = 0,          0,
conv_tr_aqchem                      = 0,          0,
biomass_burn_opt                    = 1,          0,
plumerisefire_frq                   = 30,         0,
have_bcs_chem                       = .false., .false., .false.,
aer_ra_feedback                     = 1,
aer_op_opt                          = 0,
opt_pars_out                        = 0,
diagnostic_chem                     = 0,       !0,
/


下面是rsl.error.0000文件最后几十行的输出

...

d01 2018-12-28_01:00:00 module_io.F: in wrf_write_field
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
output_wrf.b writing 2d real albedo Status =            0
Timing for wrf_ext_write_field:    0.00195 elapsed seconds
output_wrf.b writing 2d real clat
date 2018-12-28_01:00:00
ds            1           1           1
de           40          40           1
ps            1           1           1
pe           40          40           1
ms           -4          -4           1
me           46          46           1
d01 2018-12-28_01:00:00 module_io.F: in wrf_write_field
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
output_wrf.b writing 2d real clat Status =            0
Timing for wrf_ext_write_field:    0.00189 elapsed seconds
output_wrf.b writing 2d real albbck
date 2018-12-28_01:00:00
ds            1           1           1
de           40          40           1
ps            1           1           1
pe           40          40           1
ms           -4          -4           1
me           46          46           1
d01 2018-12-28_01:00:00 module_io.F: in wrf_write_field
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
output_wrf.b writing 2d real albbck Status =            0
Timing for wrf_ext_write_field:    0.00187 elapsed seconds
output_wrf.b writing 2d real emiss
date 2018-12-28_01:00:00
ds            1           1           1
de           40          40           1
ps            1           1           1
pe           40          40           1
ms           -4          -4           1
me           46          46           1
d01 2018-12-28_01:00:00 module_io.F: in wrf_write_field
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
output_wrf.b writing 2d real emiss Status =            0
Timing for wrf_ext_write_field:    0.00186 elapsed seconds
output_wrf.b writing 2d real noahres
date 2018-12-28_01:00:00
ds            1           1           1
de           40          40           1
ps            1           1           1
pe           40          40           1
ms           -4          -4           1
me           46          46           1
d01 2018-12-28_01:00:00 module_io.F: in wrf_write_field
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
d01 2018-12-28_01:00:00 module_io.F (md_calls.m4) : in wrf_put_var_ti_char_arr
output_wrf.b writing 2d real noahres Status =            0
Timing for wrf_ext_write_field:    0.00190 elapsed seconds
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source            
wrf.exe            0000000003B77A54  for__signal_handl     Unknown  Unknown
libpthread-2.17.s  00002AB3227295D0  Unknown               Unknown  Unknown
wrf.exe            00000000017DC8FC  Unknown               Unknown  Unknown
wrf.exe            0000000001788ABB  Unknown               Unknown  Unknown
wrf.exe            000000000188FB80  Unknown               Unknown  Unknown
wrf.exe            000000000187DD99  Unknown               Unknown  Unknown
wrf.exe            0000000001876B23  Unknown               Unknown  Unknown
wrf.exe            000000000051F506  Unknown               Unknown  Unknown
wrf.exe            000000000040D081  Unknown               Unknown  Unknown
wrf.exe            000000000040D03F  Unknown               Unknown  Unknown
wrf.exe            000000000040CFDE  Unknown               Unknown  Unknown
libc-2.17.so       00002AB322C5A3D5  __libc_start_main     Unknown  Unknown
wrf.exe            000000000040CEE9  Unknown               Unknown  Unknown


终于在ubuntu系统上装了中文输入法，可以正常发帖回帖了。

感谢气象家园的版主们长期维护给大家提供一个很好的学习平台！感谢各位同行朋友！Good luck！

quwj

感谢大家！用WRF3.9.1运行，总算能够正常计算完成了。

但是提高模式计算分辨率以后，计算又出现问题了，继续摸索中。。。谢谢！