WRF模式执行wrf.exe的时候报Fortran 90 ：Segmentation fault 错误（已解决）

一颗梅子在南京

之前由于模拟的区域比较小，都能正常运行，现在想要模拟整个中国区域，修改了namelist之后，执行wrf.exe的时候报如下错误：

Timing for main: time 2016-12-20_00:18:00 on domain   1:   95.24945 elapsed seconds
Timing for main: time 2016-12-20_00:19:00 on domain   2:   10.15469 elapsed seconds
Timing for main: time 2016-12-20_00:20:00 on domain   2:   12.05055 elapsed seconds
Timing for main: time 2016-12-20_00:21:00 on domain   2:   10.24895 elapsed seconds
Timing for main: time 2016-12-20_00:21:00 on domain   1:   39.96389 elapsed seconds
Timing for main: time 2016-12-20_00:22:00 on domain   2:   10.19173 elapsed seconds
Timing for main: time 2016-12-20_00:23:00 on domain   2:   10.26275 elapsed seconds
Timing for main: time 2016-12-20_00:24:00 on domain   2:   10.24037 elapsed seconds
Timing for main: time 2016-12-20_00:24:00 on domain   1:   39.03205 elapsed seconds
d01 2016-12-20_00:24:00           99  points exceeded cfl=2 in domain d01 at time 2016-12-20_00:24:00 hours
d01 2016-12-20_00:24:00  MAX AT i,j,k:           32           5          24  vert_cfl,w,d(eta)=   2.66693401      -6.97709608       2.52925158E-02
d01 2016-12-20_00:24:00         1166  points exceeded cfl=2 in domain d01 at time 2016-12-20_00:24:00 hours
d01 2016-12-20_00:24:00  MAX AT i,j,k:           31           5          24  vert_cfl,w,d(eta)=   7.77810240      -177.663757       2.52925158E-02
d01 2016-12-20_00:24:00         2228  points exceeded cfl=2 in domain d01 at time 2016-12-20_00:24:00 hours
d01 2016-12-20_00:24:00  MAX AT i,j,k:           30           5          33  vert_cfl,w,d(eta)=   86390.5391       1006058.06       8.25513527E-03
Timing for main: time 2016-12-20_00:25:00 on domain   2:   12.09914 elapsed seconds
Timing for main: time 2016-12-20_00:26:00 on domain   2:   10.30128 elapsed seconds
Timing for main: time 2016-12-20_00:27:00 on domain   2:   10.30016 elapsed seconds
Timing for main: time 2016-12-20_00:27:00 on domain   1:   40.26282 elapsed seconds
d01 2016-12-20_00:27:00          330  points exceeded cfl=2 in domain d01 at time 2016-12-20_00:27:00 hours
d01 2016-12-20_00:27:00  MAX AT i,j,k:           26           1          18  vert_cfl,w,d(eta)=   800.372620      -18192522.0       4.31730747E-02
d01 2016-12-20_00:27:00         2030  points exceeded cfl=2 in domain d01 at time 2016-12-20_00:27:00 hours
d01 2016-12-20_00:27:00  MAX AT i,j,k:          252           5          24  vert_cfl,w,d(eta)=   8.70168591      -243.581421       2.52925158E-02
d01 2016-12-20_00:27:00         3511  points exceeded cfl=2 in domain d01 at time 2016-12-20_00:27:00 hours
d01 2016-12-20_00:27:00  MAX AT i,j,k:            7           4          34  vert_cfl,w,d(eta)=   93212.8516      -523278.344       7.25505641E-03
Timing for main: time 2016-12-20_00:28:00 on domain   2:   10.31857 elapsed seconds
Timing for main: time 2016-12-20_00:29:00 on domain   2:   10.34950 elapsed seconds
Timing for main: time 2016-12-20_00:30:00 on domain   2:   12.20212 elapsed seconds
Timing for main: time 2016-12-20_00:30:00 on domain   1:   41.27123 elapsed seconds

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
#0  0x7FD9FBBE2467
#1  0x7FD9FBBE2AAE
#2  0x7FD9FB0E924F
#3  0x1610D50 in __module_ra_rrtm_MOD_taugb3
#4  0x1612C1A in __module_ra_rrtm_MOD_gasabs
#5  0x1626DBB in __module_ra_rrtm_MOD_rrtm
#6  0x162BE38 in __module_ra_rrtm_MOD_rrtmlwrad
#7  0x1257D41 in __module_radiation_driver_MOD_radiation_driver
#8  0x13518F1 in __module_first_rk_step_part1_MOD_first_rk_step_part1
#9  0xE0FC03 in solve_em_
#10  0xD2AE49 in solve_interface_
#11  0x459582 in __module_integrate_MOD_integrate
#12  0x404683 in __module

我的namelist.wps设置的区域为：

&geogrid
parent_id         =   1,   1,
parent_grid_ratio =   1,   3,
i_parent_start    =   1,  80,
j_parent_start    =   1,  50,
e_we              =  261, 361,
e_sn              =  201, 301,
geog_data_res     = 'modis_lakes','modis_lakes'
dx = 27000,
dy = 27000,
map_proj = 'lambert',
ref_lat   =  30.0,
ref_lon   = 115.0,
truelat1  =  30.0,
truelat2  =  60.0,
stand_lon = 115.0,
geog_data_path = '/usr/local/geog/geog/'
/
我的namelist.input设置如下：

&domains
time_step                           = 180,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 2,
e_we                                = 261,    361,   94,
e_sn                                = 201,    301,   91,
e_vert                              = 35,     35,    35,
p_top_requested                     = 5000,
num_metgrid_levels                  = 32,
num_metgrid_soil_levels             = 4,
dx                                  = 27000, 9000,  3000,
dy                                  = 27000, 9000,  3000,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,     80,    30,
j_parent_start                      = 1,     50,    30,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     3,     3,
feedback                            = 1,
smooth_option                       = 0
/

&physics
mp_physics                          = 4,     4,     3,
ra_lw_physics                       = 1,     1,     1,
ra_sw_physics                       = 2,     2,     1,
radt                                = 15,    15,    30,
sf_sfclay_physics                   = 1,     1,     1,
sf_surface_physics                  = 2,     2,     2,
bl_pbl_physics                      = 1,     1,     1,
bldt                                = 1,     1,     0,
cu_physics                          = 5,     5,     0,
cudt                                = 5,     5,     5,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
sf_urban_physics                    = 0,     0,     0,
num_land_cat                        = 21
类似于这样的FORTRAN溢出问题要怎么修改？

已解决：错误提示是数组溢出了，因为我修改了namelist中的积分步长，之前dx,dy都设置的是30km，然后时间积分步长time_step为180s，后来dx，dy改为27KM后，时间积分步长没改过来，所以报了上述错误。查了有关资料，有说在实际大气积分过程中，时间积分步长time_step最好为dx，dy公里数的6倍，然后将时间积分步长time_step修改为162之后就能正常运行了。

一颗梅子在南京

好像改了一下时间积分步长就好了～

andrewsoong

一颗梅子在南京 发表于 2017-1-1 11:50
好像改了一下时间积分步长就好了～

数组溢出了，减小时间步长

一颗梅子在南京

andrewsoong 发表于 2017-1-1 12:36
数组溢出了，减小时间步长

恩。。。我改了一下就好了，谢谢～

clannand

请问一下，我和你的错误一样，可是改了time_step也不行。。

薇风拂尘

楼主，同样的问题，改了time_step也还是不行，有解决的办法吗？谢谢

geomuenchen

I revised the time step. It worked just for several seconds, and then the problem came out again.
How to deal with it? Thanks a lot!

一颗梅子在南京

geomuenchen 发表于 2017-8-24 18:02
I revised the time step. It worked just for several seconds, and then the problem came out again.
H ...

这得看你的具体是什么问题了

zhaolinanshiwo

楼主，我把积分步长改为162也还是报错，只不过能模拟11个小时了，没改之前模拟1个小时就报错

sunzy025

这个问题困扰到现在啊~英文读不懂啊