python bfield_fd.py --mode single_magnet -nt 40000 -mu 100 -nx 200 -ny 200
python bfield_fd.py --mode double_magnet -nt 40000 -mu 100 -nx 200 -ny 200