grads考试例题

'set gxout line' 'set ccolor 2' ’set cmark 3' ’set cstyle 1' 'set cthick 6’ 'd h’

’draw title height'

'printim d：\\abc.png white' ;

② 利用所提供的数据文件，计算500hPa高度场的距平场(参考P114例2），并绘制1983年7月的距平场图(参考P114例3），要求〈0的值标出阴影,0线加粗，最终将图形保存。所有命令编写于.gs文件中。

’reinit’

'open d:\\data。ctl' ’set t 1 12’

'have=ave（h，t+0，t=48,12）' 'modify have seasonal’ 'set t 1 48'

’hano=h—have’ 'set time jul1983' 'set gxout shaded’ ’set cmax 0’ 'd hano’

’set gxout contour' ’d hano’

'set clevs 0' ’set cthick 10' ’d hano’

'draw title 500hpa anomalous fields for jul1983’ 'draw xlab LON'；’draw ylab LAT' ’printim d:\\abc。png white' ;

1)恰当利用所提供的数据资料，绘制一张北半球地图底图，即图形中只显示地图信息，不含其它形式的数据信息。 ’reinit’

'open d：\\model。le.ctl’ 'set cmin 9999’ ’set cmax 99999' ’set lat 0 90’ ’set lon -180 180' 'd ts'

’printim d：\\a。png white’ ;

2)绘制一张第2时刻温度垂直剖面图（纬度固定于30ºN,经度范围80ºW—0º，层次由第1层变化至第7层），要求将温度正值区域用阴影显示,负值区域用等值线显示,并加粗0值等值线,同时Y轴标注为“Z/hPa”，标题标注为“Temperature/ºC”。 'reinit'

’open d:\\model。le.ctl’ ’set lat 30' 'set lon -80 0’ 'set lev 1000 100’ 'set t 2’ 'tm=t—273'

'set gxout shaded’ 'set cmin 0' 'd tm’

’set gxout contour' 'd tm'

’set clevs 0' ’set cthick 10' ’d tm’

’draw title Temperature/`3。`1C' 'draw ylab Z/hPa'

'printim d：\\b。png white’ ;

3）利用fwrite命令将第2天西北半球500hPa高度场资料提取存放至500hPa.grd文件中。 ’reinit’

’open d:\\model.le。ctl' ’set fwrite d：\":{"h":17.534,"w":12.24,"x":565.844,"y":1023.196,"z":312},"ps":null,"t":"word。grd’ 'set gxout fwrite' 'set lat 0 90’

'set lon —180 0' 'set lev 500' ’set t 2' 'd u'

'disable fwrite'

① 利用所提供的数据文件,求区域（60-80E,20—30N）内200hPa纬向风场的区域平均值，并给出时间演变（1982.1—1985。12)图，所有命令编写于。gs文件中。 'reinit’

'open d:\\data。ctl' 'set lon 60’ 'set lat 20’ ’set lev 200’ ’set t 1 48’

’d tloop(aave(u，lon=60,lon=80，lat=20，lat=30）)'

'printim d:\\a。png white’ ;

② 利用所提供的数据文件,计算高度场 500hPa、1982—1985年逐月、沿100-120E的平均值，并画出纬度—时间剖面。所有命令编写于.gs文件中。 ’reinit'

’open d：\\data。ctl' ’set lon 100’ ’set lev 850' 'set t 1 48’

'd ave（h,t=1,t=48)'

'printim d：\\b。png white’ ；

’reinit’

’open d：\\data.ctl’ ’set lat 0 40’ ’set lon 60 150’ 'set lev 850'

’enable print d：\\height.gmf' i=1

while（i〈=48) 'set t ’i 'q time'

x1=subwrd(result，3） 'd h'

’draw title ’x1'’ ’print’ ’c'

i=i+1 endwhile

’disable print'