There is an issue when computing bgq cycling with power-of-2 or multiple species. The problem is that the kinetics is being calculated based on mass to enable including non-mobile-non-dissolved-substances.
This is not an issue with only one species in the kinetics because there will be division and then multiplication by Vol, so they cancel out.
However, in the cases with multiple species of power-of-2, we will have wrong results and have to divide the mass transfer calculated by Vol^(n-1) . I need to find a workaround...
Make sure that unit conversions, in or out, are being calculated properly.
