WRF模式的wrfbdy大小是否有限制？

870744170

最近做了一个季节模拟，从2012年1月1日~2012年5月1日，当我设置interval_seconds=21600时，real.exe只能处理到2012年的3月5日，也没有报错。当我设置interval_seconds=43200时，才可以运行到5月1日。想问问是不是WRF模式对wrfbdy有限制？还是其他什么原因？

有风吹过

没有限制，把namelist贴出来我看看

870744170

有风吹过 发表于 2015-6-1 08:05
没有限制，把namelist贴出来我看看

我用的WRF版本是3.5，run_days这些已经可以缺省了，下面是我出错的时候的namelist.input
&time_control            
run_days                 = ,
run_hours                = ,
run_minutes              = ,
run_seconds              = ,
start_year               = 2012,
start_month              = 01,
start_day                = 02,
start_hour               = 00,
start_minute             = 00,
start_second             = 00,
end_year                 = 2012,
end_month                = 05,
end_day                  = 01,
end_hour                 = 00,
end_minute               = 00,
end_second               = 00,
interval_seconds         = 21600,
input_from_file          = .true.,
history_interval         = 360,
frames_per_outfile       = 10000,
restart                  = .false.,
restart_interval         = 10000,
io_form_history          = 2,
io_form_restart          = 2,
io_form_input            = 2,
io_form_boundary         = 2,
debug_level              = 0,
/

&domains                 
time_step                = 180,
time_step_fract_num      = 0,
time_step_fract_den      = 1,
max_dom                  = 1,
s_we                     = 1,
e_we                     = 160,
s_sn                     = 1,
e_sn                     = 140,
s_vert                   = 1,
e_vert                   = 38,
num_metgrid_levels       = 38,
num_metgrid_soil_levels             = 4,

dx                       = 30000,
dy                       = 30000,
grid_id                  = 1,
parent_id                = 1,
i_parent_start           = 1,
j_parent_start           = 1,
parent_grid_ratio        = 1,
parent_time_step_ratio   = 1,
feedback                 = 1,
smooth_option            = 0,
/

&physics                 
mp_physics               = 2,        2,
ra_lw_physics            = 1,        1,
ra_sw_physics            = 1,        1,
radt                     = 30,       30,
sf_sfclay_physics        = 1,        1,
sf_surface_physics       = 1,        1,
bl_pbl_physics           = 1,        1,
bldt                     = 0,        0,
cu_physics               = 1,        1,
cudt                     = 5,        5,
isfflx                   = 1,
ifsnow                   = 0,
icloud                   = 1,
surface_input_source     = 1,
num_soil_layers          = 5,
mp_zero_out              = 0,
maxiens                  = 1,
maxens                   = 3,
maxens2                  = 3,
maxens3                  = 16,
ensdim                   = 144,
/

&fdda                    
/

&dynamics               
w_damping                = 0,
diff_opt                 = 1,
km_opt                   = 4,
base_temp                = 290.,
damp_opt                 = 0,
zdamp                    = 5000.,    5000.,
dampcoef                 = 0.01,     0.01,
khdif                    = 0,        0,
kvdif                    = 0,        0,
smdiv                    = 0.1,      0.1,
emdiv                    = 0.01,     0.01,
epssm                    = 0.1,      0.1,
time_step_sound          = 4,        4,
h_mom_adv_order          = 5,        5,
v_mom_adv_order          = 3,        3,
h_sca_adv_order          = 5,        5,
v_sca_adv_order          = 3,        3,
non_hydrostatic          = .true.,   .true.,
/

&bdy_control            
spec_bdy_width           = 5,
spec_zone                = 1,
relax_zone               = 4,
specified                = .true.,  .false.,
periodic_x               = .false.,  .false.,
symmetric_xs             = .false.,  .false.,
symmetric_xe             = .false.,  .false.,
open_xs                  = .false.,  .false.,
open_xe                  = .false.,  .false.,
periodic_y               = .false.,  .false.,
symmetric_ys             = .false.,  .false.,
symmetric_ye             = .false.,  .false.,
open_ys                  = .false.,  .false.,
open_ye                  = .false.,  .false.,
nested                   = .false.,   .true.,
/

&grib2                  
/

&namelist_quilt         
nio_tasks_per_group      = 0,
nio_groups               = 1,
/
                                                                                                                                                      

870744170

有风吹过 发表于 2015-6-1 08:05
没有限制，把namelist贴出来我看看

我用的WRF版本是3.5，run_days这些已经可以缺省了，下面是我出错的时候的namelist.input
&time_control            
run_days                 = ,
run_hours                = ,
run_minutes              = ,
run_seconds              = ,
start_year               = 2012,
start_month              = 01,
start_day                = 02,
start_hour               = 00,
start_minute             = 00,
start_second             = 00,
end_year                 = 2012,
end_month                = 05,
end_day                  = 01,
end_hour                 = 00,
end_minute               = 00,
end_second               = 00,
interval_seconds         = 21600,
input_from_file          = .true.,
history_interval         = 360,
frames_per_outfile       = 10000,
restart                  = .false.,
restart_interval         = 10000,
io_form_history          = 2,
io_form_restart          = 2,
io_form_input            = 2,
io_form_boundary         = 2,
debug_level              = 0,
/

&domains                 
time_step                = 180,
time_step_fract_num      = 0,
time_step_fract_den      = 1,
max_dom                  = 1,
s_we                     = 1,
e_we                     = 160,
s_sn                     = 1,
e_sn                     = 140,
s_vert                   = 1,
e_vert                   = 38,
num_metgrid_levels       = 38,
num_metgrid_soil_levels             = 4,

dx                       = 30000,
dy                       = 30000,
grid_id                  = 1,
parent_id                = 1,
i_parent_start           = 1,
j_parent_start           = 1,
parent_grid_ratio        = 1,
parent_time_step_ratio   = 1,
feedback                 = 1,
smooth_option            = 0,
/

&physics                 
mp_physics               = 2,        2,
ra_lw_physics            = 1,        1,
ra_sw_physics            = 1,        1,
radt                     = 30,       30,
sf_sfclay_physics        = 1,        1,
sf_surface_physics       = 1,        1,
bl_pbl_physics           = 1,        1,
bldt                     = 0,        0,
cu_physics               = 1,        1,
cudt                     = 5,        5,
isfflx                   = 1,
ifsnow                   = 0,
icloud                   = 1,
surface_input_source     = 1,
num_soil_layers          = 5,
mp_zero_out              = 0,
maxiens                  = 1,
maxens                   = 3,
maxens2                  = 3,
maxens3                  = 16,
ensdim                   = 144,
/

&fdda                    
/

&dynamics               
w_damping                = 0,
diff_opt                 = 1,
km_opt                   = 4,
base_temp                = 290.,
damp_opt                 = 0,
zdamp                    = 5000.,    5000.,
dampcoef                 = 0.01,     0.01,
khdif                    = 0,        0,
kvdif                    = 0,        0,
smdiv                    = 0.1,      0.1,
emdiv                    = 0.01,     0.01,
epssm                    = 0.1,      0.1,
time_step_sound          = 4,        4,
h_mom_adv_order          = 5,        5,
v_mom_adv_order          = 3,        3,
h_sca_adv_order          = 5,        5,
v_sca_adv_order          = 3,        3,
non_hydrostatic          = .true.,   .true.,
/

&bdy_control            
spec_bdy_width           = 5,
spec_zone                = 1,
relax_zone               = 4,
specified                = .true.,  .false.,
periodic_x               = .false.,  .false.,
symmetric_xs             = .false.,  .false.,
symmetric_xe             = .false.,  .false.,
open_xs                  = .false.,  .false.,
open_xe                  = .false.,  .false.,
periodic_y               = .false.,  .false.,
symmetric_ys             = .false.,  .false.,
symmetric_ye             = .false.,  .false.,
open_ys                  = .false.,  .false.,
open_ye                  = .false.,  .false.,
nested                   = .false.,   .true.,
/

&grib2                  
/

&namelist_quilt         
nio_tasks_per_group      = 0,
nio_groups               = 1,
/
                                                                                                                                                      

有风吹过

870744170 发表于 2015-6-1 09:19
我用的WRF版本是3.5，run_days这些已经可以缺省了，下面是我出错的时候的namelist.input
&time_control  ...

你设置maxiens选项有什么作用呢？你的模拟结果是不是出来了多个wrfout？把restart_interval设置的更大些就不用出现多个wrfout， 我暂时还看不了问题，你把namelist.wps贴出来再看看

870744170

有风吹过 发表于 2015-6-1 11:13
你设置maxiens选项有什么作用呢？你的模拟结果是不是出来了多个wrfout？把restart_interval设置的更大些 ...

maxiens选项是WRF自带的，一直就没动过，不知道是啥意思，解释解释呗！WPS的话，有两个，一个是生成SFC的，还有一个是生成PRESS的，这里只把SFC的贴出来哈，PRESS的话就是把ungrib的prefix改成PRESS
&share
wrf_core = 'ARW',
max_dom = 1,
start_date = '2012-01-02_00:00:00',
end_date   = '2012-05-01_00:00:00',
interval_seconds = 21600,
io_form_geogrid = 2,
debug_level = 0,
/

&geogrid
parent_id         = 1,
parent_grid_ratio = 1,
i_parent_start    = 1,
j_parent_start    = 1,
e_we          = 160,
e_sn          = 140,
geog_data_res = '5m',
dx = 30000,
dy = 30000,
map_proj =  'lambert',
ref_lat   = 36.5,
ref_lon   = 102.3,
truelat1  = 30,
truelat2  = 60,
stand_lon = 102.3,
geog_data_path = '/public/data/wrf/geog',
/

&ungrib
out_format = 'WPS',
prefix = 'SFC',
/

&metgrid
fg_name = 'PRESS','SFC'
io_form_metgrid = 2,

有风吹过

870744170 发表于 2015-6-1 13:54
maxiens选项是WRF自带的，一直就没动过，不知道是啥意思，解释解释呗！WPS的话，有两个，一个是生成SFC的 ...

看出来了，你用的是ERA-interim驱动资料，你看下你的3月5日附件的met_em.d01有没有问题吧

870744170

有风吹过 发表于 2015-6-1 14:37
看出来了，你用的是ERA-interim驱动资料，你看下你的3月5日附件的met_em.d01有没有问题吧

额，用的CFS资料驱动的，3月5日的met_em都查过了，和其他时刻的一样。我还把3月5日和3月6日的flxf和pgbf的驱动资料重新下载了一遍，都没有问题诶。。。

有风吹过

870744170 发表于 2015-6-1 15:13
额，用的CFS资料驱动的，3月5日的met_em都查过了，和其他时刻的一样。我还把3月5日和3月6日的flxf和pgbf ...

难道是用的资料问题？为什么用CFS资料，而不用其他的呢？这套资料有什么优点？

870744170

有风吹过 发表于 2015-6-1 15:35
难道是用的资料问题？为什么用CFS资料，而不用其他的呢？这套资料有什么优点？

主要是能自动获取吧，ECMWF的话是不能自动下载的吧？我们这边要做定时滚动预报，要是自己天天手动下载，那估计得烦死了。。。