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

chuanhao

lqouc 发表于 2014-5-21 12:20
但是这次不管怎么说，这部分应该图是对的了吧？
在输出的时候是不是可以不设置x，y维度？那样应该默认是 ...

这是所有的资料，麻烦你帮我看看

lqouc

chuanhao 发表于 2014-5-21 13:35
这是所有的资料，麻烦你帮我看看

嗯，晚上给你看看哈~

lqouc

chuanhao 发表于 2014-5-21 13:35
这是所有的资料，麻烦你帮我看看

'reinit'
'open F:\pinci2.ctl'
'set grads off'
'set grid off'
'set x 1'
'set y 1'
'set t 1 50'
'define a=tt'
'close 1'
'open F:\yuehaiwenhecheng.ctl'
'set t 1 '
'set gxout shaded'
'set cint  0.2'
'set black -0.2 0.2'
'd tcorr(a,sst,t=1,t=50)'
'cbar.gs'
'set gxout contour'
'set cint  0.2'
*'set black -0.2 0.2'
'd tcorr(a,sst,t=1,t=50)'

'enable print F:\yuehaiwenxiangguan.gmf'
'print'
'disable print'
'reinit'
;
图我画出来了，看着没问题

chuanhao

lqouc 发表于 2014-5-21 21:15
'reinit'
'open F:\pinci2.ctl'
'set grads off'

谢谢大神，我在你写的这个gs基础上加了'set gxout fwrite'
'set fwrite F:\xin\dsdb\sst\xiangguanjianyan\xishu.grd'
'set x 1 180'
'set y 1 89'
'd tcorr(a,sst,t=1,t=50)'

'disable fwrite'
得出的图终于正确了，之前图不对是怎么回事呢？

lqouc

chuanhao 发表于 2014-5-22 09:41
谢谢大神，我在你写的这个gs基础上加了'set gxout fwrite'
'set fwrite F:\xin\dsdb\sst\xiangguanjiany ...

我觉得是你的那个pinci的ctl的问题你的那个ctl里面的空间分辨率是1，而海温数据是2。
你先打开了pinci，那么grads里面就默认你之后所有的网格分辨率都是1，也就是说不论你同时打开多少ctl，grads只会以第一个的信息作为标准。那么之后的错误就不难理解了，我在你define变量之后close了file1，于是空间分辨率被重新定义，这个时候结果就对了。
当然上面说的这些都是我个人的经验，我没去看过官网的解释，也不知道这样做对不对，你可以尝试着用你原来的gs，把pinci的ctl分辨率改成2应该也是可以的。

ps：不要叫我大神了，论坛大神很多，不过我不算，我也不是很擅长grads。

chuanhao

lqouc 发表于 2014-5-22 10:31
我觉得是你的那个pinci的ctl的问题你的那个ctl里面的空间分辨率是1，而海温数据是2。
你先打开了pinci， ...

哈哈，我又学了一招，用close，非常感谢！！！