Abstarct: Magnetotelluric method is one of the most widely used subsurface resources exploration methods, which is a suitable method for detecting subsurface abnormalities at high depths. The studied area in the Moil Valley is located in the geographical coordinates and the 38ᵒ12'26.7'', 38ᵒ16'59.8'' latitude and 47ᵒ43'38.2'', 47ᵒ45'58.6'' longitude, in the south of Meshkinshahr. In this study, the most consistent model with the source of geothermal is identified by using two-dimensional modeling of Magnetotelluric data of the region, by comparing the two models. In this research, two WinGlixnk and ZondMT2D softwares which are the most well-known modeling software in the field of MT were used to compare two models to achieve the most appropriate model. Because in both softwares the different approaches are used in the term of MT methodology for modeling. In addition to the well information, the output of the two Oasis Montaj and Surfer softwares have also been used. Due to the fact that the dimensions of the region structures are often two-dimensional in the middle depths, hence for more accurate results from the location of the Anomalous zone, two-dimensional inversion was carried out by Nonlinear Conjugate gradient (NLCG) method on MT data profile. Considering the all modeling done with both approaches in the metioned softwares, the results of the two-dimensional inversion baxsed on the existence resistivities show that the geothermal source is located at the depth of 2000 meter in the approximate position below the stations 25, 21, 6, 217 and 2. As it is obvious, in the Focus modeling, it is possible to identify the source of the geothermal than the Smooth modeling. Magmatic intrusions are the main thermal source of the system, which is located in a depth of 5 to 6 kilometers, and include a rock structure from granite to gabbro-diorite.