pro pfi_totalizer, index, data, totals ndat = n_elements(index) si = size(data) npix = si(1)*si(2) totals = fltarr(ndat) for i = 0, ndat-1 do $ totals(i) = total(sxt_decomp(data(*,*,i)))/gt_expdur(index(i))*1000/npix end