Re: cd_calendar returns 12 values instead of 1

Try print (f->time) and print (TIME). My guess is that either all 12
values of the original time coordinates are identical, or else you
have sub-daily data, 12 different hours/minutes on the same calendar
date.

--Dave

On Mon, Oct 1, 2012 at 1:03 AM, Almami Johnson <almamij@gmail.com> wrote:
> Dear all,
>
> Trying to get some daily data over a specific time period. After applying
> cd_calendar, it returns 12 values instead one 1 for iStrt and for iLast. Of
> course, I got an error like this:
>
> (0) 730 4383
> (1) 761 4414
> (2) 790 4442
> (3) 821 4473
> (4) 851 4503
> (5) 882 4534
> (6) 912 4564
> (7) 943 4595
> (8) 974 4626
> (9) 1004 4656
> (10) 1035 4687
> (11) 1065 4717
> fatal:Illegal Subscript. Only scalar values are allowed in subscript ranges.
>
> Here is the portion of the code. Could you please help me on this? Thanks.
>
> undef("getPrecip")
> function getPrecip(fname[1]:string,ymStrt[1]:integer,ymLast[1]:integer,
> vname[1]:string)
> local TIME, YYYYMM, iStrt, iLast, x
>
> begin
> f = addfile(fname, "r")
> TIME = f->time
> YYYYMMDD = cd_calendar(TIME, -2) ; convert
> iStrt = ind(ymStrt.eq.YYYYMMDD) ; index of start time
> iLast = ind(ymLast.eq.YYYYMMDD) ; last time
>
> print(" " + iStrt + " " + iLast)
>
> x = f->$vname$(iStrt:iLast,{-5:25},{-10:10}) ; read specified
> time period
>
> _______________________________________________
> 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 Oct 1 19:27:12 2012