wrf-chem：wrf.exe后报错

央哥儿

做了两层domain, wrf.exe后报错如下：一开始就断了
*********************************************************************
*             PROGRAM:WRF/CHEM V3.6.1 MODEL
*                                                                   *
*    PLEASE REPORT ANY BUGS TO WRF/CHEM HELP at                     *
*                                                                   *
*              wrfchemhelp.gsd@noaa.gov                             *
*                                                                   *
*********************************************************************
  call madronich phot initialization
*************************************
Nesting domain
ids,ide,jds,jde            1          79           1          61
ims,ime,jms,jme           30          84          -4          42
ips,ipe,jps,jpe           40          79           1          30
INTERMEDIATE domain
ids,ide,jds,jde           35          66          11          36
ims,ime,jms,jme           40          71           6          32
ips,ipe,jps,jpe           50          68           9          22
*************************************
d01 2016-09-26_00:00:00  alloc_space_field: domain            2,              289538584 bytes allocated
d01 2016-09-26_00:00:00  alloc_space_field: domain            2,               10941696 bytes allocated
d01 2016-09-26_00:00:00 *** Initializing nest domain # 2 from an input file. ***
d01 2016-09-26_00:00:00 med_initialdata_input: calling input_input
INPUT LandUse = "USGS"
LANDUSE TYPE = "USGS" FOUND          33 CATEGORIES           2 SEASONS WATER CATEGORY =           16 SNOW CATEGORY =           24
forrtl: severe (174): SIGSEGV, segmentation fault occurred
调整积分步长没有解决，将max_dom调为1可以顺利跑出wrfout_d01。实在找不出解决办法，求助大家！
namelist.input:
&time_control
run_days                            = 0,
run_hours                           = 24,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2016, 2016, 2000,
start_month                         = 09,   09,   01,
start_day                           = 26,   26,   24,
start_hour                          = 00,   00,   12,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2016, 2016, 2000,
end_month                           = 09,   09,   01,
end_day                             = 27,   27,   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                  = 24, 24, 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                = 1440, 1440, 60,
auxinput7_interval_m                = 1440, 1440, 60,
auxinput8_interval_m                = 1440, 1440, 60,
auxinput13_interval_m               = 1440, 1440, 60,
io_form_auxinput5                   = 2,
io_form_auxinput6                   = 2,

io_form_auxinput7                   = 2,
io_form_auxinput8                   = 2,
io_form_auxinput12                  = 2,
io_form_auxinput13                  = 0,
debug_level                         = 0,
auxinput1_inname                    = "met_em.d<domain>.<date>",
/

&dfi_control
/

&domains
time_step                           = 270,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 2,
e_we                                = 91,   79,
e_sn                                = 73,   61,
e_vert                              = 30,    30,    28,
num_metgrid_levels                  = 32,
num_metgrid_soil_levels             = 0,
dx                                  = 45000,  15000,  3333.33,
dy                                  = 45000,  15000,  3333.33,
grid_id                             = 1,     2,     3,
parent_id                           = 1,     1,     2,
i_parent_start                      = 1,     37,    30,
j_parent_start                      = 1,     13,    30,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     3,     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                          = 0,     0,     2,
progn                               = 0,
ra_lw_physics                       = 1,     1,     1,
ra_sw_physics                       = 5,     5,     5,
radt                                = 30,    30,    30,
sf_sfclay_physics                   = 1,     1,     1,
sf_surface_physics                  = 1,     1,     1,
bl_pbl_physics                      = 7,     7,     7,
bldt                                = 0,     0,     0,
cu_physics                          = 3,     3,     3,
cu_diag                             = 1,     1,     0,
cudt                                = 5,     5,     5,
ishallow                            = 0,
isfflx                              = 1,
ifsnow                              = 0,
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                           = 0,
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,      2,      0,
scalar_adv_opt                      = 2,      2,      0,
chem_adv_opt                        = 2,      2,      0,
tke_adv_opt                         = 2,      2,      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.,.true.,.false.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

&chem
kemit                               = 1,
chem_opt                            = 301,
bioemdt                             = 30,
photdt                              = 30,
chemdt                              = 10,
io_style_emissions                  = 1,
emiss_inpt_opt                      = 1,          1,
emiss_opt                           = 5,          5,
emiss_opt_vol                       = 0,          0,
emiss_ash_hgt                       = 20000.,
chem_in_opt                         = 0,          0,
phot_opt                            = 1,          0,
gas_drydep_opt                      = 1,          1,
aer_drydep_opt                      = 1,          1,
bio_emiss_opt                       = 3,          3,
ne_area                             = 104,
depo_fact                           = 0.25,
dust_opt                            = 0,
dmsemis_opt                         = 1,
seas_opt                            = 0,
gas_bc_opt                          = 1,          1,
gas_ic_opt                          = 1,          1,
aer_bc_opt                          = 1,          1,
aer_ic_opt                          = 1,          1,
gaschem_onoff                       = 1,          1,
aerchem_onoff                       = 1,          1,
wetscav_onoff                       = 0,          0,
cldchem_onoff                       = 0,          0,
vertmix_onoff                       = 1,          1,

chem_conv_tr                        = 1,          1,
conv_tr_wetscav                     = 0,          0,
conv_tr_aqchem                      = 0,          0,
biomass_burn_opt                    = 1,          1,
plumerisefire_frq                   = 120,        120,
have_bcs_chem                       = .false., .false., .false.,
aer_ra_feedback                     = 1,
aer_op_opt                          = 1,
opt_pars_out                        = 1,
diagnostic_chem                     = 0,
/

问题少女

个人认为是你微物理参数化方案那一块设置的有问题，比如mp_physics=0，这样的，建议你仔细检查一下这一块。

央哥儿

问题少女 发表于 2017-7-27 07:25
个人认为是你微物理参数化方案那一块设置的有问题，比如mp_physics=0，这样的，建议你仔细检查一下这一块。

有考虑过调整mp_physics 0或者4 都会报同样错误

央哥儿

感觉这次报错有点让人摸不着头脑，并没有指向性的错误提示就中断了。运行记录到LANDUSE TYPE = "USGS" FOUND          33 CATEGORIES           2 SEASONS WATER CATEGORY =           16 SNOW CATEGORY =           24就停止了

绛紫

央哥儿 发表于 2017-7-27 10:17
有考虑过调整mp_physics 0或者4 都会报同样错误

domain里边num_metgrid_levels  比e_vert    值还高？不太科学~还有pyhsics里边的微物理方案设为0的貌似很少

绛紫

绛紫 发表于 2017-7-27 10:37
domain里边num_metgrid_levels  比e_vert    值还高？不太科学~还有pyhsics里边的微物理方案设为0的貌似 ...

还有积分步长可以考虑设为360

央哥儿

绛紫 发表于 2017-7-27 10:38
还有积分步长可以考虑设为360

谢谢回复~
物理参数应该是没有问题的，因为关掉d02可能正常跑d01。
问题可能出在d02嵌套上，因为convert_emiss.exe只默认识别XXX_d01，所以处理d02的清单时将nnamlist.input的column1改为d02的范围跑出wrfinput然后在去制作d02的清单，这样子其实很容易出错， 不知大家都是用什么方法做的domain2的清单？

绛紫

央哥儿 发表于 2017-7-27 11:11
谢谢回复~
物理参数应该是没有问题的，因为关掉d02可能正常跑d01。
问题可能出在d02嵌套上，因为conver ...

这个问题手册上写的很明白，每次转都是当成d01，但是real的时候跑的是3层的，convert转单独的每层的

绛紫

绛紫 发表于 2017-7-27 12:49
这个问题手册上写的很明白，每次转都是当成d01，但是real的时候跑的是3层的，convert转单独的每层的

我说的你的参数并不是说不对，是说你做模拟也好预报也好，这些物理参数一般很少有设置成0默认的，不同的区域可能物理参数不一样，但我相信你做过敏感性试验后就会发现，0和非0差别还是蛮大的

央哥儿

绛紫 发表于 2017-7-27 12:51
我说的你的参数并不是说不对，是说你做模拟也好预报也好，这些物理参数一般很少有设置成0默认的，不同的 ...

是的是的，你说的很有道理，确实是这样。我之前怀疑这块可能有问题所以就多试了几个微物理过程的参数，其中包括0，后来发现报错原因应该不是这里。
回过头检查感觉在单独转每层清单时有些地方不够明白，每次./convert_emiss.exe时namelist.input里i_parent_start直接改为domain2的参数吗，还是要乘以ratio?