Fix output to include boundary condition cells