grads垂直积分函数vint计算~

霍森布魯茨

L:常数；
g：常数；
u：u风（变量）；
q：比湿（变量）；
Q为潜热，这个公式是想算潜热的整层输送。。。。并画图


高手指点哈~谢谢啦





照着写了下，

1. 'reinit'
   
2. 'open d:\cx\q.ctl'
   
3. 'open d:\cx\u.ctl'
   
4. 'open d:\cx\v.ctl'
   
5. 'open d:\cx\ps.ctl'
   
6. 'set t 1'
   
7. 'set lev 850'
   
8. 'set lon 75 115.5'
   
9. 'set lat 15 36'
   
10. 'set grid off'
    
11. 'set xlint 5'
    
12. 'set ylint 5'
    
13. 'set grads off'
    
14. 'set map 1 1 6'
    
15. 'set mpdset cnworld'
    
16. 'draw map'
    
17. 'define a=vint(z.4/100,u.2*q.1*2510400/9.8,200)'
    
18. 'define b=vint(z.4/100,v.3*q.1*2510400/9.8,200)'
    
19. 'define div=hdivg(a,b)*100000'
    
20. 'set lon 75 115.5'
    
21. 'set lat 15 36'
    
22. 'set rgb 31 236 236 236'
    
23. 'set rgb 32 215 215 215'
    
24. 'set rgb 33 194 194 194'
    
25. 'set rgb 34 183 183 183'
    
26. 'set rgb 35 160 160 160'
    
27. 'set rgb 36 137 137 137'
    
28. 'set rgb 37 112 112 112'
    
29. 'set rgb 38 100 100 100'
    
30. 'set rgb 39 85 85 85'
    
31. 'set rgb 40 70 70 70'
    
32. 'set rgb 41 54 54 54'
    
33. 'set rgb 42 40 40 40'
    
34. 'set gxout shaded'
    
35. 'set cthick 8'
    
36. 'set clevs -6 -5 -4 -3 -2 -1 0'
    
37. 'set ccols 41 38 37 36 35 34 33 31 0'
    
38. *'set ccols 9 14 4 11 5 13 3 7'
    
39. *'set cmax 0'
    
40. 'd smth9(smth9(div))'
    
41. 'cbarn 1 2'
    
42. 'set gxout vector'
    
43. 'set arrscl 0.8 60'
    
44. 'd a;b'
    
45. 'printim e:\qrss.png white'
    
46. 'reinit'

'reinit'<br /> 'open d:\cx\q.ctl'<br /> 'open d:\cx\u.ctl'<br /> 'open d:\cx\v.ctl'<br /> 'open d:\cx\ps.ctl'<br /> 'set t 1'<br /> 'set lev 850'<br /> 'set lon 75 115.5'<br /> 'set lat 15 36'<br /> 'set grid off'<br /> 'set xlint 5'<br /> 'set ylint 5'<br /> 'set grads off'<br /> 'set map 1 1 6'<br /> 'set mpdset cnworld'<br /> 'draw map'<br /> 'define a=vint(z.4/100,u.2*q.1*2510400/9.8,200)'<br /> 'define b=vint(z.4/100,v.3*q.1*2510400/9.8,200)'<br /> 'define div=hdivg(a,b)*100000'<br /> 'set lon 75 115.5'<br /> 'set lat 15 36'<br /> 'set rgb 31 236 236 236'<br /> 'set rgb 32 215 215 215'<br /> 'set rgb 33 194 194 194'<br /> 'set rgb 34 183 183 183'<br /> 'set rgb 35 160 160 160'<br /> 'set rgb 36 137 137 137'<br /> 'set rgb 37 112 112 112'<br /> 'set rgb 38 100 100 100'<br /> 'set rgb 39 85 85 85'<br /> 'set rgb 40 70 70 70'<br /> 'set rgb 41 54 54 54'<br /> 'set rgb 42 40 40 40'<br /> 'set gxout shaded'<br /> 'set cthick 8'<br /> 'set clevs -6 -5 -4 -3 -2 -1 0'<br /> 'set ccols 41 38 37 36 35 34 33 31 0'<br /> *'set ccols 9 14 4 11 5 13 3 7'<br /> *'set cmax 0'<br /> 'd smth9(smth9(div))'<br /> 'cbarn 1 2'<br /> 'set gxout vector'<br /> 'set arrscl 0.8 60'<br /> 'd a;b'<br /> 'printim e:\qrss.png white'<br /> 'reinit'

river

霍森布魯茨 发表于 2012-12-21 12:58

                               
登录/注册后可看大图

是在联网状态下的啊。。。

那就是环境变量的问题。grads1.9以上的版本会自己添加环境变量，如果是win7的话就比较麻烦，具体的环境变量不太知道怎么设置。我上次就是grads突然闪退，怎么重装都不行。按照网上XP系统的设置了环境变量之后主程序倒是可以打开用了，可是运行gs的时候其他的东西有识别不了，比如调用的地图啊、脚本文件什么的，找不到设置这些文件的环境变量的设置方法。最后只好重装系统才解决了这个问题希望楼主不要这样了

river

坛子里有计算整层水汽通量的gs，楼主可以仿照做一下啊

mofangbao

这要怎么指点，你说你把程序大概写出来了，有些小疑问说指点还可以，这个有点大了

霍森布魯茨

river 发表于 2012-12-17 21:03

                               
登录/注册后可看大图

坛子里有计算整层水汽通量的gs，楼主可以仿照做一下啊

照着做了下，一运行grads直接闪退啊~

霍森布魯茨

mofangbao 发表于 2012-12-18 08:50

                               
登录/注册后可看大图

这要怎么指点，你说你把程序大概写出来了，有些小疑问说指点还可以，这个有点大了

照着写了下gs。。。grads闪退啊

river

霍森布魯茨 发表于 2012-12-18 13:44

                               
登录/注册后可看大图

照着做了下，一运行grads直接闪退啊~

那是grads出问题了，联网的情况下还闪退的话应该是环境变量的问题

霍森布魯茨

river 发表于 2012-12-18 21:09

                               
登录/注册后可看大图

那是grads出问题了，联网的情况下还闪退的话应该是环境变量的问题

是在联网状态下的啊。。。

霍森布魯茨

river 发表于 2012-12-21 18:42

                               
登录/注册后可看大图

那就是环境变量的问题。grads1.9以上的版本会自己添加环境变量，如果是win7的话就比较麻烦，具体的环境变 ...

亲，好几天上论坛了，工作搁置了。我之前一直在检查我的数据和程序问题呢，。。。ps说是地面压强，就用高度场就可以了吧？

river

霍森布魯茨 发表于 2012-12-24 15:33

                               
登录/注册后可看大图

亲，好几天上论坛了，工作搁置了。我之前一直在检查我的数据和程序问题呢，。。。ps说是地面压强，就用高 ...

计算整层的潜热通量肯定得从地面开始积分啊，要不怎么是整层呢。而且那个需要的是气压的值，不是等压面丫