2章 数据处理.pptVIP

第2章 数据处理;2.1 各类数据格式的特点;GRIB码数据格式：是自定义的数据格式，压缩率高，占用空间小。例如NCEP提供的数据资料。 NETCDF（.nc）等通用数据格式：自定义数据格式，精确性好，便于传输。 文件中自带描述文件。;2.2 数据文件的转换 ;2.数据存放形式 ;3. 举例（格点资料） ;编程时数组大小： X方向：150-60=90，90/2.5+1=37 Y方向：40-0=40，40/2.5+1=17 Z方向： u、v为850、 200hPa ，Z为2 sst为海表温度，只有一层，Z为1 T时次：月资料，4年，共48个月 ;Fortran程序编写如下： ;! 打开原始数据文件,如果数据文件和Fortran程序在一个文件夹下，打开数据文件时可以不写路径。 open(1,file=u.dat) open(2,file=v.dat) open(3,file=sst.dat) ! 打开目标文件，recl为记录的长度，对于直接存取文件必须指定记录长度 。 open(12,file=mhy.grd,form=unformatted, #access=direct, recl=nx*ny*4) ;;!将数据资料写入新文件中， 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;;;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;如果不考虑采用无格式直接方式存取文件的话，可以使用下面的编程形式，较为简单。;;;注意事项：;2.2 数据描述文件（.ctl） ;1. 数据描述文件的概念;2. 数据描述文件的形式 ;以下是另一个格点数据描述文件的例子： * this is an example to demonstrate the data descriptor file DSET ua.dat TITLE Upper Air Data DTYPE grid OPTION byteswapped UNDEF -9.99E33;XDEF 80 LINEAR -140.0 1.0 YDEF 50 LINEAR 20.0 1.0 ZDEF 10 LEVELS 1000 850 700 500 400 300 250 200 150 100 TDEF 4 LINEAR 0Z10apr1991 12hr ;vars 6 slp 0 99 sea level pressure z 10 99 heights t 10 99 temps td 6 99 dewpoints u 10 99 u winds v 10 99 v winds endvars;3.数据描述文件的构成 ;（2）数据描述文件中各记录的含义：;TITLE string 用字符串string简略描述数据文件的内容。即：数据的标题 UNDEF value 定义缺测值。 一般给一很大的正值或负值。GrADS在运算操作和图形操作时将忽略这些值点。;OPTIONS keywords 定义数据格式选项，keywords 有： yrevzrevsequentialbyteswapped templatebig-endiantitlee-endian 缺省时为 direct 直接记录格式（form=binary） yrev 自北