GrADS程序设计 第2章 数据处理.ppt

* * * 气象资料观测系统介绍 了解气象资料的获取方法，认识数据资料的构成。 从学科上分，气象观测属于大气科学的一个分支。它包括地面气象观测、高空气象观测、大气遥感探测和气象卫星探测等，有时统称为大气探测。由各种手段组成的气象观测系统，能观测从地面到高层，从局地到全球的大气状态及其变化。 气象观测系统 气象观测资料 数据资料处理（4D数据集） 科学试验和研究 了解大气变化特征和海-陆-气相互作用的关系 对未来进行合理预测。 第2章 数据准备 在使用GrADS绘图前，必须确保所用数据文件满足GrADS的数据格式要求。 2.1 各类数据格式的特点 Binary（二进制）数据格式：可以是格点数据或站点数据，可以用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要求。 要求编写出将这三个文件转换成一个二进制（ binary）文件的Fortran程序。 Z X Y 海平面 850hPa 200hPa SST U,V U,V 1982.1资料 。。。。。。。。。。。。。。 。。。 1985.12资料 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 打开目标文件open(12,file=‘mhy.grd,form=‘binary) c按一定格式把数据文件读入 do 100 it=1,nt do 50 iz=1,nz read(1,*) ((u(i,j,iz,it),i=1,nx),j=1,ny) 50 continue do 51 iz=1,nz read(2,*) ((v(i,j,iz,it),i=1,nx),j=1,ny) 51 continue read(3,*) ((sst(i,j,it),i=1,nx),j=1,ny) 100 continue 物理量U 物理量V 物理量SST 时次的变化 经度变化 层次变化 纬度变化 c 将数据资料写入新文件中， do 200 it=1,nt do 30 iz=1,nz write(12) ((u(i,j,iz,it),i=1,nx),j=1,ny) 30 Continue do 31 iz=1,nz write(12) ((v(i,j,iz,it),i=1,nx),j=1,ny) 31 continue write(12) ((sst(i,j,it),i=1,nx),j=1,ny) 200 continue end 按记录顺序存取方式编写的程序如下： 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) o