warning:ContourPlotInitialize: no valid values in scalar field; ContourPlot not possible:[errno=1101]

Hello

I found the composite of sst (aveX)and also have done ttest. and fount
alpha...
alpha = 100.*(1-ttest(aveX,VarX,sX,aveY,varY,sY,True,False)...I want to
overlay the values of alpha=95 as contour over the shaded plot(aveX). When I
tried by using alpha=where (alpha.ne.95.,alpha@_FillValue,alpha)..I am
getting the following waring:

   warning:ContourPlotInitialize: no valid values in scalar field;
ContourPlot not possible:[errno=1101]
Means I want to overlay the only contour corresponding to the values of
alpha=95...If I use alpha = where(alpha.lt.95.,alpha@_FillValue,alpha) it
will overlay all the contour which is having values greater than 95..I want
only outer contour

Thanks in advance

aveX=dim_avg_n_Wrap(EW_JJAS,0)
x=aveX
y=aveX
sX=28
sY=84
aveY=dim_avg_n_Wrap(JJAS_1982_02,0)
varY=dim_variance_n_Wrap(JJAS_1982_02,0)

alpha = 100.*(1. - ttest(aveX,varX,sX, aveY,varY,sY,True,False))
alpha = where(alpha.ne.95.,alpha@_FillValue,alpha)

  alpha!0="lat"
  alpha!1="lon"
  alpha&lat=a&lat
  alpha&lon=a&lon

;***********************************************************
wks=gsn_open_wks("ps","EW_JJAS_SSTA_composite_signi3")
gsn_define_colormap(wks,"BlueRed")
res=True
res90=True
;************************************************************
res_at_gsnDraw=False
res_at_gsnFrame=False
res90=res
;***********************************************************
res_at_lbLabelAutoStride=True
res_at_tmXTOn=False
res_at_tmYROn=False
res_at_cnFillOn=True
res_at_gsnSpreadColors=True
res_at_cnLinesOn=False
res_at_mpCenterLonF=-180.
;************************************************************
res90_at_cnLinesOn=True
res90_at_cnLineLabelsOn=True
res90_at_gsnContourNegLineDashPattern=2
res90_at_cnLineColor="black"
res90_at_cnLineThicknessF=1.5
base=gsn_csm_contour_map(wks,aveX,res)
sig90=gsn_csm_contour(wks,alpha,res90)
;********************************************************
overlay(base,sig90)
draw(base)
frame(wks)

end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Sep 17 2009 - 00:29:40 MDT