登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
菜鸟一枚,在论坛里搜索了相关的模拟结果对比,没有发现类似问题,只好在此开贴求助,请路过的各位大神帮忙看看问题所在。。
利用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,
/
|