求物理量垂直剖面图的绘制

gdfsy

我有物量场资料，想绘制垂直剖面图，请高手指点，谢谢！

river

你有什么资料，说清楚啊。还有就是剖面图也分好几种，你要画哪种啊？时间高度，沿某一经纬度的高度？·····

river

给你个用nc资料的时间高度剖面的例子，自己研究一下，其他的剖面也是一个意思

1. 'reinit'
   
2. 'sdfopen f:\2012\uwnd.2012.nc'
   
3. 'sdfopen f:\2012\vwnd.2012.nc'
   
4. 'sdfopen f:\2012\rhum.2012.nc'
   
5. 'sdfopen f:\2012\shum.2012.nc'
   
6. 'sdfopen f:\2012\omega.2012.nc'
   
7. 'sdfopen f:\2012\air.2012.nc'
   
8. 'set t 725 730'
   
9. *-----水汽通量散度------
   
10. 'set lev 1000 300'
    
11. 'define qdiv=hdivg(uwnd.1*shum.4*1e3/9.8,vwnd.2*shum.4*1e3/9.8)*(1e5)'
    
12. *--------GRADS设置------------
    
13. 'set grads off'
    
14. 'set grid off'
    
15. 'set csmooth on'
    
16. 'set map 15 1 9'
    
17. 'set parea 1 10.2 1.5 8'
    
18. 'set xlopts 1 6 0.14'
    
19. 'set ylopts 1 6 0.14'
    
20. 'set font 0'
    
21. 
    
22. *---------输出图形------------
    
23. 'enable print f:\2012\qdiv.gmf'
    
24. 'set lon 125'
    
25. 'set lat 35'
    
26. 'set zlog on'
    
27. 'set ylevs 1000 925 850 700 600 500 400 300 250 200 150 100'
    
28. 'set xlabs 08|12|16|20|00|04|08|12|16|20|00|04|08|12|16|20|00|04|08'
    
29. 'set gxout shaded'
    
30. 'set cint 0.3'
    
31. 'set cmax 0'
    
32. 'd qdiv'
    
33. 'run cbar.gs'
    
34. 'set ylpos 0 r'
    
35. 'set gxout contour'
    
36. 'set cint 0.3'
    
37. 'set cthick 11'
    
38. 'set clskip 1'
    
39. 'd qdiv'
    
40. 'writehz 0.4 4 气压（hPa） 1 2 1 1.3 0.6 90 11'
    
41. 'writehz 1 8.1 水汽通量散度 1 2.5 1 1.3 0.6 0 11'
    
42. 'writehz 0.7 1.0 6月26日 1 2 1 1.3 0.6 0 7'
    
43. 'writehz 2.8 1.0 6月27日 1 2 1 1.3 0.6 0 7'
    
44. 'writehz 5.9 1.0 6月28日 1 2 1 1.3 0.6 0 7'
    
45. 'writehz 9.1 1.0 6月29日 1 2 1 1.3 0.6 0 7'
    
46. 
    
47. 'print'
    
48. 'disable print'
    
49. 'printim f:/2012/qdiv.gif white'
    
50. ;

复制代码

gdfsy

river 发表于 2012-10-3 22:12

                               
登录/注册后可看大图

你有什么资料，说清楚啊。还有就是剖面图也分好几种，你要画哪种啊？时间高度，沿某一经纬度的高度？··· ...

thetse-假相当位温、vor-涡度、div-散度，想画东经100度的高度剖面和北纬28度的高度剖面图，谢谢！

abd

其实，原理跟画二维平面图有什么本质上的区别吗？

做个霸气的木头

垂直剖面图顾名思义是随垂直高度变化的图，另外还有一维是变化的，按照你的是纬度变化（经度固定）和经度变化（纬度固定），如果是grads可以读的数据应该用grads画比较方便。

river

gdfsy 发表于 2012-10-3 22:54

                               
登录/注册后可看大图

thetse-假相当位温、vor-涡度、div-散度，想画东经100度的高度剖面和北纬28度的高度剖面图，谢谢！

是GRADS可以直接读的资料吗？首先要有CTL文件，然后用GS控制一下就好了。你要画的东经100度的话就是说首先要固定经度set lon 180 ，然后让高度和纬度变化就行了set lev 1000 100 , set lat  20 60,最后d thetse什么的就行了，其他的同理。你要自己试着编gs文件啊，要不然怎么都不会成功的

易小凯

gdfsy 发表于 2012-10-3 22:54

                               
登录/注册后可看大图

thetse-假相当位温、vor-涡度、div-散度，想画东经100度的高度剖面和北纬28度的高度剖面图，谢谢！

用grads画吧

gdfsy

river 发表于 2012-10-4 18:47

                               
登录/注册后可看大图

是GRADS可以直接读的资料吗？首先要有CTL文件，然后用GS控制一下就好了。你要画的东经100度的话就是说首先 ...

我的数据是micaps第四类数据（实况物理量场数据），用grads画的话，要先把我的数据转换成二进制无格式数据，这个我不会，请赐教，谢谢！！

gdfsy

做个霸气的木头 发表于 2012-10-4 08:50

                               
登录/注册后可看大图

垂直剖面图顾名思义是随垂直高度变化的图，另外还有一维是变化的，按照你的是纬度变化（经度固定）和经度变 ...

我的数据是micaps第四类数据（实况物理量场数据），用grads画的话，要先把我的数据转换成二进制无格式数据，这个我不会，请版主赐教，谢谢！！