wrfout转换为ctl和dat之后grads画阴影图所有数据都没有

大侠MVP

grads输完命令后显示contouring：-2e+38 to 2.5e+38 interval 5e+37

arwpost转换后的ctl文件：
dset D:/wrfoutctr.dat
undef 1.e30
title  OUTPUT FROM WRF V2.2.1 MODEL
pdef  199 179 lcc  36.000  117.500  100.000   90.000  60.00000  30.00000  117.50000  10000.000  10000.000
xdef  586 linear  104.29045   0.04504504
ydef  394 linear   26.93918   0.04504504
zdef   30 linear 1 1  
tdef   25 linear 06Z18SEP2005      60MN      
VARS    2
RAINC          1  0  ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
RAINNC         1  0  ACCUMULATED TOTAL GRID SCALE PRECIPITATION (mm)
ENDVARS

大侠MVP

补充下namelist.arw:

&datetime
start_date = '2005-09-18_06:00:00',
end_date   = '2005-09-19_06:00:00',
interval_seconds = 3600,
tacc = 0,
debug_level = 0,
/

&io
io_form_input  = 2,
input_root_name = '/cmb/g5/shaochl/WRFV2/WRFV2/run/wrfout_d01_2005-09-18_06v2'
output_root_name = '/cmb/g5/shaochl/wrfoutctr'
plot = 'list'
fields = 'RAINC,RAINNC'
output_type = 'grads'
mercator_defs = .false.
/
! split_output = .false.
! 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 = '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 = .false.

! 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,

兰溪之水

options byteswapped

又是那隻貓

你可以看看里面的数据到底有没有

大侠MVP

兰溪之水 发表于 2014-9-19 11:26
options byteswapped

我用的ibm的xlf编译器，ctl里本来就没有这项。而且我也试过了把这个加上，输出还是一样

大侠MVP

又是那隻貓 发表于 2014-9-19 12:42
你可以看看里面的数据到底有没有

dat里面的数据吗？我只能用写字板打开，里面都是非零的，但是不知道具体数字是多少，有其他方法看吗？用set gxout contour画出来的图满屏幕都是很小很小孤立的圈，用shaded输出的就是我传上来这张图

denglin19901030

我也遇到了同样的问题，直接用ncl看wrfout的数据时有的，就是不知道为什么画出来就是一个一个小圈？

大侠MVP

兰溪之水 发表于 2014-9-19 11:26
options byteswapped

感谢兰溪。。按照你的提示问题已经得到解决

qixiangren34

看你的网格中心位置，你是河北的，还是山东的？