Results from the Magnetic Induction Equation

I have a Fortran program called "rmagnet.f90" which aims to solve the Magnetic Induction Equation numerically.

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.

Back