GrADS绘图实习

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码数据资料文件：月平均风场和月平均高度场分别转换成二进制无格式直接存取（GrADS数据格

2、式）文件（uv.grd和hgt.grd文件）的Fortran程序，给出相应的数据描述文件（uv.ctl和hgt.ctl文件），并利用GrADS基本操作命令显示该地区850hPa和200hpa uv风场，最后保存图形。三、实习目的：通过编写该程序，进一步熟悉数据GrADS的数据格式，熟练编写将ASCII码数据转换为二进制格式的Fortran程序以及建立相应的CTL文件。四、实习步骤：4.1 利用“写字板”程序打开*.dat数据文件，熟悉该数据资料；4.2 按要求编写Fortran程序，将所给的ASCII码数据资料文件转换成二进制无格式直接存取文件，结果保存为mhy.grd；4.3 通过“写字板”

3、或“记事本”程序编写相应的数据描述文件，保存为mhy.ctl；4.4 画出2002年1月850hpa风场图，与图1比较验证数据正确性；4.5 利用GrADS基本操作命令（open,display,print等）显示2002年7月850hpa和200hpa uv风场图形；4.5 完成实习报告1）说明所用资料 资料名称：月平均风场 (u200.dat、u850.dat、v200.dat、v850.dat)；月平均高度场 (hgt500.dat)范围：60150E，040N.层次：风场，850 hPa、200hPa；高度场，500hPa；时段：2002.12005.12共48个月。变量名：u v h

4、gt分辨率：2.5*2.5，格点数：37*17资料读取文件：data.f902）给出所编写的Fortran程序： parameter(nx=37,ny=17,nz=1,nt=48)dimension u1(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='u200.dat')open(2,file='v200.dat')open(3,file='u850.dat')open(4,file='v850.dat&#

5、39;)open(11,file='hgt500.dat')open(12,file='uv.grd',form='binary') open(13,file='hgt500.grd',form='binary')do 100 it=1,nt do 50 iz=1,nz read(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) 50 continue do 51 iz=1,nz read(2,*) (v1(i,

6、j,iz,it),i=1,nx),j=1,ny) read(4,*) (v2(i,j,iz,it),i=1,nx),j=1,ny) 51 continue do 52 iz=1,nz read(11,*) (hgt(i,j,iz,it),i=1,nx),j=1,ny) 52 continue100 continuedo 200 it=1,nt do 30 iz=1,nz write(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) 30 Continue do 31 iz=1,nz write

7、(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) 31 continue do 55 iz=1,nz write(13) (hgt(i,j,iz,it),i=1,nx),j=1,ny) 55 continue200 continueCLOSE(1)CLOSE(2)CLOSE(3)CLOSE(4)CLOSE(11)CLOSE(12)CLOSE (13)end3）给出所编写的数据描述文件uv.ctl：dset C:GrADS19uv.grdundef -9.99E+13title two lev

8、els air dataxdef 37 linear 60.0 2.5ydef 17 linear 0.0 2.5zdef 2 levels 850 200tdef 48 linear JAN2002 1movars 2u 2 99 u wind(m/s)v 2 99 v wind(m/s)endvarshgt.ctldset C:GrADS19hgt.grdundef -9.99E+13title the attitude of 500hpaxdef 37 linear 60.0 2.5ydef 17 linear 0.0 2.5zdef 1 levels 500tdef 48 linear

9、 JAN2002 1movars 1h 0 99 attitudeendvars文件*.GS'reinit''open c:datauv.ctl' *ctl路径'set grid off''set grads off''set lat 0 40'set lon 60 150'' set time jul2002''set lev 850''set gxout vector''d u;v''draw title Air data in 2

10、00201 850''printim c:imagesuv8.png white''c''set lev 200''set gxout vector''d u;v''draw title Air data in 200201 200''printim c:imagesuv2.png white'4）给出所绘制的图形2002年1月850hpa高度的风场图2002年7月850hpa 高度的风场图2002年7月200hpa高度的风场图：实习二 基本命令的上机使用（2小时）1实习资料

11、：现在“data”文件夹下有GrADS数据格式的二进制资料文件model.le.dat（一组模式输出的全球5天数值预报结果，包括了多个要素、多层，按经纬度网格存放的数据），以及对应的数据描述文件model.le.ctl。2实习要求：利用所提供的数据文件，使用基本命令编辑gs文件绘制图形一维单线图，和两维平面图。3. 实习目的：熟悉GrADS的工作环境及使用流程，掌握基本操作命令，练习编写和使用gs文件。4. 实习步骤：4.1熟悉GrADS绘图软件操作界面命令窗口和图形显示窗口。利用所提供的数据文件model.le.dat及相应的数据描述文件model.le.ctl，在命令窗口输入基本绘图指令，

12、在图形显示窗口观察相应图形例如：在文本窗口输入下列命令open <路径>model.le.ctlset lat 40set lon -180 0set lev 500set t 1d z则在图形显示窗口出现如下图形4.2 练习编写和使用gs文件。利用文本编辑器（如记事本，或者写字板）将所需要的基本命令编辑成gs文件，然后使用run命令运行该文件；Model.le.ctl：'open c:datamodel.le.ctl''set lat 40''set lon -180 0''set lev 500''set

13、t 1''d z''printim c:imagesim40500.png white'4.3 根据所学过的图形要素设置命令“set”， 编辑gs文件绘制北半球200hPa纬向风一维单线图（经度固定于120ºE）和500hPa高度场平面图，并用“enable print”命令将所画图形分别存放。Gs:'reinit''open c:datamodel.le.ctl''enable print c:imagesmodel2.gmf''set grid off''set gra

14、ds off''set lon 120''set lat 0 90''set lev 200''set t 1''set xyrev on''d u''print''c''reset''set grid off''set clip 0 11 2 6''set lat 0 90''set lev 500''set t 1''set grads off'

15、'd z''print''disable print'实习三 set命令的使用练习（2小时）1实习资料：现在“data”文件夹下有二进制资料文件uv.grd、hgt.grd以及对应的数据描述文件uv.ctl、hgt.ctl。2实习要求：1）利用所提供的数据文件，绘制出第7时刻500hPa高度场图，要求设置等值线间隔单位为40位势米，颜色，标记方式，并且给出相应标题（请包含姓名拼音与学号），最终将图形保存。所有命令编写于.gs文件中。2）利用所提供的数据文件，绘制出第7时刻850hPa风场图，要求以三种方式显示风场，并将三幅图上下排列于同一图形窗口

16、，并最终将图形保存。所有命令编写于.gs文件中。3. 实习目的：熟练使用set命令的各类参数设置，掌握.gs文件的编写格式和运行，学会保存图形文件。4. 实习步骤：4.1 在写字板或记事本程序中按要求编写*.gs文件，注意文件编写格式；4.2 启动GrADS绘图软件，运行*.gs文件；4.3完成实习报告：1）说明所用资料；资料名称：月平均风场 (uv.grd)；月平均高度场 (hgt.grdt)范围：60150E，040N.层次：风场，850 hPa、200hPa；高度场，500hPa；时段：2002.12005.12共48个月。变量名：u v h分辨率：2.5*2.5，格点数：37*172）

17、给出所编写的*.gs文件；第7时刻500hPa高度场图的gs文件hgt7.gs'reinit''open c:datahgt.ctl''set grid off''set grads off''set xlopts 3 5 0.18''set ylopts 3 5 0.18''set parea 1 10.2 0.8 7.5''set t 7''set lev 500''hgt=h/10''set cint 4''

18、;d hgt''set annot 5 8''draw title dingxiaoli 20081331001''printim c:imageshgt7.png white'第7时刻850hPa风场图的gs文件gxout.gs：3）给出所绘图形：'reinit''open c:datauv.ctl''set t 7''set lev 850''set grid off''set vpage 0 8.5 0 3.6''set pare

19、a 0.3 8.5 0.3 3.5''set gxout vector''set grads off''d u;v''set vpage 0 8.5 3.6 7.2''set grads off''set gxout barb''d u;v''set vpage 0 8.5 7.2 11''set grads off''set gxout STREAM''d u;v''printim c:imagesgx

20、out3uv850.png white'第7时刻500hPa高度场图第7时刻850hPa风场图（三种输出形式）： 从上到下依次为：grid、barb、vector实习四 各类参数设置综合练习 (2小时)1实习资料：现在“data”文件夹下有二进制资料文件uv.grd、dxgd.grd以及对应的数据描述文件uv.ctl、dxgd.ctl。2实习要求：利用所提供的数据文件，绘制出2003年7月60-150E、0-40N区域内850hPa流线图，且地形高度场超过3000米以上用黑色阴影显示，并且给出相应标题（请包含姓名拼音与学号），最终将图形保存。所有命令编写于.gs文件中。3. 实习目的：

21、掌握阴影图形与等值线图形的叠加画图方法，同时练习基础绘图指令draw命令的使用。4. 实习步骤：4.1 在写字板或记事本程序中按要求编写*.gs文件，注意文件编写格式4.2 启动GrADS绘图软件，运行*.gs文件4.3完成实习报告1）说明所用资料 资料名称：1、月平均风场 (uv.grd)；范围：60150E，040N.层次：风场，850 hPa、200hPa；时段：2002.12005.12共48个月。变量名：u v hgt分辨率：2.5*2.5，格点数：37*172、高度场（dxgd.grd）；范围: 经向0360；纬向：-90S90N层次：一个层次 地面层时段：2002.1共1个月变量

22、名：h分辨率：2.5*2.52）给出所编写的（uv073000.gs）文件'reinit''open c:datadxgd.ctl''open c:datauv.ctl''set grid off''set grads off''set map 3 1 5''set xlopts 3 5 0.18''set ylopts 3 5 0.18''set parea 0.7 10.4 0.6 7.8''set t 1''set lon

23、 60 150''set lat 0 40' 'set gxout shaded''set rbcols 1 1 1 1 1 1 1''set cmin 3000''d h''set gxout stream''set time jul2003''set grads off''set lev 850''d u.2;v.2''set annot 5 8''draw title DingXiaoli 2008

24、1331001''printim c:imagesuv200307850.png white'3）给出所绘图形实习五 各类参数设置综合练习 (2小时)1实习资料：现在“data”文件夹下有二进制资料文件uv.grd以及对应的数据描述文件uv.ctl。2实习要求：利用所提供的数据文件，画出2002年1-12月120 E、 0-40N 200hPa纬向风的纬度-时间剖面图。图型要求：（1）纬向风为西风时填色，东风绘制等值线，给出色标，0值线加粗；（2）X轴标注为“time”,Y轴标注为“lat”,标题标注为“u 100-120E”。（3）最终将图形保存为gmf格式。所有命令

25、编写于.gs文件中。3. 实习目的：掌握特殊图形处理方法，掌握阴影图形画图方法，同时学会处理特别等值线与特殊字符的方法。4. 实习步骤：4.1 在写字板或记事本程序中按要求编写*.gs文件，注意文件编写格式4.2 启动GrADS绘图软件，运行*.gs文件4.3完成实习报告1）说明所用资料月平均风场 (uv.grd)；范围：60150E，040N.层次：风场，850 hPa、200hPa；时段：2002.12005.12共48个月。变量名：u v hgt分辨率：2.5*2.5，格点数：37*172）给出所编写的*.gs文件'reinit''open c:datauv.ct

26、l''set grid off''set grads off''set xlopts 3 5 0.16''set ylopts 3 5 0.16''set parea 1 10.4 1.2 7.9''set lon 120''set lat 0 40''set lev 200''set t 1 12''set annot 5 8''enable print c:imagesuv200.gmf''set x

27、yrev on''set gxout shaded''set cmin 0''d u''cbarn 1 0 5.5 0.3''set grads off''set gxout contour''set cmax 0''set grads off''d u''set clevs 0''set ccolor 2''set cthick 10''set grads off''d

28、u''draw xlab time''draw ylab lat''draw title u 100-120E''print''disable print'3）给出所绘图形实习六 练习描述语言的使用（2小时）1实习资料：现在“data”文件夹下有全球海表温度资料sst.mnmean.nc。2实习要求：1）根据所提供的数据资料，利用fwrite命令将该范围内（纬度：-20-60ºN，经度：90ºE-200Eº，时间：1971年1月至2000年12月）海温场资料写入名为sst.gr

29、d的二进制数据文件中。注意确定写出数据的格点数。2）根据所得数据资料，利用Fortran程序将sst.grd文件转换为十进制数据文件sst.txt。3）根据所得sst.grd数据文件，编写相应的数据描述文件sst.ctl，并绘制图形。请将所绘图形与使用原来sst.mnmean.nc资料在相同维数环境下所画图形进行比较，看是否一致。3. 实习目的：掌握GrADS软件数据提取方法，学会资料处理。4. 实习步骤：4.1 在写字板或记事本程序中按要求1）编写*.gs文件，进行数据资料的提取；'reinit''sdfopen c:datasst.mnmean.nc'

30、9;set fwrite c:datasst.grd''set gxout fwrite''set time jan1971 dec2000''set lat -20 60''set lon 90 200''d sst''disable fwrite'4.2 根据实习要求2），利用Fortran程序进行二进制与十进制数据文件的转换，并查看sst.txt的数据信息；Fortran程序：parameter(nx=56,ny=41,nz=1,nt=360)dimension sst(nx,ny,n

31、z,nt)open(15,file='sst.grd',form='binary')open(16,file='sst.txt')do it=1,nt do iz=1,nz read(15) (sst(i,j,iz,it),i=1,nx),j=1,ny) enddoenddodo it=1,nt do iz=1,nz write(16,*) (sst(i,j,iz,it),i=1,nx),j=1,ny)enddoenddoclose(15)close(16)end4.3 在写字板或记事本程序中编写数据描述文件sst.ctl，并保存；Ctl：dse

32、t C:datasst.grdundef -9.99E+33title sea surface tempturexdef 56 linear 90 2ydef 41 linear -20 2tdef 360 linear jan1971 1mozdef 1 levels 0vars 1sst 0 99 surface sea temptureendvars4.4 利用sst.grd和sst.ctl文件画图，完成实习内容。实习七 变量与函数的设置（2小时）1实习资料：现在“data”文件夹下有全球海表温度资料sst.mnmean.nc。2实习要求：1）根据所提供的数据资料，计算1-12月多年平均

33、（1971年-2000年）海温场，并绘出1月海温多年平均图。2）进一步计算1971年至2000年1月海温距平场，绘制1998年1月和1999年1月海温距平场，海温正距平区填色，画色标，并分别写标题Jan1998和Jan1999。要求设置虚页，将两幅图并排放在一起。，最终将图形保存。所有命令编写于.gs文件中。3. 实习目的：掌握函数的使用，进一步练习阴影图形与等值线图形的叠加画图方法。4. 实习步骤：4.1 在写字板或记事本程序中按要求编写*.gs文件；4.2启动GrADS绘图软件，运行*.gs文件4.3完成实习报告1）说明所用资料：sst.mnmean.nc文件：名称：全球海表月平均温度资料

34、时间：1854年1月2009年9月，各月数据共1869个时间序列经度：0 36 纬度：88S 88N分辨率：2*2变量：timebnds 1层 sst 1层（海表温度）2）给出所编写的*.gs文件'reinit''sdfopen c:datasst.mnmean.nc''set grid off''set grads off''set map 1 1 5''set annot 5 8''set xlopts 3 5 0.16''set ylopts 3 5 0.16'

35、'set xlint 20''set ylint 20''set parea 0.6 10.6 0.6 8.0''set t 1 12''asst=ave(sst,t+1404,t=1764,12)''modify asst seasonal''set t 1''d asst''draw title January aveage sst''printim c:imagesavesstjanuary.png white''c

36、9;'set map 1 1 5''set xlopts 3 5 0.14''set ylopts 3 5 0.14''set xlint 40''set ylint 20''set grid off''set grads off''set parea 0.6 5.5 0.6 8.0''set t 1405 1764''nasst=sst-asst''set time JAN1998''set gxout sha

37、ded''set cmin 0''d nasst''cbarn 0.5 0 3.0 2.1''set grads off''set gxout contour''d nasst''draw title JAN1998''set grads off''set parea 6.1 10.8 0.6 8.0''set time jan1999''set gxout shaded''set cmin 0''d nasst''cbarn 0.5 0 8.4 2.1''set grads off''set gxout contour''d nasst''draw title JAN1999''printim c:imagesansstJAN1998with1999.png white'3）给出所绘图形实习八 描述语言的使用（2小时）1实习资料：现在“data”文件夹下有全球海表温度资料sst.mnmea