关于grads中aave函数的使用

jeffggg

这个是我的GS文件和图片 图片是黑色的额我不知道问题出在了哪里 大家帮忙看看吧

river

首先楼主没搞明白aave函数的意思，求完一个时刻该区域的区域平均就只剩下一个值了，楼主一共求了27个时刻的区域平均，那就有27个值，再加上楼主固定了z=1，那么就只能画出一条折线，怎么可能出填色图呢。如果要出填色图，那就必须保证高度和时间必须都是变化的才行。

   再者楼主之所以连折线图都没出，是因为aave函数使用完之后lon和lat就相当于固定的了，所以要在d之前加上set lon 1和set lat 1之类的随便一个固定值都可以。但按楼主这个gs直接加上这两句还是会提示错误，说Z T或者E是不能变化的，这就又会导致aave函数在计算过程中出错而出不了图。这是因为楼主直接使用了d，导致了维数已经固定才开始aave的计算。
   这时候就需要define了，所以楼主先定义一个define aa=aave(air,lon=110,lon=120,lat=23,lat=30)，然后加上set lon 1和set lat 1，最后再d aa就可以出图了。这个时候set lon 1和set lat 1加在define之前或者之后都是可以的，但一定要在d之前。

  所以出来折线图的话，最终的gs应该：
'reinit'
'open f:\CMAP\precippentad.nc\CMAPALONE.crl'
'set xyrev on'
*'set gxout shaded'
'set t 1 27'

'set z 1'
'set lon 1'
'set lat 1'
'define aa=aave(air,lon=110,lon=120,lat=23,lat=30)'
'set xlabs 01|03|06|09|12|15|18|21|24|27'
'set ylabs 1|2|3|4|5|6|7|8'
'd aa'
'draw title Pantadly CMAP over (110-120E,23-30N)'
*'cbarn.gs'
'printim f:\CMAP\precippentad.nc\PantadlyCMAP.gif white'
;

lqouc

尽头的尽头 发表于 2013-3-30 16:43

                               
登录/注册后可看大图

你时次设定了那么多，要么设定一个时次，要么先求时间平均才可以吧，楼主设定了试试

时间应该是没问题的吧，要是不是时间序列的话，grads只会输出一个值，就谈不上画图了。

尽头的尽头

你时次设定了那么多，要么设定一个时次，要么先求时间平均才可以吧，楼主设定了试试

lqouc

不一定能解决你的问题吧，但是set lon和lat那里不对。

尽头的尽头

楼主可以把ctl文件，还有出错的详细信息给出来，这样才能更好的发现问题

jeffggg

尽头的尽头 发表于 2013-3-30 16:48

                               
登录/注册后可看大图

楼主可以把ctl文件，还有出错的详细信息给出来，这样才能更好的发现问题

等下 我试试  ctl文件信息是这样的 1979-2008 三十年 30个第一侯的平均值到30个第27侯的平均值 就说有30个值  每个值是30年某一侯的平均值

尽头的尽头

jeffggg 发表于 2013-3-30 16:51

                               
登录/注册后可看大图

等下 我试试  ctl文件信息是这样的 1979-2008 三十年 30个第一侯的平均值到30个第27侯的平均值 就说有30个 ...

还有就是Y坐标是区域平均的值么，大小是多少，是否匹配

Aires

首先aave是用来求空间场平均的，也就是如果你时间和层固定的话，那么你lon,lat也要固定。这样求出来的才是一个值啊。另外，如果你想得到时间序列，那么你的t可以设置为变化的，但是d的时候要这样写才行：d tloop(aave(    ))，楼主试试。

jeffggg

river 发表于 2013-3-30 17:58

                               
登录/注册后可看大图

首先楼主没搞明白aave函数的意思，求完一个时刻该区域的区域平均就只剩下一个值了，楼主一共求了27个时刻的 ...

river大神！我给你跪了！收我为徒吧！！！在你的指导下我画出来了 和预期效果一样！你上面说的有2个地方我不明白希望师傅给解释下,1,'"aave函数使用完之后lon和lat就相当于固定的了"为什么？   2、“但按楼主这个gs直接加上这两........这时候就需要define了” 为什么需要define才行 不define直接用d为啥不行呢？！