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

nnlovess

---

前面发的这个帖子，第一个问题已经大部分了解了，但是提取出来的文件出了问题，希望大家能给予指教。

---

从网站上下了全球的sst数据，应该怎么截取一部分数据呢？包括经纬度范围，时间范围。下面是数据信息：是周数据，想截取lon：-60-10，lat：10-40；2000年-2009年的数据；最近才接触，看了一些教程自己写了一个，但是好像不对，请求大家指教
'reinit'
'sdfopen d:\sst.nc'
'set fwrite d:\ss.grd'
i=2000
while(i<=2009)
'set gxout fwrite'
'set lon -60 10'
'set lat 10 40'
'set z 1'
'set t 1'
i=i+1
endwhile
'disable fwrite'
;

---

通过再次学习以及大家的帮助，提取的脚本如下：
'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

昪水泽若

'set t 1'这应该是‘set t’ i
然后你中间也没有输出呀

mofangbao

好像不太对是不太负责任的提问哈，你应该说一下为啥你觉得好像不太对，不能让别人猜。
你没有看懂别人的gs，你可以试着对比一下你的和别人的，仔细对比

nnlovess

mofangbao 发表于 2013-7-8 17:36

                               
登录/注册后可看大图

好像不太对是不太负责任的提问哈，你应该说一下为啥你觉得好像不太对，不能让别人猜。
你没有看懂别人的gs ...

谢谢你的教诲啊

Rolle0524

您好，为什么我按照您的命令，输入到while(i<=1044)总是提示unknown  command？
求指点

Rolle0524

楼主您好，为什么我用你的代码提取nc的蒸发数据，大当输入到whilewhile(i<=1044)的时候，总是提示unkwown command？
求指教