> 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
This archive was generated by hypermail 2.2.0 : Mon May 08 2006 - 09:15:22 MDT