Temperature vertical profiles (a) comparing TRAPPIST-1b, c, and d (simulations T1b, T1c, and T1d), (b) several internal heat fluxes for TRAPPIST-1b (simulations T1b, T1b-Fgeo1, T1b-Fgeo5, and T1b-Fgeo25), and (c) several water partial pressures for TRAPPIST-1b (simulations T1b, T1b-3bar, T1b-1bar). The solid lines are global mean, temporal average (over 100 Earth days) temperature profiles calculated from GCM simulations. The dashed lines are calculated with a 1D inverse radiative-convective model (Turbet et al. 2019, 2020b) similar to those used in Kopparapu et al. (2013) to calculate the inner edge of the HZ.