GrADS绘图实习.docx

实习一数据文件的转换及数据描述文件的建立
一、 实习资料：
现在"data"文件夹下有ASCII码数据资料文件：月平均风场(ii2()().dar、 \ )月平均髙度场(), -，资料的水平网格范围是M X N 个网格点(M=37, N二17)，° ° ,自西向东经度60-150° E,由南至北纬度 0-40° No
二、 实习要求：
编写出将ASCII码数据资料文件：月平均风场和月平均高度场分别转换成二进制无格式直 接存取(GrADS数据格式)文件()的Fortran程序，给出相应的数据描 述文件()，并利用GrAOS基本操作命令显示该地区850hPa和200hpa uv 风场，最后保存图形。
三、 实习目的：
通过编写该程序，进一步熟悉数据GrADS的数据格式，熟练编写将ASCII码数据转换为 二进制格式的Fortran程序以及建立相应的CTL文件。
四、 实习步骤：
“写字板”程序打开乞cht数据文件，熟悉该数据资料；
，将所给的ASCII码数据资料文件转换成二进制无格式直接存取文 件， ；
“写宇板”或“记事本”程序编写相应的数据描述文件，；
，与图1比较验证数据正确性；
利用 GrADS 基本操作命令(open,display,print 等)显示 2002 年 7 月 850hpa 和 200hpa uv 风 场图形；

1）说明所用资料
资料名称：月平均风场(> . > );
月平均高度场()
范團：60-150E, 0〜40N.
层次：风场，850 hPax 200hPa；高度场，500hPa; 时段：〜。
变量名：u v hgt
分辨率：*,格点数：37*17
资料读取文件：
2）给出所编写的Fortran程序:
parameter(nx=37^iy=17^iz=l qt=48) dimension
ul (nx^iy^iz/it),vl (nx/iy^ynt) ,u2(nx/iy xQz^Qt) ,v
open(l ,file=,,)
open(2,file=,v2(X).datO
open(3,file=',)
open(4,file=,)
open(l 141«=，)
open(l 231e=',fbrm=,binary,)
open(l 3^e=^,>fbrm=tbinaiy1)
do 100it=l^it
do 50
read(l *) ((ul(L,jJzJt)4=l 4ix),j=My)
read(3,*) ((112°3期)=1皿)沪1坷)
continue
do 51 iz=la
read(2,*) ((vl(i,jjzjt)Ala),j=lqy) read(4*) ***@2①诈
continue
do 52 iz=lqz
read(l 1,*) ((hgt@,j»iz»nx),j=l qy)
continue
3)给出所编写的数据描述文件
100 continue
do 200 it=3
do 30 iz=la
write(12) ((u2Q,j血jt)j=la),j=lAy)
wnte(12) ((ulQQ期戶ia),j=Sy)
Continue do 31 iz=l^i2
write(12) ***@2码逊为=1同沪Sy)
writc(12) ***@1©立期 卫x),j=lqy)
continue
do 55 iz=l^iz
write(13) @才(12旳」=1対沪1再)
55 continue
200 continue
CLOSE(l)
CLOSE(2)
CLOSE ⑶
CLOSE(4)
CLOSE(ll)
CLOSE(12)
CLOSE (13)
end
dset C:\GrADS19\
xdef37 linear
undef・+13
ydef 17 linear
title two level