GrADS读取nc写为二进制后再画图出错

半个瞎子

向各位大神请教：
    我用GrADS将nc写为二进制文件后，自己写了ctl，然后画图，画第一层没有问题，可是画其他层次就出现了问题。
   写二进制gs文件：
'reinit'
'sdfopen ./air.2013.nc'
'set gxout fwrite'
'set fwrite ./air.2013.grd'


tt=1
while ( tt <= 365 )
zz=1
while ( zz <= 3 )
  'set t 'tt''
  'set z 'zz''
  'set x 1 144'
  'set y 1 73'
  'd air'
  zz=zz+1
endwhile
tt=tt+1
endwhile
'disable fwrite'
'reinit'


       ctl描述文件：
dset ^air.2013.grd
title NCEP data
undef 32766
xdef 144 linear 0.0 2.5
ydef 73 linear -90.0 2.5
zdef 3 levels 1000 925 850
tdef 365 linear 01JAN2012 1dy
vars 1
air 1 99 Sea Level Pressure
endvars



  画除第一层外的其他层次时出现的错误信息：
ga-> open air.2013.ctl
Scanning description file:  air.2013.ctl
Data file air.2013.grd is open as file 1
LON set to 0 360
LAT set to -90 90
LEV set to 1000 1000
Time values set: 2012:1:1:0 2012:1:1:0
E set to 1 1
ga-> set z 2
LEV set to 925 925
ga-> set t 1
Time values set: 2012:1:1:0 2012:1:1:0
ga-> d air
Data Request Warning:  Request is completely outside file limits
  Entire grid contents are set to missing data
  Grid limits of file:     X = 1 144  Y = 1 73  Z = 1 3  T = 1 365  E = 1 1
  Grid limits of request:  X = 1 145  Y = 1 73  Z = 2 2  T = 1 1  E = 1 1
  Warning issued for variable = air
Cannot contour grid - all undefined values
ga->

marclevoy

ctl写的有问题，变量名后面的1 99改成3 99就可以了

marclevoy

这里的第一个数值一是用来表示层次的，如果是单层的数据，比如地表气温，海表面气压就可以用1或者0，但是多层次的数据，如高度场，风场，就要写上对应的层数，你这个数据有3层，所以要写成3 99

朽木露琪亚

air 1 99 Sea Level Pressure
改成
air 3 99 Sea Level Pressure

爱知天气

marclevoy 发表于 2014-4-29 12:39
这里的第一个数值一是用来表示层次的，如果是单层的数据，比如地表气温，海表面气压就可以用1或者0，但是多 ...

说的不错   

半个瞎子

marclevoy 发表于 2014-4-29 12:39
这里的第一个数值一是用来表示层次的，如果是单层的数据，比如地表气温，海表面气压就可以用1或者0，但是多 ...

哦哦 谢谢了，只改了zdef，忘记改下边的了，愣是找了半天错误，惭愧啊！

半个瞎子

朽木露琪亚 发表于 2014-4-29 12:43
air 1 99 Sea Level Pressure
改成
air 3 99 Sea Level Pressure

谢谢了，原来是个小问题，自己没注意。

半个瞎子

爱知天气 发表于 2014-4-29 12:49
说的不错

恩，非常感谢大家的热心帮助。

azurewings

??

marclevoy 发表于 2014-4-29 12:39
这里的第一个数值一是用来表示层次的，如果是单层的数据，比如地表气温，海表面气压就可以用1或者0，但是多 ...

weather3699

{:5_213:}