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

Hi,

It sounds like perhaps the test for exactly 95. is causing alpha to be
filled with all _FillValue (?)

What happens if you explicitly set the contour levels to what you
want; something on the order of:

   res95_at_cnLevelSelectionMode = "ExplicitLevels"
   res95_at_cnLevels = (/ 95. /)

I hope this helps...

Rick

On Sep 17, 2009, at 12:29 AM, Sabeerali(sebi) wrote:

> 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

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