Comparison between the numerical solutions (dashed lines) and the analytical solutions of model D (using the functional form of the coefficients given in Table 1) over a wide range of irradiation temperatures for a giant planet of solar composition orbiting a sun-like star. Here we used g = 25 m/s², T_int = 100 K, and . Thick and thin lines represent convective and radiative zones, respectively.