I set the program's parameters to the following:
meshpoints=20
dt=0.001
steps in timeloop=100
I wanted to make the magnetic field B appear to be cylindrical. To do this, I defined a radius in terms of the y & z axes and defined B such that it depended upon this radius in the x-direction only. My initial condition for B was:
B=(-72*(y-0.5)exp(-36r*r),0,0)
I was limited to this form for B, since I intended to later solve this equation in terms of a vector potential A. With this in mind, I was satisfied that the modulus of my magnetic field will actually appear as 2 cylinders now.
In order to check that my results were as I hoped, I used another IDL program called, "cont.pro" to plot a contoured cross-section, shown below.
My contour plot verified that the results from my code were as I had expected.
