gs文件批量处理出错

'reinit'
var='air'
pathr='D:\'
pathw='D:\'
y=1979
while (y<=2014)
say y
'sdfopen 'pathr''var'.'y'.nc'
'set fwrite 'pathw''var'.'y'.grd'
'set gxout fwrite'
k=1
while (k<=17)
say k
'set lon 1 144'
'set lat 1 73'
'set z 'k''
j=1
while (j<=365)
'set t 'j''
'd 'var''
j=j+1
endwhile
k=k+1
endwhile
'disable fwrite'
'close 1'
y=y+1
endwhile
'reinit'
;
gs文件如上所示，到底出什么错误了，怎么处理后的数据17层中只有1层，不知出了一个什么错误，我一直没看出来，可能错误很低级，麻烦大家看一下，多谢各位。

四叶草

麻烦楼主把你grd文件的ctl文件也贴出来看看

wangdazheng

四叶草 发表于 2015-9-10 20:00
麻烦楼主把你grd文件的ctl文件也贴出来看看

我暂时并没有对其进行画图，只是在用fortran进行处理时，将层数变量设为17时，会显示“end of file during reading”,当将层数变量设置为1时，竟然能运行。

四叶草

建议你先作图看看，看看能不能出17层的图

wangdazheng

四叶草 发表于 2015-9-10 22:06
建议你先作图看看，看看能不能出17层的图

你好，我选了一年做了一下图，证明gs文件的数据处理是错误的。
ctl内容如下：
DSET D:\air.1979.grd
*OPTIONS Yrev sequential
options 365_day_calendar
UNDEF 32766
TITLE SAMPLE
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 365 LINEAR 01JAN1979 1dy
VARS 1
var 0 99  GIRD DATA
ENDVARS
做出来的图明显是错误的，请您看一下哪里出错了？

river

楼主先去看一下GrADS的基础知识，重点理解一下GrADS需要的数据存储方式是什么样的

str@nger

吃柠檬的猫

grads数据的存储方式：经度→纬度→层次→变量→时间，你的gs读取顺序好像错了

四叶草

wangdazheng 发表于 2015-9-10 23:44
你好，我选了一年做了一下图，证明gs文件的数据处理是错误的。
ctl内容如下：
DSET D:\air.1979.grd

楼主你好！
1、看了@river和@吃柠檬的猫 的回帖，再仔细看了楼主你的gs，gs中读的数据与nc数据在时间、高度层次错位，正确的读写顺序是经度→纬度→层次→变量→时间，而楼主你读数据的顺序是经度→纬度→时间→变量→层次；
2、楼主你选一年的数据做了图，ctl中ZDEF 17 Levels和var 0 99  GIRD DATA前后是相互矛盾的，按照楼主的ctl，你出图时，17层次的图都是一样的。

wangdazheng

问题已解决，抱歉打扰！多谢各位了！