The work consists of the development of a computational model for the simulation of the temperatures, using a system of differential equations in partial derivatives, which solution was found by finite differences method, using the hypothesis of the isobaric coordinates, it means that the pressure is not a function but becomes an independent variable, in the new system, z is a function, more precisely a variable closely related to the geopotential ф. The data was obtained from the NCAR's Research Data Archive Database in January 2011 at different altitudes of 2 to 4572 meters above sea level (m.a.s.l.). The grid was structured with 10 spatial points on the x-axis, 8 spatial points on the y axis, and 5 sigma’s sigma's layers at thesame altitude evaluated, representing the 80 km² covering the Pamplonita River basin in the sector between Pamplona and Cucuta. The calculations obtained by the modeling system at the time of validation of the model have a low error level. The minimum error of validating the model was 0.17 %, obtained for the results of the model at 2 m.a.s.l.