Re: month_to_annual

For some reason the 4-dimensional case had an incorrect
mask statement. A bug fix has been implemented for
the next release of NCL. See:

http://www.ncl.ucar.edu/future_release.shtml

===
I will send Matt the code offline.

D

>
>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

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Feb 16 2007 - 09:59:24 MST