advection_c.c error linking to nxyz value

[Yongqiang]

Hi Phreeqc experts,

I changed the nxyz values in the advection_c.c.

I tried nxyz = 10, 20, 30, 40, 400. All works well, except the value 10. It popped out that the following information:


WARNING: Mass of water is less than 1e-10 kilogram.
The aqueous phase may not be stable relative to given masses of minerals.
WARNING: Numerical method failed with this set of convergence parameters.

Any comments for this strange error?

Regards,
Michael

[dlparkhurst]

I think there is an assumption that nxyz is 20 or greater in the coding of the test case. Initial conditions for cell 18 and 19 (0 based index) are defined explicitly.