求助：四维nc格式数据在grads中读取并set frwrite写入二进制格式数据中报错

提拉米苏lei

想把经度、维度、高度的三维nc格式数据在grads中读取并set frwrite写入二进制格式数据中，gs文件如下：
'sdfopen geo.nc'
'sdfopen geopotential.nc'
'set fwrite geo.dat'
'set dfile 1'
'set x 1 81'
'set y 1 81'
'set t 2'
aa=1
while(aa<32)
'set z aa'
'define org1=z.1'
'set gxout fwrite'
'd org1'
aa=aa+1
endwhile
'disable fwrite'
;
'set fwrite geo-s.dat'
'set dfile 2'
'set x 1 81'
'set y 1 81'
'set t 1'
'set lev 0'
'define org2=z.2'
'set gxout fwrite'
'd org2'
'disable fwrite'
;


但是在高度层数据写入时，就是第一部分循环那里总是报错。
错误如下：
ga-> org.gs
SET error:  Missing or invalid arguments for Z option
SET error:  Missing or invalid arguments for Z option
SET error:  Missing or invalid arguments for Z option
SET error:  Missing or invalid arguments for Z option
SET error:  Missing or invalid arguments for Z option
SET error:  Missing or invalid arguments for Z option


请高手指点，我应该怎么把32层（从1000hPa到200hPa）数据写进去？

提拉米苏lei

'sdfopen geo.nc'
'sdfopen geopotential.nc'
'set fwrite geo.dat'
************************************
'set dfile 1'
'set x 1 81'
'set y 1 81'
'set t 2'
lev.1=1000;lev.2=975;lev.3=950;lev.4=925;lev.5=900;lev.6=875;lev.7=850;lev.8=825;
lev.9=800;lev.10=775;lev.11=750;lev.12=700;lev.13=650;lev.14=600;lev.15=550;lev.16=500;lev.17=450;lev.18=400;
lev.19=350;lev.20=300;lev.21=250;lev.22=225;lev.23=200
'set gxout fwrite'
aa=1
while(aa<=32)
lev=lev.aa
'set lev 'lev
'define org1=z.1'
'd org1'
aa=aa+1
endwhile
'disable fwrite'
*****************************
'set fwrite geo-s.dat'
'set gxout fwrite'
'set dfile 2'
'set x 1 81'
'set y 1 81'
'set t 1'
'set lev 0'
'define org2=z.2'
'd org2'
'disable fwrite'
;

改成这样的层数，依然报错
ga-> org.gs
SET error:  Missing or invalid arguments for LEV option
SET error:  Missing or invalid arguments for LEV option
SET error:  Missing or invalid arguments for LEV option
SET error:  Missing or invalid arguments for LEV option
SET error:  Missing or invalid arguments for LEV option
SET error:  Missing or invalid arguments for LEV option
SET error:  Missing or invalid arguments for LEV option
SET error:  Missing or invalid arguments for LEV option
SET error:  Missing or invalid arguments for LEV option

怎么办

提拉米苏lei

一点微小的错误在gs文件里怎么都看不出来，刚发完贴再看一遍就看到了，是这里把23携程32层了！！！

aa=1
while(aa<=23)

提拉米苏lei

写数据的问题解决了！！
一点微小的错误在gs文件里怎么都看不出来，刚发完贴再看一遍就看到了，是这里把23携程32层了！！！

aa=1
while(aa<=23)

我鄙视我自己

Gamile

这个错误是什么意思呀

ymy_00000

这是我的gs程序：
'renint'

'open H:/ymy/ECMWF-precipitation/06/ECMWF/grib/15_cf_EC_06.ctl'
'set gxout fwrite'
'set fwrite H:/ymy/ECMWF-precipitation/06/ECMWF/grib/15_cf_EC_06.grib'

'set lon 70 135'
'set lat 0 55'
'set z 1'
tt=1
while (tt<=147)
    'set t 'tt''
     'd tpratesfc'
      tt=tt+1
endwhile
  'disable fwrite'
  ;
  
出现的问题：set error:missing or invalid arguments for T option

请问楼主知道问题出在哪吗？
我是想要提取出一定经纬度范围的数据，存入新的文件中

提拉米苏lei

ymy_00000 发表于 2018-10-6 16:36
这是我的gs程序：
'renint'

'set t 'tt''  改成  'set t 'tt   试下

提拉米苏lei

ymy_00000 发表于 2018-10-6 16:36
这是我的gs程序：
'renint'

'set t 'tt''  改成  'set t 'tt   试下

提拉米苏lei

ymy_00000 发表于 2018-10-6 16:36
这是我的gs程序：
'renint'

'set t 'tt''  改成  'set t 'tt   试下