Error with 0D batch simulation in TOUGHREACT v4.13

Hello everyone,

I recently purchased TOUGHREACT v4.13 and am trying to reproduce the batch experiment from this paper: https://www.mdpi.com/2073-4441/14/19/3160

I am using the EOS1 module and have set up a single-element model by removing the CONNE block and defining only one element in ELEME.

I've encountered a specific issue with my input files. The simulation runs  when I use my flow(1).inp file. However, when I run the simulation with my flow.inp file (for reactive transport), it fails with an error.

If anyone knows why a single-element model would work for flow(1).inp but not for flow.inp, I would greatly appreciate any advice.

Thank you for your help,

Hibiki Iwata