CBL marginal stability condition is used instead of the heat flow boundary condition. Predicted thickness of suboceanic thermal boundary layer agrees with LAB observations. A MSC-geotherm deep segment is not affected by crustal heat production uncertainty. Model explains a xenolith-based heat flow vs. age/lithosphere thickness correlation.