数据文件的转换及数据描述文件的建立试验二

1、南京信息工程大学实验（实习）报告实验（实习）名称数据文件的转换及数据描述文件的建立日期2016.3.18指导老师马利系专业年级_班次姓名学号实习二数据文件的转换及数据描述文件的建立1 .实习资料：现有“data”文件夹ASCII码数据资料文件：月平均风场（u200.dat、u850.datv200.datv850.dat）;月平均高度场（hgt500.dat）。资料的水平网格范围是MXN个网格点（M=37N=17），分辨率为2.50X2.5,范围为：自西向东经度60-1500E,由南至北纬度0-40No时段为2002.12005.12共48个月。2 .实习要求：1）安装GrADS运行软件（1.

2、8版或2.0版），熟悉GrADS!行环境；2）根据所提供的数据文件，编写出将ASCII码数据资料文件（月平均风场和月平均高度场）分别转换成二进制无格式直接存取（GrAD缴据格式）文件（uv.grd和hgt.grd文件）的Fortran程序；3）写出相应的数据描述文件（uv.ctl和hgt.ctl文件），并利用GrADSS本操作命令（open*.ctl;d*）显示该地区850hPa和200hPauv风场及500hPa位势高度场图形，最后保存图形（printim路径*png）。3 .实习目的：通过编写该程序，进一步熟悉数据GrADS的数据格式，熟练编写将ASCII码数据转换为二进制格式的Fortr

3、an程序以及建立相应的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 利用GrADSS本操作命令（open*.ctl,d*）显示2002年7月850hPa和200hPauv风场图形及50

4、0hPa位势高度场图形，最后保存图形（printim路径*.png）；4.7 完成实习报告1）说明所用资料2）给出所编写的Fortran程序3）给出所编写的数据描述文件4）给出所绘制的图形10图12002年1月850hpa风场图5 .实验结果所用资料“data”文件夹下有ASCII码数据资料文件:月平均风场：u200.dat、u850.dat、v200.dat、v850.dat月平均高度场：hgt500.dat时间范围：2002.1-2005.12共48个月水平范围：37X17个网格点，分辨率为2.5X2.5经纬度范围：60-150E,0-40N2 .fortran程序programexam0

5、1parameter(nx=37,ny=17,nt=48)dimensionu200(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)doit=1,ntread(1,*)(u200(i,j,it),i=1,nx),j=1,ny)read(2,*)(u850(i,j,it),i=1,nx)

6、,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)enddodoit=1,ntwrite(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)enddoend3 .数据描写文件uv.ctlDSETD:datauv.ctlUNDEF-9.99E33TITLEUpperAirDataXDEF37LINEAR60.02.5YDEF17LINEAR0.02.5ZDEF2LEVELS850200TDEF48LINEARJAN20021movars2u299uwind(m/s)v299vwind(m/s)endvarshgt.ctlDESTD:datahgt.grdUNDEF-9.99E+33XDEF37linear60.02.5YDEF17linear0.02.5ZDEF1levels500