Re: warning:regCoef: 10512 array(s) contained less than 3 non-missing values

The message is telling you that one or more series
do not have enough non-missing values to perform
the calculation.

If you input 100 values and all are missing ... you
can not perform a regression.

On 4/4/12 6:53 AM, yikun liu wrote:
> Hi,
>
> I tried to draw Geopotential Height trend from 2000-2010, but there are
> two warnings:
> warning:regCoef: 10512 array(s) contained less than 3 non-missing values
> warning:ContourPlotInitialize: no valid values in scalar field;
> ContourPlot not possible:[errno=1101]
>
> I am not sure what do they mean and how they affect the map that I draw.
> And how can i correct my code so there will be no warnings?
>
> Thanks!
>
>
> Yikun
>
> P.S. Below is my code.
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>
> begin
>
> dir =
> "/media/disk/DATA/NCEP/Pressures_Level/Monthly_mean/hgt.mon.mean.nc
> <http://hgt.mon.mean.nc>"
> f1 =addfile(dir,"r")
> year=ispan(2000,2010,1)
> P1 = new((/73,144,11/),float)
>
> do i=0, 10
> T1 = f1->hgt(:,0,:,:)
> a= 254+i*12
> b= 256+i*12
> T2 = dim_avg_Wrap(T1(lat|:,lon|:,time|a:b))
> P1(:,:,i)=T2(lat|:,lon|:)+31265
> reg1 = T2
> reg1 = (/regCoef( year(:),P1(:,:,:) ) /)
>
> xwks=gsn_open_wks("pdf","Global_Seasonal_SPR_hgt_Trend")
>
> res=True
> res@cnFillOn=True ;
> res@gsnMaximize=True
> res@gsnPaperOrientation = "portrait"
> res@cnLinesOn=False
> res@cnLineLabelsOn = False
> res@gsnSpreadColors=True
> res@cnFillMode="RasterFill"
> res@lbOrientation="vertical"
> res@cnMissingValFillPattern = 0
> res@cnMissingValFillColor=0
> res@lbLabelAutoStride = True
> res@cnLevelSelectionMode= "ExplicitLevels"
> res@cnLevels =(/-10,-8,-6,-4,-2,0,2,4,6,8,10/)
> res@gsnSpreadColors= True
>
> gsn_define_colormap(xwks,"BlRe")
>
> res@gsnCenterString="SPR"
> res@gsnLeftString="Geopotential Height"
> res@gsnRightString="2000-2010 Trend "
>
> plot=gsn_csm_contour_map_ce(xwks,reg1,res)
>
>
>
> res@mpLimitMode = "LatLon"
>
> delete(plot) ; cleaning up resources used
> delete(xwks)
> delete(res)
>
> end do
>
>
> 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 Wed Apr 4 07:46:06 2012