## Problem with totalFlux calculation in adjustPhi

The variable totalFlux in cfdTools/general/adjustPhi/adjustPhi.C is calculated like this:

`scalar totalFlux = VSMALL + sum(mag(phi)).value();`

If I understand this correctly, the sum is over all internal faces. At least in my case floating2Dbox.zip, the sum is initially zero, and I get totalFlux = 1e-300 although on the boundaries I have matching in- and outflow fluxes. For the attached case, this results in crash with:

```
--> FOAM FATAL ERROR: (openfoam-2012)
Continuity error cannot be removed by adjusting the outflow.
Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow.
Total flux : 1.0000000000000000251e-300
Specified mass inflow : 1.0000000000000059952
Specified mass outflow : 1.0000000000000064393
Adjustable mass outflow : 0
[floating2Dbox.zip](/uploads/eaba426b2041929d8e0787cc7a877d0d/floating2Dbox.zip)
```

The culprit is the line:

`else if (mag(fixedMassOut - massIn)/totalFlux > 1e-8)`

triggering the error.

My suggestion is to redefine totalFlux to also include the fixed boundary flux:

```
scalar totalFlux =
VSMALL + sum(mag(phi)).value() + mag(fixedMassOut) + mag(massIn);
```

It could also be considered to replace the hardcoded 1e-8 in the error condition above by for instance the `tolerance`

of pFinal or p_rghFinal specified by the user in fvSolution.