求大神帮忙啊！grads 的 tcorr 函数出错！

chuanhao

lqouc 发表于 2014-4-28 21:58
其实没什么工作量，直接把ctl中的一个文件时间错后一个时刻，然后计算时用time设定纬度就行了。
可以自 ...

如果想把求出来的这个相关系数写入grd，怎么办呢，我写进去老出错

lqouc

chuanhao 发表于 2014-5-20 21:01
如果想把求出来的这个相关系数写入grd，怎么办呢，我写进去老出错

吧gs，ctl还有报错贴上来

chuanhao

lqouc 发表于 2014-5-20 21:28
吧gs，ctl还有报错贴上来

'reinit'

'open F:\xin\haiwen\pinci2.ctl'
'open F:\xin\dsdb\sst\67xiangguan\yuehaiwenhecheng.ctl'
'set grads off'
'set grid off'
'set x 1'
'set y 1'
'set z 1'
'set t 1 50'
'define a=t'

'set lon 0 360'
'set lat -88 88'
'set z 1'
'set t 1'
'set gxout fwrite'
'set fwrite F:\xin\dsdb\sst\xiangguanjianyan\xishu.grd'

'd tcorr(a,sst.2,t=1,t=50)'

'disable fwrite'

;

dset F:\xin\dsdb\sst\xiangguanjianyan\xishu.grd
title string
undef -9.99e+08
xdef  180 linear 0 2
ydef  89 linear  -88  2
zdef 1 linear 1 1
TDEF  1 linear   01jan1961 1yr
vars 1
s 1 0 u winds
endvars

chuanhao

lqouc 发表于 2014-5-20 21:28
吧gs，ctl还有报错贴上来

图是一个错位的，应该是ctl和grd不匹配，但是我不知道怎么改，希望大神指导，谢谢

chuanhao

chuanhao 发表于 2014-5-21 10:17
图是一个错位的，应该是ctl和grd不匹配，但是我不知道怎么改，希望大神指导，谢谢

C:\Documents and Settings\Administrator\桌面

chuanhao

I:\11.jpg

lqouc

chuanhao 发表于 2014-5-21 10:22

我觉得你新写的那个ctl维度的描述应该要和yuehaiwenhecheng.ctl里面一致。
在gs里面不要set lon/lat，改成set x/y比较好。

chuanhao

lqouc 发表于 2014-5-21 10:54
我觉得你新写的那个ctl维度的描述应该要和yuehaiwenhecheng.ctl里面一致。
在gs里面不要set lon/lat，改 ...

ctl的写法是一致的，把gs里改成set x 1 180，set y 1 89，以后我试过了，画出的图范围只有0-90度N，只有东经0-180，没有西经的部分，同时写成grd的话文件的大小也小了很多，我不明白这是怎么回事

chuanhao

这是图片。。。。。。。。。。。。。。。。

lqouc

chuanhao 发表于 2014-5-21 11:27
这是图片。。。。。。。。。。。。。。。。

但是这次不管怎么说，这部分应该图是对的了吧？
在输出的时候是不是可以不设置x，y维度？那样应该默认是整个空间的。不行就把你所有东西打个包放上来。