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

风暴侦察兵

画出的水汽通量散度图，如何将周边没有画到的部分与图边的经纬度连接起来呢？file:///C:/Users/iap-nju-qian/AppData/Roaming/Tencent/Users/12270645/QQ/WinTemp/RichOle/FKGIMKUMAO%25(2RW%60J%60EJO(O.jpg

天目神眉

图挂了，请楼主上图，没明白你的意思

river

我也没明白，楼主上图，把问题说明白点儿啊

风暴侦察兵

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

风暴侦察兵

river 发表于 2012-10-18 20:25

                               
登录/注册后可看大图

我也没明白，楼主上图，把问题说明白点儿啊

图示已发，望高手解答

风暴侦察兵

天目神眉 发表于 2012-10-18 20:04

                               
登录/注册后可看大图

图挂了，请楼主上图，没明白你的意思

图示已发，望高手解答

river

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

                               
登录/注册后可看大图

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

楼主这个是什么资料画的啊，如果是通过FORTRAN计算得来的话，应该是差分之后边界就少掉一圈。具体的我也不太清楚

风暴侦察兵

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

                               
登录/注册后可看大图

楼主这个是什么资料画的啊，如果是通过FORTRAN计算得来的话，应该是差分之后边界就少掉一圈。具体的我也不 ...

就是NCEP资料画的水汽通量散度，只是后续处理不太会，想把散度的分布覆盖到整个图面

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 16:54

                               
登录/注册后可看大图

就是NCEP资料画的水汽通量散度，只是后续处理不太会，想把散度的分布覆盖到整个图面

如果是NC资料的话是不应该出现这种问题的。你把gs贴上来吧。我怀疑是你计算的时候设置的经纬度比你要画的经纬度范围小了