grads上机实习报告.docVIP

实习一 数据文件的转换及数据描述文件的建立 1．实习资料： 现有ASCII码数据资料文件hi.dat（1973年4月29日08时我国东北、华北地区500hPa等压面位势高度场资料），资料的水平网格范围是M×N个网格点（M=20，N=16）,分辨率为3.5°×2.5°，自西向东经度从88.5°E起，由南至北纬度从32.5°N起。 2. 实习要求： 1）根据机器条件，试安装GrADS运行软件（1.8版或2.0版），熟悉GrADS运行环境； 2）根据所提供的数据文件，编写出将ASCII码数据资料文件转换成二进制无格式直接存取（GrADS数据格式）文件（.grd文件）的Fortran程序，给出相应的数据描述文件（.ctl文件），并利用GrADS基本操作命令（open *.ctl; d *）显示该地区500hPa高空形势场，最后保存图形(printim 路径*.png)。 3．实习目的： 通过编写该程序，进一步熟悉数据GrADS的数据格式，熟练编写将ASCII码数据转换为二进制格式的Fortran程序以及建立相应的CTL文件。 4. 实习步骤： 4.1 安装GrADS运行软件； 4.2 利用“写字板”程序打开hi.dat数据文件，熟悉该数据资料； 4.3 按要求编写Fortran程序，将所给的ASCII码数据资料文件转换成二进制无格式直接存取文件，结果保存为*.grd； 4.4 通过“写字板”或“记事本”程序编写相应的数据描述文件，保存为*.ctl； 4.4 利用GrADS基本操作命令（open,display,printim）显示图形，并保存图片文件*.png，存图命令参考教材P32； 4.5 完成实习报告 （2）转化为二进制文件：  parameter(ii=16,jj=20) real var(ii,jj) open(1,file=k:\grads1\hi.dat) open(7,file=k:\grads1\data.grd,form=binary) read(1,*) ((var(i,j),i=1,ii),j=1,jj) close(1) write(7)((var(i,j),i=1,ii),j=1,jj) close(7) end 数据描述文件： dset k:\grads1\data.grd undef -9.99E+33 title The Height Field Of China xdef 20 linear 88.500 3.500 ydef 16 linear 32.500 2.500 zdef 1 levels 500 tdef 1 linear APR1983 1dy vars 1 h 1 99 h500 endvars GS文件： reinit open k:\grads1\hi.ctl enable print k:\grads1\h.gmf set xlopts 1 7 0.15 set ylopts 1 7 0.15 d h draw title h in china draw xlab lon draw ylab lat print disable print ;  （4）画有关图形： 实习二 站点数据的转化 将所给有关数据转化为二进制格点数据，并配以相应的数据描述文件。 parameter(ii=37,jj=17) real var(ii,jj) open(1,file=d:\li\gd\grd\u850.dat) open(2,file=d:\li\gd\grd\v850.dat) open(3,file=d:\li\gd\grd\u200.dat) open(4,file=d:\li\gd\grd\v200.dat) open(5,file=d:\li\gd\grd\h500.dat) open(6,file=d:\li\gd\grd\tsfc.dat) open(7,file=k:\li\data.grd,form=binary) DO 200 iy=1,4 do 100 m=1,12 ccc read u850 read(1,1000) read(1,2000) ((var(i,j),I=1,ii),j=1,jj) write(7) ((var(i,j),I=1,ii),j=1,jj) ccc read u200 read(3,1000) read(3,2000) ((var(i,j),I=1,ii),j=1,jj) write(