标准差场，相关系数及相关系数场的程序示例

姗宝wsw

楼主，你好。我想问一下我对全年sst场做平均（1950-2000），然后用楼主画nino指数和冷舌指数的方法作图。我想画的是陆温和海温的比较图。
但是图出不来，报错是

想问一下原因是什么？
ctl分别是
海温：
dset E:\data\surface\sst1.grd
undef -9.99E+08
title NCEP/NCAR REANALYSIS 1950-2000
xdef 1 linear 0 2
ydef 1 linear -88 2
zdef 1 linear  0 1
tdef 51 linear 00z01jan1950 12mo
vars 1
sst 0 99 1000
endvars
陆温：
dset E:\data\surface\air1.grd
undef -9.99E+08
title NCEP/NCAR REANALYSIS 1950-2000
xdef 1 linear 0 2.5
ydef 1 linear -90 2.5
zdef 1 linear  0 1
tdef 51 linear 00z01jan1950 12mo
vars 1
air 0 99 1000
endvars

善人/jw

姗宝wsw 发表于 2013-6-21 11:19

                               
登录/注册后可看大图

楼主，你好。我想问一下我对全年sst场做平均（1950-2000），然后用楼主画nino指数和冷舌指数的方法作图。我 ...

是不是两者的分辨率不同的问题

mofangbao

姗宝wsw 发表于 2013-6-21 11:19

                               
登录/注册后可看大图

楼主，你好。我想问一下我对全年sst场做平均（1950-2000），然后用楼主画nino指数和冷舌指数的方法作图。我 ...

gs贴出来

兰溪之水

姗宝wsw 发表于 2013-6-21 11:19

                               
登录/注册后可看大图

楼主，你好。我想问一下我对全年sst场做平均（1950-2000），然后用楼主画nino指数和冷舌指数的方法作图。我 ...

xdef 1 linear 0 2
ydef 1 linear -88 2

xdef 1 linear 0 2.5
ydef 1 linear -90 2.5
应该是y的问题
把-90改成-88就OK了

姗宝wsw

mofangbao 发表于 2013-6-21 12:42

                               
登录/注册后可看大图

gs贴出来

GS是：
'reinit'
'open E:\data\test\sst.ctl'
'open E:\data\test\air.ctl'
'set grads off'
'set grid off'
'set x 1'
'set y 1'
'set z 1'
'set t 1 51'
'define sst=tloop(aave(sst,lon=72,lon=146,lat=4,lat=54))'
'define air=tloop(aave(air,lon=72,lon=146,lat=4,lat=54))'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set ccolor 1'
'set cstyle 1'
'set cmark 2'
'set cthick 8'
'd sst'
'set ccolor 1'
'set cstyle 3'
'set cmark 1'
'set cthick 8'
'd air'
'set strsiz 0.15'
'set string 1 c 8'
*'draw string 2.54 6.87 sst temp'
*'draw string 8.69 1.34 air temp'
'set strsiz 0.2'
'printim E:\data\surface\sstair.png x1000 y800 white'

姗宝wsw

兰溪之水 发表于 2013-6-21 12:46

                               
登录/注册后可看大图

xdef 1 linear 0 2
ydef 1 linear -88 2

谢谢兰溪。。我改过了，可是还是出现同样的错误

兰溪之水

姗宝wsw 发表于 2013-6-21 14:46

                               
登录/注册后可看大图

GS是：
'reinit'
'open E:\data\test\sst.ctl'

你的经纬度范围只有1个，你哪来的：
'define sst=tloop(aave(sst,lon=72,lon=146,lat=4,lat=54))'
'define air=tloop(aave(air,lon=72,lon=146,lat=4,lat=54))'

姗宝wsw

兰溪之水 发表于 2013-6-21 17:16

                               
登录/注册后可看大图

你的经纬度范围只有1个，你哪来的：
'define sst=tloop(aave(sst,lon=72,lon=146,lat=4,lat=54))'
'def ...

哦哦。。清楚了。。我先做的时间平均然后经纬都变成了1，但是再做区域平均就不对了，那我应该先做区域平均再做时间平均对吧？

兰溪之水

姗宝wsw 发表于 2013-6-21 17:53

                               
登录/注册后可看大图

哦哦。。清楚了。。我先做的时间平均然后经纬都变成了1，但是再做区域平均就不对了，那我应该先做区域平均 ...

先做时间平均经纬度怎么会变为1啊，建议你先屡清思路，然后再做吧~

姗宝wsw

善人/jw 发表于 2013-6-21 12:41

                               
登录/注册后可看大图

是不是两者的分辨率不同的问题

谢谢你。。好像不是这个问题。。