发布时间 : 星期六 文章GrADS更新完毕开始阅读
实习一 数据文件的转换及数据描述文件的建立
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文件。利用文本编辑器（如记事本，或者写字板）将所需要的基本命令编辑成gs文件，然后使用run命令运行该文件；
（3）根据所学过的图形要素设置命令“set”， 编辑gs文件绘制北半球200hPa纬向风一维单线图（经度固定于120oE）和500hPa高度场平面图，并用“enable print”命令将所画图形分别存放。
4．实习报告： 1）说明所用资料
“data”文件夹下有GrADS数据格式的二进制资料文件model.le.dat以及对应的数据描述文件model.le.ctl
（这是一组模式输出的全球5天数值预报结果，包括了多个要素、多层，按经纬度网格存放的数据，具体见数据描述文件model.le.ctl） 2）所编写的数据描述文件
A．北半球200hPa纬向风一维单线图（经度固定于120oE） 'reinit'
'open F:\\1\\model.le.ctl' 'set lon 120' 'set lat 0 90' 'set lev 200' 'd u'
'printim F:\\1\\shuang2.1.jpg' ;
B.500hPa高度场平面图'reinit'
'open F:\\1\\model.le.ctl' 'set lon 0 360' 'set lat 0 90' 'set lev 500' 'd z'
'printim F:\\1\\shuang2.2.jpg' ;
3）所绘制的图形
A．北半球200hPa纬向风一维单线图（经度固定于120oE）
B.500hPa高度场平面图
实习三 set命令的使用练习