关于在满足某一条件的区域内画符号的问题

勇敢的悲伤

这个问题困扰我了很久，终于在 river 和mofangbao 的帮助下画出了自己想要的图形。我也来分享下我的经验，主要是参照mofangbao的《GrADS绘制站点分布图-升级版》，用到了subwrd和sublin函数。

图形说明：水汽含量（比湿）正、负异常的区域分别用+、-表示。
正异常用比湿大于1e-05表示，负异常用比湿小于-1e-05表示。
经纬度：lon 30°E-180°E，lat 30°S-60°N

llon=30
while(llon<=180)
'set lon 'llon
llat=-30
while(llat<=60)
'set lat 'llat
'set t 1'
'set z 1'
'q w2xy 'llon' 'llat' '
x1=subwrd(result,3)
y1=subwrd(result,6)
'set gxout print'          *将每个经纬度对应的比湿赋值给b                           
'd shum.8'
res=result
data=sublin(res,2)
b=subwrd(data,1)
if(b>1e-05)                  比湿大于1e-05的格点上画+            
'set string 8 c 9 0'
'draw string 'x1' 'y1' +'
endif
if(b<-1e-05)                    比湿小于-1e-05的格点上画-
'set string 5 c 9 0'
'draw string 'x1' 'y1' -'
endif
llat=llat+2.5
endwhile
llon=llon+2.5
endwhile

river

呵呵，这么快就出来啊，感谢分享。配张图就更好了

mofangbao

不错啦，谢谢分享~

ric_bob

谢谢楼主分享~

qxtlyf

谢谢楼主分享~