ndown.exe执行后再执行wrf.exe结果出错

changsc

大神们好，我要做一个两层嵌套并用ndown做成one-way Nested用作模式输入结果。按照操作说明执行ndown后没问题，再执行wrf.exe结果就错了。
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:     431
  RIBX never exceeds RIC, RIB(i,kte) =             NaN THETAV(i,1) =             NaN MOL=            NaN TCONV =   0.0000000E+00 WST =             NaN KMIX =            1 UST =             NaN TST =             NaN U,V =             NaN            NaN I,J=           2           2
-------------------------------------------
[cli_0]: aborting job:
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

网上有的说改time_step，可怎么改小都是这个错误。
下面是进行该步的namelist（已按照手册用ndown生成了wrfbdy.d02,将四个*d02文件都改成了*d01）
&time_control
run_days                            = 0,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2017, 2017, 2017,
start_month                         = 09,   09,   09,
start_day                           = 01,   01,   01,
start_hour                          = 00,   00,   00,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2017, 2017, 2017,
end_month                           = 09,   09,   09,
end_day                             = 02,   02,   02,
end_hour                            = 00,   00,   00,
end_minute                          = 00,   00,   00,
end_second                          = 00,   00,   00,
interval_seconds                    = 21600
input_from_file                     = .false.,.true.,.true.,
history_interval                    = 60,  60,   60,
frames_per_outfile                  = 24, 24, 24,
restart                             = .false.,
restart_interval                    = 1440,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
io_form_auxinput2                   = 2
io_form_auxinput4                   = 2
auxinput4_inname                    = "wrflowinp_d<domain>"
auxinput4_interval                  = 360
auxinput4_end_h                     = 744
debug_level                         = 0
/

&domains
time_step                           = 60,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
e_we                                = 112,    112,    124,
e_sn                                = 106,    106,    124,
e_vert                              = 30,    30,    30,
p_top_requested                     = 5000,
eta_levels                          = 1.000, 0.993, 0.983, 0.97, 0.954,
                                       0.934, 0.909, 0.88, 0.8341715, 0.7883431,
                                       0.7425146, 0.6966861, 0.6167657, 0.5439157, 0.4776143,
                                       0.417372, 0.3627303, 0.3132596, 0.2685581, 0.2282499,
                                       0.1919836, 0.1594311, 0.1302862, 0.1042633, 0.08109632,
                                       0.06053737, 0.04235569, 0.02633654, 0.01228005, 0.
num_metgrid_levels                  = 48,
num_metgrid_soil_levels             = 4,
dx                                  = 12000, 12000,  4000,
dy                                  = 12000, 12000,  4000,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,     102,   40,
j_parent_start                      = 1,     35,    36,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     3,     3,
feedback                            = 1,
smooth_option                       = 0,
hypsometric_opt                     = 1
/

&physics
num_land_cat                        = 24,
mp_physics                          = 10,     10,     10,
progn                               = 1,     1,     1,
ra_lw_physics                       = 4,     4,     4,
ra_sw_physics                       = 4,     4,     4,
paerlev                             = 29,
levsiz                              = 59,
cam_abs_dim1                        = 4,
cam_abs_dim2                        = 36,
radt                                = 20,    20,    20,
sf_sfclay_physics                   = 7,     7,     7,
sf_surface_physics                  = 7,     7,     7,
bl_pbl_physics                      = 7,     7,     7,
bldt                                = 0,     0,     0,
cu_physics                          = 1,     1,     0,
cudt                                = 5,     5,     5,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 2,
sst_update                          = 1,
sf_urban_physics                    = 0,     1,     1,
slope_rad                           = 1,     1,     1,
topo_shading                        = 1,     1,     1,
shadlen                             = 25000.,
mp_zero_out                         = 2,
/

&fdda
grid_fdda                           = 1,     1,     1,
gfdda_inname                        = "wrffdda_d<domain>",
gfdda_interval_m                    = 360,   360,   360,
gfdda_end_h                         = 10000, 10000, 10000,
grid_sfdda                          = 0,
pxlsm_soil_nudge                    = 0,
sgfdda_inname                       = "wrfsfdda_d01",
sgfdda_interval_m                   = 180
sgfdda_interval                     = 10800
sgfdda_end_h                        = 10000
fgdt                                = 0,     0,     0,
if_no_pbl_nudging_uv                = 1,     1,     1,
if_no_pbl_nudging_t                 = 1,     1,     1,
if_no_pbl_nudging_q                 = 1,     1,     1,
if_zfac_uv                          = 0,     0,     0,
k_zfac_uv                           = 13,
if_zfac_t                           = 0,     0,     0,
k_zfac_t                            = 13,
if_zfac_q                           = 0,     0,     0,
k_zfac_q                            = 13,
guv                                 = 0.00005, 0.00005, 0.00005,
gt                                  = 0.00005, 0.00005, 0.00005,
gq                                  = 0.00001, 0.00001, 0.00001,
guv_sfc                             = 0.0000, 0.0000, 0.0000,
gt_sfc                              = 0.0000, 0.0000, 0.0000,
gq_sfc                              = 0.0000, 0.0000, 0.0000,
if_ramping                          = 1,
dtramp_min                          = 60.0,
io_form_gfdda                       = 2,
rinblw                              = 250.0,
/

&dynamics
w_damping                           = 1,
diff_opt                            = 1,
km_opt                              = 4,
diff_6th_opt                        = 2,      2,      2,
diff_6th_factor                     = 0.12,   0.12,   0.12,
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.,
moist_adv_opt                       = 2,      2,      2,
tke_adv_opt                         = 2,      2,      2,
scalar_adv_opt                      = 2,      2,      2,
use_baseparam_fr_nml                = .true.,
/

&bdy_control
spec_bdy_width                      = 5,
spec_zone                           = 1,
relax_zone                          = 4,
specified                           = .true., .false.,.false.,
nested                              = .false., .true., .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
Type  :quit<Enter>  to exit Vim              

绛紫

你这个namelist是正确的？你的时间控制那里，run的都是0？单向嵌套，feedback是0吧，积分步长是要改大不是改小

问题少女

1. run_days = 0,应该改为：run_days = 1，
2. input_from_file  = .false.,.true.,.true.,应该改为：input_from_file  = .true.,.true.,.true.,
3. io_form_auxinput2 = 2应该去掉
4. i_parent_start  = 1,     102,   40,
    j_parent_start  = 1,     35,    36,            也需要改；
5. feedback  = 1,应该改为feedback   = 0,
你先改改这些试试。

changsc

谢楼上答。据我所知使用ndown和使用feedback=0是两种不同的情况，具体是什么我也搞不清楚。现在具体细化一下问题。我认为可能是网格不匹配的问题，我在使用ndown之后生成的fine domain 的wrfinput文件与coarse domain 有两个区别可能导致此错误：
1我在使用ncdump看两个文件时发现fine domain 多了urban_wind_directions_stag = 4一项，具体如下：
fine domain--------------------------------------------
dimensions:
        Time = UNLIMITED ; // (1 currently)
        DateStrLen = 19 ;
        west_east = 111 ;
        south_north = 105 ;
        bottom_top = 29 ;
        bottom_top_stag = 30 ;
        soil_layers_stag = 2 ;
        west_east_stag = 112 ;
        south_north_stag = 106 ;
        DIM0009 = 5 ;
        land_cat_stag = 24 ;
        soil_cat_stag = 16 ;
       urban_wind_directions_stag = 4 ;
        num_ext_model_couple_dom_stag = 1 ;
coarse domain -------------------------------------------------
dimensions:
        Time = UNLIMITED ; // (1 currently)
        DateStrLen = 19 ;
        west_east = 111 ;
        south_north = 105 ;
        bottom_top = 29 ;
        bottom_top_stag = 30 ;
        soil_layers_stag = 2 ;
        west_east_stag = 112 ;
        south_north_stag = 106 ;
        DIM0009 = 5 ;
        land_cat_stag = 24 ;
        soil_cat_stag = 16 ;
        num_ext_model_couple_dom_stag = 1 ;
2.fine domain 的变量变多了几个，比如 LF_URB2D
可我在使用ndown前后并未改变除了Users' Guide里说的需要调整的网格的信息以外的信息。
请大神们解答