如何在 GrADS 中绘制正方形?

How to draw a square in GrADS?

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

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

我假设你要绘制的正方形的边界是基于lat/lon坐标,左下角在(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