求助：WRF3.6.1 segmentation fault occurred

holly

    我在使用WRF3.6.1版本运行WRF的时候，出现了段错误，如下：d01 2008-06-27_01:56:00 microphysics_driver: calling milbrandt2mom
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source            
wrf.exe            000000000168C5B8  Unknown               Unknown  Unknown
wrf.exe            000000000168B65D  Unknown               Unknown  Unknown
wrf.exe            0000000001682E3C  Unknown               Unknown  Unknown
wrf.exe            000000000166B269  Unknown               Unknown  Unknown
wrf.exe            0000000001425058  Unknown               Unknown  Unknown
wrf.exe            0000000000FD619F  Unknown               Unknown  Unknown
wrf.exe            0000000000EF815D  Unknown               Unknown  Unknown
wrf.exe            000000000050337D  Unknown               Unknown  Unknown
wrf.exe            0000000000469F99  Unknown               Unknown  Unknown
wrf.exe            0000000000469B6B  Unknown               Unknown  Unknown
wrf.exe            0000000000469AFC  Unknown               Unknown  Unknown
libc.so.6          00000039C9A1ECDD  Unknown               Unknown  Unknown
wrf.exe            00000000004699F9  Unknown               Unknown  Unknown
   
我尝试了ulimit -s unlimited ，减小积分步长，甚至不同的驱动数据，都没有用。同时，我将微物理方案9改成其它方案，模式都能正常运行。下面是我的namelist，请大家多多指教。谢谢！
&time_control
run_days                            = 0,
run_hours                           = 48,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2008, 2008,
start_month                         = 06,   06,   
start_day                           = 27,   27,   
start_hour                          = 00,   00,   
start_minute                        = 00,   00,   
start_second                        = 00,   00,   
end_year                            = 2008, 2008,
end_month                           = 06,   06,   
end_day                             = 29,   29,   
end_hour                            = 00,   00,   
end_minute                          = 00,   00,   
end_second                          = 00,   00,   
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,
history_interval                    = 180,  60,   
frames_per_outfile                  = 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                         = 1000
iofields_filename = "my_file_d01.txt","my_file_d02.txt"
ignore_iofields_warning = .true.,
/

&domains
time_step                           = 60,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
e_we                                = 203, 181,
e_sn                                = 146, 154,
e_vert                              = 38, 38,  
dx                                  = 9000, 3000,
dy                                  = 9000, 3000,
grid_id                             = 1,     2,  
parent_id                           = 0,     1,   
i_parent_start                      = 1,     100,  
j_parent_start                      = 1,     50,  
parent_grid_ratio                   = 1,     3,
parent_time_step_ratio              = 1,     3,
feedback                            = 0,
smooth_option                       = 0
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
force_sfc_in_vinterp                = 3
p_top_requested                     = 5000,
interp_type                         = 1
lowest_lev_from_sfc                 = .false.
lagrange_order                      = 1
zap_close_levels                    = 500
/

&physics
mp_physics                          = 9,     9,
ra_lw_physics                       = 4 ,    4,
ra_sw_physics                       = 7 ,    7,
radt                                = 15,    15,
sf_sfclay_physics                   = 2 ,    2,     
sf_surface_physics                  = 7 ,    7,     
bl_pbl_physics                      = 2 ,    2,     
bldt                                = 0,     0,     
cu_physics                          = 0 ,    0,     
cudt                                = 0,     0,     
isfflx                              = 1,
ifsnow                              = 0,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 2,
mp_zero_out                         = 2,
mp_zero_out_thresh                  = 1.e-8,
maxiens                             = 1,
num_land_cat                        = 24,
sf_urban_physics                    = 1,     1,
rdmaxalb                            = .false.
/

&fdda
/

&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,   0.12,   
base_temp                           = 290.
damp_opt                            = 0,
zdamp                               = 5000.,  5000.,  
dampcoef                            = 0.0,    0.0,   
khdif                               = 0,      0,      
kvdif                               = 0,      0,      
smdiv                               = 0.1,    0.1,
emdiv                               = 0.01,   0.01,
epssm                               = 0.1,    0.1,
non_hydrostatic                     = .true., .true.,
h_mom_adv_order                     = 5,      5,
v_mom_adv_order                     = 3,      3,
h_sca_adv_order                     = 3,      3,
v_sca_adv_order                     = 2,      2,
time_step_sound                     = 4,      4,
moist_adv_opt                       = 1,      1,           
scalar_adv_opt                      = 1,      1,           
iso_temp                            = 200.,
/

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

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

奔波尔灞

啊哈　不知道楼主知道原因了么　　好像是３．６．１版本用微物理方案９　　就是有问题　　所以３．７改了这个ｂｕｇ　　但是3.7.1和3.7也是有ｂｕｇ的

andrewsoong

并不是所有的方案都可以组合的

holly

andrewsoong 发表于 2014-12-21 11:42
并不是所有的方案都可以组合的

你好，你是指我的参数化方案组合有问题？如果是这样，请问我要怎么确定问题出在哪个方案上？
事实上，我还转了几个参数化方案组合，微物理方案都选的9，结果都失败了，尽管最后的出错提示不同。

andrewsoong

holly 发表于 2014-12-21 11:55
你好，你是指我的参数化方案组合有问题？如果是这样，请问我要怎么确定问题出在哪个方案上？
事实上，我 ...

我的意思是并不是所有的方案都允许你组合使用

holly

andrewsoong 发表于 2014-12-21 12:06
我的意思是并不是所有的方案都允许你组合使用

哦，我的这些组合都是user‘s Guide上没有明确说明不能组合的.

无敌泡

有这样的要求，不能随意组合，举个例子：有的物理方案前面选了1，后面某个也只能选1

无敌泡

至少前后的方案要一致

吴阿好大人

请问前辈my_file_d01.txt是怎么写的？我照着手册写但是总是报错无法打开，谢谢！

clarawho

奔波尔灞 发表于 2015-10-23 10:52
啊哈　不知道楼主知道原因了么　　好像是３．６．１版本用微物理方案９　　就是有问题　　所以３．７改了这 ...

是吗？？