求助！grads提取nc资料中遇到的问题

猫猫monty

我主要是要用grads提取1981—2010年的冬季（1，2，3月）高度场资料，用的是NCEP2的高度场。
提取的gs如下（以1981年为例）：
'reinit'
yy=1981
'set gxout fwrite '
'set fwrite E:/data/ncep_grd/hgt/hgt'yy'.grd'
'sdfopen E:/data/NCEP2/hgt/hgt.'yy'.nc'
if((math_mod(yy,4)=0&math_mod(yy,100)!=0)|math_mod(yy,400)=0)
tt=1
while(tt<=59|t>=336&t<=366)
'set t 'tt''
zz=1
while(zz<=17)
'set z 'zz''
'set x 1 144'
'set y 1 73'
'd hgt'
zz=zz+1
endwhile
tt=tt+1
endwhile
else
tt=1
while(tt<=59|t>=335&t<=365)
'set t 'tt''
zz=1
while(zz<=17)
'set z 'zz''
'set x 1 144'
'set y 1 73'
'd hgt '
zz=zz+1
endwhile   
tt=tt+1
endwhile
endif
'disable fwrite'
;
！！！此处提取1，2，3月每年90天，闰年的2月最后一天不要
分别完成1981—2010年30个年份的提取之后，鉴于之后要用一个30年完整序列，我就想对每年的头尾两天检验一下一原始的nc数据是否一致。
遇到如下问题（仍然以1981为例）：
ctl文件如下：
dset E:\data\ncep_grd\hgt\hgt1981.grd
title daily NC winter
undef -9.99e+08
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 17 levels 1000 925 850 700 600 500 400 300
250 200 150 100 70 50 30 20 10
tdef 90 linear 00Z01JAN1981 1440mn
vars 1
hgt=>hgt  17  t,z,y,x  Daily Geopotential Heights on Pressure Levels
endvars
当set t=1时，正常，可以画出和nc资料中一样的图；

当set t=90时，d hgt时，出现如下错误：

求助啊，这是什么问题？有人碰到过么？求大神帮忙！！

311wen

啊帮你顶顶吧。等大神来解答

芥末

你写入grd的文件中，t是啥？ 你的问题是写入数据的gs就不对

river

不清楚while循环里面还可以使用条件语句的啊，可能是我懂得不多···

猫猫monty

芥末 发表于 2015-11-13 23:39
你写入grd的文件中，t是啥？ 你的问题是写入数据的gs就不对

我现在知道应该是while循环选月份日期的那一句出错了，我单独分月再合成就对了

猫猫monty

river 发表于 2015-11-13 23:45
不清楚while循环里面还可以使用条件语句的啊，可能是我懂得不多···

的确是那句的错

猫猫monty

river 发表于 2015-11-13 23:45
不清楚while循环里面还可以使用条件语句的啊，可能是我懂得不多···

的确是while循环那句有问题

风之精灵

根据你gs的意图，你程序提取的应该是冬季(12,1,2月)三个月的位势高度数据，因此，1981年冬季应该被提取的就只是当年12月和次年1月、2月的数据，不能像你这样处理提取当年1月、2月和12月的数据。虽然你能提取数据成功，但你在ctl的描述文件中tdef 90 linear 00Z01JAN1981 1440mn这句这么写，就错了

风之精灵

这样处理，画图时，set t=90，时间信息是1981年3月31日，显然你的数据中没有这个时间信息。如果非要这样1、2、12月全提取当年的数据，相应的ctl文件中tdef这句话，应该像zdef一样以枚举型表示，否则一定会报你这种错误。

猫猫monty

风之精灵 发表于 2015-11-14 10:12
这样处理，画图时，set t=90，时间信息是1981年3月31日，显然你的数据中没有这个时间信息。如果非要这样1、 ...

grads中tdef的设置中不是只能写均匀增量的么？