Potential magnetic field calculator for solar physics applications using staggered grids

School of Physics and Astronomy, Queen Mary University of London, G.O. Jones Building, 327 Mile End Road, London E1 4NS, UK
e-mail: c.boocock@qmul.ac.uk

Received: 19 November 2018
Accepted: 18 March 2019

Abstract

A program has been designed to generate accurately a potential magnetic field on a staggered grid by extrapolating the magnetic field normal to the photospheric surface. The code first calculates a magnetic potential using the Green’s function method and then uses a finite differencing scheme to calculate the magnetic field from the potential. A new finite differencing formula was derived which accounts for grid staggering; it is shown that this formula gives a numerical approximation that is closest to the real potential field. It is also shown that extending the region over which normal photospheric field is specified can improve the accuracy of the potential field produced. The program is a FORTRAN 90 code that can be used to generate potential magnetic field inputs for Lare3d and other MHD solvers that use a staggered grid for magnetic field components. The program can be parallelised to run quickly over multiple computing cores. The code and supporting description are provided in the appendices.