任意两点间斜剖面做法总结【利用GRADS自带函数】

孙亦

按照大神的脚本只画了垂直速度，但是觉得剖面地形近地面处会有空白，该

1. 'reinit'
   
2. 'sdfopen E:\NCEP1\Monthly\omega.mon.mean.nc'
   
3. 'sdfopen E:\NCEP1\Monthly\uwnd.mon.mean.nc'
   
4. 'sdfopen E:\NCEP1\Monthly\vwnd.mon.mean.nc'
   
5. 'sdfopen E:\NCEP1\Monthly\hgt.mon.mean.nc'
   
6. 'open E:\NCEP1\surface_elevation\DiXing\orog0.25.ctl'
   
7. 
   
8. 'set grads off'
   
9. 'set grid off'
   
10. 'set zlog on'
    
11. *'set time may1970 sept1970'
    
12. 
    
13. 'set lev 1000 100'
    
14. 'define w6=ave(omega.1,t=6,t=818,12)'
    
15. 'define u6=ave(uwnd.2,t=6,t=818,12)'
    
16. 'define v6=ave(vwnd.3,t=6,t=818,12)'
    
17. 'define w7=ave(omega.1,t=7,t=818,12)'
    
18. 'define u7=ave(uwnd.2,t=7,t=818,12)'
    
19. 'define v7=ave(vwnd.3,t=7,t=818,12)'
    
20. 'define w8=ave(omega.1,t=8,t=818,12)'
    
21. 'define u8=ave(uwnd.2,t=8,t=818,12)'
    
22. 'define v8=ave(vwnd.3,t=8,t=818,12)'
    
23. 'define w9=ave(omega.1,t=9,t=818,12)'
    
24. 'define u9=ave(uwnd.2,t=9,t=818,12)'
    
25. 'define v9=ave(vwnd.3,t=9,t=818,12)'
    
26. 'define w=(w6+w7+w8+w9)/4'
    
27. 'define v=(v6+v7+v8+v9)/4'
    
28. 'define u=(u6+v7+u8+u9)/4'
    
29. 
    
30. 
    
31. 'define a=ave(omega.1-w,time=jun1965,time=sep1965)'
    
32. 'define b=ave(omega.1-w,time=jun1966,time=sep1966)'
    
33. 'define c=ave(omega.1-w,time=jun1972,time=sep1972)'
    
34. 'define d=ave(omega.1-w,time=jun1974,time=sep1974)'
    
35. 'define e=ave(omega.1-w,time=jun1987,time=sep1987)'
    
36. 'define f=ave(omega.1-w,time=jun2009,time=sep2009)'
    
37. 'define n=(a+b+c+d+e+f)/6'
    
38. lon1=60.0
    
39. lon2=130.0
    
40. lat1=10.0
    
41. lat2=45.0
    
42. 'define alfa=atan2('lat2-lat1','lon2-lon1')'
    
43. 'set x 1'
    
44. 'set y 1'
    
45. 'set lev 1000 100'
    
46. 
    
47. 
    
48. lon=lon1
    
49. 'collect 1 free'
    
50. 'collect 2 free'
    
51. 'collect 3 free'
    
52. 'collect 4 free'
    
53. while(lon<=lon2)
    
54. lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
    
55. 'collect 1 gr2stn(-n*100,'lon','lat')'
    
56. 'collect 3 gr2stn(orog.5(t=1),'lon','lat')'
    
57. lon=lon+1
    
58. endwhile
    
59. 'set lon 60 140'
    
60. 'set lev 1000 100'
    
61. 'set ylevs 1000 925 850 700 600 500 400 300 200 100'
    
62. 'set xlabs 10N,60E|15N,70E|20N,80E|25N,90E|30N,100E|35N,110E|40N,120E|45N,130E'
    
63. 
    
64. 'set gxout shaded'
    
65. 'set csmooth on'
    
66. 'set clab off'
    
67. 'set clevs 0'
    
68. 'set ccols 0 5 0'
    
69. 'd coll2gr(3,-u)'
    
70. 
    
71. 'set gxout stream'
    
72. 'set cint 0.2'
    
73. 'set black -0.1 0.1'
    
74. 'set clab on'
    
75. 'set dfile 1'
    
76. 'd maskout(coll2gr(1,-u),-coll2gr(3,-u))'
    
77. 
    
78. 'draw title W less along (10N,60E) to (45N,130E)'
    
79. 'cbarn 1.0 0'
    
80. 'gxprint E:\NCEP1\9.28\anypoumian5_W_less.png white'
    
81. ;

复制代码

怎么处理呢

孙亦

孙亦 发表于 2017-10-11 12:22
按照大神的脚本只画了垂直速度，但是觉得剖面地形近地面处会有空白，该
怎么处理呢

对了，不用'set gxout stream'我要画的是填色图

hwshh

地形资料很好使，谢谢

river

孙亦 发表于 2017-10-11 12:22
按照大神的脚本只画了垂直速度，但是觉得剖面地形近地面处会有空白，该
怎么处理呢

  嫌有空白，就去看我的另一篇斜剖面贴子里讲到的透明色用法代替maskout 就好了

cold_wq

我看朱禾老师的教程里面有任意3点以及3点以上剖面图，有没有人研究过？

river

cold_wq 发表于 2017-10-15 17:57
我看朱禾老师的教程里面有任意3点以及3点以上剖面图，有没有人研究过？

我倒是有那个脚本，在屏幕上单击取点，然后画垂直剖面。但是没有仔细研究过，后面有空再看看

kisscate

说的很清楚，怒赞

kisscate

说的很清楚，怒赞，顶！顶！顶！

毛毛虫106

谢谢分享，学习了很多

毛毛虫106

学到了许多，谢谢楼主分享