The iteration process is done to determine the surface temperature from Eqs. (67)–(70) by using the heat flux residual (Eqs. (71) and (72)) and its derivative with respect to temperature, D/T j as suggested by Finlayson et al. (1993) until the heat flux residual approaches zero.