Added reference comparison with GNAT and script for extracing percent differences. Also plot them.HEADmaster

1 files changed, 25 insertions, 5 deletions

diff --git a/plot.gnu b/plot.gnu
index e001571..80fd5cd 100644
--- a/plot.gnu
+++ b/plot.gnu
@@ -1,12 +1,32 @@
-#set term x11 persist
-set term cairolatex pdf size 4in,2.5in
+#set term x11 1 persist
+set term cairolatex pdf size 5in,3in
 set output "flux.tex"
 
-set key reverse top left
+set key outside top center horizontal
 set grid
 set xlabel "Distance [cm]"
 set ylabel "$\\phi$"
+set border lw 3
 
+# Draw the regions
+
+set lt 9 lw 1 lc rgb "black" dt 2
+set arrow from 2,graph 0 to 2, graph 1 nohead lt 9
+set arrow from 3,graph 0 to 3, graph 1 nohead lt 9
+set arrow from 5,graph 0 to 5, graph 1 nohead lt 9
+set arrow from 6,graph 0 to 6, graph 1 nohead lt 9
+
+
+plot 'flux.out' using 1:($2 + 3*$3):($2 - 3*$3) with filledcurves lc rgb "#CCCCFF" title "$\\pm 3\\sigma$", \
+    'flux.out' using 1:2 with lines lc rgb "red" lw 2 title "$\\overline{\\phi}_{MC}$", \
+    'test.ref' using 1:2 with lines lc rgb "green" lw 1 dt 2 title "$\\overline{\\phi}_{ref}$"
+
+
+#set term x11 2 persist
+set output "error.tex"
+set yrange [-15:15]
+set key off
+
+set ylabel "Percent Difference [\\%]"
+plot 'diff.out' using 1:2 with lines
 
-plot 'flux.out' using 1:($2 + $3):($2 - $3) with filledcurves lc rgb "#CCCCFF" title "$\\pm 1 \\sigma$", \
-    'flux.out' using 1:2 with lines lc rgb "red" title "$\\overline{\\phi}$"