《数据处理》PPT课件 (2)

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)
6
2.数据存放形式
数据排放顺序为：经度、纬度、层次， 然后是物理变量，最后是时次变化 。
Write(10,rec=1)(((temp(I,j,k),i=1,18),j=1,9),k=1,5)
精选ppt
7
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程序。
接存取文件必须指定记录长度 。 open(12,file=‘data.grd',form='unformatted', # access='direct', recl=nx*ny*4)
精选ppt
10
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 c源自文库ntinue
$ gribmap –i test.ctl ➢ 在GrADS中处理
ga->open grib.ctl
精选ppt
5
2.2 数据文件的转换
1. 转换方法 转换文件的数据存放格式，一般
是通过Power Station或Visual Fortran 等软件使用Fortran或者C语言来编程 转换。
精选ppt
open(1,file='u.dat')
open(2,file='v.dat')
open(3,file='sst.dat')
open(12,file=‘data.grd',form=‘binary')
精选ppt
13
do 100 it=1,nt
do 50 iz=1,nz
read(1,*) ((u(i,j,iz,it),i=1,nx),j=1,ny)
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
精选ppt
12
如果不考虑采用无格式直接方式存取文件的话， 可以使用下面的编程形式，较为简单。
parameter(nx=37,ny=17,nz=2,nt=48)
dimension u(nx,ny,nz,nt),
#
v(nx,ny,nz,nt), sst(nx,ny,nt)
irec=irec+1
read(3,rec=irec) ((sst(i,j,it),i=1,nx),j=1,ny)
100 continue
精选ppt
11
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)
理 ga->xdfopen absolute.ddf
精选ppt
4
GRIB格式数据
➢ T213，GFS，JRA-25，等 ➢ 常用后缀.grb ➢ 需.ctl和.idx文件一起使用 ➢ 用grib2ctl.pl生成数据说明（描述）文件.ctl
$ perl grib2ctl.pl test.grb ➢ 用gribmap –i test.ctl 生成映射文件.idx
200 continue
end
第二章 数据处理
精选ppt
1
2.1 各类数据格式的特点
一般二进制无格式直接或者顺序记录 的数据格式：可以是格点数据或站点 数据，可以用Fortran或者C语言读写。
精选ppt
2
Binary格式数据
➢ 常用后缀.dat .grd .bin , …… ➢ 需要.ctl文件一起使用 ➢ 是无格式文件(form=‘unformatted’)
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
精选ppt
14
do 200 it=1,nt
do 30 iz=1,nz
精选ppt
3
netCDF格式数据
➢ 自定义数据格式，精确性好，便于传输 ➢ ECMWF/NCEP再分析资料等 ➢ 常用后缀.nc ➢ 自带数据内容说明（自描述） ➢ 标准NC格式，GrADS可直接处理
ga->sdfopen absolute.nc ➢ 非标准NC，写.ddf描述文件（v1.8）后可处
精选ppt
8
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)
精选ppt
9
c 打开原始数据文件 open(1,file='u.dat') open(2,file='v.dat') open(3,file='sst.dat') c 打开目标文件，recl为记录的长度，对于直