Grads中IF语句的问题

Isakting

想实现一个判别的过程，在某一区域风速＞13.9m/s的时候绘制图像，否则不绘制图像，想用Grads直接画出来，但是IF语句那里出问题了，是不是我对IF的用法还不对，IF里面只能用 i=1.0 这种，而不能用 i>=1.0 呢？我的程序运行的时候还是没有筛选，还是一张一张的画图……抓狂

1. year=1990
   
2. while(year<=1990)
   
3. 
   
4. 'reinit'
   
5. 'clear'
   
6. 
   
7. 'sdfopen G:\ECwinddata\wind\'year'.nc'
   
8. 
   
9. it=1
   
10. 
    
11. while(it<=1464)
    
12. 
    
13. 
    
14. lo=105
    
15. la=0
    
16. 
    
17. while(lo<=130)
    
18. while(la<=42.5)
    
19. 
    
20. 'set lon 'lo''
    
21. 'set lat 'la''
    
22. 'define lv=mag(u10,v10)'
    
23. 'd lv'
    
24. 
    
25. 
    
26. if(lv>=13.9)then
    
27. 
    
28. 'set lon 100 150'
    
29. 'set lat 0 50'
    
30. 'set mpdset hires'
    
31. 'set poli on'
    
32. 'set map 1 1 8'
    
33. 'set t 'it''
    
34. 'q time'
    
35. time=subwrd(result,3)
    
36. 'define v=mag(u10,v10)'
    
37. 
    
38. 'H:\Data\Wind\define_colors.gs'
    
39. 'set grads off'
    
40. 'set gxout shaded'
    
41. 'set clevs 8.0 10.8 13.9 17.2 20.8 24.5 28.5 32.7 37'
    
42. 'set ccols 0 21 23 25 27 28 29 64 66 68'
    
43. 'set csmooth on'
    
44. 'set parea 0.5 7.5 0.5 10'
    
45. 'd v'
    
46. 
    
47. 'H:\Data\Wind\cbarn.gs'
    
48. 
    
49. 'set grads off'
    
50. 'set gxout contour'
    
51. 'set clevs 8.0 10.8 13.9 17.2 20.8 24.5 28.5 32.7 37'
    
52. 'set ccols 1'
    
53. 'set csmooth on'
    
54. 'd v'
    
55. 
    
56. 'set gxout barb'
    
57. 'set ccolor 1'
    
58. 'd skip(u10*2.5,10);skip(v10*2.5,10)'
    
59. 'draw title 'time''
    
60. 'printim I:\wind\1990\'year'winpic'it'.bmp white'
    
61. 'clear'
    
62. 
    
63. la=50
    
64. lo=150
    
65. 
    
66. else
    
67. la=la+0.25
    
68. lo=lo+0.25
    
69. 
    
70. endif
    
71. 
    
72. endwhile
    
73. endwhile
    
74. 
    
75. it=it+1
    
76. endwhile
    
77. 
    
78. year=year+1
    
79. endwhile
    
80. 
    

复制代码

tulalang

你也可以设置输出形式后用'set cmin 13.9'来屏蔽掉小于13.9的区域

小其其格

grads中if这样用貌似没效果的。我也试过好几次。建议你还是用楼上的'set cmin 13.9'