如何在GrADS中画一个正方形？

Question

我正在寻找一个允许在GrADS中绘制正方形的命令行。这是一个例子

我想获得一个命令行，它允许我绘制一个像我在图片示例中显示的那样的正方形。

Answer 1

我假设您要绘制的正方形的边界基于纬度/经度坐标，左下角位于（4N，74W），右上角位于（8N，72W）。以下脚本片段应该可以使用，但您必须先绘制绘图，以便设置世界坐标与页面上X，Y位置之间的缩放。

'q w2xy -74 4'  
xpos1=subwrd(result,3)
ypos1=subwrd(result,6)
'q w2xy -72 8'
xpos2=subwrd(result,3)
ypos2=subwrd(result,6)
'set line 2 1 6'  '* a thick, red line
'draw rec 'xpos1' 'ypos1' 'xpos2' 'ypos2

如果您以交互方式使用GrADS，这就是它的样子：

ga-> q w2xy  -74 4   
X = 4.77083  Y = 3.08333
ga-> q w2xy  -72 8
X = 5.74306  Y = 5.41667
ga-> set line 2 1 6
SET LINE values:  color = 2  style = 1  thickness = 6
ga-> draw rec 4.77083  3.08333 5.74306  5.41667