求助：利用Grads处理FY-2E 降水估计产品（Pre，9210格式）遇到问题

繁星一点

    之前结合本论坛中的多个帖子成功处理过FY-2E TBB产品，这次想利用类似的方法处理降水产品，可是遇到了问题，自己找不到错误原因，在此请教各位大神
————————————————————————————————————
ctl文件如下：
dset G:\AWX\20090621_30\FY2E_PRE_001_OTG_2012%m2%d2_%h200.AWX
title 2012 twpice fy2e_pre
options template
options yrev
* yrev表示y轴反向
undef 0.0
xdef 1201 linear  45 0.1
ydef 1201 linear -60 0.1
zdef 1 linear 1000 1
tdef 1 linear 06z24Aug2012 1hr
vars  1
PRE  0    -1,40,1     Geopotential Heights (m)
endvars
———————————————————————————————————————
GS文件：
'reinit'
'open G:\AWX\pre.ctl'
**************************
tt=1
while(tt<2)
*******控制坐标轴*******
'set xlint 5'
'set ylint 5'
'set xlopts 1 7 0.14'
'set ylopts 1 7 0.14'
*************************
'set t ' tt
'set grads off'
'set grid off'
'set lon 70 140'
'set lat 10 60'
*'set lon 84 91'
*'set lat 41 46'
'set rgb 21 0 0 128'
'set gxout shaded'
'set clevs  0 10 20 30 40 50 60  '
'set ccols    9 4 13 12 7 8 2 6'
*******添加地图********
'set mpdset cnworld cnhimap  '
'set map 1 1 8'
'draw map'
'd pre'
'set mpdset cnriver'
'set map 21 1 8'
'draw map'
**********************
'cbarn 1 0 5.5 8'
'draw string 2.3 0.2 time='tt''
if(rc='0')break;endif;
'enable print G:\AWX\pre_20090630.gmf '
'print'
'clear'
tt=tt+1
endwhile
pull dummy
'disable print'
'clear'
return
——————————————————————————————————————————————————
根据两个文件出图如下：

grads中的提示是：




请教各位：
  为什么画出来的图中不显示降水信息（根据grads的提示并没有什么错误）？我用micaps打开awx文件显示是有降水信息的。
我猜想出现这个错误是不是跟ctl文件中
vars 1
PRE 0 -1,40,1 Geopotential Heights (m)
这两句有关？因为在网上没有找到关于fy-2e产品的格式说明

FY2E_PRE_001_OTG_20120824_0600.AWX (983.4 KB, 下载次数: 23)

2012-8-27 09:54 上传

点击文件名下载附件

繁星一点

Forcast 发表于 2012-8-27 11:40

                               
登录/注册后可看大图

可以到国家卫星中心网站去下载相关说明，论坛提供的GS画出的TBB是不对的，你看看分析一下最大值和最小值就知 ...

国家卫星中心网站的提供的是hdf格式数据的说明 而且看不太懂  主要是不知道从哪里筛选有用的信息 以及我用grads编程时究竟哪里不对  所以觉得有点无从下手了

Forcast

这个网站是好几个帖子对AWX的读取是错误的。AWX文件有文件头，用GRADS直接读取是错位的

繁星一点

Forcast 发表于 2012-8-27 10:35

                               
登录/注册后可看大图

这个网站是好几个帖子对AWX的读取是错误的。AWX文件有文件头，用GRADS直接读取是错位的

请问能说得详细一点吗？或者您能不能提供AWX的格式说明
我也觉得应该是grads读数的格式不对  
不过我之前用论坛提供的gs文件能画TBB的图

Forcast

可以到国家卫星中心网站去下载相关说明，论坛提供的GS画出的TBB是不对的，你看看分析一下最大值和最小值就知道了

wazky

学姐您好，可不可以把画这个图里面用到的地图集发给我，最近在改论文，中国地图中的南海诸岛一直画不出来，请学姐帮帮忙

cpinging

我用的处理TBB的ctl文件好像是

vars 1
tbb 1  -1,40,1 Geopotential Heights (m)
endvars

楼主的
pre 0  -1,40,1 Geopotential Heights (m)
跟这个有关系没？

zhaolanxiai

把undef 0.0改成undef 9999,是可以作图出来，但是我对照micpas，感觉有问题

繁星一点

wazky 发表于 2012-8-28 10:04

                               
登录/注册后可看大图

学姐您好，可不可以把画这个图里面用到的地图集发给我，最近在改论文，中国地图中的南海诸岛一直画不出来， ...

这个动力论坛grads版有下的
我把我的传给你吧
map.rar (1.77 MB, 下载次数: 41)

2012-8-30 16:34 上传

点击文件名下载附件

繁星一点

cpinging 发表于 2012-8-28 16:10

                               
登录/注册后可看大图

我用的处理TBB的ctl文件好像是

vars 1

我试过用 pre 1 -1,40,1  Geopotential Heights (m)  也画不出来
而且我也不知道这个语句到底是什么意思