PRO brucedata                ;; Read in BRUCE models from FITS files

COMMON fileblock
COMMON fitsblock
COMMON procblock


fitsfile=' '
htime = 0.0

; Loop over data folders
   FOR ib = 0,nfiles-1 DO BEGIN

; Open the timing file
     timefile = folder(ib)+'/'+flabel(ib)+'0.times'
     openr,1,timefile

; Loop over data files
    FOR ifl = 0,nruns(ib)-1 DO BEGIN
     fitsname = folder(ib)+'/'+flabel(ib)+string(ifl+1,format='(i4.4)')+'.fits'
     IF nruns(ib) GT 999 THEN fitsname = folder(ib)+'/'+flabel(ib)+string(ifl+1,format='(i4.4)')+'.fits'
     fitsname = strtrim(strcompress(fitsname, /remove_all),2)
     IF ib EQ 0 THEN it = ifl
     IF ib GT 0 THEN it = total(nruns(0:ib-1))+ifl

; Read the fits file    
     fitsdata = readfits(fitsname, fitshead, /silent)

; Read the times from the timing file
      lablen = strlen( flabel(ib)+string(ifl+1,format='(i3.3)') )
      IF nruns(ib) GT 999 THEN lablen = strlen( flabel(ib)+string(ifl+1,format='(i4.4)') )
      stime = " "
      readf,1,stime
      slen = strlen ( stime )
      reads, strmid( stime, lablen+2 ), htime
;;      readf,1, htime, format='(<lablen>x,F6.3)'
      time(it) = htime / 24.0

; Copy the data
  nwave =  FXPAR(fitshead, 'NAXIS1')

  ; define the arrays if they do not exist
  if (n_elements(spec) EQ 0 OR n_elements(wave) NE nwave) THEN BEGIN
    wave   = fltarr(nwave)
    spec   = fltarr(nwave,nspec)
    specm  = fltarr(nwave)
    specf  = fltarr(nwave,nspec)
    specr  = fltarr(nwave,nspec)
  ENDIF

  IF ifl GE nbomit THEN spec(0:nwave-1,it) = fitsdata(0:nwave-1,0,0)


    ENDFOR


; Close the timing file
     close,1


   ENDFOR


; Establish wavelength data array
  w0 =  FXPAR(fitshead, 'CRVAL1')
  dw =  FXPAR(fitshead, 'CD1_1')
  IF ( abs(nint(dw*10)/10. - dw) lt 0.001 ) THEN dw = nint(dw*10)/10.
  ; this was a fudge because ndf2fits corrupted the original data
  p0 =  FXPAR(fitshead, 'CRPIX1') - 1.
  FOR iw = 0,nwave-1 DO BEGIN
   wave(iw) = w0 + (iw -p0) * dw
  ENDFOR

END