NCL 5.1.1: Issues with concatenating a subset of an array from multiple files

From: Yongxin Zhang <yongxin.fred_at_nyahnyahspammersnyahnyah>
Date: Wed Sep 22 2010 - 09:22:30 MDT

I have two netCDF files that contain only one variable, RAIN. RAIN is
dimensioned as (time, south_north, west_east) where time=25, south_north=765
and west_east=1191. I needed to read out the first 24 data (i.e., time=0 to
23) at each grid point from each file and then concatenate them to become a
new variable. I used "systemfunc" and "addfiles" but strange thing was that
when I used the following statement

obs_rain = rain_fin[:]->RAIN(0:23,:,:)

I would get obs_rain dimensioned as (24, 765, 1191) instead of (48, 765,
1191). Same thing happened even though I explicitly set
ListSetType(rain_fin,"cat"). However, if I used this statement

obs_rain00 = rain_fin[:]->RAIN(:,:,:)

I would get obs_rain00 dimensioned as (50, 765, 1191) which is right. I
tried this on two different Linux clusters at NCAR using both NCL 5.1.1 and
NCL 5.1.0 and I got the same results. I am not sure if there is a bug in NCL
5.1.1 and NCL 5.1.0 or I need to set something else. Please advise. Thanks.

Yongxin

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Sep 22 09:22:36 2010

This archive was generated by hypermail 2.1.8 : Mon Oct 04 2010 - 08:55:55 MDT