Hi will,
I think I found a bug(small) in Ec.
In the subroutine - Ecfit_strips_peaks which is called onece for each
view (U, V or W), each stack (inner or outter) and for each sector (1...6)
there is a countter 'ith', wich counts strips with energy less then defined
threshold 'ith=ith+1'. Therefore it should be initialized at the beginning
and should run from 1 to 36. But since there was not initialisation, it
grows with each call and arrea where ID_numbers of low energy strips are stored
gave error - out of bounds.
So I put one line at the beginning of the subroutine:
ith=0
It works.
Arn suggested to check it out and fix that. What you think (and Sergei).
regards stepan