area_yz = (myNi[1] - 2 * halo) * (myNi[2] - 2 * halo)