Cubadak area was located in depression zone of great Sumatra fault zone which affected by normal faults. Geothermal in this area was characterized by the appearance of hot springs type chloride-bicarbonate on Cubadak, Sawah Mudik and Talu. The manifestation temperature on this area is about 37,1° - 74,8° C. This study purpose to recontruct model of subsurface geothermal system in Cubadak area based on MT data. Model 1-D was processed by apply Bostick inversion and NLCG (Nonlinear Conjugate Gradient) for model 2-D. The result showed that lowest resistivity was on Southwest to Northeast research area. Highest resistivity spread at West to Northeast. Geothermal system was composed by caprock (2-10 ohm), reservoar (20-300 ohm) and heat source (>300 ohm).