error creating spline approximation for trXCoordPoints

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

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:34:02 2013

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