PRO Fixspect,specname

; Load star spectrum for G-type star
OpenR,lun,'\usask\893\crossect\wavelgth.txt',/Get_LUN
starwave=DblArr(1084)
ReadF,lun,starwave
Free_LUN,lun
dwl=Median(starwave-Shift(starwave,1))
;    starflux=ddread('\usask\893\crossect\g1--2v.txt)
    OpenR,lun,'\usask\893\crossect\'+StrTrim(specname,2)+'.txt',/Get_LUN
    starflux=DblArr(1084)
    ReadF,lun,starflux
    Free_LUN,lun

done=0B
starfix=starflux
REPEAT BEGIN
    Plot,starwave,starfix,xrange=[6000,9000]
    Read,'Enter wavelength bracket for bad fluxes (<0 to end): ',low,high
    IF (low+high) LT 0 THEN done=1B ELSE BEGIN
        bad=Where((low LT starwave) AND (starwave LE high) $
                                    AND (starflux LT 0.1))
        good=[Indgen(10)-10+Min(bad),Indgen(10)+1+Max(bad)]
        res=LADFit(starwave(good),starflux(good))
        starfix(bad)=res(0)+res(1)*starwave(bad)
        ENDELSE
    ENDREP UNTIL done

OpenW,lun,'\usask\893\crossect\'+StrTrim(specname,2)+'.txt',/Get_LUN
PrintF,lun,starfix,Format='(10F8.3)'
Free_LUN,lun

Return
END