Found that field_B_ also needed to be initailized otherwise if time step is odd number, y boundary has zero value
Correct this mistake
Implement Open MPI version of PDE solver
Can compile, but not tested