WRF 台风模拟结果不理想

蚊子

很长一段时间都在进行台风模拟，参考了一些文献设置，但效果一直不好。
但自己确实搞不懂问题出在哪里，因此希望各位前辈指点下。

用WRF4.0并行模拟，fnl 1°× 1°
     嵌套区域（三重-双向固定）

嵌套区域-三重双向固定

namelist.wps  如下：
&share
wrf_core = 'ARW',
max_dom = 3,
start_date = '2017-06-09_18:00:00', '2017-06-09_18:00:00', '2017-06-09_18:00:00',
end_date   = '2017-06-13_06:00:00', '2017-06-13_06:00:00', '2017-06-13_06:00:00',
interval_seconds = 21600,
io_form_geogrid = 2,
opt_output_from_geogrid_path = '/home/mode/WPS/merbok/',
debug_level = 0,
/

&geogrid
parent_id         = 1,1,2,
parent_grid_ratio = 1,3,3,
i_parent_start    = 1,24,52,
j_parent_start    = 1,24,37,
e_we          = 110,190,262,
e_sn          = 110,190,352,
geog_data_res = '10m','2m','30s',
dx = 36000,
dy = 36000,
map_proj =  'mercator',
ref_lat   = 19.15,
ref_lon   = 117.15,
truelat1  = 19.15,
truelat2  = 0,
stand_lon = 117.15,
geog_data_path = '/home/geog/WPS_GEOG',
opt_geogrid_tbl_path = '/home/mode/WPS/merbok/',
ref_x = 55.0,
ref_y = 55.0,
/

&ungrib
out_format = 'WPS',
prefix = 'FILE',
/

&metgrid
fg_name = 'FILE',
io_form_metgrid = 2,
opt_output_from_metgrid_path = '/home/mode/WPS/merbok/',
opt_metgrid_tbl_path = '/home/mode/WPS/merbok/',
/

&mod_levs
press_pa = 201300 , 200100 , 100000 ,
             95000 ,  90000 ,
             85000 ,  80000 ,
             75000 ,  70000 ,
             65000 ,  60000 ,
             55000 ,  50000 ,
             45000 ,  40000 ,
             35000 ,  30000 ,
             25000 ,  20000 ,
             15000 ,  10000 ,
              5000 ,   1000
/


&domain_wizard
grib_data_path = '/home/fnls/merbok',
grib_vtable = 'Vtable.GFS',
dwiz_name    =merbok
dwiz_desc    =merbok
dwiz_user_rect_x1 =3162
dwiz_user_rect_y1 =608
dwiz_user_rect_x2 =3524
dwiz_user_rect_y2 =986
dwiz_show_political =true
dwiz_center_over_gmt =true
dwiz_latlon_space_in_deg =10
dwiz_latlon_linecolor =-8355712
dwiz_map_scale_pct =25.0
dwiz_map_vert_scrollbar_pos =0
dwiz_map_horiz_scrollbar_pos =0
dwiz_gridpt_dist_km =36.0
dwiz_mpi_command =
dwiz_tcvitals =null
dwiz_bigmap =Y
/


namelist.input  如下
&time_control            
run_days                 = 3,
run_hours                = 12,
run_minutes              = 0,
run_seconds              = 0,
start_year               = 2017,     2017,     2017,
start_month              = 06,       06,       06,
start_day                = 09,       09,       09,
start_hour               = 18,       18,       18,
start_minute             = 00,       00,       00,
start_second             = 00,       00,       00,
end_year                 = 2017,     2017,     2017,
end_month                = 06,       06,       06,
end_day                  = 13,       13,       13,
end_hour                 = 06,       06,       06,
end_minute               = 00,       00,       00,
end_second               = 00,       00,       00,
interval_seconds         = 21600,
input_from_file          = .true.,   .true.,   .true.,
history_interval         = 180,       60,       60,
frames_per_outfile       = 1000,     1000,     1000,
restart                  = .false.,
restart_interval         = 2880,
io_form_history          = 2,
io_form_restart          = 2,
io_form_input            = 2,
io_form_boundary         = 2,
debug_level              = 0,
/

&domains                 
eta_levels               = 1.000, 0.994, 0.9881, 0.9822, 0.9763,
                0.9705, 0.9641, 0.957, 0.9491, 0.9401,
                0.9298, 0.9176, 0.9033, 0.8863, 0.8662,
                0.8425, 0.8148, 0.7827, 0.7461, 0.7053,
                0.6606, 0.6126, 0.5623, 0.5106, 0.4583,
                0.4064, 0.3557, 0.3065, 0.2595, 0.2149,
                0.1729, 0.1336, 0.097, 0.063, 0.0317,
                0.000,
time_step                = 180,
time_step_fract_num      = 0,
time_step_fract_den      = 1,
max_dom                  = 3,
e_we                     = 110,      190,      262,
e_sn                     = 110,      190,      352,
e_vert                   = 36,       36,       36,
p_top_requested          = 5000,
num_metgrid_levels       = 32,
num_metgrid_soil_levels  = 4,
dx                       = 36000,    12000,     4000,
dy                       = 36000,    12000,     4000,
grid_id                  = 1,        2,        3,
parent_id                = 1,        1,        2,
i_parent_start           = 1,       24,       52,
j_parent_start           = 1,       24,       37,
parent_grid_ratio        = 1,        3,        3,
parent_time_step_ratio   = 1,        3,        3,
feedback                 = 1,
smooth_option            = 0,
/

&physics                 
mp_physics               = 6,        6,        6,
ra_lw_physics            = 1,        1,        1,
ra_sw_physics            = 1,        1,        1,
radt                     = 30,       12,       12,
sf_sfclay_physics        = 1,        1,        1,
sf_surface_physics       = 2,        2,        2,
bl_pbl_physics           = 1,        1,        1,
bldt                     = 0,        0,        0,
cu_physics               = 6,        6,        0,
cudt                     = 5,        5,        5,
isfflx                   = 1,
ifsnow                   = 0,
icloud                   = 1,
surface_input_source     = 1,
num_soil_layers          = 4,
sf_urban_physics         = 0,        0,        0,
maxiens                  = 1,
maxens                   = 3,
maxens2                  = 3,
maxens3                  = 16,
ensdim                   = 144,
/

&fdda                    
/

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

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

&grib2                  
/

&namelist_quilt         
nio_tasks_per_group      = 0,
nio_groups               = 1,
/

模拟结果，台风路径对比

台风路径对比

模拟台风路径 1.与bst相比，距离相差太大
         2.与bst相比，台风移动延迟

islandowner

如果没啥好办法，换参数化方案试试呗

男紫汗

据说可以先用两重模拟，路径模拟好了再用三重。两重的话跑得快，能节约时间

kazami2003

TC的初始位置差蠻多的

freekiller

起始就不好，换初始场，初始化台风，同化，先把初始场搞好再谈模式

chensiqi

您好，我最近也在转台风嵌套模拟，但是我已经按照官网说的已经用dm编译并且选了3，网格也设置了移动嵌套，为什么我的输出结果rsl.out.0000里面没有ATCF路径输出阿

zzwei

你好，请问问题解决了吗？

YonghangX

chensiqi 发表于 2019-6-11 14:53
您好，我最近也在转台风嵌套模拟，但是我已经按照官网说的已经用dm编译并且选了3，网格也设置了移动嵌套， ...

同样的问题，不知道你解决了没有？