wrf后处理及原始输出资料碰到的一点问题

追风的阳光

自己刚接触wrf，想利用fnl资料作为初始场（多个时次fnl资料），运行了一个简单的降水个例运行时前面一切顺利，最后进行ARWpost后处理时，./compile后出现make：nothing to be done for 'all'.请问这是怎么一回事呢？
网上查阅说这不是错误或make clean 再编译，我没有管他，继续./wrf.exe
然后用ARWpost处理输出后的资料，用grads打开，各变量值得范围都是-3e.+38到3e.+38这么大，




上图是ARWpost处理后的U分量，时间、层次都是取的1。下图是时间取8，层次取8后出现的问题。U分量如此，其他分量也是这样。

然后我将原始输出资料用meteoinfo打开，各变量取值正常，如下图：


但是没有U、V风变量，请问上面几个错误是怎么回事呢？谢谢！
下面贴上我的namelist
namelist.wps:
[zssapr@swang01 WPS]$ vi namelist.wps

io_form_geogrid = 2,
/

&geogrid
parent_id         =   1,   1,
parent_grid_ratio =   1,   3,
i_parent_start    =   1,  46,
j_parent_start    =   1,  25,
e_we              =  121,64 ,
e_sn              =  76,64,
geog_data_res     = '5m','2m',
dx = 30000,
dy = 30000,
map_proj = 'lambert',
ref_lat   =  32.5,
ref_lon   = 120,
truelat1  =  30.0,
truelat2  =  60.0,
stand_lon = 120,
geog_data_path = '/swang02/geog'
/

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

&metgrid
fg_name = 'FILE'
io_form_metgrid = 2,
/
namelist.input:

[zssapr@swang01 em_real]$ vi namelist.input
&time_control
run_days                            = 0,
run_hours                           = 120,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2011, 2011,
start_month                         = 07,   07,
start_day                           = 10,   10,
start_hour                          = 00,   00,
start_minute                        = 00,   00,
start_second                        = 00,   00,
end_year                            = 2011, 2011,
end_month                           = 07,   07,
end_day                             = 15,   15,
end_hour                            = 00,   00,
end_minute                          = 00,   00,
end_second                          = 00,   00,
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,
history_interval                    = 180,  60,
frames_per_outfile                  = 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                             = 2,
e_we                                = 121,    64,
e_sn                                = 76,    64,
e_vert                              = 28,    28,
eta_levels = 1.000,0.990,0.978,0.964,0.946,0.922,0.894,0.860,0.817,0.766,0.707,0.644,0.576,0.507,0.444,0.380,0.324,0.273,0.228,0.188,0.152,0.121,0.093,0.069,0.048,0.029,0.014,0.000,
p_top_requested                     = 5000,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
dx                                  = 30000, 10000,
dy                                  = 30000, 10000,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,     46,
j_parent_start                      = 1,     25,
parent_grid_ratio                   = 1,     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,

namelist.ARWpost:
[zssapr@swang01 ARWpost]$ vi namelist.ARWpost

&datetime
start_date = '2011-07-10_00:00:00',
end_date   = '2011-07-15_00:00:00',
interval_seconds = 21600,
tacc = 0,
debug_level = 0,
/

&io
input_root_name = '/home/zssapr/WRF3.5.1/20110710/WRFV3/test/em_real/wrfout_d02_2011-07-10_00:00:00'
output_root_name = './rainout02'
plot = 'all_list'
fields = 'height,pressure,tk,tc'
output_type='grads'
mercator_defs = .true.
/
split_output = .true.
frames_per_outfile = 2
output_type= 'grads'
output_type= 'v5d'

plot = 'all'
plot = 'list'
plot = 'all_list'
! Below is a list of all available diagnostics
fields = 'u,v,height,geopt,theta,tc,tk,td,td2,rh,rh2,umet,vmet,pressure,u10m,v10m,wdir,wspd,wd10,ws10,slp,mcape,mcin,lcl,lfc,cape,cin,dbz,max_dbz,clfr'


&interp
interp_method = 0,

interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
/
extrapolate = .true.

interp_method = 0,     ! 0 is model levels, -1 is nice height levels, 1 is user specified pressure/height

interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
interp_levels = 0.25, 0.50, 0.75, 1.00, 2.00, 3.00, 4.00, 5.00, 6.00, 7.00, 8.00, 9.00, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 17.0, 18.0, 19.0, 20.0,

mofangbao

数据很小或者很大是缺测值问题 下面的错误应该是那个变量没那么多层次的数据

追风的阳光

mofangbao 发表于 2014-4-11 09:16
数据很小或者很大是缺测值问题 下面的错误应该是那个变量没那么多层次的数据

但是我的fnl资料U、V都是有的26层数据啊

黄小仙儿

我也遇到了这个问题，grads打开时乱码，metoinfo可以打开，果断用NCL了，处理wrf数据最方便

追风的阳光

黄小仙儿 发表于 2014-4-11 10:09
我也遇到了这个问题，grads打开时乱码，metoinfo可以打开，果断用NCL了，处理wrf数据最方便

谢谢，我要学学NCL了。。

MeteoInfo

更新了MeteoInfo来处理这种情况，参加此帖：处理读WRF模式输出的netCDF数据stagger维问题
http://bbs.06climate.com/forum.p ... 081&fromuid=106
(出处: 气象家园)