处理数据顺序不同导致出图结果不一样

Aires

处理数据顺序不同可能会导致出图结果不一样。
大家有没有作图时发现如果处理数据的顺序不同，得到的图也不一样呢？
下面以画标准差场为例：
数据说明：
dset  d:\work\oisst\%y4%m2.dat
options template
undef  -999.9
title SST monthly for 1981.12-2003.02 derived with the revised OI version2 analysis
xdef 360 linear  0.5 1
ydef 180 linear  -89.5 1
zdef 1 linear  1 1
tdef 255 linear Dec1981 1mon
vars 1
sst     1   99   clim
endvars

作图要求：画冬季（12月~2月）SST标准差场

方案一：先提取逐年冬季（12月~2月）的SST场，共66个场，再对这66个场求标准差场。
*************提取冬季66个SST场***************
'reinit'
'open d:\work\oisst\oisstmon.ctl'
'set fwrite d:\work\winter.grd'
'set gxout fwrite'
'set x 1 360'
'set y 1 180'
'set z 1'
i=1
while(i<=255)
'set t 'i
'd sst'
j=1
while(j<=2)
'set t 'i+j
'd sst'
j=j+1
endwhile
i=i+12
endwhile
'disable fwrite'
*************画出冬季SST标准差场***************
'reinit'
'open d:\work\winter.ctl'
'set grads off'
'set grid off'
'set map 1 1 10'
'set lon 85 135'
'set lat 0 25'
'set z 1'
'set t 1 66'
'define anom=sst-ave(sst,t=1,t=66)'
'set t 1'
'set gxout shaded'
'set clevs 0.65'
'set rbcols 0 15'
'd sqrt(ave(anom*anom,t=1,t=66))'
'set gxout contour'
'set clskip 2'
'set clopts 1 0 0.15'
'd sqrt(ave(anom*anom,t=1,t=66))'
'printim d:\work\1.png white'
;

方案二：先对逐年冬季SST场做平均，即把每年12月~2月SST场做平均，得到22个场，再对这22个场求标准差场
************提取冬季（12月~2月平均）逐年SST场**************
'reinit'
'open d:\work\oisst\oisstmon.ctl'
'set fwrite d:\work\sst_winter_1981_2002.grd'
'set gxout fwrite'
'set x 1 360'
'set y 1 180'
'set z 1'
'set t 1'
i=1
while(i<255)
'd ave(sst,t='i',t='i'+2)'
i=i+12
endwhile
'disable fwrite'
****************画出冬季SST标准差场******************
'reinit'
'open d:\work\sst_winter_1981_2002.ctl'
'set grads off'
'set grid off'
'set map 1 1 10'
'set lon 85 135'
'set lat 0 25'
'set z 1'
'set t 1 22'
'define anom=sst-ave(sst,t=1,t=22)'
'set t 1'
'set gxout shaded'
'set clevs 0.65'
'set rbcols 0 15'
'd sqrt(ave(anom*anom,t=1,t=22))'
'set gxout contour'
'set clskip 2'
'set clopts 1 0 0.15'
'd sqrt(ave(anom*anom,t=1,t=22))'
'printim d:\work\work4_rms.png white'
;

大家可以讨论下这种问题~

许守华

学习学习，感谢分享

小小人

谢谢楼主，我觉得可能是在grads的读取数据的默认方式不同，也就是说，第一副图，他读取的是第66个场，冬季的数据，而第二幅图，他读得是第二十二个冬季场的标准差场。我觉得作者在grads上绘图时，是否会出现一些屏幕（闪屏），然后屏幕定格在作者出来的这两幅图？

river

第二种方案里的'd ave(sst,t='i',t='i'+2)'是有点儿问题的吧，t='i'+2这个理论上是不错的，不过放在ave函数里就有点儿问题了，他识别不了代数式，也就是后面那个+2等于没用。但是可以设置一个j='i'+2然后ave(sst,t='i',t='j')这样的

Aires

river 发表于 2013-3-30 14:38

                               
登录/注册后可看大图

第二种方案里的'd ave(sst,t='i',t='i'+2)'是有点儿问题的吧，t='i'+2这个理论上是不错的，不过放在ave函数 ...

按照你的说法我试了下，这两种方法出来的图是不一样的，那哪个是正确的呢？

Aires

小小人 发表于 2013-3-30 14:29

                               
登录/注册后可看大图

谢谢楼主，我觉得可能是在grads的读取数据的默认方式不同，也就是说，第一副图，他读取的是第66个场，冬季的 ...

不是默认方式，是我用两种不同的方法处理的数据然后画的图。

river

Aires 发表于 2013-3-30 15:21

                               
登录/注册后可看大图

按照你的说法我试了下，这两种方法出来的图是不一样的，那哪个是正确的呢？

楼主方案一里提取资料循环那一块看着好别扭，我感觉有可能是那个问题

Aires

river 发表于 2013-3-30 15:50

                               
登录/注册后可看大图

楼主方案一里提取资料循环那一块看着好别扭，我感觉有可能是那个问题

带进去看这样提出来的应该是每年冬季的SST啊

river

Aires 发表于 2013-3-30 17:03

                               
登录/注册后可看大图

带进去看这样提出来的应该是每年冬季的SST啊

您这个逻辑思维太强大了，我就有点想不明白，怎么就是每年的12月~2月呢。第一次循环i=1的时候，j=1
第二次的时候i=13，j=15？这不太对啊

Aires

river 发表于 2013-3-30 18:11

                               
登录/注册后可看大图

您这个逻辑思维太强大了，我就有点想不明白，怎么就是每年的12月~2月呢。第一次循环i=1的时候，j=1
第二 ...

第一次循环的时候i=1,set t 1 ,j=1,set t i+j  就是set t 2 再循环是set t 3,然后是i=13,j=1,再循环没错啊