Numerical modeling of a desaturating geothermal reservoir