month_to_annual

hi folks, I've got a problem while using month_to_annual.
The code (below) is simple, I read in ocean model temperature data at
monthly resolution in a netcdf file which is 120 months long, then I
try to annualize it. The output data set is all fill value, which is
not right. I'm using output from the IPCC AR4 archive at pcmdi so
I'm sure that the data themselves are fine. to recreate the problem
try downloading one of the standard files like 20c3m/ocn/mo/thetao/
gfdl_cm2_0/run1/thetao_O1.195101-196012.nc, or I can send it along.

;=========================================;
; annualize.ncl
;=========================================;
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
;=========================================;
;=========================================;
   fname="thetao_O1.195101-196012.nc"
   in = addfile(fname,"r")
   y = in->thetao(:,0:1,:,:)
   lat=in->lat
   lon=in->lon
   time=in->time
   depth=in->depth(0:1)
   nlon=dimsizes(lon)
   nlat=dimsizes(lat)
   nlev=dimsizes(depth)
printVarSummary(y)
TMP=month_to_annual(y,1)
printVarSummary(TMP)
ncdf = addfile("annualized_"+fname,"c") ; open output netCDF file
ncdf->thetao=TMP
end

-------------------------------------------------------------
while running it produces the following message which seems very odd
to me.
(0) month_to_annual: some grid points have missing data: nt=0
nyr=0 num(nGood)=101216
(0) month_to_annual: some grid points have missing data: nt=12
nyr=1 num(nGood)=101216
(0) month_to_annual: some grid points have missing data: nt=24
nyr=2 num(nGood)=101216
(0) month_to_annual: some grid points have missing data: nt=36
nyr=3 num(nGood)=101216
(0) month_to_annual: some grid points have missing data: nt=48
nyr=4 num(nGood)=101216
(0) month_to_annual: some grid points have missing data: nt=60
nyr=5 num(nGood)=101216
(0) month_to_annual: some grid points have missing data: nt=72
nyr=6 num(nGood)=101216
(0) month_to_annual: some grid points have missing data: nt=84
nyr=7 num(nGood)=101216
(0) month_to_annual: some grid points have missing data: nt=96
nyr=8 num(nGood)=101216
(0) month_to_annual: some grid points have missing data: nt=108
nyr=9 num(nGood)=101216

thanks for your help,
-matthew

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Feb 15 2007 - 20:47:20 MST