In the [1]} paper the improper handling of the interfacial water also caused an issue in selecting the proper simulation time step: if larger time steps for the same timespan were used, the calculated initial temperatures required for the same serpentinization level were lower and the final temperatures higher, causing a greater overall heat production in the simulations.
The reason behind this phenomenon was that in the case of decreased time steps the serpentinization reaction started later in time.
By incorporating the new condition in the melting ice calculation, the simulations did not produce more interfacial water any more than it was possible at the given temperature, and there was no difference in the amount of heat produced as a function of time scale at the same time.
