#A demonstration package for Divergence-free Interpolation
#
# Colin P. McNally 2010,2011
# cmcnally at amnh.org
#
#  Department of Astrophysics
#  American Museum of Natural History
#  79th Street at Central Park West
#  New York, NY, USA, 10024
#
#  and
#
#  Department of Astronomy
#  Columbia University
#
# This code is not very fast, its here to show you how things can be done.
#
# This is public-domain code. This is too simple to bother with a lisence.

# The files here will do the Alfven wave problem as set up, and you can play with other things.
# The numerical method actually implemented isn't good for much research, the point it to 
# demonstrate by construction that divergence-free interpolation stencils can be used for MHD.
#
# The visualiation commands in rbfmhd.py will pop up a window the refreshes the plots on each timestep.
# However, you may need to tweak things to make it display on your matplotlib install - I leave that to the reader.
# This code runs with the MacOSX backend selected in my matplitlibrc using the line
# backend      : MacOSX

