Because the model immediately removes bubble CO2 from the domain and does not track the bubble phase explicitly, to determine the gas bubble flux out of the domain, the model is run with and without the flux term shown in equation (17), and the flux of bubbles out of the domain