We analyze two multi-chord stellar occultations by Pluto that were observed on 2012 July 18th and 2013 May 4th, and respectively monitored from five and six sites. They provide a total of fifteen light curves, 12 of which were used for a simultaneous fit that uses a unique temperature profile, assuming a clear (no haze) and pure N2 atmosphere, but allowing for a possible pressure variation between the two dates. We find a solution that satisfactorily fits (i.e., within the noise level) all of the 12 light curves, providing atmospheric constraints between ∼1190 km (pressure ∼11 μbar) and ∼1450 km (pressure ∼0.1 μbar) from Pluto's center. Our main results are: (1) the best-fitting temperature profile shows a stratosphere with a strong positive gradient between 1190 km (at 36 K, 11 μbar) and r = 1215 km (6.0 μbar), where a temperature maximum of 110 K is reached; above it is a mesosphere with a negative thermal gradient of -0.2 K km-1 up to ∼1390 km (0.25 μbar), where the mesosphere connects itself to a more isothermal upper branch around 81 K; (2) the pressure shows a small (6%) but significant increase (6σ level) between the two dates; (3) without a troposphere, Pluto's radius is found to be 5 km. Allowing for a troposphere, RP is constrained to lie between 1168 and 1195 km; and (4) the currently measured CO abundance is too small to explain the mesospheric negative thermal gradient. Cooling by HCN is possible, but only if this species is largely saturated. Alternative explanations like zonal winds or vertical compositional variations of the atmosphere are unable to explain the observed mesospheric negative thermal gradient.