diff --git a/as5013-test/hall-plots.png b/as5013-test/hall-plots.png
new file mode 100644
index 0000000000000000000000000000000000000000..b97cb0e1106a79d06331fb0e53689c03b8e5059a
Binary files /dev/null and b/as5013-test/hall-plots.png differ
diff --git a/as5013-test/hall-plots_full_view_aberrations.png b/as5013-test/hall-plots_full_view_aberrations.png
new file mode 100644
index 0000000000000000000000000000000000000000..6c0223d4e9882da3e8c05dd4969123baa211ff62
Binary files /dev/null and b/as5013-test/hall-plots_full_view_aberrations.png differ
diff --git a/as5013-test/sim/bfield_fd.py b/as5013-test/sim/bfield_fd.py
index 9086061fe1d52d53c521e5c9e0881c48484bc47f..d5e2e3944e40bc04b2c2f8014553872d587e4494 100644
--- a/as5013-test/sim/bfield_fd.py
+++ b/as5013-test/sim/bfield_fd.py
@@ -97,6 +97,7 @@ def run_single_magnet(args):
 	axs[1].set_ylabel('B_z')
 	axs[1].set_yticks([0.])
 	axs[1].set_xlim([-2*args.w,2*args.w]) #set x range based on magnet width
+	print "Section field extents: ", amin(Hy[py_sep1]), amax(Hy[py_sep1])
 
 	plt.savefig('run_single_magnet.png', bbox_inches='tight', pad_inches=0)
 	plt.show()
@@ -154,7 +155,7 @@ def run_double_magnet(args):
 	axs[1].set_ylabel('B_z')
 	axs[1].set_yticks([0.])
 	axs[1].set_xlim([-2*args.w,2*args.w]) #set x range based on magnet width
-
+	print "Section field extents: ", amin(Hy[py_sep1]), amax(Hy[py_sep1])
 	plt.savefig('run_double_magnet.png', bbox_inches='tight', pad_inches=0)
 	plt.show()
 
@@ -169,7 +170,7 @@ if __name__ == '__main__':
 	parser.add_argument("-pf","--pf", type=float, default=.75, help="plot fraction (how much of simulated space to plot)")
 	parser.add_argument("-l","--l", type=float, default=.01, help="magnet length")
 	parser.add_argument("-w","--w", type=float, default=.002, help="magnet width")
-	parser.add_argument("-sep","--sep", type=float, default=.001, help="separation between magnet bottom and graphing plane")
+	parser.add_argument("-sep","--sep", type=float, default=.0005, help="separation between magnet bottom and graphing plane")
 	parser.add_argument("-mu","--mu", type=float, default=1., help="relative permeability of magnet")
 	parser.add_argument("-Br","--Br", type=float, default=1., help="residual flux density of magnet (Teslas)")
 	parser.add_argument("-nt","--nt", type=int, default=50, help="num time steps")