数据文件的转换及数据描述文件的建立 实验二 - 图文

南京信息工程大学 实验（实习）报告

实验（实习）名称数据文件的转换及数据描述文件的建立日期2016.3.18指导老师 马利 系 专业 年级 班次姓名 学号

实习二 数据文件的转换及数据描述文件的建立

1．实习资料：

现有“data”文件夹ASCII码数据资料文件：月平均风场 (u200.dat、u850.dat、v200.dat、v850.dat)；月平均高度场 (hgt500.dat)。资料的水平网格范围是M×N个网格点（M=37，N=17）,分辨率为2.5°×2.5°，范围为：自西向东经度60-150°E，由南至北纬度0-40°N。时段为2002.1~2005.12共48个月。

2. 实习要求：

1）安装GrADS运行软件（1.8版或2.0版），熟悉GrADS运行环境；

2）根据所提供的数据文件，编写出将ASCII码数据资料文件（月平均风场和月 平均高度场）分别转换成二进制无格式直接存取（GrADS数据格式）文件 （uv.grd和hgt.grd文件）的Fortran程序；

3）写出相应的数据描述文件（uv.ctl和hgt.ctl文件），并利用GrADS基本操 作命令（open *.ctl; d *）显示该地区850hPa和200hPa uv风场及500hPa 位势高度场图形，最后保存图形(printim <路径>*.png)。。

3．实习目的：

通过编写该程序，进一步熟悉数据GrADS的数据格式，熟练编写将ASCII码数据转换为二进制格式的Fortran程序以及建立相应的CTL文件。

4. 实习步骤：

4.1 安装GrADS运行软件（1.8版或2.0版）；

4.2 利用“写字板”程序打开*.dat数据文件，熟悉该数据资料；

4.3 按要求编写Fortran程序，将所给的ASCII码数据资料文件转换成二进制无 格式直接存取文件，结果保存为uv.grd和hgt.grd；

4.4 通过“写字板”或“记事本”程序编写相应的数据描述文件：uv.ctl和

hgt.ctl；

4.5 画出2002年1月850hPa风场图，与图1比较验证数据正确性； 4.6 利用GrADS基本操作命令（open *.ctl,d *）显示2002年7月850hPa和

200hPa uv风场图形及500hPa位势高度场图形，最后保存图形(printim < 路径>*.png)； 4.7 完成实习报告 1）说明所用资料

2）给出所编写的Fortran程序 3）给出所编写的数据描述文件 4）给出所绘制的图形

图1 2002年1月850hpa风场图

5.实验结果

所用资料

“data”文件夹下有ASCII码数据资料文件： 月平均风场：u200.dat、u850.dat、

v200.dat、v850.dat 月平均高度场：hgt500.dat

时间范围：2002.1-2005.12 共48个月

水平范围：37×17个网格点，分辨率为2.5°×2.5° 经纬度范围：60-150°E， 0-40°N

2.fortran程序

program exam01

parameter(nx=37,ny=17,nt=48)

dimension u200(nx,ny,nt),u850(nx,ny,nt),v200(nx,ny,nt),v850(nx,ny,nt) open(1,file='D:\%u200.dat') open(2,file='D:\%u850.dat') open(3,file='D:\\v200.dat') open(4,file='D:\\v850.dat')

open(12,file='D:\%uv.grd',form='binary') do it=1,nt

read(1,*) ((u200(i,j,it),i=1,nx),j=1,ny) read(2,*) ((u850(i,j,it),i=1,nx),j=1,ny) read(3,*) ((v200(i,j,it),i=1,nx),j=1,ny) read(4,*) ((v850(i,j,it),i=1,nx),j=1,ny) end do do it=1,nt

write(12) ((u850(i,j,it),i=1,nx),j=1,ny) write(12) ((u200(i,j,it),i=1,nx),j=1,ny) write(12) ((v850(i,j,it),i=1,nx),j=1,ny) write(12) ((v200(i,j,it),i=1,nx),j=1,ny) end do end

3.数据描写文件 uv.ctl

DSET D:\\data\%uv.ctl UNDEF -9.99E33

TITLE Upper Air Data XDEF 37 LINEAR 60.0 2.5 YDEF 17 LINEAR 0.0 2.5 ZDEF 2 LEVELS 850 200

TDEF 48 LINEAR JAN2002 1mo vars 2

u 2 99 u wind(m/s) v 2 99 v wind(m/s) endvars

hgt.ctl

DEST D:\\data\\hgt.grd UNDEF -9.99E+33

XDEF 37 linear 60.0 2.5 YDEF 17 linear 0.0 2.5 ZDEF 1 levels 500

TDEF 48 linear JAN2002 1mo vars 1

hgt 1 99 q1 endvars 4.所绘制图形

2002年一月850hpa风场图

①??

② 2002年7月850hpa风场图