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