请大家帮个忙，看看我的图为什么这样？

漂洋过海

画的是30年平均的逐侯降水，资料已经处理过，没有缺测。
而且那三处不连续的地方，站点有资料，我读成站点grd也没错，难道是插值出错了？
下面图是35侯时的站点数据图，站点grd应该没没错
下面是ctl&gs:

dset H:\drought\1981-2010pre\pentad\30ave-pentad-sta578.grd
dtype station
stnmap H:\drought\1981-2010pre\pentad\30ave-pentad-sta578.map
undef 32766.0
title average
tdef  73 LINEAR 1jan1981 1dy
vars 1
r 0 99 eigenvalue
endvars

下面是gs：

'reinit'
'open H:\drought\pentad\prc5.ctl'
'open H:\drought\pentad\30ave-pentad-sta578.ctl'
'set gxout fwrite'
'set fwrite H:\drought\pentad\30ave-pentad-grid.grd'

'set x 1 141'
'set y 1 81'
it=1
while(it<=73)
'define a=oacres(prc.1,r.2(t='it'))'
'define b=maskout(a,prc.1-0.5)'
'd b'
say it
it=it+1
endwhile
'disable fwrite'
;

囧光光

看不懂  头疼

漂洋过海

lqouc 发表于 2013-4-9 14:56

                               
登录/注册后可看大图

你画图的gs呢，不能就断定是插值不对吧，很可能是你画图的问题，或者是差值的数据ctl有问题？就说这个插值的 ...

我后来知道是grads版本的问题导致插值出现一定的错误，1.9的版本插值出来没有问题，家园整合版的2.0就会有插值画的剖面图出现怪异现象！

Jennifer

还是不懂，我有必要自己从头学习一遍GRADS，可是论文不等人啊。愁得头发都要掉光了。

漂洋过海

river 发表于 2013-4-9 19:29

                               
登录/注册后可看大图

那就不清楚了。我个人感觉是应为你的数据都集中在国内，20°N左右已经到边界了，一共都没几个测站，无论是 ...

恩，你的感觉是对的，谢谢你啦，我最后取了110-111的经度范围，看了下还凑合画出了22°N这块....

漂洋过海

mofangbao 发表于 2013-4-9 19:29

                               
登录/注册后可看大图

能不能不要先弄那么多set什么的，直接 d oacres插值一下你的变量画个图看看 各种语句加上去再查错误不容易找 ...

恩，谢谢你的建议！

river

漂洋过海 发表于 2013-4-9 19:20

                               
登录/注册后可看大图

谢谢你了~！可是最下面的数据插值侯不应该是这样子，那块的降水按理说是有大值中心的~

那就不清楚了。我个人感觉是应为你的数据都集中在国内，20°N左右已经到边界了，一共都没几个测站，无论是内插、外插、还是中央差，边界都不会插值的很好的

mofangbao

能不能不要先弄那么多set什么的，直接 d oacres插值一下你的变量画个图看看 各种语句加上去再查错误不容易找

漂洋过海

river 发表于 2013-4-9 18:59

                               
登录/注册后可看大图

图的下面就是色标，填上颜色的都是大于2的，你说都是负的这不是赤裸裸的忽悠人呢吗。还有最下面的数据明显 ...

谢谢你了~！可是最下面的数据插值侯不应该是这样子，那块的降水按理说是有大值中心的~

river

漂洋过海 发表于 2013-4-9 18:35

                               
登录/注册后可看大图

数据是负的

图的下面就是色标，填上颜色的都是大于2的，你说都是负的这不是赤裸裸的忽悠人呢吗。还有最下面的数据明显都是不合适的，这应该是插值过程中边界条件的问题。把'set cmin 2'都改成'set cmin 0'边上的负值就画不出来了。