已经解决！绘制的700hpa水平涡度分布图，四周有一圈空白是什么情况？

tvfantasy

已经解决

1. 'reinit'
   
2. 'sdfopen D:\data\uwnd.2016.nc'
   
3. 'sdfopen D:\data\vwnd.2016.nc'
   
4. *set map*
   
5. 'set mpdset cnworld'
   
6. 
   
7. t=737
   
8. while(t<741)
   
9. <font color="red">'set lon 70 170'
   
10. 'set lat 0 60'</font>
    
11. 'set lev 700'
    
12. 'set t 't
    
13. 'define vort=hcurl(uwnd.1,vwnd.2)*10000'
    
14. <font color="red">'set lon 80 160'
    
15. 'set lat 10 50'
    
16. 'd vort'</font>
    
17. *set title with elements and time*
    
18. 'q time'
    
19. rc=sublin(result,1)
    
20. tt=subwrd(rc,3)
    
21. 'draw title 700hpa_Vorticity_'%tt%''
    
22. 'printim D:\by_pic\700hpa_Vorticity_'%tt%'.png white x1000 y800'
    
23. 'c'
    
24. t=t+1
    
25. endwhile
    
26. 'reinit'

复制代码

求助！绘制的700hpa水平涡度分布图，四周有一圈空白是什么情况？

脚本

1. 'reinit'
   
2. 'sdfopen D:\data\uwnd.2016.nc'
   
3. 'sdfopen D:\data\vwnd.2016.nc'
   
4. *set map*
   
5. 'set mpdset cnworld'
   
6. 'set lon 80 160'
   
7. 'set lat 10 50'
   
8. 
   
9. t=737
   
10. while(t<741)
    
11. 'set lev 700'
    
12. 'set t 't
    
13. 'define vort=hcurl(uwnd.1,vwnd.2)*100000'
    
14. 'd vort'
    
15. *set title with elements and time*
    
16. 'q time'
    
17. rc=sublin(result,1)
    
18. tt=subwrd(rc,3)
    
19. 'draw title 700hpa_Vorticity_'%tt%''
    
20. 'printim D:\by_pic\700hpa_Vorticity_'%tt%'.png white x1000 y800'
    
21. 'c'
    
22. t=t+1
    
23. endwhile
    
24. 'reinit'

复制代码

忽忽儿

计算的经纬度范围要取的比画图的经纬度大才行。你要画什么范围，就把计算的经纬度取大一圈试试。我没画过涡度，不过其他物理量有这个问题，你试试看。
'set lon 70 170'
'set lat 0 60'
'define vort=hcurl(uwnd.1,vwnd.2)*100000'
'set lon 80 160'
'set lat 10 50'
'd vort'

tvfantasy

忽忽儿 发表于 2017-5-26 17:57
计算的经纬度范围要取的比画图的经纬度大才行。你要画什么范围，就把计算的经纬度取大一圈试试。我没画过涡 ...

感谢指点，辛苦拉！

tvfantasy

忽忽儿 发表于 2017-5-26 17:57
计算的经纬度范围要取的比画图的经纬度大才行。你要画什么范围，就把计算的经纬度取大一圈试试。我没画过涡 ...

感谢指点，辛苦拉！

tvfantasy

tvfantasy 发表于 2017-5-27 07:42
感谢指点，辛苦拉！

你用nc资料画过地面图没，我画出来的图每天都是一样的，搜索了下也是说的不清出。

xuezhanbin

{:5_201:}{:5_201:}

hbl1021

请问楼主问题是怎么解决的啊？我也遇到了这种情况{:eb513:}{:eb513:}