Re: Missing value problem when removing Annualcycle from daily data

Change:
Dlyclm=clmDayTLL(rain, yyyyddd)
smthDlyclm=smthClmDayTLL(Dlyclm,3)
rm_anncycl=calcDayAnomTLL(rain,yyyyddd,smthDlyclm)

To:
Dlyclm=clmDayTLL(rain, yyyyddd)
rm_anncycl=calcDayAnomTLL(rain,yyyyddd,Dlyclm)

Sabeerali(sebi) wrote:
> I want to remove annual cycle from daily data...My data is having
> missing values...So it shows the following error. What I can do for
> removing Annualcycle?
>
> (0) smthClmDay_TLL: No missing values allowed: ezfftf does not allow
> (0) smthClmDay_TLL: nFill=292068
>
> Any help would be appreciated.
>
>
> 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
> f=addfile("/home/sebi/data/IMD/daily/rf_1901_03.nc
> <http://rf_1901_03.nc>","r")
> rain=f->rain(17897:37619,0,:,:)
> t=f->time(17897:37619)
> rain!1="lat"
> rain!2="lon"
> TIME=ut_calendar(t,0)
> year=floattointeger(TIME(:,0))
> month=floattointeger(TIME(:,1))
> day=floattointeger(TIME(:,2))
> ddd=day_of_year(year,month,day)
> yyyyddd=year*1000+ddd
> Dlyclm=clmDayTLL(rain, yyyyddd)
> smthDlyclm=smthClmDayTLL(Dlyclm,3)
> rm_anncycl=calcDayAnomTLL(rain,yyyyddd,smthDlyclm)
> printVarSummary(rm_anncycl)
> end
>
>
> --
> **********************************
> Sabeerali
> Climate and Global Modeling Division
> Indian Institute of Tropical Meteorology
> Pashan, Pune, 411 008
> ****************************************
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> 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 Mar 18 08:32:03 2010