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
"
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
Received on Wed Apr 4 06:53:19 2012
This archive was generated by hypermail 2.1.8 : Mon Apr 09 2012 - 13:43:03 MDT