Updated

    :lowres  => 350.,

    :highres => 150.)

job = :all

job = :highres

chosen_epochs = dict_chosen_epochs[job]

Nloop = 6

(all_scale, all_fwhm) = deserialize("output/scale_fwhm.dat")

for loop in 1:Nloop

    if !isfile("output/results_$(job)_$(loop).dat")

    file = "output/results_$(job)_$(loop).dat"

    if !isfile(file)

        source = QSO{q1927p654}("1ES 1927+654", Z, ebv=EBV);

        source.options[:min_spectral_coverage][:OIII_5007] = 0.5

        @gp :zoom "set grid" :-

            OIII_norm[id] = bestfit2[id][:OIII_5007].norm.val

            all_scale[chosen_epochs[id]] *= bestfit2[id][:OIII_5007].norm.val

        end

        serialize("output/results_$(job)_$(loop).dat", (res, all_scale, OIII_norm))

        serialize(file, (res, all_scale, OIII_norm))

    else

        (res, all_scale, OIII_norm) = deserialize("output/results_$(job)_$(loop).dat")

        (res, all_scale, OIII_norm) = deserialize(file)

    end

end

(scale_evol, _) = deserialize("output/scale_fwhm.dat")

(all_scale_evol, _) = deserialize("output/scale_fwhm.dat")

OIII_norm_evol = fill(NaN, length(chosen_epochs))

for loop in 1:Nloop

    ss = deserialize("output/scale_$(job)_$(loop).dat")

    scale_evol = hcat(scale_evol, ss)

    file = "output/results_$(job)_$(loop).dat"

    (res, all_scale, OIII_norm) = deserialize(file);

    all_scale_evol = hcat(all_scale_evol, all_scale)

    OIII_norm_evol = hcat(OIII_norm_evol, OIII_norm)

end

scale_evol = 2 .*

    (scale_evol[:, 2:end] .- scale_evol[:, 1:end-1]) ./

    (scale_evol[:, 2:end] .+ scale_evol[:, 1:end-1])

all_scale_evol = 2 .*

    (all_scale_evol[:, 2:end] .- all_scale_evol[:, 1:end-1]) ./

    (all_scale_evol[:, 2:end] .+ all_scale_evol[:, 1:end-1])

@gp "set grid" :-

for epoch in chosen_epochs

    @gp :- 1:Nloop scale_evol[epoch, :] "w lp t '$epoch'"

    @gp :- 1:Nloop all_scale_evol[epoch, :] "w lp t '$epoch'"

end

OIII_norm_evol = OIII_norm_evol[:,2:end]

@gp :aa "set grid" :-

for id in 1:length(chosen_epochs)

    @gp :- :aa 1:Nloop OIII_norm_evol[epoch, :] "w lp t '$(chosen_epochs[id])'"

end