GrADS使用中的实用技巧之一 修改描述文件格式达到简化操作的目的

/xin村儿/

师兄，为什么我这样改写后，grads能打开ctl，但是不识别变量，不能画图，不知道错在哪里？
dset  E:/My-Data/NCEP-II/hgt.mon.mean.nc
title Mounthly NCEP/DOE Reanalysis 2
dtype netcdf
undef 32766
unpack scale_factor add_offset 1 32066
options yrev
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 35 linear Jan1979 1yr
vars 12
hgt1 17 99 Geopotential height(gpm)
hgt2 17 99 Geopotential height(gpm)
hgt3 17 99 Geopotential height(gpm)
hgt4 17 99 Geopotential height(gpm)
hgt5 17 99 Geopotential height(gpm)
hgt6 17 99 Geopotential height(gpm)
hgt7 17 99 Geopotential height(gpm)
hgt8 17 99 Geopotential height(gpm)
hgt9 17 99 Geopotential height(gpm)
hgt10 17 99 Geopotential height(gpm)
hgt11 17 99 Geopotential height(gpm)
hgt12 17 99 Geopotential height(gpm)
endvars

四叶草

这个变化确实用的比较巧，不过还真没有尝试过，先收藏学习下。

river

/xin村儿/ 发表于 2014-10-29 13:07
师兄，为什么我这样改写后，grads能打开ctl，但是不识别变量，不能画图，不知道错在哪里？
dset  E:/My-Da ...

你是给nc资料编写ctl ，那就肯定要和别的资料的ctl 有不同之处，你看我另一个帖子没看仔细，变量的表示方法怎么能和其他的一样呢。举个例子：
hgt=>hgt 17  t,z,y,x  Geopotential height(gpm)

/xin村儿/

river 发表于 2014-10-29 14:59
你是给nc资料编写ctl ，那就肯定要和别的资料的ctl 有不同之处，你看我另一个帖子没看仔细，变量的表示方 ...

我改为了hgt01=>hgt01 17  t,z,y,x  Geopotential height(gpm)，还是出错“Error for variable 'hgt01'
师兄生成 【dset  hgt.grd】用的是fortran还是gs，能否发给我参考一下（798839427@qq.com)

river

/xin村儿/ 发表于 2014-10-29 20:04
我改为了hgt01=>hgt01 17  t,z,y,x  Geopotential height(gpm)，还是出错“Error for variable 'hgt01' ...

对了，NC资料是一种特定的带有自描述文件的格式，不能随便自己改造。你转换成普通的二进制文件吧。利用grads命令set fwrite就可以。

道狼52203

{:eb502:}{:eb502:}

modly25

好东西必须顶一发！