Geothermal model breaking after increasing injection / injection duration

[Mai MRST]

Hello MRST team!

I'm using the geothermal model and have been trying to run a few models (based on HTATES example) and I am running into some problems: it breaks when I increase injection rate and injection duration.

Relevant model properties: full sandstone model, cells are 20 m in x, y, z (8000 m3 cubes)

Model breaking when extended time from 4 months injection to 6 months injection:

Solving timestep 014/216: 160 Days -> 180 Days

Solver did not converge in 25 iterations for timestep of length 20 Days. Cutting timestep.

!!! Simulation resulted in fatal error !!!
Exception thrown: Specify the coordinates as vectors or matrices of the same size, or as a vector and a matrix that share the same length in at least one dimension.

Error using plot
Specify the coordinates as vectors or matrices of the same size, or as a vector and a matrix that share the same
length in at least one dimension.

Model breaking after increasing injection rate from 1000 m3/day to 5000 m3/day:

Solving timestep 006/192: 10 Days -> 20 Days

Solver did not converge in 25 iterations for timestep of length 10 Days. Cutting timestep.

Solver did not converge in 25 iterations for timestep of length 5 Days. Cutting timestep.

!!! Simulation resulted in fatal error !!!

What is the explanation for the model breaking and is there a way around it?

Lastly, I have been using initializeGeothermalEquilibrium module and was wondering, does it also account for salinity in the geothermal model as well?

Many thanks in advance for your support and answers!