Added the two D simulation code with halo boundary condition.