1] If the steady state temperature field, the diffusion temperature time step, or the initial lithostatic pressure are activated, they are performed also after restarting the simulation:
a.) This implies that the vector of pressure from the restart database is overwritten by the guess of lithostatic pressure (=BAD).
b.) This implies that the current temperature field coming from the restart database will be modified.
2] This PR introduces a solution:
-> TSSolve is an argument of these function: if the timestep is different than 0, these function are not applied.
1] If the steady state temperature field, the diffusion temperature time step, or the initial lithostatic pressure are activated, they are performed also after restarting the simulation: a.) This implies that the vector of pressure from the restart database is overwritten by the guess of lithostatic pressure (=BAD). b.) This implies that the current temperature field coming from the restart database will be modified. 2] This PR introduces a solution: -> TSSolve is an argument of these function: if the timestep is different than 0, these function are not applied.