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

繁星一点

Crazy_Merlin 发表于 2012-8-31 10:33

                               
登录/注册后可看大图

个人感觉有两种可能：ctl问题，数据问题。不知道降水估计产品的分辨率为多少？我用meteoinfo作图时和你传的 ...

我也觉得奇怪 个人觉得是ctl的问题  因为感觉画出来的图错位了  

Crazy_Merlin

能不能把您的代码贴出来看看？

繁星一点

Crazy_Merlin 发表于 2012-8-31 10:50

                               
登录/注册后可看大图

能不能把您的代码贴出来看看？

ctl文件：
dset G:\AWX\20090621_30\FY2E_PRE_001_OTG_20120824_0600.AWX
title 2012 twpice fy2e_pre
options yrev
* yrev表示y轴反向
undef 99.9
xdef 1001 linear  55 0.1
ydef 501 linear 0 0.1
zdef 1 levels 1000 1
tdef 1 linear 06z24Aug2012 1hr
vars  1
PRE  0    -1,40,1      FY-2E_PRE (mm)
endvars
————————————————————————————————
gs文件：
'reinit'
'open G:\AWX\pre.ctl'

**************************
tt=1
*******控制坐标轴*******
'set xlint 10'
'set ylint 10'
'set xlopts 1 7 0.14'
'set ylopts 1 7 0.14'
*************************
'set t ' tt
'set grads off'
'set grid off'
'set lon 55 155'
'set lat 0 50'
'set rgb 21 0 0 128'
'set gxout shaded'
'set clevs  0 5 15 30 70 140   '
'set ccols    0 9 4 13 7 12 8  '
*******添加地图********
'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;
'printim G:\AWX\pre_201208240600.png white '
'print'
'clear'
'disable print'
'clear'
return

Crazy_Merlin

繁星一点 发表于 2012-8-31 10:53

                               
登录/注册后可看大图

ctl文件：
dset G:\AWX\20090621_30\FY2E_PRE_001_OTG_20120824_0600.AWX
title 2012 twpice fy2e_pre
...

我试了，但是结果和你的差不多~郁闷了~会不会是因为undef 的问题！不知道降水估计里面无效值是怎么定义的？你可以找高手问问，我也是初学的！

繁星一点

Crazy_Merlin 发表于 2012-8-31 11:32

                               
登录/注册后可看大图

我试了，但是结果和你的差不多~郁闷了~会不会是因为undef 的问题！不知道降水估计里面无效值是怎么定义的 ...

好的  我再琢磨琢磨到底怎么回事
太谢谢你啦

Samadhi

有包含钓鱼岛的底图吗。。

wazky

繁星一点 发表于 2012-8-30 16:34

                               
登录/注册后可看大图

这个动力论坛grads版有下的
我把我的传给你吧

谢谢学姐还想问学姐个问题，您知道怎么在2.5*2英寸（虚页面）中画南海底图吗？

knut26

不知道最后解决了没有？以来每个卫星产品都不能通用ctl啊！

txxrk

有解决的吗？

qazgsy90

我也出现了这样的问题，我想错误可能出现在这句话上。pre   1    -1,40,1    Geopotential Heights (m)