PRO mktable, lbl, lat, lon, height, width

lbl = strcompress(string(lbl), /remove_all)
lat1 = strcompress(string(lat), /remove_all)
lon1 = strcompress(string(lon), /remove_all)
height1 = strcompress(string(height), /remove_all)
width1 = strcompress(string(width), /remove_all)

cd, 'fits'
spawn, 'ls -I *.fits >> wholelist'
spawn, 'ls -I *data >> bdatcheck'
readcol, 'bdatcheck', ch, format = 'a', /silent
IF n_elements(ch) NE 0 THEN BEGIN
  thereisbaddata = 1
  cd, 'bad_data' & spawn, 'ls -I *.fits >> ../wholelist' & spawn, 'ls -I *.fits >> badlist'
  spawn, 'cp badlist ../badlist'  & spawn, 'cp badlist ../../badlist'
  cd, '..'
ENDIF ELSE thereisbaddata = 0
spawn, 'sort wholelist >> wholelist2'

openw, 90, '../data_table', width = 200

;read in tables for identifiers and variables
readcol, '../'+lbl+'_total_flux.dat', fn, tflux, format = 'a,a', /silent
readcol, '../'+lbl+'_WLSG.dat', fn, WLSG, format = 'a,a', /silent
readcol, '../'+lbl+'_av_NL_grad.dat', fn, avNL, format = 'a,a', /silent
readcol, '../'+lbl+'_eff_dist.dat', fn, effdist, format = 'a,a', /silent
readcol, '../'+lbl+'_NL_length.dat', fn, NLlen, format = 'a,a', /silent
IF thereisbaddata EQ 1 THEN readcol, 'bad_data/badlist', blist, format = 'a', /silent
readcol, 'wholelist2', fits, format = 'a', /silent

date_tabl = strarr(n_elements(fits))
time_tabl = strarr(n_elements(fits))
neglect = intarr(n_elements(fits))
neg1 = [-1]

FOR i1 = 0, n_elements(fits) - 1 DO BEGIN
  IF thereisbaddata EQ 1 THEN BEGIN
    FOR i2 = 0, n_elements(blist) - 1 DO BEGIN
      IF fits(i1) EQ blist(i2) THEN neglect(i1) = 1 ;the fitsfile names for all and bad match
    ENDFOR ; i2 loop
  ENDIF
ENDFOR ;i1 loop


j1 = 0

FOR i1 = 0, n_elements(fits) -1 DO BEGIN
  fitsdata = readfits(fits(i1),hd, /silent)
  IF fitsdata EQ neg1 THEN BEGIN ; was in the bad data section, so have to reread to keep order
     fits(i1) = string('bad_data/' + fits(i1))
     fitsdata = readfits(fits(i1), hd, /silent)
  ENDIF
  date_tabl = strmid(fxpar(hd, 'T_OBS'),0,10)
  time_tabl = strmid(fxpar(hd,'T_OBS'),11,8)
  fits_id = strmid(fxpar(hd,'DATAFILE'),13,9)

  IF neglect(i1) EQ 1 THEN BEGIN
    printf, 90, lbl, ',', date_tabl, ',', time_tabl, ',', fits_id, ',', lat1, ',', lon1, ',', $ 
      height1, ',', width1, ',', 'N/A', ',', 'N/A', ',', 'N/A', ',', 'N/A', ',', 'N/A'
  ENDIF ELSE BEGIN
    printf, 90, lbl, ',', date_tabl, ',', time_tabl, ',', fits_id, ',', lat1, ',', lon1, ',', $ 
     height1, ',', width1, ',', tflux(j1), ',', WLSG(j1), ',', avNL(j1), ',', effdist(j1), ',',$
     NLlen(j1)
  ENDELSE
IF j1 LT (n_elements(WLSG) -1) and neglect(i1) NE 1 THEN j1 = j1 +1 
ENDFOR ;i1 loop

close, 90
END