Updated

    :lowres   => [8,9,10,11,12,14,15,16,18,21,22,23,24,25],

    :highres  => [5, 7, 13, 17, 20],

    :limited  => [1,2,3,13],

    :all      => collect(5:26))

    :all      => collect(1:26))

deleteat!(dict_chosen_epochs[:all], 4)  #insufficient coverage in na_Hg

# Can't use resolution smaller than 150 km / s otherwise some line

# will be neglected because of spectral coverage

    :highres => 150.)

job = :highres

job = :all

Nloop = 6

for loop in 1:Nloop

    chosen_epochs = dict_chosen_epochs[job]

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

@gp "set grid" :-

for loop in 1:Nloop

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

    scale_evol = hcat(scale_evol, ss)

scale_evol = 2 .*

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

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

@gp "set grid" ylog=true :-

for epoch in 1:size(scale_evol)[1]

    @gp :- 1:Nloop scale_evol[epoch, :] "w lp notit"

@gp "set grid" :-

for epoch in chosen_epochs

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

end