The retrieval of surface composition from IR measurements of airless bodies requires a model capable of computing the significant thermal gradients present in the top few hundred microns of the regolith. In this study we introduce a model which reproduces most of the features found in controlled experiments made in the simulated lunar environment emission chamber (SLEEC). Although the model presented here is forced by a lower boundary held at a fixed temperature, we conclude that a similar al...