ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/JSOC/proj/globalhs/scripts/idl/arrsumm.pro
Revision: 1.1
Committed: Fri Apr 23 09:15:53 2010 UTC (13 years, 5 months ago) by tplarson
Branch: MAIN
CVS Tags: Ver_6-0, Ver_6-1, Ver_6-2, Ver_6-3, Ver_6-4, Ver_9-1, Ver_5-14, Ver_5-13, Ver_5-12, Ver_5-11, Ver_5-10, Ver_LATEST, Ver_9-3, Ver_9-41, Ver_9-2, Ver_5-9, Ver_8-8, Ver_8-2, Ver_8-3, Ver_8-0, Ver_8-1, Ver_8-6, Ver_8-7, Ver_8-4, Ver_8-5, Ver_7-1, Ver_7-0, Ver_9-5, Ver_9-4, Ver_8-10, Ver_8-11, Ver_8-12, Ver_9-0, HEAD
Log Message:
idl functions, procedures, scripts, and input files needed by doretidl, dogfidl, and dopkgn

File Contents

# Content
1 FUNCTION ARRSUMM, arr, skip=skip, nsecs=nsecs
2
3 if (n_elements(skip) eq 0) then skip=1
4
5 na=n_elements(arr)
6 arr=long(arr)
7 arr=arr[sort(arr)]
8 arr0=[arr,0]-[0,arr]
9 diff=arr0[1:na-1]
10
11 w1=where(diff gt skip,nw1)
12 nsecs=nw1+1
13 if (nw1 eq 0) then begin
14 return,[arr[0],arr[na-1]]
15 endif
16
17 secs=lonarr(2,nw1+1)
18 secs[*,0]=[arr[0],arr[w1[0]]]
19 for i=0,nw1-2 do begin
20 secs[*,i+1]=[arr[w1[i]+1],arr[w1[i+1]]]
21 endfor
22 secs[*,nw1]=[arr[w1[nw1-1]+1],arr[na-1]]
23
24 return,secs
25
26 END