grads提取nc数据 特定范围内多层次多时次的数据

ploung

想要截取25~40N，75~105E范围内8个垂直层35年1月份的数据，，，
自己编的了下面的gs文件，但是while好像不能嵌套循环，
试了很多种方法，可以截取特定的500或400hPa等，就是不能一下全部写完
垂直层好像不能循环额？都只能写好第一层，第二层的数据就出问题
是什么问题啊？应该怎么办呢？

'reinit'
'sdfopen F:\temperature.nc'
'set fwrite F:\t.tp.grd'
'set gxout fwrite'
j=1
while(j<=420)
'set t 'j''
'set x 76 106'
'set y 26 41'
i=1
while(i<=10)
'set z 'i''
'd t'
i=i+1
endwhile
j=j+12
endwhile
'disable fwrite'
;

一个不小心

首先吧我认为你的经纬度是不是有问题？直接写25~40N，75~105E就行了吧
然后吧我觉得应该把高度层的循环放到外面去，这样数据不会混乱

river

gs没有太大问题，你说是8个高度层，设置的怎么是10层？请确定原始数据里的变量名是不是t？还有你提取完画图有问题，不一定代表你提取的数据有问题，还有可能是你写的ctl 有问题

Figo

试试用CDO做

ploung

river 发表于 2014-11-12 08:31
gs没有太大问题，你说是8个高度层，设置的怎么是10层？请确定原始数据里的变量名是不是t？还有你提取完画图 ...

非常感谢，一直掉到gs有问题的坑里面了，结果是ctl的问题

river

ploung 发表于 2014-11-12 14:24
非常感谢，一直掉到gs有问题的坑里面了，结果是ctl的问题

所以有了问题，一定要把所有过程中可能出现的问题都想一想，不能钻牛角尖

ploung

想算一下每层的区域平均又出问题了，，，
所有层的图全部是一样的，，又麻烦大家帮帮忙，，有没有人知道是怎么回事啊
我的gs
'reinit'
'sdfopen F:\temperature.nc'
'set fwrite F:\t.tp.mo.mean.grd'
'set gxout fwrite'
i=1
while(i<=10)
'set z 'i''
'set t 1 420'
'd aave(t,lon=75,lon=105,lat=25,lat=40)'
i=i+1
endwhile
'disable fwrite'
;
ctl文件
dset F:\t.tp.mo.mean.grd
title reanalysis data
undef -9.99E+33
xdef 1 linear 1 1
ydef 1 linear 1 1
zdef 10 levels 600 500 400 300 250 200 150 100 70 50
tdef 420 linear Jan1979 1mo
vars 1
t 0 99 temperature  
endvars

ploung

Figo 发表于 2014-11-12 09:31
试试用CDO做

cdo？？grads可以用吗？