四维变分同化运行报错

ljf

各位家园的前辈好！！！我是一名WRF新手。
我在对嵌套域d02运行4D-var时老是报错：da_med_latbound_in: error opening wrfbdy_d01 for reading. 查了一圈资料（论坛、WRF user guide、Google）都没找到病灶。难道四维变分同化的嵌套域d02运行时，需要把d01的wrfbdy_d01和wrfinput_d01都链接到工作目录下？按我的理解只需要wrfinput_d02就可以啊，这报错到底是为什么？麻烦前辈们指点。

报错截图如下所示：

                               
登录/注册后可看大图

附件分别为real和wrfda使用的namelist.input文件。
namelist.input_real (4.21 KB, 下载次数: 0)

2023-11-8 19:51 上传

点击文件名下载附件
real使用的namelist
namelist.input_da (4.19 KB, 下载次数: 0)

2023-11-8 19:51 上传

点击文件名下载附件
同化使用的namelist

ljf

重新贴一下real的namelist.input：
&time_control
start_year                          = 2023, 2023,
start_month                         = 01,   01,
start_day                           = 10,   10,
start_hour                          = 00,   00,
start_minute                        = 00,   00,
start_second                        = 00,   00,
end_year                            = 2023, 2023,
end_month                           = 01,  01,
end_day                             = 10,  10,
end_hour                            = 06,  06,
end_minute                          = 00,  00,
end_second                          = 00,  00,
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,
fine_input_stream                   = 0, 0
history_interval                    = 180,  180,
frames_per_outfile                  = 1, 1,
restart                             = .false.,
restart_interval                    = 50000,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
debug_level                         = 0
inputout_interval                   = 360, 360,
inputout_begin_h                    = 6, 6,
inputout_end_h                      = 6, 6,
/
&dfi_control
/
&domains
time_step                           = 120,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 2,
s_we                                = 1,    1,
e_we                                = 61,  61,
s_sn                                = 1,    1,
e_sn                                = 61,  61,
s_vert                              = 1,     1,
e_vert                              = 45,   45,
num_metgrid_levels                  = 34,
dx                                  = 27000, 9000,
dy                                  = 27000, 9000,
grid_id                             = 1,     2,   
parent_id                           = 1,     1,
i_parent_start                      = 1,    11,
j_parent_start                      = 1,    11,
parent_grid_ratio                   = 1,     3,
parent_time_step_ratio              = 1,     3,
feedback                            = 1,
smooth_option                       = 0
p_top_requested                     = 5000
/

&physics
mp_physics                          = 3,     3,
ra_lw_physics                       = 1,     1,
ra_sw_physics                       = 1,     1,
radt                                = 30,    30,
sf_sfclay_physics                   = 1,     1,
sf_surface_physics                  = 2,     2,
bl_pbl_physics                      = 1,     1,
bldt                                = 0,     0,
cu_physics                          = 1,     1,
cudt                                = 5,     5,
isfflx                              = 1,
ifsnow                              = 0,
icloud                              = 1,
surface_input_source                = 3,
num_soil_layers                     = 4,
mp_zero_out                         = 0,
maxiens                             = 1,
maxens                              = 3,
maxens2                             = 3,
maxens3                             = 16,
ensdim                              = 144,
num_land_cat                        = 21,
/
&fdda
/
&dynamics
hybrid_opt = 0,
use_theta_m = 0,
w_damping                           = 0,
diff_opt                            = 1,
km_opt                              = 1,
diff_6th_opt                        = 0,     0,
diff_6th_factor                     = 0.12,  0.12,
damp_opt                            = 0,
base_temp                           = 290.,
iso_temp                            = 200.,
zdamp                               = 5000.,  5000.,
dampcoef                            = 0.2,    0.2,
khdif                               = 0,      0,
kvdif                               = 0,      0,
/
&bdy_control
spec_bdy_width                      = 5,
spec_zone                           = 1,
relax_zone                          = 4,
specified                           = .true., .false.,
nested                              = .false., .true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

ljf

下面是wrfda的namelist.input:
&wrfvar1
var4d=true,
var4d_lbc=false,
var4d_bin=3600,
/
&wrfvar2
/
&wrfvar3
ob_format=2
/
&wrfvar4
USE_SYNOPOBS =  T,
USE_SHIPSOBS =  T,
USE_METAROBS =  T,
USE_SOUNDOBS =  T,
USE_MTGIRSOBS =  T,
USE_TAMDAROBS =  T,
USE_PILOTOBS =  T,
USE_AIREPOBS =  T,
USE_GEOAMVOBS =  T,
USE_POLARAMVOBS =  T,
USE_BOGUSOBS =  T,
USE_BUOYOBS =  T,
USE_PROFILEROBS =  T,
USE_SATEMOBS =  T,
USE_GPSZTDOBS =  F,
USE_GPSPWOBS =  F,
USE_GPSREFOBS =  F,
USE_QSCATOBS =  F,
USE_AIRSRETOBS =  F,
use_ssmiretrievalobs=false
use_amsuaobs = F,
use_amsubobs = F,
use_mhsobs   = F,
use_airsobs  = F,
use_eos_amsuaobs = F,
USE_OBS_ERRFAC=F
/
&wrfvar5
check_max_iv=true,
/
&wrfvar6
max_ext_its=1,
ntmax=100,
eps=0.01,
orthonorm_gradient=.true.,
/
&wrfvar7
cv_options=3,
je_factor=1.0,
/
&wrfvar8
/
&wrfvar9
/
&wrfvar10
test_transforms=false,
test_gradient=false,
/
&wrfvar11
check_rh=0,
/
&wrfvar12
/
&wrfvar13
/
&wrfvar14
/
&wrfvar15
/
&wrfvar16
alphacv_method=2,
ensdim_alpha=0,
alpha_truncation=0,
alpha_corr_type=3,
alpha_corr_scale=1500,
alpha_std_dev=1.0,
/
&wrfvar17
analysis_type="3D-VAR",
/
&wrfvar18
analysis_date= 2023-01-10_00:00:00,
/
&wrfvar19
/
&wrfvar20
/
&wrfvar21
time_window_min="2023-01-10_00:00:00",
/
&wrfvar22
time_window_max="2023-01-10_06:00:00",
/
&wrfvar23
/
&time_control
start_year                          = 2023
start_month                         = 01
start_day                           = 10
start_hour                          = 00
start_minute                        = 00
start_second                        = 00
end_year                            = 2023
end_month                           = 01
end_day                             = 10
end_hour                            = 06
end_minute                          = 00
end_second                          = 00
run_hours=6,
interval_seconds=21600,
/
&domains
s_we                                = 1
e_we                                = 61
s_sn                                = 1
e_sn                                = 61
s_vert                              = 1
e_vert                              = 45
dx                                  = 9000
dy                                  = 9000
time_step = 120,
/
&dfi_control
/
&tc
/
&physics
mp_physics                          =  3,   3,
ra_lw_physics                       =  1,   1,
ra_sw_physics                       =  1,   1,
radt                                =  30,  30,
sf_sfclay_physics                   =  1,   1,
sf_surface_physics                  =  2,   2,
bl_pbl_physics                      =  1,   1,
bldt                                =  0,   0,
cu_physics                          =  1,   1,
cudt                                =  5,   5,
isfflx                              = 1,
ifsnow                              = 0,
icloud                              = 1,
surface_input_source                = 3,
num_soil_layers                     = 4,
mp_zero_out                         = 0,
maxiens                             = 1,
maxens                              = 3,
maxens2                             = 3,
maxens3                             = 16,
ensdim                              = 144,
num_land_cat                        = 21,
co2tf=0,
/
&scm
/
&fdda
/
&dynamics
hybrid_opt                          = 0,
use_theta_m                         = 0,
diff_opt                            = 1,
km_opt                              = 1,
w_damping                           = 0,
diff_6th_opt                        = 0,     0,
diff_6th_factor                     = 0.12,  0.12,
damp_opt                            = 0,
base_temp                           = 290.,
iso_temp                            = 200.,
zdamp                               = 5000.,  5000.,
dampcoef                            = 0.2,    0.2,
khdif                               = 0,      0,
kvdif                               = 0,      0,
time_step_sound                     = 0,      0,
/
&bdy_control
specified = .true., .false.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
&perturbation
trajectory_io=true,
enable_identity=false,
jcdfi_use=false,
jcdfi_diag=1,
jcdfi_penalty=1000.0,
/

ljchen1989

4dvar现在已经能嵌套了嘛？

ljf

ljchen1989 发表于 2023-11-9 13:59
4dvar现在已经能嵌套了嘛？

应该可以吧，我看有人文章里做WRF 4Dvar用的嵌套域：“two nested domains are configured with ...”，还是说他用的其他什么方法

zzmmzz

同学，你这个问题解决了么？