grads批量输出，每年的值为什么都一样呢？求大侠们指导

isoasis · 发表于 2018-7-7 12:12:49

'reinit'
'set gxout fwrite'
'set fwrite E:\data\ncep_data\daily\xun\121\hgt_121_1953-2018.grd'
yy=1953
while(yy<=2018)
'sdfopen E:\data\ncep_data\daily\hgt.'yy'.nc'
'set x 1 144'
'set y 1 73'
'set z 6'
'define m121'yy'=ave(hgt,t=335,t=344)'
'd m121'yy''

yy=yy+1
endwhile
'disable fwrite'

river · 发表于 2018-7-7 21:05:32

你自己写的描述文件呢？
从你的提问来看你想的是每年有一组数据，但是你的合并脚本合并出来的意思是一个资料里有那么多个变量

isoasis · 发表于 2018-7-7 22:35:29

river 发表于 2018-7-7 21:05
你自己写的描述文件呢？
从你的提问来看你想的是每年有一组数据，但是你的合并脚本合并出来的意思是一个资 ...

谢谢回复

只有1个变量'set z 6'
DSET E:\data\ncep_data\daily\xun\data_xun\hgt_121-1953-2017.grd
TITLE qu
UNDEF -9.962921e+36.
XDEF 144 LINEAR 0 2.5
YDEF 73 LINEAR  -90 2.5
ZDEF 1 LEVELS  500
TDEF 65 LINEAR 00Z01DEC1953  1yr
VARS 1
hgt 0 99 qu shum uwnd component
ENDVARS

isoasis · 发表于 2018-7-7 22:39:53

isoasis 发表于 2018-7-7 22:35
谢谢回复只有1个变量'set z 6'
DSET E:\data\ncep_data\daily\xun\data_xun\hgt_121-1953-2017 ...

最后在循环内加了'close 1'
可以出结果了，但提示
data request waring:request is completely outside file limits

river · 发表于 2018-7-8 08:34:33

isoasis 发表于 2018-7-7 22:35
谢谢回复只有1个变量'set z 6'
DSET E:\data\ncep_data\daily\xun\data_xun\hgt_121-1953-2017 ...

只有一个变量？那你define的变量后面加了'yy' 干什么？你加了'yy' ，循环一次就变一次，怎么可能是同一个变量。你用set z 6 来证明只有一个变量是什么逻辑？层次和变量有什么关系？

按照你的描述文件，定义变量后面的'yy' 是完全没有必要的

qxsdjt · 发表于 2018-7-10 18:23:29

定义yy没有影响，可以不定义。

qxsdjt · 发表于 2018-7-10 18:27:29

2018年才过了一半，你dayly的数据设置t335肯定是超过数据的limit了。

isoasis · 发表于 2018-7-11 10:26:05

qxsdjt 发表于 2018-7-10 18:27
2018年才过了一半，你dayly的数据设置t335肯定是超过数据的limit了。

多谢大侠热心回复

		自动登录	找回密码
密码			立即注册

[图形美化] grads批量输出，每年的值为什么都一样呢？求大侠们指导