Grads 条件IF语句的问题

勇敢的悲伤

我想判断每个格点上变量（shum）数值大于1e-05，并且在这个格点上画+，用if语句判断，结果并没有判断，而是所有的格点都画了+，这是为什么，是IF用的不对吗？
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)
if(shum>1e-05）
raw mark 1 'x1' 'y1' 0.1'
'set line 3 1 6'
endif
llat=llat+2.5
endwhile
llon=llon+2.5
endwhile

river

grads的判断好像不能简单地等同于fortran里面的吧，你这个虽然判断了变量的值，但是这个变量没有和经纬度建立联系，也就是你这个判断语句仅仅作用在了变量上，而对于经纬度却没有限制作用，所以就都化成那样的了。给你个清风帖子的连接，去看一下，可能会有收获的。

楼主完成之后可以来发帖分享一下，谢谢

勇敢的悲伤

river 发表于 2013-10-31 12:25

                               
登录/注册后可看大图

grads的判断好像不能简单地等同于fortran里面的吧，你这个虽然判断了变量的值，但是这个变量没有和经纬度建 ...

恩恩，好的。我试一试~~

勇敢的悲伤

river 发表于 2013-10-31 12:25

                               
登录/注册后可看大图

grads的判断好像不能简单地等同于fortran里面的吧，你这个虽然判断了变量的值，但是这个变量没有和经纬度建 ...

根据那个帖子我在Grads里一步一步运行了一下，到上面第三句的时候，显示
Syntax Error: Invalid Operand，这个是为什么，result不是不用定义吗？

river

勇敢的悲伤 发表于 2013-10-31 17:30

                               
登录/注册后可看大图

根据那个帖子我在Grads里一步一步运行了一下，到上面第三句的时候，显示
Syntax Error: Invalid Operand ...

这个只能在gs里面用吧

勇敢的悲伤

river 发表于 2013-10-31 17:53

                               
登录/注册后可看大图

这个只能在gs里面用吧

哦哦，好的~

勇敢的悲伤

river 发表于 2013-10-31 17:53

                               
登录/注册后可看大图

这个只能在gs里面用吧

现在我已经画出来了大于某个数值的区域，但是小于某个数值的区域仍然出现刚开始的情况（就是每个给点都会画上符号），这个是为什么呢？用subwrb函数得出的应该是字符串吧，这个怎么能直接和数字比较呢？

llon=30
while(llon<=180)
'set lon 'llon
llat=-30
while(llat<=0)
'set lat 'llat
'set t 1'
'set z 1'
'q w2xy 'llon' 'llat' '
x1=subwrd(result,3)
y1=subwrd(result,6)
'set gxout print'
'd shum.8'
res=result
ln=sublin(res,2)
i=1
while(i<=2257)
data=sublin(res,2)
a=subwrd(data,i)    → 此时a应该是字符串，怎么能和1e-05比较呢？
if(a>1e-05)            →  如果改为if（a<-1e-05） 画出来的图就没有判断直接全部画出符号
'draw mark 1 'x1' 'y1' 0.05'
'set line 8 1 7'
endif
i=i+1
endwhile
llat=llat+2.5
endwhile
llon=llon+2.5
endwhile

river

勇敢的悲伤 发表于 2013-11-1 10:36

                               
登录/注册后可看大图

现在我已经画出来了大于某个数值的区域，但是小于某个数值的区域仍然出现刚开始的情况（就是每个给点都会 ...

那是不是你的数据没有小于那个值的呢。我也不太会弄这几个函数命令，我帮你呼叫清风吧@mofangbao

勇敢的悲伤

river 发表于 2013-11-1 10:58

                               
登录/注册后可看大图

那是不是你的数据没有小于那个值的呢。我也不太会弄这几个函数命令，我帮你呼叫清风吧@mofangbao

有的，我画过这个变量的阴影图，是有比这个值小的数的。谢谢啦。。。多亏你的帮忙，现在成功一半了~~

mofangbao

river 发表于 2013-11-1 10:58

                               
登录/注册后可看大图

那是不是你的数据没有小于那个值的呢。我也不太会弄这几个函数命令，我帮你呼叫清风吧@mofangbao

你需要给出ctl的信息。
这种错误应该比较好查找，只要你在判断之前say一下变量的结果基本就能看出问题所在了。
我不知道你这里为什么是三重循环，如果说外面两层是经纬度的循环，那里面这个是什么，为什么每个点上都要进行2257次循环，而且这么多次用的是同一个res，也就是截取出来的a结果是一样的？
还有外层经纬度循环最好用x y的形式，而不是lon，lat，除非你明白lon、lat这种需要和数据点一一对应。