[ncl-talk] Question regarding array indexing

[Prashanth Bhalachandran]

Dave, 
This is a very elegant solution. My sincere appreciation. 

Prashanth 
> On Oct 17, 2017, at 2:10 PM, Dave Allured - NOAA Affiliate <dave.allured at noaa.gov> wrote:
> 
> Prashanth,
> 
> I like using a modulo calculation and vector subscripting for cyclic indexing.  This reduces possible errors from hard coded indexing.  Something like this:
> 
>     theta = 20
>     lons = 10.0 * ispan (0,35,1)
>     offset = mod (lons + 360 - theta, 360)
> 
>     iq1 = ind (offset .le. 90)
>     iq2 = ind (offset .ge. 90 .and. offset .le. 180)
>     iq3 = ind (offset .ge. 180 .and. offset .le. 270)
>     iq4 = ind (offset .ge. 270)
> 
> Note that this can create overlapping index points in some cases.  For no overlap, change .le. to .lt. in 3 places.
> 
> In the mod expression, it is important to add an extra 360 degrees, to prevent negative values in the first mod argument.  What this really does is convert the NCL mod function to a true mathematical modulo function for cyclic applications.  Understand the difference here:
> 
>     https://en.wikipedia.org/wiki/Modulo_operation <https://en.wikipedia.org/wiki/Modulo_operation>
> 
> --Dave
> 
> 
> On Tue, Oct 17, 2017 at 1:19 PM, Prashanth Bhalachandran <prashanth.bhalachandran at gmail.com <mailto:prashanth.bhalachandran at gmail.com>> wrote:
> Hello all, 
> I have a relatively simple question on array indexing. 
> 
> Say, I have an array that is cyclic that contains values from 0 to 360 in steps of 5 degrees (dimension [72]). And I am trying to divide into four quadrants starting from a given theta. Say that the index of the given theta is 13 (out of 73), I define my quadrants that are each 90 degrees as follows : 
> 
> Array : theta [73]
> Quadrant 1 : theta(13:31)
> Quadrant 2: theta(31:49)
> Quadrant 3: theta(49:67)
> Quadrant 4: theta(67:13)
> 
> Now, as you’ll notice, having an index of 67:13 doesn’t imply 67:72 and 0:13, it simply implies 13:67 backwards. How do I index the array such that it takes it in a cyclic fashion? For example (what I don’t want is), in a simpler array of just five dimensions : Variable: a
> Type: integer
> Total Size: 20 bytes
>             5 values
> Number of Dimensions: 1
> Dimensions and sizes:	[5]
> Coordinates: 
> (0)	1
> (1)	2
> (2)	3
> (3)	4
> (4)	5
> 
> ncl 3> print(a(4:2)) 
> 
> 
> Variable: a (subsection)
> Type: integer
> Total Size: 12 bytes
>             3 values
> Number of Dimensions: 1
> Dimensions and sizes:	[3]
> Coordinates: 
> (0)	5
> (1)	4
> (2)	3
> 
> What I would have liked it is for it to display 5 1 2 3. 
> 
> Please let me know how this is possible. 
> 
> Thank you all, 
> Prashanth

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171017/5b426a46/attachment.html>