Temperature inversion layers originate based on the combined forcing of local- and large-scale synoptic meteorology. A numerical procedure based on a linear interpolation function of variable length that minimizes an error function set a priori is proposed to extract thermodynamic information of the multilayered thermal structure. The method is demonstrated to detect surface-based inversion and multilayered elevated inversions present often in high-latitude atmospheres.