zero gas phase in reaction-transport model

[GeeqC]

In the meantime, my reaction-transport model with decoupled methane gas runs fairly well, although a warning appears:

"WARNING: Element Mtg is contained in gas Mtg(g) (which has 0.0 mass),
but is not in solution or other phases."

It seems this warning appears as in some cells the aqueous solution becomes undersaturated with respect to CH4 gas so that the gas content goes to zero. Then also no gas phase should be present at all.

While the program seems to run fine, I wonder if this warning has some unexpected side effects and whether it would be possible to turn it off.

[dlparkhurst]

WARNINGs have no effect on the simulation; however, ERRORs cause the simulation to fail.

In this case, you have defined Mtg(g), but there are no moles of Mtg in the solution or reactants of the cell. Presumably, Mtg will be introduced in an infilling solution or produced by a KINETIC reaction, but at least initially, Mtg(g) cannot react in any way, which is why the warning was generated.

You can remove or limit warning messages if the workers run a file or string that includes "PRINT; -warnings n", where n is an integer. Each worker of PhreeqcRM will emit at most n warning messages. If n is zero, no warning messages will be produced.

[GeeqC]

Indeed, I added the following line after the SOLUTION 1 data block:

oss << "PRINT; -warnings 0" << "\n";

This works!
No more warnings.