关于 水汽通量散度，各位帮忙看看

bolv6666

1. 
   
2. 'reinit'
   
3. 'open e:/sr/diagnose/fnl20120729/fnl20120729.ctl'
   
4. 'enable print E:\SR\Diagnose\FNL20120729\mpvnew\PIC\QFLUX&QFDIV.gmf'
   
5. 'set mpdset cnworld'
   
6. 'set mpdraw on'
   
7. 'set poli oN'
   
8. 'set map 1 1 0.5'
   
9. 'set csmooth on'
   
10. 'set grid off'
    
11. 'set grads off'
    
12. 'set xlopts 1 7 0.15'
    
13. 'set ylopts 1 7 0.15'
    
14. 
    
15. 'set lon 100 130'
    
16. 'set lat 25 45'
    
17. 'set lev 850'
    
18. 'set t 12'
    
19. 
    
20. 'define t=TMPprs'
    
21. 'define rh=RHprs'
    
22. 'define u=UGRDprs'
    
23. 'define v=VGRDprs'
    
24. 'define prs=lev'
    
25. 'define g=9.8'
    
26. #饱和水汽压的计算，单位：hPa
    
27. 'define es=6.112*exp(17.67*(t-273.15)/(t-29.65))'
    
28. #饱和比湿的计算，单位：克/千克
    
29. 'define qs=1000*0.622*es/(prs-0.378*es)'
    
30. #求实际比湿
    
31. 'define q=rh*qs/100'
    
32. #计算水汽通量 |V|*q/g
    
33. 'define q1=q*u/g'
    
34. 'define q2=q*v/g'
    
35. 
    
36. #显示水汽通量散度
    
37. 'set gxout shaded'
    
38. 'set cmax 0'
    
39. 'define qfdiv=hdivg(q1,q2)*1e6'
    
40. 'd qfdiv'
    
41. 'cbarn 1 1'
    

复制代码

下面两幅图上边是micaps实况水汽通量散度，下边是上述程序计算出来的水汽通量散度，可见两者辐合中心有严重偏差，到底哪里错了呢？

实况

river

'set t 12'
看看是不是时间设置有问题啊，世界时和咱们的时间差了8小时哦。还有就是fnl资料里应该有比湿的资料shum，看看直接用比湿计算水汽通量散度的结果怎么样

bolv6666

fnl 里面没有shum资料啊，只有一层spfh。况且我这个程序也是拷贝的论坛里的，好像大家都没发现什么问题，奇怪了。时间肯定错不了，是第十二个时次的数据，因为别的假相当位温什么的都计算正确。哪位大侠再帮帮忙哦

xiaoxiao2013

我的也是这样的，自己做出来的图和实况相差很远。求解释

xiaoxiao2013

看你的图可以出省界，我的怎么不出来啊，我的现在连中国地图都出不来，我只画出来过一次，后来怎么画都不出来。

seaou

云图怎么加色标

seaou

哪位能告诉，云图怎么加色标

bolv6666

xiaoxiao2013 发表于 2013-4-7 12:28

                               
登录/注册后可看大图

我的也是这样的，自己做出来的图和实况相差很远。求解释

是吗，那可真奇了怪了，以前没有注意过这个问题。不过我看了下700hPa的图和实况岁仍然有差距，但能大体说的过去，只是850差的太多了

xiaoxiao2013

seaou 发表于 2013-4-8 10:30

                               
登录/注册后可看大图

云图怎么加色标

作图的其他项都设好了就加上‘cbarn’嘿嘿刚学的。就用了！

xiaoxiao2013

seaou 发表于 2013-4-8 10:30

                               
登录/注册后可看大图

云图怎么加色标

作图的其他项都设好了就加上‘cbarn’嘿嘿刚学的。就用了！