This paper presents a new three-dimensional geothermal reservoir simulator. The solution is based on finite difference techniques and uses a pressure semi-implicit, enthalpy explicit scheme. The equations are solved for each cell directly, while a new procedure is introduced for treating transitions from one to two phase conditions (and vice versa). A stability analysis is included proving the unconditional stability of the scheme. The model is applied to two cases for geothermal simulator testing.