Simulation of geothermal subsidence