gnuplot-examples/pl-hf2s6.gp at main · ac-dias/gnuplot-examples · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132


set terminal pngcairo color size 1000,500 enhanced


set palette rgbformulae 33,13,10 negative


set style line 1 \
    linecolor rgb 'blue' \
    linetype 1 linewidth 3 \
    pointtype 7 pointsize 1.5

set style line 2  \
    linecolor rgb 'red' \
    linetype 1 linewidth 2 \
    pointtype 7 pointsize 1.5

set style line 3 dt "-" \
    linecolor rgb 'black' \
    linetype 1 linewidth 3 \
    pointtype 7 pointsize 1.5

set style line 4  dt "-" \
    linecolor rgb 'dark-green' \
    linetype 1 linewidth 2 \
    pointtype 7 pointsize 1.5


#set ylabel "Absorption Intensity (a.u)" font "Times-Roman,20"
#set xlabel "Photon Energy (eV)" font "Times-Roman,20"


#set yrange[-6:6]

set xrange [0:4]
set xtics 1


#set ytics 0.2

set tics font "Times-Roman,20"
set lmargin 10
set bmargin 4


set output 'pl_hf2s6.png'

#unset key

#normalize data

stats "./Hf2S6/out/bse_absorption_coef.dat" using 2 nooutput

m1 = STATS_max

stats "./Hf2S6/out/bse_absorption_coef.dat" using 3 nooutput

m2 = STATS_max

stats "./Hf2S6/out/ipa_absorption_coef.dat" using 2 nooutput

m3 = STATS_max

stats "./Hf2S6/out/ipa_absorption_coef.dat" using 3 nooutput

m4 = STATS_max


stats "./Hf2S6/out/bse_PL.dat" using 2 nooutput

n1 = STATS_max

stats "./Hf2S6/out/bse_PL.dat" using 3 nooutput

n2 = STATS_max

stats "./Hf2S6/out/ipa_PL.dat" using 2 nooutput

n3 = STATS_max

stats "./Hf2S6/out/ipa_PL.dat" using 3 nooutput

n4 = STATS_max


set multiplot layout 1,2

set yrange [0:1.1]
############################################################################

set xlabel "Photon Energy (eV)" font "Times-Roman,20"
set ylabel "Normalized Intensity (a.u)" font "Times-Roman,20" offset -2,0
set format y "%.1f"      # sempre 2 casas decimais
set label "(a)" at graph 0.02,0.95 front font "Times-Roman,20"
set key at graph 0.43,0.88 font "Times-Roman,15"

plot './Hf2S6/out/bse_absorption_coef.dat' u 1:($2/m1) w lines linestyle 1 title "BSE Ab" \
,'./Hf2S6/out/bse_PL.dat' u 1:($2/n1) w lines linestyle 2 title "BSE PL"\
,'./Hf2S6/out/ipa_absorption_coef.dat' u 1:($2/m3) w lines linestyle 3 title "IPA Ab"\
,'./Hf2S6/out/ipa_PL.dat' u 1:($2/n3) w lines linestyle 4 title "IPA PL"\

unset ylabel
unset label
unset format y
unset key

set format y "%.1f"      # sempre 2 casas decimais
set label "(b)" at graph 0.02,0.95 front font "Times-Roman,20"
set key at graph 0.43,0.88 font "Times-Roman,15"


plot './Hf2S6/out/bse_absorption_coef.dat' u 1:($3/m2) w lines linestyle 1 title "BSE Ab" \
,'./Hf2S6/out/bse_PL.dat' u 1:($3/n2) w lines linestyle 2 title "BSE PL"\
,'./Hf2S6/out/ipa_absorption_coef.dat' u 1:($3/m4) w lines linestyle 3 title "IPA Ab"\
,'./Hf2S6/out/ipa_PL.dat' u 1:($3/n4) w lines linestyle 4 title "IPA PL"\


############################################################################


unset multiplot

q