Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.naic.edu/~tghosh/software/coraccum-original.pro Дата изменения: Fri Feb 21 17:53:30 2003 Дата индексирования: Sun Apr 10 23:17:37 2016 Кодировка:

;+
;NAME:
;coraccum - accumulate a record in a summary rec
;SYNTAX: coraccum,binp,badd,new=new,scl=scl,array=array,brd=brd
;ARGS: binp[]: {corget} input data
; badd: {corget} accumulate data here
;KEYWORDS:
; new : if keyword set then this is the first call, alloc badd.
; scl : float if provided, then scale the brec data by scl before
; adding. This can be used to weight data by g/t.
; array : if set then treat badd as an array and add element wise
; the array binp[n] to badd[n]. Use with the new keyword.
; By default if badd is an array then binp will be added
; element wise.
; brd[] : int if specified then only process the specified boards.
; numbering is 1,2,3,4.
;DESCRIPTION:
; Accumulate 1 or more records worth of data into badd. If keyword
;/new is set then allocate badd before adding. The header values in
;badd will come from the first record added into badd.
; The badd.b1.accum variable will be incremented by numrecs(binp)*scl for
;each call. When corplot is called, it will scale the data by 1./badd.b1.accum.
; When calling coraccum with the new keyword, you can include the
;array keyword. This will allocate badd to be the same dimension as
;binp. All future calls using badd will add binp element wise to badd.
;This can be used when accumulating multiple maps.
; Accumulated data must be of the same type (numlags, numbsbc, bw,etc..).
;If you have observations to accumulate with only partial overlap of the
;data types, you can use the brd keyword to specify which boards to accum. The
;accumlated data must still be of the same type.
;
;Example:
;
; print,corget(lun,b)
; coraccum,b,badd,/new
; print,corget(lun,b)
; coraccum,b,badd
; corplot,badd
;
;; Add n scans together element wise:
; for i=0,n-1 do begin
; print,corinpscan(lun,b)
; coraccum,b,bsum,new=(i eq 0),/array
; endfor
;;
;; input an entire scan and then plot the average of the records
;; (this can also be done directly by corinpscan).
; print,corinpscan(lun,b,scan=scan)
; coraccum,b,bsum,/new
; corplot,bsum
;;
;; let scan 1 be:4 brds,2pol, 1024 lags
;; let scan 2 be:2 brds,2 pol,512 lags followed by 2 brds,2 pol, 1024 lags.
;; To accumulate brds 1,2 of scan 1 with brds 3,4 of scan 2:
;
; print,corinpscan(lun,b,scan=scan1)
; coraccum,b,badd,/new,brd=[1,2]
; print,corinpscan(lun,b,scan=scan2)
; coraccum,b,badd,brd=[3,4]
;
; help,badd.b1,/st
;** Structure <3957c8>, 4 tags, length=9200, refs=2:
; H STRUCT -> HDR Array[1]
; P INT Array[2]
; ACCUM DOUBLE 2.00000
; D FLOAT Array[1024, 2]
;-
;history:
; 12aug02 - added brd option
pro coraccum,b,baccum,new=new,scl=scl,array=array,brd=brd
;
on_error,2
; print,"New coraccum"
if n_elements(scl) eq 0 then scl=1.
mult=1.
usebrd=keyword_set(brd)
if keyword_set(new) then begin
if usebrd then begin
if keyword_set(array) then begin
baccum=corsubset(b,brd)
endif else begin
baccum=corsubset(b[0],brd)
endelse
endif else begin
if keyword_set(array) then begin
baccum=b
endif else begin
baccum=b[0]
endelse
endelse
mult=0. ; to zero it out
for i=0,n_tags(baccum[0])-1 do begin
baccum.(i).accum=0.
endfor
endif
;
; check that some things match here
;
numrecsinp=(size(b))[1]
numrecsout=n_elements(baccum)
outarray=0L
if (numrecsout gt 1) then begin
if (numrecsinp ne numrecsout) then $
message,'array accumlation requires matching inp,accum dimensions'
outarray=1L
endif
if usebrd then begin
numbrdsinp=n_elements(brd)
brdind=brd-1
endif else begin
numbrdsinp=n_tags(b)
brdind=indgen(numbrdsinp)
endelse
numbrdsout=n_tags(baccum)
if numbrdsinp ne numbrdsout then $
message,'number of boards does not match, input and accum bufs'
;
; verify that the data types are the same
;
for j=0,numbrdsout-1 do begin
jj=brdind[j]
if baccum[0].(j).h.cor.bwnum ne b[0].(jj).h.cor.bwnum then begin
lin=string(format=$
'("Bandwidth mismatch inpbrd:",i2," accum brd:",i2)',j+1,jj+1)
message,lin
endif
if baccum[0].(j).h.cor.numsbcout ne b[0].(jj).h.cor.numsbcout then begin
lin=string(format=$
'("Numsbc mismatch inpbrd:",i2," accum brd:",i2)',j+1,jj+1)
message,lin
endif
if baccum[0].(j).h.cor.lagsbcout ne b[0].(jj).h.cor.lagsbcout then begin
lin=string(format=$
'("Numlags mismatch inpbrd:",i2," accum brd:",i2)',j+1,jj+1)
message,lin
endif
if baccum[0].(j).h.cor.lagconfig ne b[0].(jj).h.cor.lagconfig then begin
lin=string(format=$
'("lagconfig mismatch inpbrd:",i2," accum brd:",i2)',j+1,jj+1)
message,lin
endif
endfor

for i=0,numrecsinp-1 do begin
for j=0,numbrdsout-1 do begin
jj=brdind[j]
baccum[i*outarray].(j).d=baccum[i*outarray].(j).d*mult + $
b[i].(jj).d*scl
endfor
if outarray eq 0 then mult=1. ;
endfor
if outarray then begin
inc=scl
endif else begin
inc=scl*numrecsinp
endelse
for i=0,n_tags(baccum)-1 do baccum.(i).accum=baccum.(i).accum + inc
return
end