怎么画T-logp图？

from1214

资料是NCEP2.5*2.5的日均数据，时间是冬季的，用 grads怎么画 T-logp图（只需画层结和状态曲线）？我看了自带的gs文件有点不太懂，请各位大侠赐教啊！

topmad

传说中的谁原来分享过一个画图的代码你可以去看看

尽头的尽头

多用论坛搜索功能

流觞曲水

哪一年月日哪个时次哪个站点的？

from1214

1. 'reinit'
   
2. 'sdfopen D:\nc\air.1984.nc'
   
3. 'sdfopen D:\nc\rhum.1984.nc'
   
4. 'sdfopen D:\nc\uwnd.1984.nc'
   
5. 'sdfopen D:\nc\vwnd.1984.nc'
   
6. 'enable print e:\2013gs\ndgmf\TLOGP.gmf'
   
7. 
   
8. 'set lon 107.5'
   
9. 'set lat 27.5'
   
10. *'set t 9'
    
11. 'set lev 1000 100'
    
12. 'define tt=ave(air.1-273.16,t=16,t=37)'
    
13. 
    
14. 
    
15. 'rh=ave(rhum.2,t=16,t=37)'
    
16. 'td=tt-((14.55+0.114*tt)*(1-0.01*rh) + pow((2.5+0.007*tt)*(1-0.01*rh),3) + (15.9+0.37*tt)*pow((1-0.01*rh),14))'
    
17. 
    
18. 'define u=ave(uwnd.3,t=16,t=37)'
    
19. 'define v=ave(vwnd.4,t=16,t=37)'
    
20. 'define uv=mag(u,v)'
    
21. 'define pi=3.14159'
    
22. if (u>=0.)
    
23. if (v>=0.)
    
24. 'define dir=180+atan2(u,v)*180/pi'
    
25. endif
    
26. endif
    
27. if (v<=0.)
    
28. if (u>=0.)
    
29. 'define dir=270+atan2(u,v)*180/pi'
    
30. endif
    
31. endif
    
32. if (u<=0.)
    
33. if (v>=0.)
    
34. 'define dir=180-atan2(u,v)*180/pi'
    
35. endif
    
36. endif
    
37. if (v<=0.)
    
38. if (u<=0.)
    
39. 'define dir=270-atan2(u,v)*180/pi'
    
40. endif
    
41. endif
    
42. *if(ugrdprs > 0.);'define dir=90-180/pi*atan2(vgrdprs,ugrdprs)';endif
    
43. *if(ugrdprs < 0.);'define dir=270-180/pi*atan2(vgrdprs,ugrdprs)';endif
    
44. 'rc=plotskew(tt,td,uv,dir)'
    
45. 
    
46. 
    
47. *c
    
48. *'reinit'

复制代码

我参照了一下帖子写的，代码省略了plotskew.gs里面的内容，但是运行出错，不知道原因出在哪，求助啊

一线天

from1214 发表于 2013-6-9 20:18
我参照了一下帖子写的，代码省略了plotskew.gs里面的内容，但是运行出错，不知道原因出在哪，求助啊

可以参考这个，我用这个画出了图http://bbs.06climate.com/forum.php?mod=viewthread&tid=34626

一线天

from1214 发表于 2013-6-9 20:18
我参照了一下帖子写的，代码省略了plotskew.gs里面的内容，但是运行出错，不知道原因出在哪，求助啊

可以参考这个，我用这个画出了图http://bbs.06climate.com/forum.php?mod=viewthread&tid=34626