再看看ncl，可以说ncl支持绝大多数各种数据的读写，包括netcdf，hdf,以及二进制数据甚至ASCII码如果说是十进制数据或许你会更熟悉，前两种数据一般都有头文件，不能用C语言或者Fortran读取，都需要插件才可以读取，ncl可以直接读取，matlab中也可以直接读取NETCDF格式.nc的数据。grads中可以读取按照说明存取的NC数据，这种数据必须又正确的时间说明，也就是说时间必须是真实的，有些模式模拟出来的数据grads的sd命令是打不开的，因为一般模式都是nonleap run，都是平年，没有设定闰年，造成了时间说明不真实，grads就会报错。那么grads如何使用NC数据呢所以建议使用ncl转换数据，将NC数据，hdf数据或者十进制数据转换成grd数据，供grads使用。
这样说明是在是太空洞了，那么下面我举个例子吧。
eg1 使用grads将netcdf数据转换成grd数据 reinit var.1air;var.2hgt;var.3uwnd var.4vwnd;var.5omega;var.6shum j6 while j6 set fwrite/disk3/users/Rao_Jian/ERA-Interim-daily/entropy/var.j.daily.1979-2010.grd set gxout fwrite i1979 while i2010 sd/disk3/users/lbq/ERA-Interim-daily/pressure/var.j.interim.i.nc tt1 ifi1980|i1984|i1988|i1992|i1996|i2000|i2004|i2008 while tt366 set t tt zz1 while zz37 set x 1 240 set y 1 121 set z zz d var.j zzzz1 endwhile tttt1 endwhile else while tt365 set t tt zz1 zz1 while zzsst9601679,, ; load 50 year data duringfrom 1950 to 2009 注意这类文件的后缀名一般为.nc /.hdf/ .grb/.hdfeos/.ccm 2.grads中的描述文件与ncl中的什么对应 描述文件（.ctl是一个纯文本文件，我们有数据grads还是不能出图，需要一个描述文件来指定他的存储数据个数，维度（时间长度、层数和经纬度信息）。ncl是没有这种文件的，那么ncl在绘图时是通过什么信息来获得维度信息的呢 其实，在ncl里又一个概念叫做coordinate,其实就是维度信息。在ctl中，是通过xdef、ydef、zdef、tdef来声明数据的经纬度信息的。也就是说，在一个ctl文件中声明的所有数据，其相同维度的信息必须相同维数可以不同）。那么我们还是举个例子好了 eg4.1ctl描述文件 dset/disk3/users/Rao_Jian/ERA-DATA/entropy/Var3D.grd undef -999 title monthly zonally mean uwnd,vwnd,temperaturein THETA-PVLAT coordinate options yrev xdef1 linear 0.000 2.500 ydef200 linear 0.000 0.75 zdef16levels 270 280 290 300 315 330 350 370 400 450 500 550 600650 700 750 tdef389 linear jan1979 1mo vars 3 ulat2 16 -99 means zonal mean u wind in theta-pvcoordinate vlat2 16 -99 means zonal mean v wind in theta-pvcoordinate tlat2 16 -99 means zonal mean temperature intheta-pv coordinate endvars ; 上面的例子中，声明了三个变量，每个变量都是16层，说明都是空间三维的（可以不是的，此时16改为0），黄色的信息是维度信息，也就是说这三个变量ulat2、vlat2、tlat2相同的经度、维度、层数、时间。你要说了，我的数据哪那么巧都是一样的维度，要是不一样怎么办 对不起，必须是相同维度的信息完全一样的数据才可以放到一个描述文件中，否则要分开在不同的文件中描述，这也就是为什么在一些使用者的批处理脚本中.gs中会打开多个文件的原因之一，因为变量的信息可能不同。
eg4.2ncl中维度声明形式举例 ;下面再来看看ncl中是如何声明维度信息的 dudP00“time“ dudP01“lat“ dudP02“lon“ dudP03“level“ dudP0 isentropic levels pvlvllatulatitude pvlvllonulongitude pvlvllatunits “degrees_north“ pvlvllonunits “degrees_east“ pvlvllevunits“K“ pvlvltimeunits “monthssince 1900-01-01“ 注意 pvlvllatunits “degrees_north“是接在pvlvllat后面的，也就是说，degrees_north属性是变量pvlvl的坐标lat的属性 3.grads中的批处理命令文件与ncl中的什么对应 这个不需要解释了吧，gs文件自然就相当与ncl脚本文件了，好好理解一下吧。说的更具体一点，是相当于ncl中graphics部分，因为ncl不仅具有绘图功能，还具有强大的计算功能，尤其是巨大的函数库。ncl的缺点是循环比较慢，优点就是有许多优化函数，不需要使用者自己编写子程序（如果你愿意，也可以自己编写，将脚本存在够得着的地方，用的时候load。
4.grads中的图形类型设置与ncl的图形接口csm函数对应 说明 grads图形类型设置set gxout grahics-style ncl图形类型设置gsn_csm_choice[_choise[_choice[]]] 一般等值线图 set gxout contour gsn_csm_contour无地图底图）gsn_csm_map_ce,之前都要设置cnLinesOnTrue 等值线 robinson投影 （cylindrical equidistant map） set gxout contour set mproj robinson gsn_csm_contour_map_ce 等值线极射赤面投影 （Polar Stereographic Projections） set gxout contour set mproj nps/sps gsn_csm_contour_map_polar之前设置 gsnPolar “NH“/gsnPolar “SH“ 填色等值线图 set gxout shaded gsn_csm_contour,之前设置 gsn_define_colormapwks,colors,cnFillOnTrue 填色等值线等值线 set gxout shaded d .. set gxout contour d .. 注意顺序不能反 同上，之前