added hall plots

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")