WRF使用4层嵌套模拟，运行一段时间后报错，仅某个rsl.error报错如下

钟小继

WRF使用4层嵌套模拟，运行一段时间后报错，仅某个rsl.error报错如下：

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.
Backtrace for this error:
#0  0x7F18D705A697
#1  0x7F18D705ACDE
#2  0x7F18D65562EF
#3  0x2671A8D in __module_sf_sfclayrev_MOD_psim_stable
#4  0x2671B23 in __module_sf_sfclayrev_MOD_zolri2
#5  0x2671DBA in __module_sf_sfclayrev_MOD_zolri
#6  0x2675F59 in __module_sf_sfclayrev_MOD_sfclayrev1d
#7  0x267AA4D in __module_sf_sfclayrev_MOD_sfclayrev
#8  0x1DD0923 in __module_surface_driver_MOD_surface_driver
#9  0x16BBFC3 in __module_first_rk_step_part1_MOD_first_rk_step_part1
#10  0x10DE0D6 in solve_em_
#11  0xFE8079 in solve_interface_
#12  0x46A0DA in __module_integrate_MOD_integrate
#13  0x46A6BB in __module_integrate_MOD_integrate
#14  0x46A6BB in __module_integrate_MOD_integrate
#15  0x46A6BB in __module_integrate_MOD_integrate
#16  0x405E83 in __module_wrf_top_MOD_wrf_run

然后rsl.out文件中末尾都是如下表示：
allocated
d03 2017-08-22_12:00:00 *** Initializing nest domain # 4 from an input file. ***
d03 2017-08-22_12:00:00 med_initialdata_input: calling input_input
INPUT LandUse = "MODIFIED_IGBP_MODIS_NOAH"
LANDUSE TYPE = "MODIFIED_IGBP_MODIS_NOAH" FOUND          33  CATEGORIES           2  SEASONS WATER CATEGORY =           17  SNOW CATEGORY =           15
INITIALIZE THREE Noah LSM RELATED TABLES
d03 2017-08-22_12:00:00  Initializing OML with HML0 =    200.000000
INPUT LandUse = "MODIFIED_IGBP_MODIS_NOAH"
LANDUSE TYPE = "MODIFIED_IGBP_MODIS_NOAH" FOUND          33  CATEGORIES           2  SEASONS WATER CATEGORY =           17  SNOW CATEGORY =           15
INITIALIZE THREE Noah LSM RELATED TABLES
d03 2017-08-22_12:00:00  Initializing OML with HML0 =    200.000000
INPUT LandUse = "MODIFIED_IGBP_MODIS_NOAH"
LANDUSE TYPE = "MODIFIED_IGBP_MODIS_NOAH" FOUND          33  CATEGORIES           2  SEASONS WATER CATEGORY =           17  SNOW CATEGORY =           15
INITIALIZE THREE Noah LSM RELATED TABLES
d03 2017-08-22_12:00:00  Initializing OML with HML0 =    200.000000
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE   1 IS    124 IE    247 JS    212 JE    247
WRF NUMBER OF TILES =   1
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE   1 IS    213 IE    424 JS    364 JE    424
WRF NUMBER OF TILES =   1



最内层的时间步未0.5s，初步猜测不是时间步的问题导致的数值不稳定，在国外网站看到说是可能物理方案不兼容。。但是都运行起来了很长时间才报错的。不知道是不是namelist.input的错误。贴下namelist的设置，希望大家帮忙看看哪里出了问题，我已经卡在这里很多天了，一直在尝试，总是失败。。。。。
&time_control
run_days                            = 0,
run_hours                           = 24,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2017, 2017, 2017,2017,
start_month                         = 08,   08,   08,  08,
start_day                           = 22,   22,   22,  22,
start_hour                          = 12,   12,   12,  12,
start_minute                        = 00,   00,   00,  00,
start_second                        = 00,   00,   00,  00,
end_year                            = 2017, 2017, 2017,2017,
end_month                           = 08,   08,   08,  08,
end_day                             = 23,   23,   23,  23,
end_hour                            = 12,   12,   12,  12,
end_minute                          = 00,   00,   00,  00,
end_second                          = 00,   00,   00,  00,
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,.true.,.true.
history_interval                    = 180,   60,   60,  60,
frames_per_outfile                  = 1,  1000, 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                         = 0
/

&domains
time_step                           = 45,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 4,
e_we                                = 178,   211,   247,  424,
e_sn                                = 178,   211,   247,  424,
e_vert                              = 30,    30,    30,   30,
p_top_requested                     = 5000,
num_metgrid_levels                  = 32,
num_metgrid_soil_levels             = 4,
dx                                  = 12000, 4000,  1333.33,  444.444458,
dy                                  = 12000, 4000,  1333.33,  444.444458,
grid_id                             = 1,     2,     3,     4,
parent_id                           = 1,     1,     2,            3,
i_parent_start                      = 1,     45,   60,   53,
j_parent_start                      = 1,     55,   75,   53,
parent_grid_ratio                   = 1,     3,     3,   3,
parent_time_step_ratio              = 1,     3,     3,   10,
feedback                            = 1,
smooth_option                       = 0
/

&physics
mp_physics                          = 3,     3,     3,     3,
ra_lw_physics                       = 5,     5,     5,     5,
ra_sw_physics                       = 2,     2,     2,     2,
radt                                = 30,    30,    30,    30,
sf_sfclay_physics                   = 91,    91,    91,    91,
sf_surface_physics                  = 2,     2,     2,     2,
bl_pbl_physics                      = 1,     1,     1,     1,
bldt                                = 0,     0,     0,     0,
cu_physics                          = 1,     1,     0,     0,
cudt                                = 5,     5,     5,     5,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
sf_urban_physics                    = 0,     0,     0,  0,
isftcflx                            = 1,
sf_ocean_physics                    = 1,
oml_hml0                            = 200,
oml_gamma                           = 0.14,
omdt                                = 1,
/

&fdda
grid_fdda                           = 1,     0,     0,
gfdda_inname                        = "wrffdda_d<domain>",
gfdda_end_h                         = 24,    23,    23,
gfdda_interval_m                    = 360,   360,   360,
fgdt                                = 0,     0,     0,
fgdtzero                            = 0,     0,     0,
if_no_pbl_nudging_uv                = 1,     0,     0,
if_no_pbl_nudging_t                 = 1,     0,     0,
if_no_pbl_nudging_ph                = 1,     0,     0,
if_zfac_uv                          = 1,     0,     0,
  k_zfac_uv                          = 10,   10,    10,
if_zfac_t                           = 1,     0,     0,
  k_zfac_t                           = 10,   10,    10,
if_zfac_ph                          = 1,     0,     0,
  k_zfac_ph                          = 10,   10,    10,
dk_zfac_uv                          = 1,     1,     1,
dk_zfac_t                           = 1,     1,     1,
dk_zfac_ph                          = 1,     1,     1,
guv                                 = 0.0003,     0.0003,     0.0003,
gt                                  = 0.00,     0.0003,     0.0003,
gph                                 = 0.00,     0.0003,     0.0003,
xwavenum                            = 2,
ywavenum                            = 2,
if_ramping                          = 1,
dtramp_min                          = 60.0,
io_form_gfdda                       = 2,
/

&dynamics
w_damping                           = 0,
diff_opt                            = 1,      1,      1,   1,
km_opt                              = 4,      4,      4,   4,
diff_6th_opt                        = 0,      0,      0,   0,
diff_6th_factor                     = 0.12,   0.12,   0.12,   0.12,
base_temp                           = 290.
damp_opt                            = 0,
zdamp                               = 5000.,  5000.,  5000.,  5000.,
dampcoef                            = 0.2,    0.2,    0.2   0.2,
khdif                               = 0,      0,      0,   0,
kvdif                               = 0,      0,      0,   0,
non_hydrostatic                     = .true., .true., .true.,  .true.,
moist_adv_opt                       = 1,      1,      1,     1,   
scalar_adv_opt                      = 1,      1,      1,     1,   
/

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

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

yuzhang110

算是板凳吗？还是地板了？

陈闯

可能是内存不足导致的。不要满核心运算试试。

沿江西

你的dt和分辨率设置不大合理。建议你调大一点dx和dt

志鹏

解决了吗？你的DY,DX设置感觉有点问题，你嵌套前是否用了WRFdomainwizad设置？另外4重嵌套，请查看并行计算资源够不够，不够或者积分溢出都会导致未完成跳出RUNNING

crazytiy

应该是溢出了

shubabu

感觉应该是溢出吐核了把

向南V5

估计是溢出了

cjy1999

计算量太大了？？

michellexi

陈闯 发表于 2018-9-4 10:24
可能是内存不足导致的。不要满核心运算试试。

您好，请问如何调整运算使用的核数呢