第2章数据处理.pptVIP

气象资料观测系统介绍;从学科上分，气象观测属于大气科学的一个分支。它包括地面气象观测、高空气象观测、大气遥感探测和气象卫星探测等，有时统称为大气探测。由各种手段组成的气象观测系统，能观测从地面到高层，从局地到全球的大气状态及其变化。;气象观测系统气象观测资料数据资料处理（4D数据集）科学试验和研究了解大气变化特征和海-陆-气相互作用的关系对未来进行合理预测。;第2章数据准备;2.2数据文件的转换;2.数据存放形式;3.举例（格点资料）;;Fortran程序编写如下：;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按一定格式把数据文件读入

do100it=1,nt

do50iz=1,nz

read(1,*)((u(i,j,iz,it),i=1,nx),j=1,ny)

50continue

do51iz=1,nz

read(2,*)((v(i,j,iz,it),i=1,nx),j=1,ny)

51continue

read(3,*)((sst(i,j,it),i=1,nx),j=1,ny)

100continue;c将数据资料写入新文件中，

do200it=1,nt

do30iz=1,nz

write(12)((u(i,j,iz,it),i=1,nx),j=1,ny)

30Continue

do31iz=1,nz

write(12)((v(i,j,iz,it),i=1,nx),j=1,ny)

31continue

write(12)((sst(i,j,it),i=1,nx),j=1,ny)

200continue

end;按记录顺序存取方式编写的程序如下：;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

do100it=1,nt

do50iz=1,nz

irec=irec+1

read(1,rec=irec)((u(i,j,iz,it),i=1,nx),j=1,ny)

50continue

do51iz=1,nz

irec=irec+1

read(2,rec=irec)((v(i,j,iz,it),i=1,nx),j=1,ny)

51continue

irec=irec+1

read(3,rec=irec)((sst(i,j,it),i=1,nx),j=1,ny)

100continue;c将数据资料写入新文件中，irec重新进行记录累

irec=0

do200it=1,nt

do30iz=1,nz

irec=irec+1

write(12,rec=irec)((u(i,j,iz,it),i=1,nx),j=1,ny)

30Continue

do31iz=1,nz

irec=irec+1

write(12,rec=irec)((v(i,j,iz,it),i=1,nx),j=1,ny)

31continue

irec=irec+1

write(12,rec=irec)((sst(i,j,it),i=1,nx),j=1,ny)

200continue

end;数据处理;2.3数据描述文件;1.数据描述文件的概念;;2.数据描述文件的形式;XDEF80LINEAR-140.01.0

YDEF50LINEAR20.01.0

ZDEF10LEVELS1000850700500400300250200150100

TDEF4LINEAR0Z10apr199112hr;