zonal图对数坐标轴出现的错误，求助

扬风三郎

zonal t随纬度和高度的变化，画等高线图，垂直坐标设成气压对数坐标，理想的是1000，,100,10,1,10-1,10-2,10-3,10-4,10-5。
66层气压。  画图出错，出不来气压对数。报错：
warning:_NhlCreateSplineCoordApprox: Attempt to create spline approximation for Y axis failed: consider adjusting trYTensionF value
warning:IrTransInitialize: error creating spline approximation for trYCoordPoints; defaulting to linear
warning:IrTransInitialize: LogAxis requires all positive data extent: defaulting trYAxisType to LinearAxis

程序如下：
16   zave = dim_avg_Wrap(T)                             ; calculate zonal ave
     17 ;************************************************
     18   p=in->lev_p                                    
     19 ;************************************************
     20 ; create plot
     21 ;************************************************
     22   wks  = gsn_open_wks("ps","zonalt")               ; Opens a ps file
     23   res                          = True             ; plot mods wanted
     24   res@sfYArray                 = p            ; put eta on y-axis
     25   res@sfXArray                 = zave&lat         ; put lat on x-axis
     26   res@tiXAxisString            = "latitude"       ; x-axis label
     27   res@tiYAxisString            = "p"     ; y-axis label
     28 res@cnInfoLabelOn=False        
     29 res@trXReverse                 =True
     30 res@cnLevelSpacingF=10         
     31   res@gsnYAxisIrregular2Log = True             ; convert y-axis to linear
     32 res@cnLevelSpacingF=10
     33   plot = gsn_csm_contour(wks,zave,res)         ; creates plot
     34 end

尽头的尽头

NCL有专门的对数压力图的函数http://www.ncl.ucar.edu/Applications/height_lat.shtml

扬风三郎

尽头的尽头 发表于 2012-12-12 13:18

                               
登录/注册后可看大图

NCL有专门的对数压力图的函数http://www.ncl.ucar.edu/Applications/height_lat.shtml

试过了，错误一样的。。不知道错误是什么意思。。我看的函数是http://www.ncl.ucar.edu/Applications/axes.shtml

hitwjk

求助同样的问题

和桥豆

res@gsnAddCyclic = False  ？