We present an algorithm to produce the necessary geometric information for finite
volume calculations in the context of Cartesian grids with embedded boundaries.
Given an order of accuracy for the overall calculation, we show what accuracy is
required for each of the geometric quantities and we demonstrate how to calculate
the moments using the divergence theorem. We demonstrate that, for a known
flux, these moments can be used to create a flux divergence of the expected
order.