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
Merge branch 'jet' into jet_engin

Conflicts:
	boundary_parker.f90