求助！wrf 3.4.1运行./real.exe出现问题

咫尺天涯BLSL

大家好，我在使用GFS数据跑WRFV3.4.1，在运行到./real.exe时出现了以下问题

--- NOTE: sst_update is 0, setting io_form_auxinput4 = 0 and auxinput4_interval = 0 for all domains
--- NOTE: grid_fdda is 0 for domain      1, setting gfdda interval and ending time to 0 for that domain.
--- NOTE: both grid_sfdda and pxlsm_soil_nudge are 0 for domain      1, setting sgfdda interval and ending time to 0 for that domain.
--- NOTE: obs_nudge_opt is 0 for domain      1, setting obs nudging interval and ending time to 0 for that domain.
bl_pbl_physics /= 4, implies mfshconv must be 0, resetting
--- NOTE: num_soil_layers has been set to      4
REAL_EM V3.4.1 PREPROCESSOR
  *************************************
  Parent domain
  ids,ide,jds,jde             1           74            1           61
  ims,ime,jms,jme            -4           79           -4           66
  ips,ipe,jps,jpe             1           74            1           61
  *************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
    alloc_space_field: domain             1 ,                 130679256  bytes allocated

Time period #   1 to process = 2017-04-22_00:00:00.
Time period #   2 to process = 2017-04-22_03:00:00.
Time period #   3 to process = 2017-04-22_06:00:00.
Time period #   4 to process = 2017-04-22_09:00:00.
Time period #   5 to process = 2017-04-22_12:00:00.
Time period #   6 to process = 2017-04-22_15:00:00.
Time period #   7 to process = 2017-04-22_18:00:00.
Time period #   8 to process = 2017-04-22_21:00:00.
Time period #   9 to process = 2017-04-23_00:00:00.
Time period #  10 to process = 2017-04-23_03:00:00.
Time period #  11 to process = 2017-04-23_06:00:00.
Time period #  12 to process = 2017-04-23_09:00:00.
Time period #  13 to process = 2017-04-23_12:00:00.
Time period #  14 to process = 2017-04-23_15:00:00.
Time period #  15 to process = 2017-04-23_18:00:00.
Time period #  16 to process = 2017-04-23_21:00:00.
Time period #  17 to process = 2017-04-24_00:00:00.
Time period #  18 to process = 2017-04-24_03:00:00.
Time period #  19 to process = 2017-04-24_06:00:00.
Time period #  20 to process = 2017-04-24_09:00:00.
Time period #  21 to process = 2017-04-24_12:00:00.
Time period #  22 to process = 2017-04-24_15:00:00.
Time period #  23 to process = 2017-04-24_18:00:00.

Time period #  24 to process = 2017-04-24_21:00:00.
Time period #  25 to process = 2017-04-25_00:00:00.
Total analysis times to input =   25.

-----------------------------------------------------------------------------

Domain  1: Current date being processed: 2017-04-22_00:00:00.0000, which is loop #   1 out of   25
configflags%julyr, %julday, %gmt:         2017          112    0.000000
  metgrid input_wrf.F first_date_input = 2017-04-22_00:00:00
  metgrid input_wrf.F first_date_nml = 2017-04-22_00:00:00
d01 2017-04-22_00:00:00 Timing for input          0 s.
d01 2017-04-22_00:00:00          flag_soil_layers read from met_em file is  1
Max map factor in domain 1 =  0.99. Scale the dt in the model accordingly.
d01 2017-04-22_00:00:00  t(i,j,k) was 0 at             1            1            1 , setting Qv to 0
d01 2017-04-22_00:00:00  t(i,j,k) was 0 at             2            1            1 , setting Qv to 0


………………
………………
………………
d01 2017-04-22_00:00:00  t(i,j,k) was 0 at            71           60           34 , setting Qv to 0
d01 2017-04-22_00:00:00  t(i,j,k) was 0 at            72           60           34 , setting Qv to 0
d01 2017-04-22_00:00:00  t(i,j,k) was 0 at            73           60           34 , setting Qv to 0
Using sfcprs3 to compute psfc
Converged znw(kte) should be about 0.0 =   -1.7695129E-08
all_dim =             1
order =             2
i,j =             1            1
p array =              NaN
f array =     0.000000
p target=              NaN             NaN             NaN             NaN
             NaN             NaN             NaN             NaN
             NaN             NaN             NaN             NaN
             NaN             NaN             NaN             NaN
             NaN             NaN             NaN             NaN
             NaN             NaN             NaN             NaN
             NaN             NaN             NaN             NaN
             NaN             NaN
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    4142
troubles, the interpolating order is too large for this few input values
-------------------------------------------
FORTRAN STOP


--------------------------------------------------------------------------------------------------------------------------------
--------------------------------------------------------------------------------------------------------------------------------
以下是我的namelist.input
&time_control
run_days                             = 0,
run_hours                            = 12,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                            = 2017, 2000, 2000,
start_month                         = 04,   01,   01,
start_day                             = 22,   24,   24,
start_hour                            = 00,   12,   12,
start_minute                         = 00,   00,   00,
start_second                         = 00,   00,   00,
end_year                              = 2017, 2000, 2000,
end_month                           = 04,   01,   01,
end_day                               = 25,   25,   25,
end_hour                              = 00,   12,   12,
end_minute                           = 00,   00,   00,
end_second                           = 00,   00,   00,
interval_seconds                     = 10800
input_from_file                       = .true.,.true.,.true.,
history_interval                       = 180,  60,   60,
frames_per_outfile                  = 1000, 1000, 1000,
restart                                   = .false.,
restart_interval                       = 5000,
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,
e_we                                = 74,    112,   94,
e_sn                                 = 61,    97,    91,
e_vert                               = 28,    28,    28,
p_top_requested                     = 5000,
num_metgrid_levels                 = 34,
num_metgrid_soil_levels           = 4,
dx                                  = 30000, 10000,  3333.33,
dy                                  = 30000, 10000,  3333.33,
grid_id                                = 1,     2,     3,
parent_id                            = 1,     1,     2,
i_parent_start                      = 1,     31,    30,
j_parent_start                      = 1,     17,    30,
parent_grid_ratio                  = 1,     3,     3,
parent_time_step_ratio          = 1,     3,     3,
feedback                              = 1,
smooth_option                      = 0
/

&physics
mp_physics                          = 3,     3,     3,
ra_lw_physics                       = 1,     1,     1,
ra_sw_physics                       = 1,     1,     1,
radt                                     = 30,    30,    30,
sf_sfclay_physics                   = 1,     1,     1,
sf_surface_physics                  = 2,     2,     2,
bl_pbl_physics                       = 1,     1,     1,
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                     = 4,
sf_urban_physics                    = 0,     0,     0,
/

&fdda
/

&dynamics
w_damping                        = 0,
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.,
nested                              = .false., .true., .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

-------------------------------------------------------------------------------------------------------
-------------------------------------------------------------------------------------------------------
说明：
1.使用的数据是2017年4月22日00时的gfs.t00z.pgrb2.0p50.anl，gfs.t00z.pgrb2.0p50.f000，gfs.t00z.pgrb2.0p50.f003，……gfs.t00z.pgrb2.0p50.f072
2.我使用了官网中的范例数据运行了一遍，可以跑通real.exe过程，也生成了wrfbdy等文件。
3.我已根据官网的提示更新了rd_grib2.F以及Vtable.GFS文件



综上，我考虑是否是我met_em文件有问题还是其他哪里出了问题？困扰几日了，实在是没办法解决，烦请各位帮忙看下，谢谢大家了！

寂寞的木偶

文件里的设置有三层，有一层 要一致

寂寞的木偶

寂寞的木偶 发表于 2017-4-27 09:47
文件里的设置有三层，有一层 要一致

时间也不对 12h

咫尺天涯BLSL

我的max_dom设置的是1层啊，后面的改不改不是不会影响么？时间这个不太懂，这个应该怎么改呢？谢谢你~

1luochen1

输出met文件的气压看看，好多NaN

PLAjim

请问楼主的问题解决了吗，我也遇到了同样的问题，怀疑是met文件出错，但是不太确定。。。

咫尺天涯BLSL

PLAjim 发表于 2017-5-25 20:46
请问楼主的问题解决了吗，我也遇到了同样的问题，怀疑是met文件出错，但是不太确定。。。

好像是gfs数据的问题，我下载的数据里没有地温，后来用fnl的数据后（也没有地温），然后把namelist里面的sf_surface_physics过程取消（数字选0），然后就可以跑了。gfs还没来得及测试是不是这个问题