资源描述
实习8一实习要求绘制第二个时次850hPa位势高度，范围lon：-80到-20,lat：15到60,并通过grads函数客观标记北大西洋低压在第一个时次和第二个时次的中心（用字符D标记）：低压中心位置范围在lat：35到55，lon：-55到-30。并将中心用直线相连（drawline）。实习步骤：1打开model.ctl2设置经纬度、时间、高度范围3画第一个时次的位势高度（hgt）4判断低压中心的位置（使用上面说的低压中心大致范围）使用min和minloc函数4.1计算x方向的低压中心格点编号dminloc（min（变量名，lat二？，lat二？），lon二？，lon二？）这个会输出很多行内容，全部都在result里面，x方向的低压中心格点编号在最后一行，设计一个循环,读取最后一行的低压中心格点编号，使用subwrd和sublin函数，中间需要判断if。（自己在ga-后面试试再写gs）4.2计算y方向的低压中心格点编号，方法同上，display语句改为dminloc（min（变量名，lon二？,lon=？）,lat二？,lat二？）4.3使用qgr2xy得到低压中心在图中的坐标（画布上的坐标）并用subwrd和sublin（如果用得到，自己在ga-后面试试再写gs）通过赋值语句保存这个坐标5清除上面的图（临时绘制）6画第二个时次的位势高度（hgt）7判断这个时次的低压中心的位置，方法同上8通过drawstring在两个时次的低压中心写D，字型Roman，大小0.29通过drawline画线，连接两个D10保存图片二代码reinitopend:dyfmodel.ctlsetlon-80-20setlat1560setlev850sett1dhgtdminloc(min(hgt(t=1),lat=35,lat=55),lon=-55,lon=-30)j=0i=1while(j=0)str1=sublin(result,i)str2=subwrd(str1,1)if(str2=Result)xpos=subwrd(str1,4)j=1endifi=i+1endwhiledminloc(min(hgt(t=1),lon=-55,lon=-30),lat=35,lat=55)j=0i=1while(j=0)str1=sublin(result,i)str2=subwrd(str1,1)if(str2=Result)ypos=subwrd(str1,4)j=1endifi=i+1endwhileqgr2xyxposyposx1=subwrd(result,3)y1=subwrd(result,6)setfont1setstrsiz0.2drawstringx1y1Dcsett2dhgtdminloc(min(hgt(t=2),lat=35,lat=55),lon=-55,lon=-30)j=0i=1while(j=0)str1=sublin(result,i)str2=subwrd(str1,1)if(str2=Result)xpos=subwrd(str1,4)j=1endifi=i+1endwhiledminloc(min(hgt(t=2),lon=-55,lon=-30),lat=35,lat=55)j=0i=1while(j=0)str1=sublin(result,i)str2=subwrd(str1,1)if(str2=Result)ypos=subwrd(str1,4)j=1endifi=i+1endwhileqgr2xyxposyposx2=subwrd(result,3)y2=subwrd(result,6)setfont1setstrsiz0.2drawstringx2y2Ddrawstringx1y1Ddrawlinex1y1x2y2printimd:dyfa.gif三运行结果
展开阅读全文