只想画副高脊线

469482462

这是我画副高脊线的gs文件，但是画出了很多脊线，我只想要副高脊线，怎么修改呢？

1. 'reinit'
   
2. iyear=1958
   
3. 'sdfopen **\hgt.'iyear'.nc'
   
4. 'sdfopen **\uwnd.'iyear'.nc'
   
5. 'set grads off'
   
6. 'set grid off'
   
7. 'set map 1 1 8'
   
8. 'set mpdset hires'
   
9. 'set mpdset cnworld'
   
10. 'set vpage 0 11 0 8.5'
    
11. 'set parea 1.0 10 1 8'
    
12. 'set lat 0 90'
    
13. 'set lon 40 180'
    
14. 'set lev 500'
    
15. 'set t 1'
    
16. k=152
    
17. while(k<=273)
    
18. 'set t ' k
    
19. 'set xlopts 1 1 0.16'
    
20. 'set ylopts 1 1 0.16'
    
21. *'set xlint 5'
    
22. *'set ylint 5'
    
23. 'set grads off'
    
24. 'set grid off'
    
25. 'set cint 2'
    
26. 'set gxout contour'
    
27. 'set cthick 5'
    
28. 'set cterp on'
    
29. 'set csmooth on'
    
30. *'set clevs 0'
    
31. *'set clab off'
    
32. 'd hgt/10.0'
    
33. 
    
34. 'set gxout contour'
    
35. 'set clevs 0'
    
36. 'define ms=cdiff(uwnd.2,y)'
    
37. 'set cthick 8'
    
38. 'set ccolor 2'
    
39. 'd maskout(uwnd.2,ms)'
    

'reinit'<br /> iyear=1958<br /> 'sdfopen **\hgt.'iyear'.nc'<br /> 'sdfopen **\uwnd.'iyear'.nc'<br /> 'set grads off'<br /> 'set grid off'<br /> 'set map 1 1 8'<br /> 'set mpdset hires'<br /> 'set mpdset cnworld'<br /> 'set vpage 0 11 0 8.5'<br /> 'set parea 1.0 10 1 8'<br /> 'set lat 0 90'<br /> 'set lon 40 180'<br /> 'set lev 500'<br /> 'set t 1'<br /> k=152<br /> while(k<=273)<br /> 'set t ' k<br /> 'set xlopts 1 1 0.16'<br /> 'set ylopts 1 1 0.16'<br /> *'set xlint 5'<br /> *'set ylint 5'<br /> 'set grads off'<br /> 'set grid off'<br /> 'set cint 2'<br /> 'set gxout contour'<br /> 'set cthick 5'<br /> 'set cterp on'<br /> 'set csmooth on'<br /> *'set clevs 0'<br /> *'set clab off'<br /> 'd hgt/10.0'<br /> <br /> 'set gxout contour'<br /> 'set clevs 0'<br /> 'define ms=cdiff(uwnd.2,y)'<br /> 'set cthick 8'<br /> 'set ccolor 2'<br /> 'd maskout(uwnd.2,ms)'<br />

mofangbao

楼主理解这里面的每一句话吗？不理解的话就一句句的翻手册，弄懂了你就会画了

469482462

mofangbao 发表于 2012-4-25 10:09

                               
登录/注册后可看大图

楼主理解这里面的每一句话吗？不理解的话就一句句的翻手册，弄懂了你就会画了

我理解，每句话我也知道是什么意思，就是揣摩不出这种现象应该怎么消除~才想请大家帮助~

mofangbao

有一个方法你试试，还是使用maskout功能，把hgt<586的地方mask掉，也就是

1. 'd maskout(maskout(uwnd.2,ms),(hgt/10.0)-586)'
   

'd maskout(maskout(uwnd.2,ms),(hgt/10.0)-586)'<br />

用这个方法需要满足你的两个文件的维数是一致的才行（格点分辨率），否则的话先把数据保存起来，然后编写个程序把其他地区手动赋值为缺测。
如果楼主通过这个方法做出来了，别忘了贴出你做出来的图哦

469482462

mofangbao 发表于 2012-4-25 10:48

                               
登录/注册后可看大图

有一个方法你试试，还是使用maskout功能，把hgt

谢谢楼主，你太给力了，我怎么都没想到可以用这个方法~已出图

mofangbao

469482462 发表于 2012-4-25 11:00

                               
登录/注册后可看大图

谢谢楼主，你太给力了，我怎么都没想到可以用这个方法~已出图

副高当然就得根据副高的特点来想方法啦，嘿嘿 欢迎多来支持！有好东西别忘了来分享啊

469482462

mofangbao 发表于 2012-4-25 11:09

                               
登录/注册后可看大图

副高当然就得根据副高的特点来想方法啦，嘿嘿 欢迎多来支持！有好东西别忘了来分享啊

恩恩，会的~~~画图中发现，用这个方法'd maskout(maskout(uwnd.2,ms),(hgt/10.0)-586)'不能用在垂直流场经向剖面图中，因为各个层副高中心的值不一样，请教这个应该怎么解决呢？

mofangbao

469482462 发表于 2012-4-26 09:56

                               
登录/注册后可看大图

恩恩，会的~~~画图中发现，用这个方法'd maskout(maskout(uwnd.2,ms),(hgt/10.0)-586)'不能用在垂直流场经 ...

没做过相关的，这种自己想想啦

陌小默的

这个学习了~~

shyongwei

学习了，幸苦了