INTEGRAL_FLUX
  ! 2 x 2 x - with the connection between 1 and 5 carved out: 3 total
  ! 3 m^2 \perp Z
  ! @ 1 mm/s for 2 s = -6 kg total, -3 kg/s
  NAME Z3
  INVERT_DIRECTION
  POLYGON
    1. 0. 1.
    2. 0. 1.
    2. 2. 1.
    0. 2. 1.
    0. 1. 1.
    1. 1. 1.
  /
END

INTEGRAL_FLUX
  ! 2 x 1 x - top
  ! 2 m^2 \perp Z
  ! @ 1 mm/s for 2 s = 4 kg total, 2 kg/s
  NAME top
  POLYGON
    0. 0. 2.
    2. 1. 2.
  /
END

INTEGRAL_FLUX
  NAME entire_top
  ! 4 m^2 \perp Z
  ! @ 1 mm/s for 2 s = 8 kg total, 4 kg/s
  PLANE
    0. 0. 2.
    1. 0. 2.
    0. 1. 2.
  /
END

# coordinates and directions
INTEGRAL_FLUX
  ! 1 m^2 \perp X, Y, Z = 3 m^2 total
  ! @ 1 mm/s for 2 s = -2 kg total, -1 kg/s
  ! negative due to -0.5 direction in z
  NAME top_north_east_corner_z
  COORDINATES_AND_DIRECTIONS
    1.5 1.5 2. 0. 0. -0.5  ! top face of cell 8
  /
END

# =============================== Z plane at 0 ================================
INTEGRAL_FLUX
  ! @ 1 mm/s for 2 s = 8 kg total, 4 kg/s
  NAME Z_at_0
  PLANE
    1. 0. 0.
    0. 1. 0.
    0. 0. 0.
  /
END

! the four below should sum to the above
INTEGRAL_FLUX Z_at_0_1
  POLYGON
    0. 0. 0.
    1. 1. 0.
  /
END

INTEGRAL_FLUX Z_at_0_2
  POLYGON
    1. 0. 0.
    2. 1. 0.
  /
END

INTEGRAL_FLUX Z_at_0_3
  POLYGON
    0. 1. 0.
    1. 2. 0.
  /
END

INTEGRAL_FLUX Z_at_0_4
  POLYGON
    1. 1. 0.
    2. 2. 0.
  /
END

# =============================== Z plane at 2 ================================
INTEGRAL_FLUX
  ! @ 1 mm/s for 2 s = 8 kg total, 4 kg/s
  NAME Z_at_2
  PLANE
    1. 0. 2.
    0. 1. 2.
    0. 0. 2.
  /
END

! the four below should sum to the above
INTEGRAL_FLUX Z_at_2_5
  POLYGON
    0. 0. 2.
    1. 1. 2.
  /
END

INTEGRAL_FLUX Z_at_2_6
  POLYGON
    1. 0. 2.
    2. 1. 2.
  /
END

INTEGRAL_FLUX Z_at_0_7
  POLYGON
    0. 1. 2.
    1. 2. 2.
  /
END

INTEGRAL_FLUX Z_at_2_8
  POLYGON
    1. 1. 2.
    2. 2. 2.
  /
END

