【已解决】对于nc文件的处理

nnlovess

前面发的这个帖子，第一个问题已经大部分了解了，但是提取出来的文件出了问题，希望大家能给予指教。'reinit'
'sdfopen d:\sst.nc'
'set fwrite d:\ss.grd'
'set gxout fwrite'
i=523
while(i<=1044)
'set t 'i''
'set lon -60 10'
'set lat 10 40'
'set z 1'
'sstave=ave(sst,t='i+0',t='i+3')'
'd sstave'
i=i+4
endwhile
'disable fwrite'
;


ctl：
dset d:\ss.grd
undef  32800
xdef 71 linear  -60  1.0
ydef 31 linear  10  1.0
zdef  1 linear 0 1
tdef  131 linear 00z01jan2000 1mo
vars 1
sst 0 99
endvars

已经能够做出结果了，但是发现结果是不对的 ，我通过直接用nc文件
set t 523
ave（sst,t+0,t+3）
做出来第一个时间的结果是
通过提取的文件ss.grd,显示的第一个结果是

现在不知道问题出在了哪里，希望大家能给我解惑

---

---

这几日没有更新，还是river比较厉害啊，一下就看到问题了，用set x y确实解决了问题，现在把正确的结果贴上来，改动的地方是黄色高亮部分，同时也很感谢大家都帮助，同时也感受到写程序真的需要很细心！
'reinit'
'sdfopen d:\sst\sst.nc'
'set fwrite d:\sst\ss1.grd'
'set gxout fwrite'
'set x -60 10'
'set y 100 130'
'set z 1'
i=523
while(i<=1044)
'set t 'i''
'sstave=ave(sst,t='i',t='i+3')'
'd sstave'
i=i+4
endwhile
'disable fwrite'
;


ctl：
dset d:\ss.grd
undef  -9.99e33
xdef 71 linear  -60  1.0
ydef 31 linear  10  1.0
zdef  1 linear 0 1
tdef  131 linear 00z02jan2000 1mo
vars 1
sst 0 99
endvars

kongfeng0824

我感觉第一个直接处理后的应该接近结果吧。处理后的grd是

kongfeng0824

不知道你的grd的是否有问题？

river

我觉得吧，你这个图还是很像的，只要把左边那个角拿下来接到图的最右面，就和直接用nc画出来的一样了。你还是把这两句设置'set lon -60 10' 'set lat 10 40' 换成对应的set x 和set y 来试一下吧

nnlovess

kongfeng0824 发表于 2013-7-9 18:55

                               
登录/注册后可看大图

不知道你的grd的是否有问题？

你好，我怎么才能知道生成出来的grd是不是有问题呢？我检查了处理的程序也没发现哪里不对，真是纠结啊

nnlovess

river 发表于 2013-7-10 08:26

                               
登录/注册后可看大图

我觉得吧，你这个图还是很像的，只要把左边那个角拿下来接到图的最右面，就和直接用nc画出来的一样了。你还 ...

谢谢你，问题已经解决了！

river

nnlovess 发表于 2013-7-15 21:38

                               
登录/注册后可看大图

谢谢你，问题已经解决了！

顺便在标题上加上【已解决】吧