./wrf.exe 默认的CONUS参数化组合能运行，自己设置的参数化方案报段错误

muyuqingye

我按照试验需求，设计了一个控制性试验和两个敏感性试验，控制性试验和敏感性试验1都转通了，但替换关键区内海平面气压的敏感性试验2没有跑通，奇怪的是，我把敏感性试验2的参数化方案调成CONUS 方案组合是可以跑通的，但是按照控制性试验的参数化方案就会在12:30:04出现段错误，想请教一下这是什么原因导致的？有什么解决方法吗？下面是我的namelist.input文件 &time_control
run_days                            = 0,
run_hours                           = 23,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 1979, 1979, 1979, 1979, 1979, 1979, 1979,
start_month                         = 06,   06,   06,   06,   06,   06,   06,
start_day                           = 01,   01,   01,   01,   01,   01,   01,
start_hour                          = 00,   00,   00,   00,   00,   01,   01,
end_year                            = 1979, 1979, 1979, 1979, 1979, 1979, 1979,
end_month                           = 06,   06,   06,   06,   06,   06,   06,
end_day                             = 01,   01,   01,   01,   01,   01,   01,
end_hour                            = 23,   23,   23,   23,   23,   23,   23,
interval_seconds                    = 3600,
input_from_file                     = .true.,.true.,.true.,.true.,true.,true.,true.,
history_interval                    = 60,   60,   60,   60,   60,   60,   60,
frames_per_outfile                  = 1,    1,    1,    1,    1,    1,    1,
restart                             = .false.,
restart_interval                    = 7200,
io_form_history                     = 2,
io_form_restart                     = 2,
io_form_input                       = 2,
io_form_boundary                    = 2,
/

&domains
time_step                           = 125,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 7,
e_we                                = 73, 171, 56, 121, 51,  71,  71,
e_sn                                = 57, 121, 56,  56, 61,  81,  71,
e_vert                              = 40,    40,   40,  40,  40,  40,  40,
p_top_requested                     = 5000,
num_metgrid_levels                  = 38,
num_metgrid_soil_levels             = 4,
dx                                  = 125000,25000,5000,5000,5000,5000,5000,
dy                                  = 125000,25000,5000,5000,5000,5000,5000,
grid_id                             = 1,     2,    3,    4,    5,    6,    7,
parent_id                           = 1,     1,    2,    2,    2,    2,    2,
i_parent_start                      = 1,   20,  54,  62,  95, 101,   95,
j_parent_start                      = 1,   17,  12,  24,  34,  56,   83,
parent_grid_ratio                   = 1,    5,   5,    5,    5,    5,    5,
parent_time_step_ratio              = 1,    5,   5,    5,    5,    5,    5,
feedback                            = 1,
smooth_option                       = 0,
smooth_cg_topo                      =.true.,
/

&physics
mp_physics                          =  3,     3,    3,    3,    3,    3,    3,
cu_physics                          =  1,     1,    0,    0,    0,    0,    0,
ra_lw_physics                       =  1,     1,    1,    1,    1,    1,    1,
ra_sw_physics                       =  1,     1,    1,    1,    1,    1,    1,
bl_pbl_physics                      =  1,     1,    1,    1,    1,    1,    1,
sf_sfclay_physics                   =  1,     1,    1,    1,    1,    1,    1,
sf_surface_physics                  =  2,     2,    2,    2,    2,    2,    2,
radt                                = 125, 125,   125, 125,  125,  125,  125,

bldt                                = 0,     0,     0,   0,    0,    0,    0,
cudt                                = 125, 125,   125, 125,  125,  125,  125,
icloud                              = 1,
num_land_cat                        = 21,
sf_urban_physics                   = 0,     0,     0,    0,    0,    0,    0,
/

&fdda
/

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

&bdy_control
spec_bdy_width                      = 5,
specified                           = .true.
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

2417404759

我的想法是可不可以把smooth打开？因为之前我遇到过好像是我模拟的时候会出现某个点的数值出现过大的情况导致段错误，或者检查一下rsl文件是不是有cfl错误？

muyuqingye

我刚刚检查了，确实有cfl错误，我添加smooth_cg_topo = .true.和w_damping = 1，在试着调epssm参数，但我不知道是应该增大epssm还是减小epssm

ayzqs

减少时间步长试一试

muyuqingye

ayzqs 发表于 2024-6-2 18:39
减少时间步长试一试

我的积分步长已经是一倍网格距了，是不是不能再小了呀？

muyuqingye

2417404759 发表于 2024-6-2 10:27
我的想法是可不可以把smooth打开？因为之前我遇到过好像是我模拟的时候会出现某个点的数值出现过大的情况导 ...

您好！我设定了smooth_cg_topo= .true. 和smooth_option=1，我搜索WRF手册，貌似没有其他平滑的选项了。我还根据WRF论坛中的建议把epssm参数调到了最大的0.9，还是在相同时刻出现段错误，也仍然有cfl错误发生，由于我的模拟区域包括了青藏高原东部，参考论坛上的其他帖子，我觉得是地形引起的积分中断，您还有其他的建议吗？

2417404759

我觉得那可能会是需要斟酌一下参数化方案了，这个我对高原部分的不是很了解，或者你搜一下在高原地区的WRF模拟参数化方案的设置？

muyuqingye

2417404759 发表于 2024-6-5 09:51
我觉得那可能会是需要斟酌一下参数化方案了，这个我对高原部分的不是很了解，或者你搜一下在高原地区的WRF ...

非常感谢您的建议！我调换一下参数化方案，谢谢！