1 | ps = 0; eps = 1; |
---|
2 | |
---|
3 | if (!ps) set terminal x11; |
---|
4 | if (ps && !eps) \ |
---|
5 | set terminal postscript enhanced color solid lw 2 "Times-Roman" 20; |
---|
6 | if (ps && eps) \ |
---|
7 | set terminal postscript eps enhanced color solid lw 2 "Times-Roman" 20; |
---|
8 | |
---|
9 | set grid; |
---|
10 | |
---|
11 | set style line 1 lt 1 lw 1 lc rgb "blue"; |
---|
12 | set style line 2 lt 1 lw 1 lc rgb "green"; |
---|
13 | set style line 3 lt 1 lw 1 lc rgb "red"; |
---|
14 | |
---|
15 | if (ps) set output "beam_envelope_1.ps" |
---|
16 | set title "{/Symbol e}^*_x"; |
---|
17 | set xlabel "s [m]"; set ylabel "[nm{/Symbol \327}rad]"; |
---|
18 | set yrange [0:]; |
---|
19 | set y2range [-1.5:20]; |
---|
20 | |
---|
21 | plot "linlat.out" using 3:4 axis x1y2 notitle with fsteps lt 1 lw 1 \ |
---|
22 | lc rgb "black", \ |
---|
23 | "beam_envelope.out" using 3:(1e9*sqrt($4*$5-$6*$6)) \ |
---|
24 | notitle with lines ls 1; |
---|
25 | |
---|
26 | if (!ps) pause -1; |
---|
27 | |
---|
28 | if (ps) set output "beam_envelope_2.ps" |
---|
29 | set title "{/Symbol e}^*_y"; |
---|
30 | set xlabel "s [m]"; set ylabel "[nm{/Symbol \327}rad]"; |
---|
31 | set yrange [0:]; |
---|
32 | set y2range [-1.5:20]; |
---|
33 | |
---|
34 | plot "linlat.out" using 3:4 axis x1y2 notitle with fsteps lt 1 lw 1 \ |
---|
35 | lc rgb "black", \ |
---|
36 | "beam_envelope.out" using 3:(1e9*sqrt($7*$8-$9*$9)) notitle \ |
---|
37 | with lines ls 1; |
---|
38 | |
---|
39 | if (!ps) pause -1; |
---|
40 | |
---|
41 | if (ps) set output "beam_envelope_3.ps" |
---|
42 | set title "{/Symbol s}_x"; |
---|
43 | set xlabel "s [m]"; set ylabel "[{/Symbol m}m]"; |
---|
44 | set yrange [0:]; |
---|
45 | set y2range [-1.5:20]; |
---|
46 | |
---|
47 | plot "linlat.out" using 3:4 axis x1y2 notitle with fsteps lt 1 lw 1 \ |
---|
48 | lc rgb "black", \ |
---|
49 | "beam_envelope.out" using 3:(1e6*sqrt($4)) notitle with lines ls 1; |
---|
50 | |
---|
51 | if (!ps) pause -1; |
---|
52 | |
---|
53 | if (ps) set output "beam_envelope_4.ps" |
---|
54 | set title "{/Symbol s}_{p_x}"; |
---|
55 | set xlabel "s [m]"; set ylabel "[{/Symbol m}rad]"; |
---|
56 | set yrange [0:]; |
---|
57 | set y2range [-1.5:20]; |
---|
58 | |
---|
59 | plot "linlat.out" using 3:4 axis x1y2 notitle with fsteps lt 1 lw 1 \ |
---|
60 | lc rgb "black", \ |
---|
61 | "beam_envelope.out" using 3:(1e6*sqrt($5)) notitle with lines ls 1; |
---|
62 | |
---|
63 | if (!ps) pause -1; |
---|
64 | |
---|
65 | if (ps) set output "beam_envelope_5.ps" |
---|
66 | set title "{/Symbol s}_y"; |
---|
67 | set xlabel "s [m]"; set ylabel "[{/Symbol m}m]"; |
---|
68 | set yrange [0:]; |
---|
69 | set y2range [-1.5:20]; |
---|
70 | |
---|
71 | plot "linlat.out" using 3:4 axis x1y2 notitle with fsteps lt 1 lw 1 \ |
---|
72 | lc rgb "black", \ |
---|
73 | "beam_envelope.out" using 3:(1e6*sqrt($7)) notitle with lines ls 3; |
---|
74 | |
---|
75 | if (!ps) pause -1; |
---|
76 | |
---|
77 | if (ps) set output "beam_envelope_6.ps" |
---|
78 | set title "{/Symbol s}_{p_y}"; |
---|
79 | set xlabel "s [m]"; set ylabel "[{/Symbol m}rad]"; |
---|
80 | set yrange [0:]; |
---|
81 | set y2range [-1.5:20]; |
---|
82 | |
---|
83 | plot "linlat.out" using 3:4 axis x1y2 notitle with fsteps lt 1 lw 1 \ |
---|
84 | lc rgb "black", \ |
---|
85 | "beam_envelope.out" using 3:(1e6*sqrt($8)) notitle with lines ls 3; |
---|
86 | |
---|
87 | if (!ps) pause -1; |
---|
88 | |
---|
89 | if (ps) set output "beam_envelope_7.ps" |
---|
90 | set title "Transverse Coupling Angle"; |
---|
91 | set xlabel "s [m]"; set ylabel "[{/Symbol \260}]"; |
---|
92 | set autoscale y; |
---|
93 | set y2range [-1.5:20]; |
---|
94 | |
---|
95 | plot "linlat.out" using 3:4 axis x1y2 notitle with fsteps lt 1 lw 1 \ |
---|
96 | lc rgb "black", \ |
---|
97 | "beam_envelope.out" using 3:10 notitle with lines ls 1; |
---|
98 | |
---|
99 | if (!ps) pause -1; |
---|