实习8
1.实习要求
绘制第二个时次850hPa位势高度,范围lon:-80到-20, lat:15到60,并通过grads函数客观标记北大西洋低压在第一个时次和第二个时次的中心(用字符D标记):低压中心位置范围在lat:35到55,lon:-55到-30。并将中心用直线相连(draw line)。 实习步骤:
1打开l
2设置经纬度、时间、高度范围
3画第一个时次的位势高度(hgt)
4判断低压中心的位置(使用上面说的低压中心大致范围)
使用min和minloc函数
'd minloc(min(变量名,lat=??,lat=??),lon=??,lon=??)'
这个会输出很多行内容,全部都在result里面,x方向的低压中心格点编号在最后一行,设计一个循环,读取最后一行的低压中心格点编号,使用subwrd和sublin函数,中间需要判断if。(自己在ga->后面试试再写gs) 4.2计算y方向的低压中心格点编号,方法同上,display语句改为
'd minloc(min(变量名, lon=??,lon=??), lat=??,lat=??)'
4.3使用q gr2xy得到低压中心在图中的坐标(画布上的坐标)并用subwrd和sublin(如果用得到,自己在ga->后面试试再写gs)通过赋值语句保存这个坐标 5清除上面的图(临时绘制)
6画第二个时次的位势高度(hgt)
7判断这个时次的低压中心的位置,方法同上
8通过draw string在两个时次的低压中心写D,字型Roman,大小0.2
9通过draw line画线,连接两个D
10保存图片
2.代码
'reinit'
'open d:\l'
'set lon -80 -20'
'set lat 15 60'
'set lev 850'
'set t 1'
'd hgt'
'd minloc(min(hgt(t=1),lat=35,lat=55),lon=-55,lon=-30)'
j=0
i=1
while(j=0)
str1=sublin(result,i)
str2=subwrd(str1,1)
if(str2='Result')
xpos=subwrd(str1,4)
j=1
endif
i=i+1
endwhile
'd minloc(min(hgt(t=1),lon=-55,lon=-30),lat=35,lat=55)'
j=0
i=1
while(j=0)
北大语料库 str1=sublin(result,i)
str2=subwrd(str1,1)
if(str2='Result')
ypos=subwrd(str1,4)
j=1
endif
i=i+1
endwhile
'q gr2xy 'xpos' 'ypos''
x1=subwrd(result,3)
y1=subwrd(result,6)
'set font 1'
'set strsiz 0.2'
'draw string 'x1' 'y1' D'
'c'
'set t 2'
'd hgt'
'd minloc(min(hgt(t=2),lat=35,lat=55),lon=-55,lon=-30)'
j=0
没有革命的理论就没有革命的运动i=1
while(j=0)
str1=sublin(result,i)
str2=subwrd(str1,1)
if(str2='Result')
xpos=subwrd(str1,4)
促销策略分析 j=1
endif
i=i+1
endwhile
'd minloc(min(hgt(t=2),lon=-55,lon=-30),lat=35,lat=55)'
j=0
i=1
while(j=0)
str1=sublin(result,i)
str2=subwrd(str1,1)
传送带设计 if(str2='Result')
ypos=subwrd(str1,4)
j=1
endif
i=i+1
清教运动
endwhile
'q gr2xy 'xpos' 'ypos''
x2=subwrd(result,3)
y2=subwrd(result,6)
'set font 1'
'set strsiz 0.2'
'draw string 'x2' 'y2' D'国家干部
'draw string 'x1' 'y1' D'
'draw line 'x1' 'y1' 'x2' 'y2''
'printim d:\dyf\a.gif'
3.运行结果