make variables for resisitivity
debbuged constant boundary condition
a
modification of simulation box
copy from initial_parker.f90
make initial condition for yokoyama and shibata 1996
a
change the number of mesh in x-direction
initial condition for yokoyama et al.1996
fixed a bug of inverse gamma
position of half grid
apply boundary condition also for magnetic potential
debugged perturbation