*** PLACE TO STORE BITS REMOVED FROM FF_SUMMARY ; CHUNK 1 head = headfits(cfile[j]) history = sxpar(head, 'HISTORY') if ((flg = hist_find('Saturation removal', history))) ne -1 $ then f1 = (['-', 'S'])[flg gt 0] $ else if sxispar(head, 'BIAS') then $ f1 = (['-', 'S'])[sxpar(head, 'BIAS') gt 0] $ else f1 = '.' if ((flg = hist_find('Large-scale flatfields', history))) ne $ -1 then $ f1 = f1+(['-', 'L'])[flg gt 0] $ else if sxispar(head, 'LSFF') then $ f1 = f1+(['-', 'L'])[sxpar(head, 'LSFF') gt 0] $ else f1 = f1+'.' if ((flg = hist_find('Dark offset subtraction', history))) ne $ -1 then $ f1 = f1+(['-', 'D'])[flg gt 0] $ else if sxispar(head, 'DARK') then $ f1 = f1+(['-', 'D'])[sxpar(head, 'DARK') ne 0] $ else f1 = f1+'.' if ((flg = hist_find('Spike removal', history))) ne $ -1 then $ f1 = f1+(['-', 'C'])[flg gt 0] $ else if sxispar(head, 'CRAYS') then $ f1 = f1+(['-', 'C'])[sxpar(head, 'CRAYS') ne 0] $ else f1 = f1+'.' if ((flg = hist_find('Moon and planet removal', history))) $ ne -1 then $ f1 = f1+(['-', 'B'])[flg gt 0] $ else if sxispar(head, 'BRIGHT') then $ f1 = f1+(['-', 'B'])[sxpar(head, 'BRIGHT') gt 0]$ else f1 = f1+'.' if ((flg = hist_find('Camera equalization', history))) $ ne -1 then $ f1 = f1+(['-', 'E'])[flg gt 0] $ else if sxispar(head, 'GAIN') then $ f1 = f1+(['-', 'G'])[sxpar(head, 'GAIN') gt 0]$ else f1 = f1+'.' if ((flg = hist_find('Pointing', history))) ne $ -1 then $ f1 = f1+(['-', 'L'])[flg gt 1] $ else if sxispar(head, 'LPOINT') then $ f1 = f1+(['-', 'L'])[sxpar(head, 'LPOINT') gt 0] $ else f1 = f1+'.' if ((flg = hist_find('Pointing', history))) ne $ -1 then $ f1 = f1+(['-', 'G'])[flg gt 0] $ else if sxispar(head, 'GPOINT') then $ f1 = f1+(['-', 'G'])[sxpar(head, 'GPOINT') gt 0] $ else f1 = f1+'.' if ((flg = hist_find('Zodiacal subtraction', history))) ne $ -1 then $ f1 = f1+(['-', 'Z'])[flg gt 0] $ else if sxispar(head, 'ZODY') then $ f1 = f1+(['-', 'Z'])[sxpar(head, 'ZODY') gt 0] $ else f1 = f1+'.' if ((flg = hist_find('Star subtraction', history))) ne $ -1 then $ f1 = f1+(['-', 'S'])[flg gt 0] $ else if sxispar(head, 'STELLAR') then $ f1 = f1+(['-', 'S'])[sxpar(head, 'STELLAR') gt 0] $ else if sxispar(head, 'BSTARS') then $ f1 = f1+(['-', 'B'])[sxpar(head, 'BSTARS') gt 0] $ else f1 = f1+'.' flg = -1 for icam = 1, 3 do begin mstr = 'Camera '+string(icam, format = '(I1)')+' bad pixel ' + $ 'mask' flg = flg > hist_find(mstr, history, /bool, /line_break) endfor if flg ne -1 then f1 = f1+(['-', 'H'])[flg gt 0] $ else if sxispar(head, 'FLIPPERS') then $ f1 = f1+(['-', 'H'])[sxpar(head, 'FLIPPERS') gt 0] $ else f1 = f1+'.' ; This messy stuff is because early SMEI Aitoffs have several ; different ways of storing their dates. if sxispar(head, 'START_UT') then begin jd0 = epoch2jd(sxpar(head, 'START_UT')) date = jd2doy(jd0) h0 = date[2] m0 = date[3] s0 = date[4] d = 0 yr = date[0] day = date[1] if sxispar(head, 'STOP_UT') then begin jd1 = epoch2jd(sxpar(head, 'STOP_UT')) d = (jd1-jd0)*1440. endif else d = 0 ; if d lt 0 then d = d+1440. endif else if sxispar(head, 'MJD-OBS') then begin jd0 = 2400000.d0 + sxpar(head, 'MJD-OBS') date = jd2doy(jd0) h0 = date[2] m0 = date[3] s0 = date[4] d = 0 yr = date[0] day = date[1] endif else if sxispar(head, 'OBS_MJD') then begin jd0 = 2400000.d0 + sxpar(head, 'OBS_MJD') date = jd2doy(jd0) h0 = date[2] m0 = date[3] s0 = date[4] d = 0 yr = date[0] day = date[1] endif else begin h0 = 0 m0 = 0 s0 = 0 d = 0 yr = 0 day = 0 endelse if sxpar(head, 'NAXIS') eq 3 then npl = sxpar(head, 'NAXIS3') $ else npl = 1 case strtrim(sxpar(head, 'CTYPE1'), 2) of 'ELON-AIT': project = 'Aitoff ' 'ELON-ARC': project = 'Fisheye' 'PA-RECT': project = 'Rectang' else: project = 'Unknown' endcase desc[j] = string(yr, day, h0, m0, d, $ sxpar(head, 'NAXIS1'), sxpar(head, 'NAXIS2'), $ npl, f1, project, $ format = $ "(I4.4,' ',I3.3,' ',I2.2,':',I2.2," + $ "' ',F6.1,' ',I4,'x',I4,' ',I1,' ',A,' ',A)") ; CHUNK 2 yr = 0l & day = 0l h0 = 0l & m0 = 0l & s0 = 0l d = 0. & nax1 = 0l & nax2 = 0l & nax3 = 0 & f1p = '' reads, hdesc[locs[j]], yr, day, h0, m0, d, nax1, nax2, nax3, $ f1p, format = $ "(I5,I4,I3,1x,I2,F7.0,I5,1x,I4,1x,i1,1x,A)" f1ps = strsplit(f1p, ' ', /extract, count = nstr) ; help, f1p, f1ps ; stop f1 = f1ps[0] if nstr eq 2 then project = f1ps[1] $ else begin smei_msg, /warn, $ "Obsolete img_hdr.txt found, remake it to get " + $ "projection information" project = " " endelse