
../gsvit2d aperture.par ; ../postprocess2dnfff xnff.txt -c -80 -d 80 -e 361 -z 5500 -Z 7500 -s 3 -t 1 -p 0; cp output.txt output_300_7ksteps.txt


plot "output_300_4ksteps.txt" w l t "dtmult 0.5, height 300, 4k steps", "output_500_4ksteps.txt" w l t "dtmult 0.5, height 500, 4k steps", "output_300_4ksteps_dt08.txt" u 1:($2/sqrt(0.8/0.5)) w l t "dtmult 0.8, height 300, 4k steps (/sqrt(0.8/0.5))", "output_300_7ksteps.txt" w l t "dtmult 0.5, height 300, 7k steps"

plot "output_300_7ksteps_p0.txt" w l t "single slit", "output_300_7ksteps_p1.txt" u 1:($2/9) w l t "2 slits", "output_300_7ksteps_p2.txt" u 1:($2/25) w l t "5 slits", "output_300_7ksteps_p3.txt" u 1:($2/49) w l t "7 slits"


../gsvit2d aperture.par ; ../postprocess2dnfff xnff.txt -c -90 -d 90 -e 721 -z 9000 -Z 10000 -s 3 -t 1 -p 4 ; cp output.txt output_p4.txt

../gsvit2d 3aperture.par ; ../postprocess2dnfff xnff.txt -c -90 -d 90 -e 361 -z 8000 -Z 9000 -r 0 -R 160  -s 3 -t 1 -p 0 ; cp output.txt output_3s.txt



plot "xb_1" u 3:2 w l t "single aperture, analytical", "xb_3" u 3:2 w l t "3 apertures, analytical",  "xb_9" u 3:2 w l t "9 apertures, analytical", "output_p1.txt" u ($1/2-90):($2/42.58/9) w l t "3 apertures, periodic", "output_p4b.txt" u ($1/4-90):($2/21.55/91) w l t "9 apertures, periodic", "output_3s.txt" u ($1/2-90):($2/360) w l t "3 apertures, manual"


