求问在WRF模拟风速时下雨天的模拟风速明显比实测风速大怎么解决

opeveu

各位大佬，还是以前的问题。过了两个月还是没有解决“WRF模拟风速时下雨天的模拟风速明显比实测风速大，且大很多”。在这里再想问问各位同学和老师是否有知道怎么解决这类问题的（在namelist.input里面设置了top_wind=1和nudging）。万分感谢！！


！！！还有一种思路，来跟大家请教一下。在日常生活中，我们知道，大部分情况下，下雨前的风速会比较大，当雨一开始下的时候，风就会停止。这种现象在WRF模拟中会表现出来吗？根据wrfout文件提取出来的结果，感觉提取出来的风速都是雨前的风速，从而使模拟的风速明显比实测大。
！！！！而且，我们的天气预报中对风速的预报会出现平均风力和阵风的预报。通过对研究地点的实测风速和当地的天气预报风速进行对比，我发现这两者是比较吻合的，反而WRF模拟的风速比天气预报和自己的实测都很大。比如针对一次暴风雪天气，实测风速（10min平均）是在9m/s左右，天气预报是5级风（8.0~10.7m/s）,而模拟的风速达到16m/s。所以我想在模拟降水天气（降雨和降雪）时是不是有什么特殊的设置，比如在namelist.input里面有什么选项开关。
针对以上问题，请求家园里的大佬进行解答，或者有什么链接和资料。万分感谢
以下是我namelist.input的设置。
&time_control
run_days                            = 13,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2020, 2020, 2020, 2020
start_month                         = 09,   09,   09,   09
start_day                           = 29,   29,   29,  29
start_hour                          = 00,   00,   00,   00
end_year                            = 2020, 2020, 2020, 2020
end_month                           = 10,   10,   10,   10
end_day                             = 12,   12,   12,   12
end_hour                            = 00,   00,   00,   00
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,.true.,.true.,
history_interval                    = 180,  180,   180,   60,
frames_per_outfile                  = 1000, 1000, 1000,  1000,
restart                             = .false.,
restart_interval                    = 7200,
io_form_history                     = 2,
io_form_restart                     = 2,
io_form_input                       = 2,
io_form_boundary                    = 2,
debug_level                           = 0,
/

&domains
time_step                           = 60,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 4,
e_we                                = 101,   109,  109,   109,
e_sn                                = 101,   109,   109,   109,
e_vert                              = 58,    58,    58,    58,
p_top_requested                     = 5000,
eta_levels=1,0.9965,0.993,0.989666667,0.986333333,0.983,0.9765,0.97,0.962,0.954,0.944,0.934,0.9215,0.909,0.8945,0.88,0.8547147,0.8294294,0.8041441,0.7788588,0.7317745,0.6868948,0.6441336,0.6034079,0.5646371,0.5277435,0.492652,0.45929,0.4275872,0.3974758,0.3688904,0.3417678,0.3160468,0.2916686,0.2685761,0.2467147,0.2260311,0.2064745,0.1879954,0.1705464,0.1540817,0.1385571,0.1239301,0.1102429,0.09761681,0.08596946,0.07522503,0.0653135,0.05617033,0.04773594,0.03995539,0.03277799,0.02615699,0.020004926,0.014415,0.009217507,0.004422931,0,
num_metgrid_levels                  = 34,
num_metgrid_soil_levels             = 4,
dx                                  = 18000, 6000, 2000, 666.667,
dy                                  = 18000, 6000, 2000, 666.667,
grid_id                             = 1,     2,     3,    4
parent_id                           = 0,     1,     2,    3
i_parent_start                      = 1,     32,    36,    36,
j_parent_start                      = 1,     32,    36,    36,
parent_grid_ratio                   = 1,     3,     3,     3,
parent_time_step_ratio              = 1,     3,     3,     3,
feedback                            = 1,
smooth_option                       = 0,
smooth_cg_topo                      = .true.,
/

&physics
mp_physics                          = 6,    6,    6,    6,
cu_physics                          = 10,    10,    10,    10,
ra_lw_physics                       = 4,    4,    4,    4,
ra_sw_physics                       = 4,    4,    4,    4,
bl_pbl_physics                      =1,    1,   1,    1,
sf_sfclay_physics                       =1,    1,   1,    1,
sf_surface_physics                    = 4,    4,    4,    4,
topo_wind                             =1,    1,   1,    1,
ua_phys                               =.true.,
radt                                = 27,    9,    3,    1,
bldt                                = 0,     0,     0,     0,
cudt                                = 5,     5,     5,     5,
icloud                              = 1,
num_land_cat                        = 21,
sf_urban_physics                    = 0,     0,     0,     0,
/

&fdda
grid_fdda                           = 1,   1,   1,   1
gfdda_inname                        = "wrffdda_d<domain>"
gfdda_interval_m                    = 360,   360,   360,   360
gfdda_end_h                         = 6,   6,   6,   6
io_form_gfdda                       = 2
fgdt                                = 0,   0,   0,   0
if_no_pbl_nudging_uv                = 0,   0,   0,   0
if_no_pbl_nudging_t                 = 1,   1,   1,   1
if_no_pbl_nudging_q                 = 1,   1,   1,   1
guv                        = 0.0001,   0.0001,   0.0001,   0.0001
gt                         = 0.0001,   0.0001,   0.0001,   0.0001
gq                         = 0.0001,   0.0001,   0.0001,   0.0001
if_ramping                = 0
dtramp_min              =60.0
/

&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                            = 3,
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.
nested                              = .false., .true., .true., .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

&scm
scm_wind_adv                   =.true.
scm_force_wind_largescale                   =.true.
/

astiny

你nudging只nudge了t和q啊。。。
风速大没办法的。改改参数化方案吧。

opeveu

astiny 发表于 2021-3-9 09:45
你nudging只nudge了t和q啊。。。
风速大没办法的。改改参数化方案吧。

谢谢回复。参数化方案试过很多了，现在用的方案是比较之后选出来的。其实，在不同的参数化方案情况下，模拟的结果并没有很大的变化，结果差值有但是比较小。另外针对您说的nudging，我在使用手册上看到是当if_no_pbl_nudging_uv=0时是nudging in the PBL，=1时是no nudging of u and v in the PBL。所以在nudging设置中应该是没错的。您可以看一看。我们可以继续交流。

astiny

opeveu 发表于 2021-3-9 11:01
谢谢回复。参数化方案试过很多了，现在用的方案是比较之后选出来的。其实，在不同的参数化方案情况下，模 ...

我的意思是，可以强行nudging u和v。这样效果会变好，但是好像不太符合物理规律。。。

opeveu

astiny 发表于 2021-3-9 09:45
你nudging只nudge了t和q啊。。。
风速大没办法的。改改参数化方案吧。

或者您知道在预报中是否针对这种雷雨天气有什么特殊设置，看别的链接与论坛，好像预报与模拟的namelist.input并没有什么不同。

opeveu

astiny 发表于 2021-3-9 11:02
我的意思是，可以强行nudging u和v。这样效果会变好，但是好像不太符合物理规律。。。

请问，强行nudging u和v是怎么设置的吗？有例子或者文本吗？设置Grid Nudging或者设置Spectral Nudging？谢谢您的及时回复。

astiny

opeveu 发表于 2021-3-9 11:06
请问，强行nudging u和v是怎么设置的吗？有例子或者文本吗？设置Grid Nudging或者设置Spectral Nudging？ ...

就是namelist里nudge_uv参数设1. 不要管别人说行不行。。。

astiny

opeveu 发表于 2021-3-9 11:03
或者您知道在预报中是否针对这种雷雨天气有什么特殊设置，看别的链接与论坛，好像预报与模拟的namelist.i ...

不太清楚降水的风速模拟，你可以尝试看看别人文章里的参数化设置。你可以先试试增加垂直层数看看，然后其他参数可以试试改改边界层、路面这些参数化方案。
根本上来说，某种参数化方案模拟结果偏大，解决方案应该是在现有的方案基础上做敏感实验，测不同种类参数化方案的排列组合。一般风速可能还是和边界层、陆面影响大一些。辐射也有可能有影响。肯定不会是微物理参数化的问题，也不太可能是积云参数化的问题。无脑加层数应该也是一种改进的路径，以我的经验来说至少不会有负面效果。除此之外还可以尝试改wps的地形、换再分析资料/预报资料等方法。
还是等多试。

opeveu

astiny 发表于 2021-3-9 11:47
不太清楚降水的风速模拟，你可以尝试看看别人文章里的参数化设置。你可以先试试增加垂直层数看看，然后其 ...

好的。谢谢

opeveu

astiny 发表于 2021-3-9 11:40
就是namelist里nudge_uv参数设1. 不要管别人说行不行。。。

nudge_uv参数设1是哪个选项啊能具体说一下不？