Re: daily climatology

> I want to derive a daily climatology from a monthly climatology using linear
> interpolation. My climatology variable (SST) has dimensions 12,nlat,nlon. Is
> there any function that will do that?
>
> At the end of the day I want to calculate weekly anomalies so I am looking at
> building my daily climatology to make my weekly anomalies.
________________________________________

http://www.ncl.ucar.edu/Document/Functions/Built-in/linint1.shtml
or, better
http://www.ncl.ucar.edu/Document/Functions/Contributed/linint1_Wrap.shtml

If SST(time ,lat,lon) with sizes (12,nlat,nlon)

   midDayMon = (/15.5,45,74.5,105,135.5,166,196.5,227.5,258,288.5,319,349.5/)
   midDayMon_at_long_name = "middle day of month"

   day = ispan(1,365,1)
   day!0 = "day"
   day_at_long_name = "day of year: no leap"
   day_at_units = "1=Jan 1, 32=Feb 1, ..., 365-Dec 31"

                            ; reorder and select cylic conditions
   sstClmDay = linint1_Wrap (midDayMon, SST(lat|:,lon|:,time|:), True, day, 0)
                                                                 ^^^^
                                                                 important
   printVarSummary(sstClmDay) ; (lat,lon,365)

   sstDayClm = sstClmDay(day|:,lat|:,lon|:)) ; (365,lat,lon)
   printVarSummary(sstDayClm) ; (lat,lon,365)
; delete(sstClmDay)

Good Luck
D

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri May 05 2006 - 17:30:20 MDT