IDL等绘图命令对应关系1览表.docVIP

IDL、NCL、GRADS、MATLAB绘图命令对应关系一览表 1.gradsncl的什么文件对应？grads只支持按照一定顺序存储的二进制数据文件，后缀名以.grd或者.dat或者.grb结束。这种数据必须以时间为最外层，然后每个变量按照向量形式存储，每个变量由外向内的存储顺序是高度（或等压面）==》纬度==》经度。这必须注意，否则画图容易出现一堆一堆的乱线条，这就说明你的数据没有按照grads的要求存储。再看看ncl，可以说ncl支持绝大多数各种数据的读写，包括netcdf，hdf,以及二进制数据甚至ASCII码(如果说是十进制数据或许你会更熟悉)，前两种数据一般都有头文件，不能用C语言或者Fortran读取，都需要插件才可以读取，ncl可以直接读取，matlab中也可以直接读取NETCDF格式(.nc)的数据。grads中可以读取按照说明存取的NC数据，这种数据必须又正确的时间说明，也就是说时间必须是真实的，有些模式模拟出来的数据grads的sdfopen命令是打不开的，因为一般模式都是nonleap run，都是平年，没有设定闰年，造成了时间说明不真实，grads就会报错。那么grads如何使用NC数据呢？所以建议使用ncl转换数据，将NC数据，hdf数据或者十进制数据转换成grd数据，供grads使用。这样说明是在是太空洞了，那么下面我举个例子吧。eg1) 使用grads将netcdf数据转换成grd数据reinitvar.1=air;var.2=hgt;var.3=uwndvar.4=vwnd;var.5=omega;var.6=shumj=6while (j=6)set fwrite/disk3/users/Rao_Jian/ERA-Interim-daily/entropy/var.j.daily.1979-2010.grdset gxout fwritei=1979while (i=2010)sdfopen/disk3/users/lbq/ERA-Interim-daily/pressure/var.j.interim.i.nctt=1if(i=1980|i=1984|i=1988|i=1992|i=1996|i=2000|i=2004|i=2008)while (tt=366)set t ttzz=1while (zz=37)set x 1 240set y 1 121set z zzd var.jzz=zz+1endwhilett=tt+1endwhileelsewhile (tt=365)set t ttzz=1zz=1while (zz=37)set x 1 240set y 1 121set z zzd var.jzz=zz+1endwhilett=tt+1endwhileendifi=i+1name=/disk3/users/lbq/ERA-Interim-daily/pressure/var.j.interim.i.ncclose 1endwhiledisable fwritej=j+1endwhileeg2)使用ncl将数据输出成二进制数据，grads可以使用(截取部分)patho? ?? ?????= /disk3/users/Rao_Jian/Hadley/system(rm +patho+ev_ts.grd)system(rm +patho+ev_ts2.grd)system(rm +patho+ev_ts.txt)system(rm+patho+ev_ts2.txt)do nt=0,719fbindirwrite(patho+ev_ts.grd,ev_ts(:,nt));写成二进制数据，供grads使用end dofbindirwrite(patho+ev_ts2.grd,ev_ts(time|:,evn|:))asciiwrite(patho+ev_ts.txt,ev_ts(time|:,evn|:))；写成十进制数据，可以贴到EXCEL中使用asciiwrite(patho+ev_ts2.txt,ev_ts);print(ev_ts(0,:))此外ncl中还有其他的读写函数，如fbinrecread,fbinrecwrite,fbinwrite,fbinreadncl读取netcdf3/4、hdf4、grib-1/2也是通过文件操作的eg3)ncl支持直接读取格式文件的读法fi? ?? ?? ?????= addfile(pathi+HadISST_sst.nc,r)sst0? ?? ?? ? =fi-sst(960:1679,:,:)