PRO fitsdata                ;; Read in WHT driftmode data from FITS files

COMMON fileblock
COMMON fitsblock


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

; Loop over data files
    FOR ifl = 0,nruns(ib)-1 DO BEGIN
     fitsname = folder(ib)+'/w'+folder(ib)+'_'+string(ifl+1,format='(i3)')+'_1d.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)

; Extract the date
     time(it) = FXPAR(fitshead, 'HJD')

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


    ENDFOR
   ENDFOR


; Establish wavelength data array
  w0 =  FXPAR(fitshead, 'CRVAL1')
  dw =  FXPAR(fitshead, 'CD1_1')
  FOR iw = 0,nwave-1 DO BEGIN
   wave(iw) = w0 + iw * dw
  ENDFOR


END