【已解决】数据处理

午后阳光

请教各位大神，我的数据范围为0-90N，0-360，分辨率2.5*2.5，时间尺度81，高度100-1hPa（11层），ctl如下：
dset i:\sswnew\ozone\ssw\sswa1.dat
title  1-100hPa Ozone mass mixing ratio
undef -9.99e+33
xdef 144 linear -180 2.5
ydef 37 linear 0 2.5
zdef 11 levels 100 70 50 30 20 10 7 5 3 2 1
tdef 81 linear 01NOV1979 1dy
vars 1
o3 11 99 ozone mass mixing ratio[ppb]
endvars

再用grads提取其中60-90N，某一层（如10hPa）上的数据，其他不变。提取数据的gs如下：
'reinit'
'open  I:\sswnew\ozone\ssw\sswa1.ctl'
'set lon 0 360'
'set lat 60 90'
'set z 6'
'set gxout fwrite'
'set fwrite I:\sswnew\ozone\ssw\eofq\jp\10hpa\sswjp.dat'
i=1
while(i<=81)
    'd o3'
i=i+1
endwhile
'disable fwrite'

提取后数据维数应该是(144,13,81)；但在MATLAB中做eof分析时，提示数组维数错误，
我检查后发现MATLAB显示总的数组个数为152685，而144*13*81=151632
我不知道哪里出了问题，会出现1053个数据的差误，请教位大神指点，不甚感激！

lqouc

set lon 的时候多算了一次0度。改成set lon 0 357.5

午后阳光

问题解决了，把提取数据的gs中经度的提取改为如下：
*'set lon 0 360'
'set x 1 144'
就可以了，我在grads中试过，当set x 145仍能绘图，且与set t 1图形重合，猜测应该是经度的问题，虽然暂时还不能解释其原理，但这样改过后数组维数就正确了

午后阳光

lqouc 发表于 2015-7-27 21:36
set lon 的时候多算了一次0度。改成set lon 0 357.5

嗯，很谢谢你，已经解决了，你的方法也是对的。

lqouc

午后阳光 发表于 2015-7-27 21:40
嗯，很谢谢你，已经解决了，你的方法也是对的。

问题是都解决了，不过你的处理方法会造成结果与你原来的设计有所不同，你原来的设计是经度0～357.5排列，改用set x的话就是-180～177.5的排列，这个后续绘图和计算需要注意了。

wsjqxjy

{:eb502:}

午后阳光

lqouc 发表于 2015-7-27 22:14
问题是都解决了，不过你的处理方法会造成结果与你原来的设计有所不同，你原来的设计是经度0～357.5排列， ...

嗯，没注意到这个细节，太谢谢了，又学习了！

yinan

{:5_275:}{:5_275:}