PRO merge_ucles,  file, nw, bad, wm,fm, wb=wb,fn=fn,fs=fs,iraf=iraf


; Examples


; merge_ucles,'JL87_UCLES_99_echomop.asc',1024, [12,14,15], wm,fm
; merge_ucles,'JL87_UCLES_02_s.cr.dat',2048, [],wm,fm,/iraf


  If KeyWord_Set(iraf) Then BEGIN
     ucles_read,     file,nw,w,f,/iraf
  ENDIF ELSE BEGIN
     ucles_read,     file,nw,w,f
  ENDELSE

  ; only use top layer in 3d echelle spectrum
  ; f = f(*,*,0)

  wc = w
  IF n_elements(bad) GT 0 THEN BEGIN
     ucles_calibcor,nw,w,bad,wc
  ENDIF

  ucles_norm,     nw,wc,f,fn,fs

  ucles_resample, nw,wc,fn,fs, wb,wi

  ucles_merge,    nw,wb,wi, fn,fs, wm,fm



  xs=100.
  window,1
    FOR i = wm(0),wm(n_elements(wm)-1)-xs,2 DO plot,wm,fm,xrange=[i*1.,i*1.+xs],xstyle=1,yrange=[0,1.2]

  window,0
  w0 = wm(0) & wd = ( wm(n_elements(wm)-1) - wm(0) ) / 10
  FOR n=0,9 DO BEGIN 
  ucles_exam, nw, w,f, wm,fm, wb,fs, [w0+n*wd,w0+(n+1)*wd], max(fs)  
  wait, 3
  ENDFOR



END