nc高度场提取500hPa数据，生成grd文件，编写ctl后，画图出错，而且没有588线出现

字俣丞

大神帮帮忙，nc高度场提取500hPa数据，生成grd文件，编写ctl后，画图出错，而且没有588线出现


GRADS转grd的脚本


'reinit'
'sdfopen e:\as\NCEP\h.nc'
'set gxout fwrite'
'set fwrite e:\as\NCEP\h2.grd'
'set x 1 144'
'set y 1 73'
'set z 6'
t0=30
while (t0<=378)
'set t 't0''
'd hgt/10'
t0=t0+12
endwhile      
'disable fwrite'



ctl文件

  DSET  E:\as\NCEP\h2.grd
  title Geopotential
  undef -32767
  TITLE Geopotential data sample
  xdef 37 linear 90 2.5                                                                                                                              
  ydef 29 linear 10 2.5
  ZDEF 1 linear 0 1
  TDEF 30 linear jun1981 12mo
  VARS 1
  hgt 0 99 Geopotential  data
  ENDVARS

字俣丞

我提取的是1981-2010年6月份的值

弱水三千

'set x 1 144'
'set y 1 73'
和ctl文件不一致，
xdef 37 linear 90 2.5                                                                                                                              
ydef 29 linear 10 2.5

字俣丞

弱水三千 发表于 2018-1-24 11:37
'set x 1 144'
'set y 1 73'
和ctl文件不一致，

非常感谢，已经出来了

字俣丞

弱水三千 发表于 2018-1-24 11:37
'set x 1 144'
'set y 1 73'
和ctl文件不一致，

非常感谢，已经出来了

字俣丞

弱水三千 发表于 2018-1-24 11:37
'set x 1 144'
'set y 1 73'
和ctl文件不一致，

非常感谢，已经正常了