Comparison of 1D temperature profiles of the hydrostatic models at the end of the simulation for two different radiation methods — flux-limited diffusion (FLD, used in this work) and moment transfer (M1, Melon Fuksman & Klahr 2022). Overall we see very good agreement with temperatures with FLD and M1, with the optically thin atmosphere being captured quite well. The dotted black line shows the analytical solution for the temperature profile in the optically thin atmosphere (Chiang & Goldreich 1997). The midplane temperatures are also in good agreement, with some minor deviations. The M1 method shows higher amplitude bumps than the FLD method. We also include the RADMC-3D calculations corresponding to our FLD model, these correspond to the green lines.