Re: subscript time coordinate of OLR data

As noted you must use the time units on the file not
any arbitrary units. Something like:

f = addfile("...", "r) ; or f = addfiles(fils, "r”)
                                   ; ALL times on file
    TIME = f->time ; TIME = f[:]->time
    YYYYMM = cd_calendar(TIME, -1) ; convert
    ymStrt = 190801 ; year-month start
    ymLast = 200712 ; year-month last
    iStrt = ind(ymStrt.eq.YYYYMM) ; index of start time
    iLast = ind(ymStrt.eq.YYYYMM) ; last time
    x = f->X(iStrt:iLast,...) ; read only specified times

    xAvg = dim_avg_n (x, 0) ; dim_avg_n_Wrap

On 6/25/12 6:29 AM, Carl Schreck wrote:
> Hi Jouky,
>
> If you run ncl_filedump on olr.day.mean.nc <http://olr.day.mean.nc>, it
> should show you the correct timeUnits:
>
> ---
> double time ( time )
> * units :hours since 1-1-1 00:00:0.0*
> long_name :Time
> actual_range :( 17298624, 17631720 )
> delta_t :0000-00-01 00:00:00
> standard_name :time
> axis :T
> ---
>
> So in this case you will have to modify the other snippet of code to
> account for the times being in hours instead of days:
>
> nt = 1 + round( *24 ** ( maxDate - minDate ), 3 )
>
> Hope this helps!
>
> Carl
>
> On Mon, Jun 25, 2012 at 3:21 AM, juki juki <juky_emc2@yahoo.com
> <mailto:juky_emc2@yahoo.com>> wrote:
>
> Dear All ncl users,
>
> I need your help for subscripting the time coordinate of OLR data.
> For example, I want to subscript the following time:
>
> timeUnits = "days since 1800-01-01 00:00:00"
> inStart = ut_inv_calendar( 1980, 01, 01, 00, 0, 0, timeUnits, 0 )
> inEnd = ut_inv_calendar( 2000, 12, 31, 00, 0, 0, timeUnits, 0 )
>
> varName = "olr"
>
> pathIn = "olr.day.mean.nc <http://olr.day.mean.nc>"
>
> ;; Open the input files
> fin = addfile( pathIn, "r" )
> time = fin->time({inStart:inEnd})
> lat = fin->lat
> lon = fin->lon
>
> data = short2flt(fin->$varName$({inStart:inEnd},:,:))
>
> -----------------
> With this subscripting technique, I got error "subcript out of
> range. The error is due to value of inStart and inEnd being not the
> same as one of value of time. I try to use another timeunits for
> example, timeUnits = "hours since 1-1-1 00:00:0.0", I can read the
> data but when making another process for example:
>
> maxDate =inEnd
> minDate =inStart
>
> nx = 144
> nt = 1 + round( maxDate - minDate, 3 )
> nVar = 1; suppose only 1
>
> datax = new( (/ nVar, nt, nx /), float )
> datax!0 = "var"
> datax!1 = "time"
> datax!2 = "lon"
>
> datax(0,:,:) = dim_avg_n(data, 1 )
>
> I got error again because number of row of datax and data is not the
> same in which value of maxDate - minDate with timeUnits = "hours
> since 1-1-1 00:00:0.0" is much bigger than total number of row of
> OLR data. THUS, how to subscript the time coordinate of OLR data ?
> What time unit must be used ?
>
> Thanks for help,
>
> Regards,
> Jouky
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
>
> --
> Carl J. Schreck III, PhD
> Postdoctoral Research Associate
> Cooperative Institute for Climate and Satellites (CICS-NC)
> North Carolina State University and NOAA's National Climatic Data Center
> 151 Patton Avenue
> Asheville, NC 28801
> Tel: 828-257-3140
> carl@cicsnc.org <mailto:carl@cicsnc.org>
> http://monitor.cicsnc.org/mjo/
> http://www.atmos.albany.edu/student/carl/
>
>
>
> _______________________________________________
> 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 Mon Jun 25 07:21:29 2012