第2章数据处理课件.ppt

第2章 数据处理 在使用GrADS绘图前，必须使所用数据文件满足GrADS的数据格式要求。 2.1 各类数据格式的特点 一般二进制无格式直接或者顺序记录的数据格式：可以是格点数据或站点数据，可以用Fortran或者C语言读写。 GRIB码数据格式：是自定义的数据格式，压缩率高，占用空间小。如NCEP提供的数据资料。 NETCDF等通用数据格式：自定义数据格式，精确性好，便于传输。 2.2 数据文件的转换 1. 转换方法 转换文件的数据存放格式，一般是通过Power Station或Visual Fortran等软件使用Fortran或者C语言来编程转换。 2.数据存放形式 数据排放顺序为：经度、纬度、层次，然后是物理变量，最后是时次变化 。 3. 举例（格点资料） 现有ASCII码（十进制存放格式）数据资料文件u.dat、v.dat和SST.dat，其空间范围：60-150°E，0-40°N；层次：u、v为850、 200hPa；时段：1982.1-1985.12；分辨率：2.5*2.5。要求编写出将这三个文件转换成二进制无格式直接存取（Grads格式.grd）文件的Fortran程序。 Fortran程序编写如下： C 定义一个X,Y,Z方向的格点数以及总时次nt parameter(nx=37,ny=17,nz=2,nt=48) c 定义数组 dimension u(nx,ny,nz,nt), v(nx,ny,nz,nt), sst(nx,ny,nt) c 打开原始数据文件 open(1,file=u.dat) open(2,file=v.dat) open(3,file=sst.dat) c 打开目标文件，recl为记录的长度，对于直接存取文件必须指定记录长度 。 open(12,file=‘mhy.grd,form=unformatted, # access=direct, recl=nx*ny*4) c按一定格式把数据文件读入，irec为记录累加器 irec=0 do 100 it=1,nt do 50 iz=1,nz irec=irec+1 read(1,rec=irec) ((u(i,j,iz,it),i=1,nx),j=1,ny) 50 continue do 51 iz=1,nz irec=irec+1 read(2,rec=irec) ((v(i,j,iz,it),i=1,nx),j=1,ny) 51 continue irec=irec+1 read(3,rec=irec) ((sst(i,j,it),i=1,nx),j=1,ny) 100 continue c 将数据资料写入新文件中， irec重新进行记录累 irec=0 do 200 it=1,nt do 30 iz=1,nz irec=irec+1 write(12,rec=irec) ((u(i,j,iz,it),i=1,nx),j=1,ny) 30 Continue do 31 iz=1,nz irec=irec+1 write(12,rec=irec) ((v(i,j,iz,it),i=1,nx),j=1,ny) 31 continue irec=irec+1 write(12,rec=irec) ((sst(i,j,it),i=1,nx),j=1,ny) 200 continue end 如果不考虑采用无格式直接方式存取文件的话，可以使用下面的编程形式，较为简单。 2.2 数据描述文件 GrADS绘图不能直接使用“数据文件”，而是通过“数据描述文件”间接使用“数据文件”。 1. 数据描述文件的概念 数据描述文件是原始数据文件的描述文件。 2. 数据描述文件的形式 以下是一个格点数据描述文件的例子： * this is an example to demonstrate the data descriptor file DSET d:\grads\ua.dat TITLE Upper Air Data DTYPE grid OPTION byteswapped UNDEF -9.99E33 XDEF 80 LINEAR -140.0 1.0 X轴说明80格点数,-140.0起始值,1.0间隔值.数据为实型 YDEF 50 LINEA