Re: Linear interpolation in time

> I want to interpolate a 1x1 grid pentad data to daily data. I tried using
>the function dspnt3 but it says need to have same number of x,y values. In my
>data the number of lat lon values are unequal.
> Is there any function in NCL that would interpolate the data in time. Any
>help would be appreciated.
>_______________________________________________

I suggest linear interpolation. The following is
completely untested. You will have to debug. I
am very busy.

consider x5(t5,lat,lon) where t5 is the pentad time dimension

          dim_x5 = dimsizes( x5 )
          ntim5 = dim_x5(0)
          nlat = dim_x5(1)
          mlon = dim_x5(2)
          
          ntim = 5*ntim5 ; 5x the number of pentads = # days
          
          time5 = ispan(3,ntim,5)
          time5!0 = "time5"
          time5_at_units = "pentad day"
          time5&time5 = time5
          
          time = ispan(1,ntim,1) ; create time/day coordinate variable
          time!0 = "time"
          time_at_units = "day"
          time&time = time
          
          xTemp = linint1_Wrap (time5, x5(lat|:,lon|:,time5), False, time, 0)
          printVarSummary(xTemp)
          
          xDay = xTemp(time|:,lat|:,lon|:)
          delete(xTemp)
          
          copy_VarAtts(x5, xDay)
          
          xDay!0 = "time"
          xDay!1 = "lat"
          xDay!2 = "lon"
          xDay&time= time
          xDay&lat = x5&lat
          xDay&lon = x5&lon
          
          printVarSummary( xDay )
          printminMax( xDay, True )
          
          printtime+" "+(xDay(;,{45},{90})) ; time series at one opint

The use linint2_Wrap to interpolate the 2.5 daily data to 1x1

good luck

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jan 24 2007 - 16:23:24 MST