GrADS绘图学习技巧与实例.docx

实用标准文案 以下技巧总结都是笔者从学习实践过程中总结出来的， 基本的问 题。不求全面，希望对读者学习有用，如果有问题，敬请留言指正， 以促进交流学习！ 1、软件综述 ： grads 软件是一款绘图软件除了绘制图形，还可以提取数据，主 要应用是在大气科学中，当然只要是数据处理成 grads 能够读取的数据文件就可以进行相关绘图。 软件版本问题，软件本身不是很大， 我接触到 1.8、1.9、 2.0 版本的， 1.8 版本的安装很多情况还要修改环境变量、 1.9 版本的不识别 ‘ sdfopen’命令，最稳定的版本是 2.0 版本，所以笔者推荐学习者安装 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 数据的 fortran 程序： parameter(nx=56,ny=41,nz=1,nt=360) dimension sst(nx,ny,nz,nt) real sst open(15,file=sst.grd,form=binary) ！固定的用 form=‘ binary ’就是二进制数据 open(16,file=sst.txt) ！新建 txt 文件 do it=1,nt do iz=1,nz read(15) ((sst(i,j,iz,it),i=1,nx),j=1,ny) ！ read 后只有文件号，数据是无格式的 enddo enddo do it=1,nt do iz=1,nz write(16,*) ((sst(i,j,iz,it),i=1,nx),j=1,ny) ！输出时是 txt 文件可直接看的数据， 有格式输出，有 * enddo enddo 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，大气的数据要根据数据的 次确定几 。 4、 grd 、ctl、、gs、nc 详述 grd 文件：grd 数据不可直接看， 二 制无格式数据， 的 只有 1 和 0，而且数据 没有 隔， grads grd 是根据 ctl 行划分的，根据 ctl 中的 度、 度、 次、 ，精度 行数据分 。 ctl 例： dset C:\data\sst.grd undef -9.99E+33 title sea surface tempture xdef 56 linear 90 2 ydef 41 linear -20 2 tdef 360 linear jan1971 1mo zdef 1 levels 0 vars 1 sst 0 99 surface sea tempture endvars 上例数据的数据 序是（以下是数据 的 度） 度 度（ t=1） -20 90~200 的 度 20S 的从 90E~200 的 56 个数据 -18 90~200 的 度 18S 的从 90E~200 的 56 个数据 . . . 90~200 的 度 60N 的从 90E~200 的 56 个数据 以上 一 的数据，接下来是 t=2， t=3?? t=360 的数据，每个 点的每一 是如上格式， 程 取和写如的数据循 序依次是： 、 次、 度、 度， 者参照 中的 fortran 程序加以理解。 ctl 文件： 具体其他指 上都有，我 里 的是 sst 后面的 0 表示一 ，如果是两 以上 是 2,3??， 1 是固定用法， sst 后面的 99 是默 置； undef -9.99E+33 此 的数 决定了 件将文件中的那些 定 不 制 的数据，所以 个 一定要与数据 。 gs 文件，批命令文件，与 ctl 一 是用 事本 写，另存 .gs文件，文件都是命令，方便大段的命令 写、修改，很常用，尤其