wrfoutput剖面地形疊加問題

ml11178

各位先進、大大們，大家好:

小弟目前正在學習繪製任意兩點間的剖面，目前已經有初步的成果，

但是在繪製地形疊加上面遇到相當大的瓶頸，附上初步成果以及gs檔，希望有高人能不吝指導






這個gs檔是試過板上很多方法失敗後，想說先單獨畫出地形的剖面看看，但還是不行....

gs檔裡面，如果我將前面的*號打開，後面關掉會報錯，報錯圖如上:

那如果我將前面關掉，後面那段打開，不會報錯，但是也沒有任何東西畫出來...

1. 'reinit'
   
2. 'open tsdata/201208120300.ctl'
   
3. 'set grads off'
   
4. 'set grid off'
   
5. 'set zlog on'
   
6. 'set csmooth on'
   
7. lon1 = 121.407994
   
8. lon2 = 121.568592
   
9. lat1 = 25.244876
   
10. lat2 = 24.852246
    
11. *'define alfa=atan2('lat2-lat1','lon2-lon1')'
    
12. *'set x 1'
    
13. *'set y 1'
    
14. *'set z 1'
    
15. *lon = lon1
    
16. *'collect 1 free'
    
17. *while (lon <= lon2)
    
18. * lat = lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
    
19. * 'collect 1 gr2stn(hgt,'lon','lat')'
    
20. * lon = lon + 0.0035
    
21. *endwhile
    
22. *'set lon 121.407994 121.568592'
    
23. *'d coll2gr(1,-u)'
    
24. 
    
25. *---------------------------------------------
    
26. ln=lon1
    
27. i=1
    
28. while(ln<=lon2)
    
29. *------read terrain height--------------------
    
30. lt=lat1+(lat2-lat1)*(ln-lon1)/(lon2-lon1)
    
31. 'set lon 'ln
    
32. 'set lat 'lt
    
33. 'set z 1'
    
34. 'd HGT*0.001'
    
35. res1=sublin(result,2)
    
36. top=subwrd(res1,4)
    
37. 'q w2xy 'ln' 0.04858'
    
38. res2=sublin(result,1)
    
39. tx=subwrd(res2,3)
    
40. ty=subwrd(res2,6)
    
41. 'q w2xy 'ln' 'top
    
42. res3=sublin(result,1)
    
43. mt=subwrd(res3,6)
    
44. if(mt>ty)
    
45. 'set line 1 1 6'
    
46. 'draw line 'tx' 'ty' 'tx' 'mt
    
47. endif
    
48. ln=ln+0.0035
    
49. i=i+1
    
50. endwhile
    
51. *---------------------------------------------------------------

复制代码