资源描述：
《GrADS绘图实习.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、.实习一数据文件的转换及数据描述文件的建立一、实习资料：现在“data”文件夹下有ASCII码数据资料文件：月平均风场(u200.dat、u850.dat、v200.dat、v850.dat)、月平均高度场(hgt500.dat)，2002.1-2005.12共48个月，资料的水平网格范围是M×N个网格点（M=37，N=17）,分辨率为2.5°×2.5°，自西向东经度60-150°E，由南至北纬度0-40°N。二、实习要求：编写出将ASCII码数据资料文件：月平均风场和月平均高度场分别转换成二进制无格式直接
2、存取（GrADS数据格式）文件（uv.grd和hgt.grd文件）的Fortran程序，给出相应的数据描述文件（uv.ctl和hgt.ctl文件），并利用GrADS基本操作命令显示该地区850hPa和200hpauv风场，最后保存图形。三、实习目的：通过编写该程序，进一步熟悉数据GrADS的数据格式，熟练编写将ASCII码数据转换为二进制格式的Fortran程序以及建立相应的CTL文件。四、实习步骤：4.1利用“写字板”程序打开*.dat数据文件，熟悉该数据资料；4.2按要求编写Fortran程序，将所给的
3、ASCII码数据资料文件转换成二进制无格式直接存取文件，结果保存为mhy.grd；4.3通过“写字板”或“记事本”程序编写相应的数据描述文件，保存为mhy.ctl；4.4画出2002年1月850hpa风场图，与图1比较验证数据正确性；4.5利用GrADS基本操作命令（open,display,print等）显示2002年7月850hpa和200hpauv风场图形；4.5完成实习报告Word资料.1）说明所用资料资料名称：月平均风场(u200.dat、u850.dat、v200.dat、v850.dat)；月
4、平均高度场(hgt500.dat)范围：60~150E，0~40N.层次：风场，850hPa、200hPa；高度场，500hPa；时段：2002.1~2005.12共48个月。变量名：uvhgt分辨率：2.5*2.5，格点数：37*17资料读取文件：data.f902）给出所编写的Fortran程序：Word资料.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
5、,nz,nt),hgt(nx,ny,nz,nt)open(1,file='u200.dat')open(2,file='v200.dat')open(3,file='u850.dat')open(4,file='v850.dat')open(11,file='hgt500.dat')open(12,file='uv.grd',form='binary')open(13,file='hgt500.grd',form='binary')do100it=1,ntdo50iz=1,nzread(1,*)((u1(i
6、,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)Word资料.51continuedo52iz=1,nzread(11,*)((hgt(i,j,iz,it),i=1,nx),j=1,ny)52continue100continu
7、edo200it=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,
8、nx),j=1,ny)55continue200continueCLOSE(1)CLOSE(2)CLOSE(3)CLOSE(4)CLOSE(11)CLOSE(12)CLOSE(13)endWord资料.3）给出所编写的数据描述文件Word资料.uv.ctl：dsetC:GrADS19uv.grdundef-9.99E+13titletwolevelsairdataxdef37linear60.02.5ydef
