- pass the spacial increment to the reconstruction subroutine since
some interpolation methods require it; in addition, move obtaining
the spacial interval and its inversion to the subroutine update();
- remove subroutines magtocen(), expand_mag(), expand_mag_bnd() since
these subroutines are not used anymore; they are required in the case
of use of the magnetic field described on the staggered mesh;
- implement the MPI version of subroutine boundary_variables(); right
now we are sending full variable arrays between blocks, but this is
not optimal; better version should exchange only subdomain needed for
boundary update;
- add a global variable in the configuration module to determine if the
boundaries are periodic along each direction; this flag is useful if
we can skip some work which is only for non-periodic boundaries, such
as applying specific boundary conditions;
- with all variables cell centered we do not need to perform additional
update between blocks at the same level since now everything can be
done at once; the only restriction is that the boundary update has to
be done along each direction separately;
- subroutine numerical_flux() is not used in the GLM-MHD approach, so
remove it; later, during the implementation of GALERKIN methods, it
may be reintroduced;
- add parameters for the reconnection problem to the config module;
- properly initialize the antiparallel magnetic field component and the
perpendicular component initial perturbation;
- rewrite the function check_ref() using the subroutine to convert the
conserved variables to primitive ones;
- if the MHD is used take the gradients of magnetic field components
into account in determining the refinement indicator;
- remove subroutine advance(); this subroutine is good for GALERKIN
approach; it's gonna be reintroduced during this method
implementation;
- remove subroutine update_flux(); this subroutine is called from
advance() and is not used in the GLM-MHD approach;
- make GNU Fortran happy bu adding () to the subroutine calls;
- in subroutine prolong_block() prolongate the magnetic field
components and scalar potential during the block refinements if the
GLM-MHD equations are used;
- in the GLM-MHD set of equations variables Bx and Psi are couples,
thus find their values and substitute to the left and right states
before entering the approximate Riemann solver;