转换nc 192 x 94 Gaussian Longitude/Latitude 数据出问题

nuistzhou

rt，转换nc 192 x 94 Gaussian Longitude/Latitude    1.875° x ~1.904° from 0E to 358.125E and 88.542N to 88.542S    数据出问题    。注：数据是1948.1---至今的。我需要的是1957.12-----2011.2月.于是仅仅提取了这个时间区间内的东亚数据
于是 有如下gs：
'reinit'
'sdfopen e:\bs\shtfl.sfc.mon.mean.nc'
'set fwrite e:\bs\st.grd'
'set gxout fwrite'
i=119
while(i<759)
'set t 'i''
'set lat 69.375  140.625'
'set lon 0.952368  59.9986'
'st=(shtfl(t='i')+shtfl(t='i+1')+shtfl(t='i+2'))/3'
'd st'
i=i+12
endwhile
'disable fwrite'
;




编出的grd经fortran读出检验，第一年的时候就已经出问题了，不知道为何。
求高人解答


刚刚改了点东西，'reinit'
'sdfopen e:\bs\shtfl.sfc.mon.mean.nc'
'set fwrite e:\bs\st.grd'
'set gxout fwrite'
i=120
while(i<=756)
'set t 'i''
'set lat 69.375  140.625'
'set lon 0.952368  59.9986'
'st=(shtfl(t='i')+shtfl(t='i+1')+shtfl(t='i+2'))/3'
'd st'
i=i+12
endwhile
'disable fwrite'
;

ctl：
dset e:\bs\st.grd
undef -9.96921e+36
xdef 39 linear  69.375  1.875
ydef 32 linear  0.952368  1.904
zdef  1 linear 0 1
tdef  54 linear 00dec1957 1yr
vars 1
st 0 99
endvars


出了一张图：

只有下面一部分的。。。。很郁闷
已解决！
是自己gs里面lon和lat 设置反了。

cuixindong

{:eb317:}{:eb317:}{:eb317:}{:eb317:}

cuixindong

{:eb317:}{:eb317:}{:eb317:}{:eb317:}

易小凯

lat和lon的set是不是有问题啊？

nuistzhou

易小凯 发表于 2012-3-31 12:14

                               
登录/注册后可看大图

lat和lon的set是不是有问题啊？

其实刚刚发现了。。。。解决了