Select Git revision
Code owners
Assign users and groups as approvers for specific file changes. Learn more.
dnu.plt 7.04 KiB
ps = 0; eps = 0; action_angle = 0; pert = 0;
#N = 1; N_x = 33; N_y = 16;
N = 3; N_x = 11; N_y = 5;
# APS
#N = 40; N_x = 0; N_y = 0;
#N = 15; N_x = 2; N_y = 1;
#N = 3; N_x = 6; N_y = 2;
# MAX-IV
#N = 20; N_x = 2; N_y = 0;
font_size = 30; line_width = 2;
if (!ps) set terminal x11;
if (ps && !eps) \
set terminal postscript enhanced color solid \
lw line_width "Times-Roman" font_size;
if (ps && eps) \
set terminal postscript eps enhanced color solid \
lw line_width "Times-Roman" font_size;
#set multiplot;
# left adjusted labels
set key Left;
set grid;
set style line 1 lt 1 lw line_width lc rgb "blue";
set style line 2 lt 1 lw line_width lc rgb "dark-green";
set style line 3 lt 1 lw line_width lc rgb "red";
set style line 4 lt 1 lw line_width lc rgb "dark-orange";
set clabel "%5.2f"; set key left;
set palette rgbformulae 22, 13, -31 negative;
if (ps) set output "dnu_1.ps"
#set size 0.5, 0.5; set origin 0.0, 0.5;
set title "{/Symbol n}_x vs. A_{x,y}";
if (!action_angle) \
set xlabel "A_{x,y} [mm]"; set ylabel "{/Symbol n}_x"; \
if (!pert) \
plot "dnu_dAx.out" using 1:(N*(N_x+$5)) title "A_x" \
with linespoints ls 1, \
"dnu_dAy.out" using 2:(N*(N_x+$5)) title "A_y" with linespoints ls 3;
if (pert) \
plot "dnu_dAx.out" using 1:(N*(N_x+$5)) title "A_x" \
with linespoints ls 1, \
"dnu_dAy.out" using 2:(N*(N_x+$5)) title "A_y" \
with linespoints ls 3, \
"dnu_dAx_pert.out" using 1:(N*(N_x+$3)) title "A_x (pert)" \
with linespoints ls 2, \
"dnu_dAy_pert.out" using 2:(N*(N_x+$3)) title "A_y (pert)" \
with linespoints ls 4;
if (action_angle) \
set xlabel "A_{x,y} [mm]"; set ylabel "{/Symbol n}_x"; \
set ytics nomirror; set y2tics; \
plot "dnu_dAx.out" using 3:5 title "J_x" with linespoints ls 1, \
"dnu_dAy.out" using 4:5 title "J_y" with linespoints ls 3;
if (!ps) pause -1;
if (ps) set output "dnu_2.ps"
#set size 0.5, 0.5; set origin 0.5, 0.5;
set title "{/Symbol n}_y vs. A_{x,y}";
if (!action_angle) \
set xlabel "A_{x,y} [mm]"; set ylabel "{/Symbol n}_y"; \
if (!pert) \
plot "dnu_dAx.out" using 1:(N*(N_y+$6)) title "A_x" \
with linespoints ls 1, \
"dnu_dAy.out" using 2:(N*(N_y+$6)) title "A_y" \