Theory_tutorial/wf.gnu at master · AlessandroPastore81/Theory_tutorial · 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
#Print the first n eigenstates for proton and neutron
#Number of states ns
ns = 5
set terminal pdf
set output 'Neutron_wave_function.pdf'
set xlabel "r[fm]"
do for [j=0:ns]{
s=j+1
set multiplot
set ylabel "Neutron Wave-function n=".s.''
#Uncomment 'wfs_test' part if IW
plot 'neutron_wfs.dat' i j pt 2 lt rgb "red" title 'Neutron Wave function n='.s.'',#'wfs_test.dat' i j pt 1 lt rgb "purple" title 'Test function n='.s.'', '' u 1:(0) w l
unset multiplot}
reset

set terminal pdf
set output 'Proton_wave_function.pdf'
set xlabel "r[fm]"
do for [j=0:ns]{
s=j+1
set multiplot
set ylabel "Proton Wave-function n=".s.''
plot "protons_wfs.dat" i j pt 2 lt rgb "red" title 'Proton Wave function n='.s.''
unset multiplot}
reset

set terminal pdf
set output 'DensityN.pdf'
set xlabel "r[fm]"
set multiplot
set ylabel "Neutron Particle density [fm^3]"
plot "DensityN.dat" u 1:2 pt 2 lt rgb "red" title 'Neutron Particle Density'
unset multiplot
set multiplot
set ylabel "Neutron Kinetic density [fm^-5]"
plot "DensityN.dat" u 1:3 pt 2 lt rgb "red" title 'Neutron Kinetic Density'
unset multiplot
set multiplot
set ylabel "Neutron Spin-orbit density [fm^-4]"
plot "DensityN.dat" u 1:4 pt 2 lt rgb "red" title 'Neutron Spin-orbit Density'
unset multiplot
reset


set terminal pdf
set output 'DensityP.pdf'
set xlabel "r[fm]"
set multiplot
set ylabel "Proton Particle density [fm^-3]"
plot "DensityP.dat" u 1:2 pt 2 lt rgb "red" title 'Proton Particle Density'
unset multiplot
set multiplot
set ylabel "Proton Kinetic density [fm^-5]"
plot "DensityP.dat" u 1:3 pt 2 lt rgb "red" title 'Proton Kinetic Density'
unset multiplot
set multiplot
set ylabel "Proton Spin-orbit density [fm^-4]"
plot "DensityP.dat" u 1:4 pt 2 lt rgb "red" title 'Proton Spin-orbit Density'
unset multiplot
reset


set terminal pdf
set output 'PotentialN.pdf'
set xlabel "r[fm]"
set multiplot
set ylabel "Neutron V(r) Potential"
plot "neutron_potential.dat" u 1:2 pt 2 lt rgb "red" title 'V(r)'
unset multiplot
set multiplot
set ylabel "Neutron Spin-orbit Potential"
plot "neutron_potential.dat" u 1:3 pt 2 lt rgb "red" title 'Vso(r)'
unset multiplot
set multiplot
set ylabel "Neutron Total Potential"
plot "neutron_potential.dat" u 1:($2+$3) pt 2 lt rgb "red" title 'Vtot(r)'
unset multiplot
reset

set terminal pdf
set output 'PotentialP.pdf'
set xlabel "r[fm]"
set multiplot
set ylabel "Proton V(r) Potential"
plot "proton_potential.dat" u 1:2 pt 2 lt rgb "red" title 'V(r)'
unset multiplot
set multiplot
set ylabel "Proton Spin-orbit Potential"
plot "proton_potential.dat" u 1:3 pt 2 lt rgb "red" title 'Vso(r)'
unset multiplot
set multiplot
set ylabel "Proton Coulomb Potential"
plot "proton_potential.dat" u 1:4 pt 2 lt rgb "red" title 'Vc(r)'
unset multiplot
set multiplot
set ylabel "Proton Total Potential"
plot "proton_potential.dat" u 1:($2+$3+$4) pt 2 lt rgb "red" title 'Vtot(r)'
unset multiplot
reset

set terminal pdf
set output 'eigenn.pdf'
set multiplot
set key off
unset xtics
set ylabel "Energies [MeV]"
plot 'neutron_singleparticles.dat' using (-1):($1):(2):(0) with vectors nohead lw 4 lt rgb "red"
unset multiplot
reset

set terminal pdf
set output 'eigenp.pdf'
set multiplot
set key off
unset xtics
set ylabel "Energies [MeV]"
plot 'proton_singleparticles.dat' using (-1):($1):(2):(0) with vectors nohead lw 4 lt rgb "red"
unset multiplot
reset