[gschaider]

In ReactingParcelIO.C there is an error that makes it fail if reading data from a previous run where some of the parcels have for instance a Y_liquid of 0 (division by zero). This patch fixes this:

ReactingParcelIO.patch

BTW: I haven't tested it, but it occurred to me during fixing this that if there was for instance a species H2O in the liquid and the solid phase of the parcel they both would get written to lagrangian/YH2O. I'd propose calling them YlH2O and YsH2O just to avoid such clashes (however unlikely they might be)