grads处理季度降水值

canny_maomao

模式转出来的6,7,8三个月的数据，想做一个夏季降水平均值，但是自己编的gs画不出，求助
'reinit'
'open F:\resultdata\ATM.2009060100.ctl'
'open F:\resultdata\ATM.2009070100.ctl'
'open F:\resultdata\ATM.2009080100.ctl'
'set gxout shaded'
'set lat 41.4 43.5'
'set lon 82.7 89'
'set z 1'
'define a1=ave(tpr.1,t=1,t=120)'
'define a2=ave(tpr.2,t=1,t=124)'
'define a3=ave(tpr.3,t=1,t=124)'
'define a4=((a1+a2+a3)/3)'
'd a4'

以下是在grads运行gs之后显示的：

mofangbao

gs中的经纬度超出了数据（ctl）中的范围

canny_maomao

mofangbao 发表于 2011-12-7 19:50

                               
登录/注册后可看大图

gs中的经纬度超出了数据（ctl）中的范围

啊 不应该会超出啊 …… 是哪里出了问题 请问怎么修改？

mofangbao

canny_maomao 发表于 2011-12-7 20:00

                               
登录/注册后可看大图

啊 不应该会超出啊 …… 是哪里出了问题 请问怎么修改？

那时间维呢，仔细检查一下吧，不好说怎么改，只有你检查出来了才知道怎么改，你看grads提示的很明确了

冬日情愫+

可以试试对默认的经纬度范围进行画图，看是否能画出图，如果不能说明问题可能是出在时间维上。

猎人

你这个t怎么到了120多？

canny_maomao

猎人 发表于 2011-12-8 08:54

                               
登录/注册后可看大图

你这个t怎么到了120多？

每6小时一次，一天4次，6月是120个时次，78月是124个时次……

canny_maomao

冬日情愫+ 发表于 2011-12-8 01:02

                               
登录/注册后可看大图

可以试试对默认的经纬度范围进行画图，看是否能画出图，如果不能说明问题可能是出在时间维上。

单独画每个月的图是没有问题的，但是3个一起画就有问题了……不知道为什么

冬日情愫+

那就是应该是时间维的问题了。从你的gs文件可以看出，你默认的时间应该是第一个打开文件的t=1，而此时间对于第二、三个文件时超出范围的。不妨在定义第二个和第三个文件之前将该文件设为当前文件'set dfile 2'   'set dfile 3' 多试试吧！

canny_maomao

冬日情愫+ 发表于 2011-12-8 14:47

                               
登录/注册后可看大图

那就是应该是时间维的问题了。从你的gs文件可以看出，你默认的时间应该是第一个打开文件的t=1，而此时间对于 ...

真的哎 加上'set dfile'就OK了~~
谢谢~