Updated

               filename="output/results_$(job)_$(loop)_rebin4.html", rebin=4)

        # Find best normalization for [OIII]

        model2 = deepcopy(res.model);

        multi2 = deepcopy(res.multi);

        for id in 1:length(chosen_epochs)

            for cname in collect(keys(model2[id]))

            for cname in collect(keys(multi2[id]))

                if cname == :OIII_5007

                    model2[id][cname].norm.fixed = false

                    multi2[id][cname].norm.fixed = false

                else

                    freeze(model2[id], cname)

                    freeze(multi2[id], cname)

                end

            end

        end

        mzer = GFit.cmpfit()

        mzer.config.ftol = mzer.config.gtol = mzer.config.xtol = 1.e-6

        bestfit2 = fit!(model2, source.data, minimizer=mzer);

        bestfit2 = fit!(multi2, source.data, minimizer=mzer);

        OIII_norm = fill(0., length(chosen_epochs))

        for id in 1:length(chosen_epochs)