Re: Time, let it regular

Hello

Possibly you did not name a dimension.
NCL functions that are _Wrap expect variables
to be netCDF variable ... named dimensions.

Normally, you start with

    time(time) ; original unequal spaced time
    x(time,lat,lon)
    ntim = dimsizes(time) ; 365

you want

    TIME(TIME) ; equally spaced time
    X(TIMR,lat,lon)
    NTIM = ntim

    TIME = fspan(1,365,365)
or, maybe
    TIME = fspan(time(0),time(ntim-1),365)
or, ?????

Add

     TIME!0 = "TIME"
     TIME_at_units = time_at_units ; assuming 'time since ...'

xNew = linint1_Wrap(time, x(lat|:,lon|:,time|:), False, TIME,0)

=======

If you are new to NCL, please read the mini-language manual

http://www.ncl.ucar.edu/Document/Manuals/

Good luck

Gwenaelle wrote:
> Hi Dennis !
>
> Thanks very much for your answer.
>
> When I do the interpolation (see bellow) following your advice, I have a
> warning:
>
> warning:Dimension (0) has not been defined
>
> What does this mean?
>
> Thanks!
>
> On Jun 18, 2009, at 12:09 PM, Dennis Shea wrote:
>
>> esp365 = fspan(1,365,365) ; euqally spaced time
>>
>> xNew = linint1_Wrap(usp365, x(lat|:,lon|:,time|:), False, esp365,0)
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Jun 18 2009 - 11:56:07 MDT