pro PATT, l0, klassatt, pat

;+
; NAME:
;       PATT
; PURPOSE:
;       Searching for invalid points during attractor periods
; EXPLANATION:
;
;
;
; CALLING SEQUENCE:
;       PATT, level_0, klassatt, pat
; INPUTS:
;       level_0: level 0 time series in the format (n element vector)
;       klassatt: output of findattract, containing location and strength of attractor values ((3, m) vector)
; OUTPUTS:
;       pat: n vector: valid points 1.0, invalid points !values.f_nan
;
; PROCEDURES USED:
;       None
;
;
; MODIFICATION HISTORY:
;       Package 'SPM level 1.5' written between 2000 and 2003, Richard Wachter, PMOD/WRC
;		was a problem with the histogram needed the /nan in the histogram and not in the max
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

 na=(size(klassatt))[2]

 l=n_elements(l0)

 pat=replicate(1B, l)

 for i=0, na-1 do begin

  if klassatt[2,i] ge 400 then begin

    seq=l0[klassatt[0,i]:klassatt[1,i]]
    h=max(histogram(seq, omin=m, binsize=1L,/nan), ma, /nan)
    maxim=m+ma
    seq0=seq-maxim

    pat[klassatt[0,i]+where(abs(seq0) le 5)]=0B

    endif

  endfor

end