Re: error creating spline approximation for trXCoordPoints

From: Xi Chang <xi.chang01_at_nyahnyahspammersnyahnyah>
Date: Tue Dec 03 2013 - 02:37:12 MST

I solved this,
I forgot to switch off this resource : "resa@gsnAddCyclic = True"

Thank you.

On Tue, Dec 3, 2013 at 10:33 AM, Xi Chang <xi.chang01@gmail.com> wrote:

> Hallo,
> When I wanted to overlay my plot with a significant plot, I just simply
> got this error.. anybody can tell me why?
>
>
> ;============================
>
> N1 = equiv_sample_size (m1(lev|:,lat|:,time|:), 0.5,0)
> N2 = equiv_sample_size (m2(lev|:,lat|:,time|:), 0.5,0)
> N3 = equiv_sample_size (m3(lev|:,lat|:,time|:), 0.5,0)
>
> prob = ttest(mm1,var1,N1,mm3,var3,N3, True, True) ; t-test
> prob!0="var"
> prob!1="lev"
> prob!2="lat"
> prob&lev=fin_1->lev
> prob&lat=fin_1->lat
>
> ;============================
>
> wks = gsn_open_wks ("x11", "ly_season" )
> plot = new(1,graphic)
>
> res = True ; plot mods desired
> res@gsnDraw = False
> res@gsnFrame = False
>
> res@cnLineLabelsOn = True ; turn on line labels
> res@cnFillOn = False ; turn on color fill
> res@gsnLeftString = " "
> res@gsnRightString = " U (ms~S~-1~N~)"
> res@gsnCenterString = " "
> res@gsnContourZeroLineThicknessF = 3
> res@cnInfoLabelOn = False
>
> res@gsnContourNegLineDashPattern = 1
>
> res@vpWidthF = 0.50
> res@vpHeightF = 0.65
> p1 = (/ 840,100,10,1,0.1/)
> tdz1 = stdatmus_p2tdz(p1)
> Z1 = tdz1(2,:)/1000
> res@tiYAxisString = "Pressure [hPa]"
> res@tmYLLabelsOn = True
> res@tmYLOn = True
> res@tmYLMode = "Explicit"
> res@tmYLLabels =
> (/"10~S~3~N","10~S~2~N","10~S~1~N","10~S~0~N","10~S~-1~N"/)
> res@tmYLValues = Z1
>
> res@tmXBMode = "Explicit"
> res@tmXBLabels = (/ 10, 20, 30, 40, 50, 60, 70, 80, 90 /)
> res@tmXBValues = res@tmXBLabels
>
> res@tmXMajorGrid = True ; implement x grid
> res@tmXMajorGridThicknessF = .5 ; 2.0 is default
> res@tmXMajorGridLineDashPattern = 2 ; select short dash
> lines
>
> res@tmYMajorGrid = True ; implement y grid
> res@tmYMajorGridThicknessF = .5 ; 2.0 is default
> res@tmYMajorGridLineDashPattern = 2 ; select short dash
> lines
>
>
> res@cnLevelSpacingF = 5.
> res@gsnLeftString="a). 1955-2005 (JFM) "
> plota1 = gsn_csm_contour(wks,dif2({0.1:54},:),res)
>
> ;=============
>
> resa = True ; res2 probability plots
> resa@gsnDraw = False ; Do not draw plot
> resa@gsnFrame = False ; Do not advance frome
> resa@gsnAddCyclic = True
> resa@vpWidthF = 0.50
> resa@vpHeightF = 0.65
> resa@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
> resa@cnMinLevelValF = 0.00 ; set min contour level
> resa@cnMaxLevelValF = 0.15 ; set max contour level
> resa@cnLevelSpacingF = 0.05 ; set contour spacing
> resa@cnInfoLabelOn = False
> resa@cnLinesOn = False ; do not draw contour lines
> resa@cnLineLabelsOn = False ; do not draw contour labels
> resa@cnFillScaleF = 0.9 ; add extra density
> resa@tmXBMode = "Explicit"
> resa@tmXBLabels = (/ 10, 20, 30, 40, 50, 60, 70, 80, 90 /)
> resa@tmXBValues = resa@tmXBLabels
>
>
> plota2 = gsn_csm_contour(wks,prob(0,{0.1:54},:), resa)
> plota2 = ShadeLtContour(plota2, 0.5, 1) ; shade signif <0.1
> overlay (plota1, plota2)
> plot(0)=plota1
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Dec 3 02:37:21 2013

This archive was generated by hypermail 2.1.8 : Wed Dec 04 2013 - 20:42:38 MST