WRF-Chem V4.1 风速模拟低估

小鸣REN

使用WRF-Chem V4.1 对北京地区进行模拟，但是风速始终是低估10%——25%左右。我尝试了不同的物理参数化方案，以及2015、2019、2020年不同月份，都是低估。我感到很困惑，wrf对风速的模拟不是一向高估吗？！怎么我模拟的都是低估。
以下为我的物理参数化方案，气象数据用的是fnl，不同的pbl方案也试过了，都是低估。想请大佬们给点意见和建议，谢谢！

&physics
mp_physics                       = 2,     2,     6,     
progn                               = 1,
ra_lw_physics                    = 4,     4,     1,     
ra_sw_physics                   = 4,     4,     1,     
radt                                  = 30,    30,    4,     
sf_sfclay_physics                = 1,     1,     1,     
sf_surface_physics              = 2,     2,     2,     
bl_pbl_physics                    = 1,     1,     1,     
bldt                                   = 0,     0,     0,     
cu_physics                         = 5,     5,     1,     
cu_diag                              = 1,     1,
cudt                                   = 5,     5,     5,     
isfflx                                  = 1,
ifsnow                                = 1,
icloud                                 = 1,
surface_input_source            = 3,
num_soil_layers                    = 4,
sf_urban_physics                   = 0,
num_land_cat                       = 21,
num_soil_cat                        = 16,
slope_rad                             = 1,     1,     1,
topo_shading                        = 1,     1,     1,
shadlen                                = 30000.,
mp_zero_out                         = 2,
mp_zero_out_thresh               = 1.e-8
maxiens                                = 1,
maxens                                 = 3,
maxens2                               = 3,
maxens3                               = 16,
ensdim                                  = 144,
sst_update                              = 0,
cu_rad_feedback                     = .true.
/

&noah_mp
opt_crs                             = 1,
opt_sfc                             = 1,
opt_btr                             = 1,
opt_run                             = 3,
opt_crop                            = 1,
/


&fdda
grid_fdda                           = 1,     0,     0,
gfdda_inname                        = "wrffdda_d<domain>",
gfdda_end_h                         = 24000, 36, 36,
gfdda_interval_m                    = 360,   360,   360,
fgdt                                = 0,     0,     0,
if_no_pbl_nudging_uv                = 0,     1,     1,
if_no_pbl_nudging_t                 = 0,     1,     1,
if_no_pbl_nudging_q                 = 0,     1,     1,
if_zfac_uv                          = 0,     0,     1,
  k_zfac_uv                          = 10,   10,     1,
if_zfac_t                           = 0,     0,     1,
  k_zfac_t                           = 10,   10,     1,
if_zfac_q                           = 0,     0,     1,
  k_zfac_q                           = 10,   10,     1,
guv                                 = 0.0003,     0.0003,     0.0003,
gt                                  = 0.0003,     0.0003,     0.0003,
gq                                  = 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,
sfs_opt                             = 0,      0,     0,   
diff_6th_opt                        = 2,2,
diff_6th_factor                     = 0.10,
base_temp                           = 290.,
damp_opt                            = 3,
zdamp                               = 5000., 5000., 5000.,
dampcoef                            = 0.05,  0.05,  0.05,  
khdif                               = 0,     0,     0,     
kvdif                               = 0,     0,     0,     
non_hydrostatic                     = .true.,.true.,.true.,
mix_full_fields                     = .true.,.true.,.true.,
h_mom_adv_order                     = 5,
v_mom_adv_order                     = 3,
h_sca_adv_order                     = 5,
v_sca_adv_order                     = 3,
moist_adv_opt                       = 1,     1,     1,   
scalar_adv_opt                      = 1,     1,     1,   
use_baseparam_fr_nml                = .true.
gwd_opt                             = 0,
/

qq531884069

建议看看我们的这个文章
https://www.sciencedirect.com/sc ... i/S1309104221002713

小鸣REN

qq531884069 发表于 2021-12-3 06:38
建议看看我们的这个文章
https://www.sciencedirect.com/science/article/pii/S1309104221002713

谢谢，去拜读一下