simulationofgeothermalsubsidence