Commit 29aab2d6 authored by Giorgio Calderone's avatar Giorgio Calderone
Browse files

Minor changes

parent 06baedf2

run.jl

+19 −2
Original line number Diff line number Diff line
@@ -79,7 +79,7 @@ save(:allepochs, term=term_hw, output="$(opt.path)/all_epochs_zoom.png") 
opt.subsets[:lowres]  = [8,9,10,11,12,14,15,16,18,21,22,23,24,25]

opt.subsets[:highres] = [5, 7, 13, 17, 20]

opt.subsets[:test]    = [7, 13]

opt.subsets[:floyds]  = [6, 8, 9, 10, 11, 12, 14, 15, 16, 18, 19, 21, 22, 23]

opt.subsets[:floyds]  = [6, 8, 9, 10, 11, 12, 14, 15, 16, 18, 19, 21, 22, 23] #, 24, 25, 26]



multi_epoch_analysis!(opt, epochs, :test)

multi_epoch_analysis!(opt, epochs, :highres)
@@ -179,6 +179,10 @@ save(term=term_fw, output="$(opt.path)/evolution_oiii_norm_zoom.png") 
@gp :- tab.contslope tab.Hb             "w l notit   lc rgb 'gray'"

@gp :- tab.contslope tab.Ha 1:nrow(tab) "w p t 'H{/Symbol a} int. lum.'  pt 1 ps 3 lc var" ylog=true

@gp :- tab.contslope tab.Hb 1:nrow(tab) "w p t 'H{/Symbol b} int. lum.'  pt 2 ps 3 lc var"

i = findfirst(tab.date .== "20180727")

@gp :- tab.contslope[i] tab.Ha[i]       "w p notit pt 6 ps 4 lc rgb 'black'"

@gp :- tab.contslope[i] tab.Hb[i]       "w p notit pt 6 ps 4 lc rgb 'black'"

@gp:- yr=[4, 50]

save(term=term_fw, output="$(opt.path)/slope_vs_lines.png")
@@ -187,7 +191,10 @@ save(term=term_fw, output="$(opt.path)/slope_vs_lines.png") 
@gp :- tab.contslope tab.c5100             "w l notit dt 2  lc rgb 'gray'"

@gp :- tab.contslope tab.c3500 1:nrow(tab) "w p t 'Cont. {/Symbol l}L_{/Symbol l} at 3500A'  pt 1 ps 3 lc var"

@gp :- tab.contslope tab.c5100 1:nrow(tab) "w p t 'Cont. {/Symbol l}L_{/Symbol l} at 5100A'  pt 2 ps 3 lc var"

@gp :- yr=[500, 3000]

@gp :- yr=[300, 3000]

i = findfirst(tab.date .== "20180727")

@gp :- tab.contslope[i] tab.c3500[i]       "w p notit pt 6 ps 4 lc rgb 'black'"

@gp :- tab.contslope[i] tab.c5100[i]       "w p notit pt 6 ps 4 lc rgb 'black'"

save(term=term_fw, output="$(opt.path)/slope_vs_cont.png")
@@ -196,6 +203,10 @@ save(term=term_fw, output="$(opt.path)/slope_vs_cont.png") 
@gp :- tab.c5100 tab.Ha             "w l notit dt 2  lc rgb 'gray'"

@gp :- tab.c3500 tab.Ha 1:nrow(tab) "w p t 'Cont. {/Symbol l}L_{/Symbol l} at 3500A' pt 1 ps 3 lc var"

@gp :- tab.c5100 tab.Ha 1:nrow(tab) "w p t 'Cont. {/Symbol l}L_{/Symbol l} at 5100A' pt 2 ps 3 lc var"

@gp :- yr=[10,50] xr=[400, 2.5e3]

i = findfirst(tab.date .== "20180727")

@gp :- tab.c3500[i] tab.Ha[i]       "w p notit pt 6 ps 4 lc rgb 'black'"

@gp :- tab.c5100[i] tab.Ha[i]       "w p notit pt 6 ps 4 lc rgb 'black'"

save(term=term_fw, output="$(opt.path)/cont_vs_lines.png")
@@ -208,6 +219,12 @@ save(term=term_fw, output="$(opt.path)/cont_vs_lines.png") 
@gp :- tab.day tab.c5100 1:nrow(tab) "w p t 'Cont. {/Symbol l}L_{/Symbol l} at 5100A' pt 2 ps 3 lc var"

@gp :- tab.day tab.Ha    1:nrow(tab) "w p t 'H{/Symbol a} int. lum.' pt 3 ps 3 lc var"

@gp :- tab.day tab.Hb    1:nrow(tab) "w p t 'H{/Symbol b} int. lum.' pt 4 ps 3 lc var"

@gp :- yr=[3, 1e4]

i = findfirst(tab.date .== "20180727")

@gp :- tab.day[[i,i]] tab.c3500[[i,i]]       "w p notit pt 6 ps 4 lc rgb 'black'"

@gp :- tab.day[[i,i]] tab.c5100[[i,i]]       "w p notit pt 6 ps 4 lc rgb 'black'"

@gp :- tab.day[[i,i]] tab.Ha[[i,i]]          "w p notit pt 6 ps 4 lc rgb 'black'"

@gp :- tab.day[[i,i]] tab.Hb[[i,i]]          "w p notit pt 6 ps 4 lc rgb 'black'"

save(term=term_fw, output="$(opt.path)/evolution_day.png")