A computer program is described that can be used to design a fracture stimulation treatment for a geothermal reservoir. The program uses state-of-the-art methods to calculate the temperature of the fracture fluid as a function of time and distance in the fracture. This information is used to determine the temperature dependent properties of the fracture fluid. These fluid properties are utilized to calculate the fracture geometry as a function of time. The fracture geometry and temperature distribution of the fracture fluid are coupled so the subroutines that calculate these distributions have been made interactive