GrADS..docVIP

实习一 数据文件的转换及数据描述文件的建立 1．实习目的： (1)通过编写程序，熟悉数据GrADS的数据格式； (2)熟练编写将ASCII码数据转换为二进制格式的Fortran程序； (3)熟练建立相应的CTL文件。 2. 实习要求： (1)编写Fortran程序，将ASCII码数据资料文件转换成二进制文件：月平均风场（uv.grd） 月平均高度场（hgt.grd） (2)写出二进制文件相应的数据描述文件： uv.ctl，hgt.ctl (3)利用GrADS基本操作命令显示该地区850hPa和200hPa uv风场，最后保存图nm形。 3 实习步骤： (1) 利用“写字板”程序打开*.dat数据文件，熟悉该数据资料； (2 )按要求编写Fortran程序，将所给的ASCII码数据资料文件转换成二进制无格式直接存取文件，结果保存为**.grd； (3) 通过“写字板”或“记事本”程序编写相应的数据描述文件，保存为**.ctl； (4) 画出2002年1月850hPa风场图，与图1比较验证数据正确性； 4.实习报告 1）所用资料 “data”文件夹下有ASCII码数据资料文件： 月平均风场：u200.dat、u850.dat、 v200.dat、v850.dat 月平均高度场：hgt500.dat 经纬度范围：60-150°E， 0-40°N 时间范围：2002.1-2005.12 共48个月 水平范围：37×17个网格点，分辨率为2.5°×2.5° 2）所编写的Fortran程序 parameter(nx=37,ny=17,nt=48) dimension u200(nx,ny,nt),v200(nx,ny,nt),h500(nx,ny,nt), u850(nx,ny,nt),v850(nx,ny,nt) open(1,file=e:\dwj\u200.dat) open(2,file=e:\dwj\v200.dat) open(3,file=e:\dwj\u850.dat) open(4,file=e:\dwj\v850.dat) open(5,file=e:\dwj\hgt500.dat) do it=1,nt read(1,*)((u200(i,j,it),i=1,nx),j=1,ny) read(2,*)((v200(i,j,it),i=1,nx),j=1,ny) read(3,*)((u850(i,j,it),i=1,nx),j=1,ny) read(4,*)((v850(i,j,it),i=1,nx),j=1,ny) read(5,*)((h500(i,j,it),i=1,nx),j=1,ny) enddo close(1);close(2);close(3);close(4);close(5) open(12,file=e:\dwj\uv0.grd,form=binary) 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) enddo close(12) open(12,file=e:\dwj\hgt0.grd,form=binary) do it=1,nt write(12) ((h500(i,j,it),i=1,nx),j=1,ny) enddo close(12) end 3）所编写的数据描述文件 reinit open F:\1\uv.ctl set lon 60 150 set lat 0 40 set t 1 d u;v printim F:\1\shuang.jpg ; 4）所绘制的图形 实习二 基本命令的上机使用 1. 实习目的： 熟悉GrADS的工作环境及使用流程，掌握基本操作命令，练习编写和使用gs文件。 2．实习要求： 利用所提供的数据文件，使用基本命令编辑gs文件绘制图形一维单线图和两维平面图。 3. 实习步骤： （1）熟悉GrADS绘图软件操作界面——命令窗口和图形显示窗口。利用所提供的数据文件model.le.dat及相应的数据描述文件model.le.ctl，在命令窗口输入基本绘图指令，在图形显示窗口观察相应图形 （2）练习编写和使用gs文件。利用文本编辑器（如记事本