求教 关于用gfs数据进行wrf预报

aspen

从ftpprd.ncep.noaa.gov下载了gfs的数据，数据名为gfs.t00z.pgrb2.1p00.f0**
修改了namelist.wps和namelist.input   
运行了  link_grib.csh  
geogrid.exe
ungrib.exe
metgrid.exe
这些看上去都正常生成文件  
但运行real.exe 显示如下
再运行wrf.exe 也报错。  
求指教可能是哪儿的问题
多谢。



Ntasks in X            1 , ntasks in Y            1
--- WARNING: traj_opt is zero, but num_traj is not zero; setting num_traj to zero.
--- 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.
--- NOTE: bl_pbl_physics /= 4, implies mfshconv must be 0, resetting
Need MYNN PBL for icloud_bl = 1, resetting to 0
--- NOTE: RRTMG radiation is not used, setting:  o3input=0 to avoid data pre-processing
--- NOTE: num_soil_layers has been set to      4
REAL_EM V3.8.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 ,             144493204  bytes allocated
Time period #   1 to process = 2017-02-13_00:00:00.
Time period #   2 to process = 2017-02-13_06:00:00.
Time period #   3 to process = 2017-02-13_12:00:00.
Time period #   4 to process = 2017-02-13_18:00:00.
Time period #   5 to process = 2017-02-14_00:00:00.
Time period #   6 to process = 2017-02-14_06:00:00.
Time period #   7 to process = 2017-02-14_12:00:00.
Time period #   8 to process = 2017-02-14_18:00:00.
Time period #   9 to process = 2017-02-15_00:00:00.
Total analysis times to input =    9.

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

Domain  1: Current date being processed: 2017-02-13_00:00:00.0000, which is loop #   1 out of    9
configflags%julyr, %julday, %gmt:        2017          44   0.00000000   
metgrid input_wrf.F first_date_input = 2017-02-13_00:00:00
metgrid input_wrf.F first_date_nml = 2017-02-13_00:00:00
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:     864
input_wrf.F: SIZE MISMATCH:  namelist ide,jde,num_metgrid_levels=          74          61          27 ; input data ide,jde,num_metgrid_levels=          74          61          32
-------------------------------------------
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:     864
input_wrf.F: SIZE MISMATCH:  namelist ide,jde,num_metgrid_levels=          74          61          27 ; input data ide,jde,num_metgrid_levels=          74          61          32
-------------------------------------------

愛、星晴

你把namelist.input 里面的num_metgrid_levels ＝32 就可以了

aspen

愛、星晴 发表于 2017-2-14 00:12
你把namelist.input 里面的num_metgrid_levels ＝32 就可以了

果然  谢谢  
不过这个数是从哪儿得到的啊？  
谢谢。  

aspen

愛、星晴 发表于 2017-2-14 00:12
你把namelist.input 里面的num_metgrid_levels ＝32 就可以了

看到了  报错里说不匹配的。  

loftydongshi

去年5月12日开始gfs的边界层数据变成32层了，之前的是27层

wqowen

楼主，您好。我最近也在学习WRF预报，渴望交流学习。我也下载了一样的数据，请问namelist.wps 和namelist.input是如何设定的？能够发一个casegei学习一下。350052646@qq.com 不胜感激

aspen

wqowen 发表于 2017-9-10 22:30
楼主，您好。我最近也在学习WRF预报，渴望交流学习。我也下载了一样的数据，请问namelist.wps 和namelist.i ...

namelist.wps


&share,
wrf_core = 'ARW',
max_dom = 1,
start_date = '2017-09-03_00:00:00','2017-08-16_00:00:00'
end_date   = '2017-09-10_00:00:00','2017-08-18_00:00:00'
interval_seconds = 10800,
io_form_geogrid = 2,
/,
,
&geogrid,
parent_id         =   1,1
parent_grid_ratio =   1,3
i_parent_start    =   1,36
j_parent_start    =   1,29
e_we              =  100,88
e_sn              =  68,70
geog_data_res     = '10m','2m'
Dx = 76500,
Dy = 76500,
map_proj = 'lambert',
ref_lat   =  40.938,
ref_lon   = 128.414,
Truelat1  =  40.938,
Truelat2  =  40.938,
stand_lon = 128.414,
geog_data_path = '/root/WPS/geog',
/,
,
&ungrib,
out_format = 'WPS',
prefix = 'FILE',
/,
,
&metgrid,
fg_name = 'FILE',
io_form_metgrid = 2,
/,

namelist.input


&time_control
run_days                            = 7,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2017, 2017, 2000,
start_month                         = 09,   08,   01,
start_day                           = 03,   16,   24,
start_hour                          = 00,   00,   12,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2017, 2017, 2000,
end_month                           = 09,   08,   01,
end_day                             = 10,   23,   25,
end_hour                            = 00,   00,   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                                = 100,    88,   94,
e_sn                                = 68,    70,    91,
e_vert                              = 32,    32,    30,
p_top_requested                     = 5000,
num_metgrid_levels                  = 32,
num_metgrid_soil_levels             = 4,
dx                                  = 76500, 7600,  3333.33,
dy                                  = 76500, 7600,  3333.33,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,     36,    30,
j_parent_start                      = 1,     29,    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,      1,      1,
km_opt                              = 4,      4,      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,  
do_avgflx_em                        = 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,
/

wqowen

aspen 发表于 2017-9-13 14:44
namelist.wps

非常感谢，楼主分享。

wqowen

aspen 发表于 2017-9-13 14:44
namelist.wps

非常感谢，楼主分享。

金玉姻缘

楼主，本人刚开始做GFS预报，之前一直是模拟，没预报，能提供下怎么修改namelist吗