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

meimeideqiqi

{:eb302:}{:eb302:}

cibaouc

楼主你好！地形垂直方向的插值，怎么弄呢

river

cibaouc 发表于 2017-11-26 21:44
楼主你好！地形垂直方向的插值，怎么弄呢

不会，你自己研究一下吧。但我已经提供了不需要垂直插值直接叠加地形的方法，自己去看我另一篇帖子（这个贴子里有链接）

yuhuiqing6

第一次接触，求问大神们set x 14 16 是什么意思呀

river

yuhuiqing6 发表于 2017-11-27 17:18
第一次接触，求问大神们set x 14 16 是什么意思呀

这个是设置经度范围。建议你还是先去看看grads的基础知识

fujianbeige

好~~~~~~~~~~~~~~~~~~

黄叁3

大神你好，我是在画高度场斜剖面时出现了问题，以下是我的脚本：

1. 'reinit'
   
2. 'set xlopts 1 5 0.13'
   
3. 'set ylopts 1 5 0.16'
   
4. 'set clopts -1 -1 0.1'
   
5. 'open f:/a2017/wrfout2/hs612xd01.ctl'
   
6. 'set parea 1 7 3 8'
   
7. 'set mproj scaled'
   
8. 'set grads off'
   
9. 'set grid off'
   
10. 'set mpdset hires'
    
11. 'set map 1 1 7'
    
12. 'set font 1'
    
13. 'set cthick 6'
    
14. **********************
    
15. 'set x 1'
    
16. 'set y 1'
    
17. lon1=120
    
18. lon2=121
    
19. lat1=37
    
20. lat2=36
    
21. 'set z 1 10'
    
22. 'set t 7'
    
23. lon=lon1
    
24. 'collect 1 free'
    
25. 'collect 2 free'
    
26. 'collect 3 free'
    
27. 'collect 4 free'
    
28. 'collect 5 free'
    
29. 'collect 6 free'
    
30. 'collect 7 free'
    
31. while(lon<=lon2)
    
32. lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
    
33. 'collect 1 gr2stn(tc,'lon','lat')'
    
34. 'collect 2 gr2stn(w*100,'lon','lat')'
    
35. 'collect 3 gr2stn(u,'lon','lat')'
    
36. 'collect 4 gr2stn(v,'lon','lat')'
    
37. 'collect 5 gr2stn(w,'lon','lat')'
    
38. 'collect 6 gr2stn(tsk-273.15,'lon','lat')'
    
39. 'collect 7 gr2stn(pblh+hgt,'lon','lat')'
    
40. lon=lon+0.045
    
41. endwhile
    
42. **********************
    
43. 'set lon 120 121'
    
44. 'set z 1 10'
    
45. 'set t 7'
    
46. 'set gxout shaded'
    
47. 'define_colors'
    
48. 'set clevs  -0.007 -0.006 -0.005 -0.004 -0.003 -0.002 -0.001 0 0.001 0.002 0.003 0.004'
    
49. 'set ccols  59 58 57 56 55 54 53 52 22 24 25 26 27'
    
50. 'd coll2gr(5,10)'
    
51. 'cbar'
    
52. 
    
53. 'gxprint f:/a2017/ctl/xiepoumian702.png white x1000 y800'
    
54. 'reinit'

复制代码

然后就是每次运行画图来的图都不太一样，而且后面我画温度场每次数值也都不一样，我想问这是lon=lon+分辨率这步出了问题么

黄叁3

画温度场我就写了这两句
'set gxout contour'
'd coll2gr(1,10)'
结果只出了一根0线，想知道原因谢谢解答

river

黄叁3 发表于 2017-12-17 17:14
大神你好，我是在画高度场斜剖面时出现了问题，以下是我的脚本：

这个我也不清楚，但从两张图我倒是没看出有那么不一样。至于经度的增量一般选和资料分辨率比较接近的值就可以，这个值的大小不会改变出图的大趋势，只是改变线条的疏密什么的。

黄叁3

river 发表于 2017-12-17 22:17
这个我也不清楚，但从两张图我倒是没看出有那么不一样。至于经度的增量一般选和资料分辨率比较接近的值就 ...

感谢回复，还是想问为什么温度只会输出0线，我的是wrfout的资料