有关wrf自带tc bogus结果的讨论

栗子

用WRFV3.4添加一个bogus到初始场中，对wrfinput文件进行检验，发现风场正常且达到我的要的强度，但是slp却很诡异，与风场极不匹配，且有时会冒出一个高值中心与低值中心并存。请问你们的bogus转出来存在这个情况吗？我应该要怎么修正？第一张是原始分析场；第二张是添加bogus之后的，可以看到气压减少的非常少，而且中心偏了。越往高层，到500hpa就会出现一个高位势中心。如果在namelist里把lattc和lontc改一下，添加vortex的中心气压值就会有很大变化

这是添加bogus之后的初始场slp和U

随缘

那你的namelist什么怎么设置的啊？

栗子

随缘 发表于 2013-9-4 12:29

                               
登录/注册后可看大图

那你的namelist什么怎么设置的啊？

这是我的namelist.input
&time_control
run_days                            = 0,
run_hours                           = 120,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2010, 2000, 2000,
start_month                         = 10,   01,   01,
start_day                           = 17,   24,   24,
start_hour                          = 00,   12,   12,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2010, 2000, 2000,
end_month                           = 10,   01,   01,
end_day                             = 22,   25,   25,
end_hour                            = 00,   12,   12,
end_minute                          = 00,   00,   00,
end_second                          = 00,   00,   00,
interval_seconds                    = 10800
input_from_file                     = .true.,.true.,.true.,
history_interval                    = 180,  60,   60,
frames_per_outfile                  = 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                           = 90,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
e_we                                = 241,    112,   94,
e_sn                                = 226,    97,    91,
e_vert                              = 28,    28,    28,
p_top_requested                     = 5000,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
dx                                  = 30000, 10000,  3333.33,
dy                                  = 30000, 10000,  3333.33,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,     31,    30,
j_parent_start                      = 1,     17,    30,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     3,     3,
feedback                            = 1,
smooth_option                       = 0
/

&physics
mp_physics                          = 8,     3,     3,
ra_lw_physics                       = 1,     1,     1,
ra_sw_physics                       = 1,     1,     1,
radt                                = 30,    30,    30,
sf_sfclay_physics                   = 1,     1,     1,
sf_surface_physics                  = 2,     2,     2,
bl_pbl_physics                      = 1,     1,     1,
bldt                                = 0,     0,     0,
cu_physics                          = 1,     1,     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,
/

&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.,
/

&tc
insert_bogus_storm                  = .true.,
remove_storm                        = .false.,
num_storm                           = 1,
latc_loc                            = 18.26,
lonc_loc                            = 127.79,
vmax_meters_per_second              = 60,
rmax                                = 90000,
vmax_ratio                          = 0.8,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

栗子

我发现如果我将latc_loc和lontc_loc改一下，得到的bogus的slp就会强了很多，如图，同时，那个高压也增强了

随缘

栗子 发表于 2013-9-4 15:25

                               
登录/注册后可看大图

我发现如果我将latc_loc和lontc_loc改一下，得到的bogus的slp就会强了很多，如图，同时，那个高压也增强了
...

bogus的涡旋强度主要是下面三个参数决定的吧
vmax_meters_per_second              = 60,
rmax                                = 90000,
vmax_ratio                          = 0.8,

栗子

随缘 发表于 2013-9-4 18:14

                               
登录/注册后可看大图

bogus的涡旋强度主要是下面三个参数决定的吧
vmax_meters_per_second              = 60,
rmax         ...

是的，所以才觉得比较怪呢。位置变了，强度竟然也跟着变。你用bogus的时候画出的初始场有这个问题吗

随缘

栗子 发表于 2013-9-4 19:50

                               
登录/注册后可看大图

是的，所以才觉得比较怪呢。位置变了，强度竟然也跟着变。你用bogus的时候画出的初始场有这个问题吗

最大风速和slp不匹配是有的

栗子

随缘 发表于 2013-9-4 20:31

                               
登录/注册后可看大图

最大风速和slp不匹配是有的

这个，好吧。。。

东风急流

请问楼主这个rmax值是怎么确定的呢

吴阿好大人

请问楼主你的虚假高压是怎么去掉的？非常感谢！！！！