forked from AlessandroPastore81/Theory_tutorial
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathwf.gnu
109 lines (100 loc) · 3.01 KB
/
wf.gnu
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
#Print the first n eigenstates for proton and neutron
#Number of states ns
ns = 5
set terminal pdf
set output 'wfn.pdf'
set xlabel "r[fm]"
do for [j=0:ns]{
s=j+1
set multiplot
set ylabel "Neutron Wave-function n=".s.''
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 'wfp.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 'spe.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