Use refined (not unrefined) density discretization to compute area elements (Closes #66 on Gitlab)