Skip to content
Snippets Groups Projects
Select Git revision
  • a8d1bdb4d73570d67039cbc6f7fdaa46bf16a0a7
  • main default protected
  • 2.12.1
  • 2.6.0
  • 2.4.1
  • 2.0.17
  • 2.0.13
  • 2.0.10
  • 2.0.8
  • 2.0.7
  • 2.0.5
  • 2.0.4
  • 2.0.3
  • 2.0.2
  • 2.0.1
  • 1.59.0
  • 1.57.0
  • 1.55.0
  • 1.54.0
19 results

.gitlab-ci.yml

Blame
  • 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" \
             with linespoints ls 3;
      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" \
             with linespoints ls 3, \
             "dnu_dAx_pert.out" using 1:(N*(N_y+$4)) title "A_x (pert)" \
             with linespoints ls 2, \
             "dnu_dAy_pert.out" using 2:(N*(N_y+$4)) title "A_y (pert)" \
             with linespoints ls 4;
    if (action_angle) \
      set xlabel "A_{x,y} [mm]"; set ylabel "{/Symbol n}_y"; \
      plot "dnu_dAx.out" using 3:6 title "J_x" with linespoints ls 1, \
           "dnu_dAy.out" using 4:6 title "J_y" with linespoints ls 3;
    if (!ps) pause -1;
    
    if (ps) set output "dnu_3.ps"
    #set size 0.5, 0.5; set origin 0.0, 0.0;
    set title "Chromaticity";
    set xlabel "{/Symbol d} [%]"; set ylabel "{/Symbol n}_x";
    set y2label "{/Symbol n}_y";
    set ytics nomirror; set y2tics;
    if (!pert) \
      plot "chrom2.out" using 1:(N*$2) title "{/Symbol n}_x" with lines ls 1, \
           "chrom2.out" using 1:(N*$3) axis x1y2 title "{/Symbol n}_y" \
           with lines ls 3;
    if (pert) \
      plot "chrom2.out" using 1:(N*$2) title "{/Symbol n}_x" with lines ls 1, \
           "chrom2.out" using 1:(N*$3) axis x1y2 title "{/Symbol n}_y" \
           with lines ls 3, \
           "chrom2_pert.out" using 1:(N*(N_x+$2)) axis x1y1 \
           title "{/Symbol n}_x (pert)" with lines ls 2, \
           "chrom2_pert.out" using 1:(N*(N_y+$3)) axis x1y2 \
           title "{/Symbol n}_y (pert)" with lines ls 4;
    if (!ps) pause -1;
    
    if (ps) set output "dnu_4.ps"
    #set size 0.5, 0.5; set origin 0.0, 0.0;
    set title "Rec. Chrom.: Quadradic Deviation";
    set xlabel "{/Symbol d} [%]"; set ylabel "{/Symbol n}_x";
    set y2label "{/Symbol n}_y";
    set ytics nomirror; set y2tics;
    plot "chrom2.out" using 1:((N*$2-33.134)**2) title "{/Symbol n}_x" with lines ls 1, \
         "chrom2.out" using 1:((N*$3-16.180)**2) axis x1y2 title "{/Symbol n}_y" \
         with lines ls 3;
    if (!ps) pause -1;
    
    fract(x) = x - int(x);
    
    #exit();
    
    set noztics; unset colorbox; set key left;
    set view 0, 0, 1, 1;
    
    if (ps) set output "dnu_5.ps"
    set title "Distance from \"3rd\" Order Resonances";
    set xlabel "{/Symbol d} [%]"; set ylabel "{/Symbol Dn}"; unset y2label;
    set ytics nomirror; unset y2tics;
    set yrange [0:1];
    splot "chrom2.out" using 1:(fract(N*($2))):(0.0) \
          title "{/Symbol n}_x" with lines palette z, \
          "chrom2.out" using 1:(fract(N*($3))):(0.1) \
          title "{/Symbol n}_y" with lines palette z, \
          "chrom2.out" using 1:(fract(N*2*($2))):(0.2) \
          title "2{/Symbol n}_x" with lines palette z, \
          "chrom2.out" using 1:(fract(N*2*($3))):(0.3) \
          title "2{/Symbol n}_y" with lines palette z, \
          "chrom2.out" using 1:(fract(N*(($2)+($3)))):(0.6) \
          title "{/Symbol n}_x+{/Symbol n}_y" with lines palette z, \
          "chrom2.out" using 1:(fract(N*(($2)-($3)))):(0.7) \
          title "{/Symbol n}_x-{/Symbol n}_y" with lines palette z, \
          "chrom2.out" using 1:(fract(N*3*($2))):(0.8) \
          title "3{/Symbol n}_x" with lines palette z, \
          "chrom2.out" using 1:(fract(N*(($2)+2*($3)))):(0.9) \
          title "{/Symbol n}_x+2{/Symbol n}_y" with lines palette z, \
          "chrom2.out" using 1:(fract(N*(($2)-2*($3)))):(1.0) \
          title "{/Symbol n}_x-2{/Symbol n}_y" with lines palette z;
    if (!ps) pause -1;
    
    if (ps) set output "dnu_6.ps"
    set title "Distance from \"4th\" Order Sextupolar Resonances";
    set xlabel "{/Symbol d} [%]"; set ylabel "{/Symbol Dn}"; unset y2label;
    set yrange [0:1];
    set ytics nomirror; unset y2tics;
    splot "chrom2.out" using 1:(fract(N*4*($2))):(0.0) \
          title "4{/Symbol n}_x" with lines palette z, \
          "chrom2.out" using 1:(fract(N*4*($3))):(0.2) \
          title "4{/Symbol n}_y" with lines palette z, \
          "chrom2.out" using 1:(fract(N*(2*($2)+2*($3)))):(0.8) \
          title "2{/Symbol n}_x+2{/Symbol n}_y" with lines palette z, \
          "chrom2.out" using 1:(fract(N*(2*($2)-2*($3)))):(1.0) \
          title "2{/Symbol n}_x-2{/Symbol n}_y" with lines palette z;
    if (!ps) pause -1;
    
    if (ps) set output "dnu_7.ps"
    set title "Distance from Resonances \"5th\" Order Sextupolar Resonances";
    set xlabel "{/Symbol d} [%]"; set ylabel "{/Symbol Dn}";
    set ytics nomirror; unset y2tics;
    splot "chrom2.out" using 1:(fract(N*5*($2))):(0.0) \
          title "5{/Symbol n}_x" with lines palette z, \
          "chrom2.out" using 1:(fract(N*5*($3))):(0.1) \
          title "5{/Symbol n}_y" with lines palette z, \
          "chrom2.out" using 1:(fract(N*(($2)+4*($3)))):(0.2) \
          title "{/Symbol n}_x+4{/Symbol n}_y" with lines palette z, \
          "chrom2.out" using 1:(fract(N*(($2)-4*($3)))):(0.8) \
          title "{/Symbol n}_x-4{/Symbol n}_y" with lines palette z, \
          "chrom2.out" using 1:(fract(N*(3*($2)+2*($3)))):(0.9) \
          title "3{/Symbol n}_x+2{/Symbol n}_y" with lines palette z, \
          "chrom2.out" using 1:(fract(N*(3*($2)-2*($3)))):(1.0) \
          title "3{/Symbol n}_x-2{/Symbol n}_y" with lines palette z;
    if (!ps) pause -1;