GrADS绘图实习.doc

实习一数据文件的转换及数据描述文件的建立一、实习资料:现在“data”文件夹下有ASCII码数据资料文件:月平均风场(、、、)、月平均高度场(),-,资料的水平网格范围是M×N个网格点(M=37,N=17),°×°,自西向东经度60-150°E,由南至北纬度0-40°N。二、实习要求:编写出将ASCII码数据资料文件:月平均风场和月平均高度场分别转换成二进制无格式直接存取(GrADS数据格式)文件()的Fortran程序,给出相应的数据描述文件(),并利用GrADS基本操作命令显示该地区850hPa和200hpauv风场,最后保存图形。三、实习目的:通过编写该程序,进一步熟悉数据GrADS的数据格式,熟练编写将ASCII码数据转换为二进制格式的Fortran程序以及建立相应的CTL文件。四、实习步骤:“写字板”程序打开*.dat数据文件,熟悉该数据资料;,将所给的ASCII码数据资料文件转换成二进制无格式直接存取文件,;“写字板”或“记事本”程序编写相应的数据描述文件,;,与图1比较验证数据正确性;(open,display,print等)显示2002年7月850hpa和200hpauv风场图形;)说明所用资料资料名称:月平均风场(、、、);月平均高度场()范围:60~150E,0~:风场,850hPa、200hPa;高度场,500hPa;时段:~。变量名:uvhgt分辨率:*,格点数:37*17资料读取文件:)给出所编写的Fortran程序:parameter(nx=37,ny=17,nz=1,nt=48)dimensionu1(nx,ny,nz,nt),v1(nx,ny,nz,nt),u2(nx,ny,nz,nt),v2(nx,ny,nz,nt),hgt(nx,ny,nz,nt)open(1,file='')open(2,file='')open(3,file='')open(4,file='')open(11,file='')open(12,file='',form='binary')open(13,file='',form='binary')do100it=1,ntdo50iz=1,nzread(1,*)((u1(i,j,iz,it),i=1,nx),j=1,ny) read(3,*)((u2(i,j,iz,it),i=1,nx),j=1,ny)50continuedo51iz=1,nzread(2,*)((v1(i,j,iz,it),i=1,nx),j=1,ny) read(4,*)((v2(i,j,iz,it),i=1,nx),j=1,ny)51continuedo52iz=1,nzread(11,*)((hgt(i,j,iz,it),i=1,nx),j=1,ny)52continue100continuedo200it=1,ntdo30iz=1,nzwrite(12)((u2(i,j,iz,it),i=1,nx),j=1,ny)write(12)((u1(i,j,iz,it),i=1,nx),j=1,ny) 30Continuedo31iz=1,nzwrite(12)((v2(i,j,iz,it),i=1,nx),j=1,ny)write(12)((v1(i,j,iz,it),i=1,nx),j=1,ny)31continuedo55iz=1,nzwrite(13)((hgt(i,j,iz,it),i=1,nx),j=1,ny)55continue200continueCLOSE(1)CLOSE(2)CLOSE(3)CLOSE(4)CLOSE(11)CLOSE(12)CLOSE(13)end3):dsetC:\GrADS19\-+