grads画站点图

弱水三千

我的资料是江淮地区92个站资料，为什么画图时提示插值错误，说less than two stations,grid will be all missing values....
不理解，求解释，将万分感谢！！！我的ctl见附件

弱水三千

dset    f:\1\grid.grd
undef  -999.0
title     sample
xdef  13 linear 110.0 1
ydef  7 linear 28.0 1
zdef  1  linear 1000 1
tdef  1 linear 01jan1956 1mo
vars 1
g 0 99
endvars

dset f:\1\corsummer3new.grd
dtype station
stnmap f:\1\cor.map
undef -999.0
title correlation
tdef 1 linear JUN1951 1mo
vars 1
c 0 99 cor data
endvars

gs文件如下：
'reinit'
'open f:\1\grid.ctl'
'open f:\1\corsummer3.ctl'

'set map 1 1 6';
'set mpdset hires cnworld'
'set grid off'
'set lon 110 122'
'set lat 28 34'
'set z 1'
'set t 1'

'define a=oacres(g,c.2)'
'define a1=maskout(a,g-0.5)'
'define aa=smth9(a1)'
'set gxout shaded'
'set grads off'
'd aa'
'set gxout contour'
*'cnbasemap aa'
'set grads off'
'd aa'
'cbarn -0.5 -0.3 0 0.3 0.5'
'printim f:\1\1.png x1000 y800 white'
*pull dummp
'c'
;

river

看看是不是生成grd的fortran程序经纬度写反了，应该是纬度在前经度在后

弱水三千

river 发表于 2013-5-7 10:19

                               
登录/注册后可看大图

看看是不是生成grd的fortran程序经纬度写反了，应该是纬度在前经度在后

我看了一下，经纬度顺序没有写反啊，实在是不知道哪出问题了啊。。。。

river

弱水三千 发表于 2013-5-7 12:46

                               
登录/注册后可看大图

我看了一下，经纬度顺序没有写反啊，实在是不知道哪出问题了啊。。。。

两个时间设置不一致的缘故吧，试试把时间设置改成一样的
tdef  1 linear 01jan1956 1mo
tdef 1 linear JUN1951 1mo

弱水三千

river 发表于 2013-5-7 14:58

                               
登录/注册后可看大图

两个时间设置不一致的缘故吧，试试吧时间设置改成一样的
tdef  1 linear 01jan1956 1mo
tdef 1 linear  ...

时间改成一致后就可以出图了，我一直以为一个时间层设置就可以随便设置的呢，现在明白啦，学习了，谢谢啦！！！

river

弱水三千 发表于 2013-5-7 17:29

                               
登录/注册后可看大图

时间改成一致后就可以出图了，我一直以为一个时间层设置就可以随便设置的呢，现在明白啦，学习了，谢谢啦 ...

你这个帖子霉运好重啊，害我损失4金钱呢

弱水三千

river 发表于 2013-5-7 17:48

                               
登录/注册后可看大图

你这个帖子霉运好重啊，害我损失4金钱呢

啊？不清楚怎么回事啊，真是不好意思啊。。。

river

弱水三千 发表于 2013-5-7 20:42

                               
登录/注册后可看大图

啊？不清楚怎么回事啊，真是不好意思啊。。。

这是论坛的一个随机机制，挺好玩儿的。我就是说说，不怪你哈，是我人品差

弱水三千

river 发表于 2013-5-7 20:47

                               
登录/注册后可看大图

这是论坛的一个随机机制，挺好玩儿的。我就是说说，不怪你哈，是我人品差

好吧，出门得注意了。。。