[ncl-talk] select year from the model ouput

Dennis Shea shea at ucar.edu
Wed Dec 17 07:40:24 MST 2014


Another alternative

TIME = cd_calendar(time, 0)
  i2010 = ind(TIME(:,0).eq.2010)

  x2010= x(i2010,...)

On Sat, Dec 13, 2014 at 2:07 PM, Xi Chang <xi.chang01 at gmail.com> wrote:
>
> sorry, I solved it.
>
> iStrt    = ind(ymd.eq.ymdStrt)              ; index start
> iLast   = ind(ymd.eq.ymdLast)             ; index last
>
> On Sat, Dec 13, 2014 at 7:14 PM, Xi Chang <xi.chang01 at gmail.com> wrote:
>>
>> Hallo NCL,
>>
>> I have a very simple question regaring my CAM output.
>> I'd like to extract data only for year 2010 (for instance) from the data
>> which has
>> time dimension as follows:
>>
>> lat = 96 ;
>> lon = 144 ;
>> time = UNLIMITED ; // (8833 currently)
>> variables:
>> float lat(lat) ;
>> lat:long_name = "latitude" ;
>> lat:units = "degrees_north" ;
>> float lon(lon) ;
>> lon:long_name = "longitude" ;
>> lon:units = "degrees_east" ;
>> * float time(time) ;*
>> * time:long_name = "time" ;*
>> * time:units = "days since 1765-01-01 00:00:00" ;*
>> * time:calendar = "noleap" ;*
>> * time:bounds = "time_bnds" ;*
>> * int date(time) ;*
>> * date:long_name = "current date (YYYYMMDD)" ;*
>>
>> Any solution how to do that?
>> Thank you
>>
>> Chang.
>>
>>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20141217/718b084d/attachment.html 


More information about the ncl-talk mailing list