关于NCL画垂直剖面图

lleoiu

2017.9.14:
在NCL的6.2之后的版本中，gsn_csm_pres_hgt_vector函数里默认的vcMapDirection=False
“This resource now defaults to False for this function in NCL V6.2.0.”
http://www.ncl.ucar.edu/Document ... es_hgt_vector.shtml




用NCL画垂直剖面图的时候，总是感觉水平方向上的矢量要比垂直方向上的大（垂直方向上已经乘上了－100），并且画出来的垂直方向总是与GRADS画的方向相反。如图1（NCL），图2（grads)。

fig1

这是NCL默认的设置。
The problem is that by default NCL maps the vector direction into the underlying coordinate space. This is good when drawing on a map projection, but for vertical plots where the two coordinate axes have different units that have very different numerical spacing between the coordinate values it leads to distortion. There is a simple fix, which is to set the resource vcMapDirection to False.


加上res@vcMapDirection   = False这一设置之后，画的图与GRADS一样，垂直方向上只需乘上100即可。如图3.

1. begin
   
2.          
   
3.     f2 = addfile("/data/NCEP-2.5/monthly/pressure/vwnd.mon.mean.nc","r")
   
4.     f3 = addfile("/data/NCEP-2.5/monthly/pressure/omega.mon.mean.nc","r")
   
5. 
   
6.    vwnd1 = short2flt(f2->vwnd(7,0:11,{0:90},{95:110}))
   
7.    vwnd11 = dim_avg_n_Wrap(vwnd1,2)
   
8. 
   
9.    omega1 = short2flt(f3->omega(7,0:11,{0:90},{95:110})) ;{45:20}
   
10.    omega11 = dim_avg_n_Wrap(omega1,2)
    
11.    w = omega11*100.
    
12.          
    
13.    wks  = gsn_open_wks("eps","vv5")
    
14.    gsn_define_colormap(wks,"gui_default")
    
15. 
    
16.    res                  = True  
    
17.    res@cnLinesOn        = False
    
18.    res@cnLineLabelsOn   = False
    
19.    res@cnInfoLabelOn    = False
    
20.    res@tmXTOn           = False   ; and right tickmarks.
    
21.    res@tmYROn           = False
    
22.    res@vcRefMagnitudeF  = 5.                 
    
23.    res@vcRefAnnoString1On = True
    
24.    res@vcRefAnnoString1   = "5m/s"
    
25.    res@vcRefAnnoOrthogonalPosF   = -1.095  
    
26.    res@vcRefAnnoString2On = False     
    
27.    res@vcRefLengthF     = 0.027                 
    
28.    res@vcGlyphStyle     = "LineArrow
    
29.    res@vcMinDistanceF   = 0.008
    
30.    res@vcMapDirection   = False
    
31. 
    
32.    res@vpWidthF         = 0.75
    
33.    res@vpHeightF        = 0.6
    
34.    res@tiMainFontHeightF = 0.02
    
35.    res@gsnLeftString   = ""
    
36.    res@gsnRightString  = ""            
    
37.    res@tiYAxisString   = "p/hPa"
    
38.    res@tmYRMode        = "Automatic"
    
39.    res@tmXBMode   = "Explicit"
    
40. 
    
41.    plot = gsn_csm_pres_hgt_vector(wks,vwnd11,vwnd11,w,res)   
    
42. end   
    
43. 
    
44. 
    

复制代码

freekiller

收藏了，一定会用到。谢谢

kongfeng0824

谢谢分享            

小拱猪Solo

这个好厉害  收藏了！

guigui

非常实用！

lemondaodao

看看，谢谢

吕丽

刚好试试看

宁小才

收藏了，一定会用到。谢谢

小胖子

请问楼主有没有试过用NCL画垂直剖面图时，加地形进去？？

lleoiu

小胖子 发表于 2014-1-4 16:10

                               
登录/注册后可看大图

请问楼主有没有试过用NCL画垂直剖面图时，加地形进去？？

没有诶。还没有画过加地形的剖面图