求助：wrf模拟风速较观测值偏大

lvxiaott

    菜鸟一枚，在论坛里搜索了相关的模拟结果对比，没有发现类似问题，只好在此开贴求助，请路过的各位大神帮忙看看问题所在。。
利用wrf模拟了2013年12月的气象场，用ARWpost抽取其中的风速风向后，发现与实际观测值差距较大，其中风速总体趋势尚可，但数值几乎为观测值的2倍左右。。模式风速处理是通过抽取U10、V10得到，风速计算公式为sqrt(U10*U10+V10*V10)。
    总共做了3层网格嵌套，所有数据均提取自最内层，初步考虑可能是参数化方案选取不当导致下垫面与实际情况不符合，请各位大神帮忙捉捉虫。。namelist如下：

&time_control
run_days                            = 32,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2013, 2013, 2013, 2013,
start_month                         = 12,   12,   12,   12,
start_day                           = 01,   01,   01,   01,
start_hour                          = 00,   00,   00,   00,
start_minute                        = 00,   00,   00,   00,
start_second                        = 00,   00,   00,   00,
end_year                            = 2014, 2014, 2014, 2014
end_month                           = 01,   01,   01,   01,
end_day                             = 01,   01,   01,   01,
end_hour                            = 00,   00,   00,   00,
end_minute                          = 00,   00,   00,   00,
end_second                          = 00,   00,   00,   00,
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,.true.,.true.,
history_interval                    = 60,  60,   60,   60,
frames_per_outfile                  = 24, 24, 24, 24,
restart                             = .false.,
restart_interval                    = 0,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
debug_level                         = 0
auxinput11_interval_s               = 60, 60, 60, 60, 60,
auxinput11_end_h                    = 24, 24, 24, 24, 24,
/
&domains
time_step                           = 60,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 3,
e_we                                = 134,    235,   154,    109,
e_sn                                = 134,    247,   208,    163,
e_vert                              = 28,     28,    28,     28,
p_top_requested                     = 10000,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
dx                                  = 27000, 9000,  3000,  1000,
dy                                  = 27000, 9000,  3000,  1000,
grid_id                             = 1,     2,     3,     4,
parent_id                           = 0,     1,     2,     3,
i_parent_start                      = 1,     28,    103,   10,
j_parent_start                      = 1,     31,    82,    118,
parent_grid_ratio                   = 1,     3,     3,     3,
parent_time_step_ratio              = 1,     3,     3,     3,
feedback                            = 1,
smooth_option                       = 0
/
&physics
mp_physics                          = 2,     2,     2,     2,
ra_lw_physics                       = 1,     1,     1,     1,
ra_sw_physics                       = 1,     1,     1,     1,
radt                                = 27,    9,     3,     1,
sf_sfclay_physics                   = 1,     1,     1,     1,
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                              = 0,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 5,
sf_urban_physics                    = 1,     1,     1,     1,
maxiens                             = 1,
maxens                              = 3,
maxens2                             = 3,
maxens3                             = 16,
ensdim                              = 144,
/
!&fdda
! obs_nudging
!obs_nudge_opt                       = 1,1,1,0,0
!max_obs                             = 150000,
!fdda_start                          =     0.,     0.,     0.,     0.,     0.
!fdda_end                            = 99999., 99999., 99999., 99999., 99999.
!obs_nudge_wind                      = 1,1,1,1,1
!obs_coef_wind                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
!obs_nudge_temp                      = 1,1,1,1,1
!obs_coef_temp                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
!obs_nudge_mois                      = 1,1,1,1,1
!obs_coef_mois                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
!obs_rinxy                           = 240.,240.,180.,180,180
!obs_rinsig                          = 0.1,
!obs_twindo                          = 0.6666667,0.6666667,0.6666667,0.6666667,0.6666667,
!obs_npfi                            = 10,
!obs_ionf                            = 1, 1, 1, 1, 1,
!obs_idynin                          = 0,
!obs_dtramp                          = 40.,
!obs_prt_freq                        = 10, 10, 10, 10, 10,
!obs_prt_max                         = 10
!obs_ipf_errob                       = .true.
!obs_ipf_nudob                       = .true.
!obs_ipf_in4dob                      = .true.
!obs_ipf_init                        = .true.

! grid_nudging
! grid_fdda                           = 1,     1,     1     1,
! gfdda_inname                        = "wrffdda_d<domain>",
! gfdda_end_h                         = 96,    96,   96,    96,
! gfdda_interval_m                    = 360,   360,   360,  360,
fgdt                                = 0,     0,     0,    0,
if_no_pbl_nudging_uv                = 0,     0,     0,    0,
if_no_pbl_nudging_t                 = 0,     0,     0,    0,
if_no_pbl_nudging_q                 = 0,     0,     0,    0,
if_zfac_uv                          = 0,     0,     0,    0,
  k_zfac_uv                          = 10,   10,    10,    10,
if_zfac_t                           = 0,     0,     0,    0,
  k_zfac_t                           = 10,   10,    10,    10,
if_zfac_q                           = 0,     0,     0,    0,
  k_zfac_q                           = 10,   10,    10,    10,
guv                                 = 0.0003,     0.0003,     0.0003,     0.0003,
gt                                  = 0.0003,     0.0003,     0.0003,     0.0003,
gq                                  = 0.0003,     0.0003,     0.0003,     0.0003,
if_ramping                          = 1,   
dtramp_min                          = 60.0,
io_form_gfdda                       = 2,
/
&dynamics
w_damping                           = 1,
diff_opt                            = 1,
km_opt                              = 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,
/
&dfi_control
  /
&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,
/

沿江西xi

风和雨一样，要做好，一直都比较难。不过风速比风向，相对容易一点。因为你没有给出scale，所以不清楚你到底模拟大了多少。但总体来看，如果你这个是城市站点，建议你采用BEP的城市冠层方案；如非，同时你用的是USGS的LULC，建议你用MODIS的LULC，并考虑适当加大某些类型的Z。另外，如果你一定要用YSU的PBL，建议你开启TOPO WIND方案
简单的经验建议，希望对你有用

zhangxiao

auxinput11_interval_s这个是什么意思呢 一直没有看懂说明文档的解释

周亚端

想问问楼主对应站点的模拟值是如何提取的

imgrus

这个搞清楚可以自己发文章你
楼主你还是自己查文献吧
论坛找答案比较困难

大气杨娃娃

沿江西xi 发表于 2014-11-14 08:23
风和雨一样，要做好，一直都比较难。不过风速比风向，相对容易一点。因为你没有给出scale，所以不清楚你到 ...

根据本人模拟经验，在有明显强天气过程的时次，wrf模拟的近地层风速与实际观测值有明显的偏差，其风速值明显较观测值大很多。如何优化WRF模拟边界层参数过程是个值得研究的问题。我一般用YSU。不知道你在这方面有没有过什么研究？

沿江西xi

大气杨娃娃 发表于 2014-11-14 10:55
根据本人模拟经验，在有明显强天气过程的时次，wrf模拟的近地层风速与实际观测值有明显的偏差，其风速值 ...

如果是单纯的wrf模拟，我也一般偏向YSU；如果要加上chem，我会用MYJ

中尺度风速模拟

你好 请问怎么用ARWpost来提取风速风向呢？我之前用的是ncl

中尺度风速模拟

也许你可以试着用read_wrf_nc 修改下垫面

chongzika

中尺度风速模拟 发表于 2015-4-27 10:27
你好 请问怎么用ARWpost来提取风速风向呢？我之前用的是ncl

ncl能直接提取风速么