关于wrfda同化后更新边界的问题

林子木

我做完变分同化后更新边界条件 然后向前预报29个小时后断掉，查rsl文件发现在边界上有好几个点垂直速度过大。
下面是我的parame.in文件
&control_param
da_file = './wrfvar_output'
wrf_bdy_file = './wrfbdy_d01'
wrf_input = 'wrfinput_d01'
domain_id = 1
cycling = .true.
debug = .true.
low_bdy_only = .false.
update_lsm = .false.

var4d_lbc = .false.
三层的domain ,同化的第一层。
以下是namelist
  &time_control
run_days                            = 0,
run_hours                           = $1,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = $2,   $2, $2,
start_month                         = $3,   $3, $3,
start_day                           = $4,   $4, $4,
start_hour                          = $5,   $5, $5,
start_minute                        = 00,   00, 00,
start_second                        = 00,   00, 00,
end_year                            = $6,   $6, $6,
end_month                           = $7,   $7, $7,
end_day                             = $8,   $8, $8,
end_hour                            = $9,   $9, $9,
end_minute                          = 00,   00, 00,
end_second                          = 00,   00, 00,
interval_seconds                    = 21600
input_from_file                     = .true.,.false.,.false.,
history_interval                    = 180,  180, 180,
frames_per_outfile                  = 1, 1, 1,
restart                             = .false.,
restart_interval                    = 3600,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
debug_level                         = 0
/
&domains
time_step                           = 120,  
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 3,
e_we                                = 260, 376, 190,
e_sn                                = 200, 310, 190,
e_vert                              = 45,45,45,
dx                                  = 27000,9000,3000,
dy                                  = 27000,9000,3000,
grid_id                             = 1,2,3
parent_id                           = 0,1,2
i_parent_start                      = 1,85,211,     
j_parent_start                      = 1,40,108,     
parent_grid_ratio                   = 1,3,3  
parent_time_step_ratio              = 1,3,3
feedback                            = 1,
smooth_option                       = 1,
vortex_interval                     = 15, 15,15
max_vortex_speed                    = 40, 40,40
corral_dist                         = 28, 18, 18
track_level                         = 85000,
time_to_move                        = 720,180,0
p_top_requested                     = 5000,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
eta_levels                          = 1.000, 0.998, 0.996, 0.995, 0.993,
                                       0.990, 0.988, 0.985, 0.983, 0.980,
                                       0.975, 0.970, 0.960, 0.954, 0.945,
                                       0.934, 0.909, 0.880, 0.845, 0.807,
                                       0.765, 0.719, 0.672, 0.622, 0.571,
                                       0.520, 0.468, 0.420, 0.376, 0.335,
                                       0.298, 0.263, 0.231, 0.202, 0.175,
                                       0.150, 0.127, 0.106, 0.088, 0.070,
                                       0.055, 0.040, 0.026, 0.013, 0.000
/

&physics
mp_physics                          = 6,6,6,
ra_lw_physics                       = 4,4,4
ra_sw_physics                       = 4,4,4
radt                                = 10,
sf_sfclay_physics                   = 1,1,1,
sf_surface_physics                  = 2,2,2,
bl_pbl_physics                      = 1,1,1,
bldt                                = 0,0,0,
cu_physics                          = 6,6,6,   
cudt                                = 0,0,0,
isfflx                              = 2,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
sf_urban_physics                    = 0,     0,     0,
/

&fdda
/

&dynamics
w_damping                           = 1,
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.,
periodic_x                          = .false.,.false.,.false.,
symmetric_xs                        = .false.,.false.,.false.,
symmetric_xe                        = .false.,.false.,.false.,
open_xs                             = .false.,.false.,.false.,
open_xe                             = .false.,.false.,.false.,
periodic_y                          = .false.,.false.,.false.,
symmetric_ys                        = .false.,.false.,.false.,
symmetric_ye                        = .false.,.false.,.false.,
open_ys                             = .false.,.false.,.false.,
open_ye                             = .false.,.false.,.false.,
nested                              = .false., .true., .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/



[img]file:///C:\Documents and Settings\Administrator\Application Data\Tencent\Users\767429827\QQ\WinTemp\RichOle\$_}56NQX)7LA2WN[6EJO[VY.png[/img]

林子木

解决了 ，根本不是这个地方的问题

白兔糖0217

林子木 发表于 2015-9-2 14:13
解决了 ，根本不是这个地方的问题

楼主，是什么问题？？？

林子木

白兔糖0217 发表于 2015-10-14 15:27
楼主，是什么问题？？？

是因为我用的移动网格，同化的场和预报的场不是同一个

杨雨轩

请问楼主有没有遇到过边界无法更新的情况，我提交了作业，跑了一晚上都没有进展

林子木

杨雨轩 发表于 2016-10-26 11:03
请问楼主有没有遇到过边界无法更新的情况，我提交了作业，跑了一晚上都没有进展

没有哎，是不是时间不一致，比如你同化的06时的资料，但wrfinput是00时刻的

杨雨轩

林子木 发表于 2016-10-26 18:47
没有哎，是不是时间不一致，比如你同化的06时的资料，但wrfinput是00时刻的

时间是一致的，我也考虑到了时间不匹配的问题，所以后来我尝试不考虑spin-up，所有时间段都是一样的，但同样出现了问题

林子木

杨雨轩 发表于 2016-10-26 19:10
时间是一致的，我也考虑到了时间不匹配的问题，所以后来我尝试不考虑spin-up，所有时间段都是一样的，但 ...

有没有可能是你的exe没编译好，你链接一个以前用过的绝对没问题的exe试试

今朝几朵白云生

林子木 发表于 2016-10-26 18:47
没有哎，是不是时间不一致，比如你同化的06时的资料，但wrfinput是00时刻的

楼主，遇到这种问题该怎么办？我同化了06是的资料，更新侧边界条件后，wrfbdy的起始时刻还是00时。更新之后不应该也是06时吗？

林子木

今朝几朵白云生 发表于 2018-9-30 10:19
楼主，遇到这种问题该怎么办？我同化了06是的资料，更新侧边界条件后，wrfbdy的起始时刻还是00时。更新之 ...

工作后完全跟以前的东西没关系，所以都忘了好难过