求助水汽通量散度的画图问题

river

风暴侦察兵 发表于 2012-10-20 18:49

                               
登录/注册后可看大图

'reinit'
'sdfopen e:\data\ncep\uwnd.2012.nc'
'sdfopen e:\data\ncep\vwnd.2012.nc'

计算通量散度什么的要先在大范围里计算，再在小范围里画图，要不然就会出现你这种情况。这个倒不是错误，式计算的时候需要差分之类的（具体没研究过）就把边界计算没了。所以你要把你计算过程放在一个大范围，画图的时候再设置成你要的范围。改了一下你的gs，画出了几张图：

修改过的gs（注意我设置的时间和你的不同）
如下：

1. 'reinit'
   
2. 'sdfopen f:\2012\uwnd.2012.nc'
   
3. 'sdfopen f:\2012\vwnd.2012.nc'
   
4. 'sdfopen f:\2012\shum.2012.nc'
   
5. 'set grid off'
   
6. 'set mpdset cnworld cnriver'
   
7. 'set map 15 1 2'
   
8. 'set xlopts 1 4 0.2'
   
9. 'set ylopts 1 4 0.2'
   
10. 'set parea 1.5 10 1 8'
    
11. 'set csmooth on'
    
12. 'set lev 850'
    
13. 'set t 1 4'
    
14. 'define qdiv=hdivg(shum.3*1e3*uwnd.1/9.8,shum.3*1e3*vwnd.2/9.8)*1e5'
    
15. nt=1
    
16. while(nt<=4)
    
17. 'set t 'nt''
    
18. 'set ylint 5'
    
19. 'set xlint 5'
    
20. 'set lat 20 45'
    
21. 'set lon 100 135'
    
22. 'set grads off'
    
23. 'set gxout shaded'
    
24. 'd qdiv'
    
25. 'cbarn 0.85 0 5.75 0.35'
    
26. 'set gxout contour'
    
27. 'd qdiv''cbarn 0.85 0 5.75 0.35'
    
28. 'q time'
    
29. rc=sublin(result,1)
    
30. tt=subwrd(rc,3)
    
31. 'draw title 850hPa Divergence of Moisture Flux '
    
32. 'printim f:\2012\1\850hPa-divergence-of-moisture-flux'tt'.gif white'
    
33. 'c'
    
34. nt=nt+1
    
35. endwhile ;

'reinit' <br /> 'sdfopen f:\2012\uwnd.2012.nc' <br /> 'sdfopen f:\2012\vwnd.2012.nc' <br /> 'sdfopen f:\2012\shum.2012.nc'<br /> 'set grid off' <br /> 'set mpdset cnworld cnriver' <br /> 'set map 15 1 2'<br /> 'set xlopts 1 4 0.2' <br /> 'set ylopts 1 4 0.2' <br /> 'set parea 1.5 10 1 8'<br /> 'set csmooth on'<br /> 'set lev 850' <br /> 'set t 1 4'<br /> 'define qdiv=hdivg(shum.3*1e3*uwnd.1/9.8,shum.3*1e3*vwnd.2/9.8)*1e5' <br /> nt=1 <br /> while(nt<=4)<br /> 'set t 'nt'' <br /> 'set ylint 5'<br /> 'set xlint 5' <br /> 'set lat 20 45'<br /> 'set lon 100 135' <br /> 'set grads off'<br /> 'set gxout shaded' <br /> 'd qdiv'<br /> 'cbarn 0.85 0 5.75 0.35' <br /> 'set gxout contour' <br /> 'd qdiv''cbarn 0.85 0 5.75 0.35' <br /> 'q time' <br /> rc=sublin(result,1)<br /> tt=subwrd(rc,3) <br /> 'draw title 850hPa Divergence of Moisture Flux ' <br /> 'printim f:\2012\1\850hPa-divergence-of-moisture-flux'tt'.gif white' <br /> 'c' <br /> nt=nt+1 <br /> endwhile ;

另外你原来计算的量级太小，影响colorbar的美观，所以我把shum的量级扩大了10的3次方，最后水汽通量又扩大了10的5次方，你自己看你出来的单位应该是什么吧。有问题再交流。

风暴侦察兵

river 发表于 2012-10-20 19:36

                               
登录/注册后可看大图

计算通量散度什么的要先在大范围里计算，再在小范围里画图，要不然就会出现你这种情况。这个倒不是错误， ...

高手啊！厉害，以后有啥问题多向你请教哈，我想问一个fnl资料的里的rhprs是不是对应着NCEP的rhum？

river

风暴侦察兵 发表于 2012-10-20 19:58

                               
登录/注册后可看大图

高手啊！厉害，以后有啥问题多向你请教哈，我想问一个fnl资料的里的rhprs是不是对应着NCEP的rhum？

是的。你应该仔细看你下载资料的网站的说明，都代表了什么物理量，单位是什么，他都会有注释。fnl的ctl里就有注释

juliekxkl

风暴侦察兵 发表于 2012-10-19 09:29

                               
登录/注册后可看大图

图是这样的，有什么办法能把水汽通量散度覆盖到整个要画的图面？

循环之前可以设置‘set lon 95 140'
                            'set lat 15 50'后面循环里再继续设置你的'set lat 20 45'
'set lon 100 135'

caojie

把经纬度均多设置大一个单位

xiaoyuyly

我画水汽通量散度则碰到这样的问题，批量出图时，第一张图是满屏的，第二张开始就是两面有空的地方，找了好长时间的原因，一直没找到，请哪位大侠指点，谢谢

落然

收到了，学学习

皮卡丘皮皮虾

学习学习，谢谢分享

大气小白

这个是差分计算必然存在的问题，差分计算散度的时候要用到第i-1，i+1,j-1,j+1个格点，就意味着4个边的值是算不出来的