gnuplot-examples/optics-hf2se6.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

set terminal pngcairo color size 1500,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 'blue' \
    linetype 1 linewidth 3 \
    pointtype 7 pointsize 1.5

set style line 4 dt '-'  \
    linecolor rgb 'red' \
    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 tics font "Times-Roman,20"
set lmargin 10
set bmargin 4

set autoscale y

set output 'optics_hf2se6.png'

#unset key


set multiplot layout 1,3


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

set xlabel "Photon Energy (eV)" font "Times-Roman,20"
set ylabel "Abs. Coeff. (10^{4} cm^{-1})" font "Times-Roman,20" offset -1,0

set label "(a)" at graph 0.02,0.95 front font "Times-Roman,20"
set key at graph 0.50,0.88 font "Times-Roman,15"

plot './Hf2Se6/out/bse_absorption_coef.dat' u 1:($2/10000) w lines linestyle 1 title "BSE A_{x,x}" \
,'./Hf2Se6/out/bse_absorption_coef.dat' u 1:($3/10000) w lines linestyle 2 title "BSE A_{y,y}"\
,'./Hf2Se6/out/ipa_absorption_coef.dat' u 1:($2/10000) w lines linestyle 3 title "IPA A_{x,x}"\
,'./Hf2Se6/out/ipa_absorption_coef.dat' u 1:($3/10000) w lines linestyle 4 title "IPA A_{y,y}"\

unset label
unset format y
unset key

set ylabel "Refractive Index" font "Times-Roman,20" offset -3,0
set format y "%.2f"      # sempre 2 casas decimais
set label "(b)" at graph 0.02,0.95 front font "Times-Roman,20"
set key at graph 0.50,0.88 font "Times-Roman,15"

set ytics 1.0

plot './Hf2Se6/out/bse_refractive_index.dat' u 1:2 w lines linestyle 1 title "BSE {/Symbol h}_{x,x}"\
,'./Hf2Se6/out/bse_refractive_index.dat' u 1:3 w lines linestyle 2 title "BSE {/Symbol h}_{y,y}"\
,'./Hf2Se6/out/ipa_refractive_index.dat' u 1:2 w lines linestyle 3 title "IPA {/Symbol h}_{x,x}"\
,'./Hf2Se6/out/ipa_refractive_index.dat' u 1:3 w lines linestyle 4 title "IPA {/Symbol h}_{y,y}"\

unset label
unset key
unset ytics

set ylabel "Reflectivity" font "Times-Roman,20" offset -3,0
set format y "%.2f"      # sempre 2 casas decimais
set label "(c)" at graph 0.02,0.95 front font "Times-Roman,20"
set key at graph 0.50,0.88 font "Times-Roman,15"

set ytics font "Times-Roman,20"
set ytics 0.2

plot './Hf2Se6/out/bse_reflectibility.dat' u 1:2 w lines linestyle 1 title "BSE R_{x,x}"\
,'./Hf2Se6/out/bse_reflectibility.dat' u 1:3 w lines linestyle 2 title "BSE R_{y,y}"\
,'./Hf2Se6/out/ipa_reflectibility.dat' u 1:2 w lines linestyle 3 title "IPA R_{x,x}"\
,'./Hf2Se6/out/ipa_reflectibility.dat' u 1:3 w lines linestyle 4 title "IPA R_{y,y}"\


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


unset multiplot

q