GRADS实习作业.doc

实习一1）编写简单的gs文件，分别绘制一个简单的一维图形（曲线）、二维图形（等值线），并生成图形文件（三种格式），学习用gv32.exe查看并转换gmf图形文件。reinitopen d:gradsdata.ctlenable print d:gradshomeworkex01_1.gmfset gxout lineset x 1 37set y 1set z 1set t 1d uprintdisable print*printim d:gradshomeworkex01_1.png; reinitopen d:gradsdata.ctlenable print d:gradshomeworkex01_2.gmfset gxout contourset x 1 37set y 1 17set z 1set t 1d uprintdisable print*printim d:gradshomeworkex01_2.png;2）熟悉三类坐标的用法，掌握查询命令Q，特别掌握其坐标转换命令的用法，并在grads命令提示符下运行 “Q w2xy 120 20”,看看其转换结果。 实习三 1）掌握linebarcontourshaded形式的绘图；学会fwrite的用法。注：line与contour参考实习一 reinitopen d:gradsdata.ctlenable print d:gradshomeworkex03_1.gmfset gxout barset x 1 37set y 1set z 1set t 1d uprintdisable print*printim d:gradshomeworkex03_1.png; reinitopen d:gradsdata.ctlenable print d:gradshomeworkex03_2.gmfset gxout shadedset x 1 37set y 1 17set z 1set t 1d uprintdisable print*printim d:gradshomeworkex03_2.png; 2）掌握阴影图和等值线一起绘图时的顺序；reinitopen d:gradsdata.ctlenable print d:gradshomeworkex03_2.gmfset gxout shadedset x 1 37set y 1 17set z 1set t 1d uset gxout contourd uprintdisable print*printim d:gradshomeworkex03_2.png; 3）自定义绘制一张linebar图（同样资料叠加）。reinitopen d:gradsdata.ctlenable print d:gradshomeworkex03_4.gmfset gxout barset x 1 37set y 1set z 1set t 1d uset gxout lined uprintdisable print*printim d:gradshomeworkex03_4.png; 自定义绘制一张contourshaded图（同样资料叠加）注：参考上一问。 4）从原始资料中写出20-30N，80-100E范围内资料，并配以数据描述文件，简单绘图验证正确与否。 reinitopen d:gradsdata.ctlset fwrite d:gradshomeworkex03_5.grdset gxout fwriteset lon 80 100set lat 20 30set t 1 48d udisable fwrite;相应的ctl文件dset d:gradshomeworkex03_5.grdundef -9.99E+33title NCEP/NCAR REANALYSIS PROJECTxdef 5 linear 80 2.500ydef 5 linear 20 2.500zdef 2 levels 850 200tdef 48 linear JAN1982 1movars 4U 2 99 u wind (m/s)V 2 99 v wind (m/s)H 1 99 H500TSFC 1 99 TSFC DATAendvars reinitopen d:gradshomeworkex03_5.ctlset gxout shadedset lat 20 30set lon 80 100set t 1d uset gxout contourd uprintim d:gradshomeworkex03_5.png;实习四1）求出500hPa高度场的气候场、距平场，并给出图形。reinitopen d:gradsdata.ctlset t 1 12have=ave(h,t+0,t=48,12)modify have seasonalset t 1 48hano=h-haved hanoprintim d:gradshomeworkex04_1.png;2）求出高度场 500hPa、1982-1985年逐月、沿100-120E的平均值，并画出纬度-时间剖面。（掌握ave函数用法）reinitopen d:gradsdata.ctlset lon 110tt=1while(tt=48)aa=ave(u(t=tt),lon=100,lon=120)d aactt=tt+1endwhile;reinitopen d:gradshomeworkex04_2.ctlenable print d:gradshomeworkex04_2.gmfset t 1 48set xyrev ond aaprintdisable print;3） 掌握aave函数用法，熟悉其它函数用法。 求区域（60-80E，20-30N）内200hPa纬向风场的区域平均值。 reinit open d:gradsdata.ctl set x 1 set y 1 set t 1 48 aa=aave(u(lev=200),lon=60,lon=80,lat=20,lat=30)d aa;实习五画出1982年12月地表海表温度场，大于270C打阴影，其余画等直线。reinitopen d:gradsdata.ctlenable print d:gradshomeworkex05.gmfset gxout shadedset cmin 300set x 1 37set y 1 17set z 1set t 12d TSFCset gxout contourd TSFCprintdisable print;实习六画出u850的沿100E的时间-纬度剖面，并在高中心处标黄色字符“G”。reinitopen d:gradsdata.ctlset lon 100set lat 0 40set t 1 48set z 1set xyrev ond uq w2xy AUG1982 12x1=subwrd(result,3)y1=subwrd(result,6)set string 7 c 8 0set strsiz 0.2draw string x1 y1 Gq w2xy JUN1984 10x2=subwrd(result,3)y2=subwrd(result,6)set string 7 c 8 0set strsiz 0.2draw string x2 y2 Gq w2xy AUG1984 10x3=subwrd(result,3)y3=subwrd(result,6)set string 7 c 8 0set strsiz 0.2draw string x3 y3 Gq w2xy JUN1985 11x4=subwrd(result,3)y4=subwrd(result,6)set string 7 c 8 0set strsiz 0.2draw string x4 y4 Gprintim d:gradshomeworkex06.png;实习七画任一张等值线图,要求：1）在其上实现0线加粗（或者指定等值线加粗）；2）等值线上标值调大；3）坐标标值调大。4）给出加上标题。X轴和y轴加说明。reinitopen d:gradsdata.ctlset parea 1 10 0 8.5set lon 60 150set lat 0 40set time Dec1982set gxout shadedset xlopts 1 2 0.15set ylopts 1 2 0.15d tsfcset gxout contourset c