Updated

    :lowres  => 350.,

    :highres => 150.)

job = :highres

job = :all

chosen_epochs = dict_chosen_epochs[job]

Nloop = 6

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

end

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

OIII_norm_evol = fill(NaN, length(chosen_epochs))

for loop in 1:Nloop

    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)

    println(res.bestfit.cost / res.bestfit.dof)

end

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 all_scale_evol[epoch, :] "w lp t '$epoch'"

end

OIII_norm_evol = OIII_norm_evol[:,2:end]

@gp :aa "set grid" :-

@gp "set grid" :-

for id in 1:length(chosen_epochs)

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

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

end