资源描述：
《grads绘图学习技巧与实例》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、Nanjinguniversityofinformationscience&technology阿木以下技巧总结都是笔者从学习实践过程中总结出来的，基本的问题。不求全面，希望对读者学习有用，如果有问题，敬请留言指正，以促进交流学习！1、软件综述：grads软件是一款绘图软件除了绘制图形，还可以提取数据，主要应用是在大气科学中，当然只要是数据处理成grads能够读取的数据文件就可以进行相关绘图。软件版本问题，软件本身不是很大，我接触到1.8、1.9、2.0版本的，1.8版本的安装很多情况还要修改环境变量、1.9版本的不识别‘sdfopen’命令，最稳定的版本是2.0版本，所以笔者推
2、荐学习者安装2.0版本，选择默认安装路径就可以。2、文件类型简述：grads处理的是网格数据，可以处理的数据类型有：grd、grib、nc（海洋常用的数据），cdf（雷达卫星数据），其中nc、cdf数据都是自带描述文件，不需要ctl，grib数据要通过命令生成ctl、index数据才可以调用，常用的是grd数据，需要ctl。3、数据文件转换：grads软件识别的数据是二进制无格式数据，文件类型是‘binary’，写入和生成时是不需要格式的如read(20)sst(i,j,iz,it),20为文件号，通常是十进制数据与grd数据间转换，这里给一个grd转换成txt数据的fortra
3、n程序：parameter(nx=56,ny=41,nz=1,nt=360)dimensionsst(nx,ny,nz,nt)realsstopen(15,file='sst.grd',form='binary')！固定的用form=‘binary’就是二进制数据open(16,file='sst.txt')！新建txt文件doit=1,ntdoiz=1,nzread(15)((sst(i,j,iz,it),i=1,nx),j=1,ny)！read后只有文件号，数据是无格式的enddoenddodoit=1,ntdoiz=1,nzwrite(16,*)((sst(i,j,iz,i
4、t),i=1,nx),j=1,ny)！输出时是txt文件可直接看的数据，有格式输出，有*enddoenddo11Nanjinguniversityofinformationscience&technology阿木close(15)close(16)end写程序时：注意格点数要与数据对应，如：上程序对应的数据是经度90~200，纬度-20~60，时间：1971.01~2000.12共360个月的海面温度数据，数据格点精度2*2，nx=(200-90)/2+1,ny=(60-(-20))/2+1,nt=360，nz=1，大气的数据要根据数据的层次确定几层。1、grd、ctl、、gs、
5、nc详述grd文件：grd数据不可直接看，为二进制无格式数据，简单的说只有1和0，而且数据间没有间隔，grads识别grd是根据ctl进行划分的，根据ctl中的经度、纬度、层次、时间，精度进行数据分块。ctl实例：dsetC:datasst.grdundef-9.99E+33titleseasurfacetempturexdef56linear902ydef41linear-202tdef360linearjan19711mozdef1levels0vars1sst099surfaceseatemptureendvars上例数据的数据顺序是（以下是数据对应的经纬度）纬度经度（
6、t=1）-2090~200的纬度20S的从90E~200的56个数据-1890~200的纬度18S的从90E~200的56个数据...6090~200的纬度60N的从90E~200的56个数据以上为一层的数据，接下来是t=2，t=3……t=360的数据，每个时间点的每一层是如上格式，编程时读取和写如的数据循环顺序依次是：时间、层次、纬度、经度，读者参照3中的fortran程序加以理解。ctl文件：具体其他指导书上都有，我这里强调的是sst后面的0表示一层，如果是两层以上则是2,3……，1层是固定用法，sst后面的99是默认设置；undef-9.99E+33此处的数值决定了软件将文
7、件中的那些值认定为不绘制的数据，所以这个值一定要与数据对应。11Nanjinguniversityofinformationscience&technology阿木gs文件，批命令文件，与ctl一样是用记事本编写，另存为.gs文件，文件都是命令，方便大段的命令编写、修改，很常用，尤其时绘制的图要求比较多时必须用，免得在命令窗口重复输入命令浪费时间，如下例子：'reinit''openc:datauv.ctl''sett7''setlev850''setgridoff''set
