To control the convergence of a simulation, it would be interesting to print the time evolution of the usual variables of the energy budget: incident stellar radiation (ISR), absorbed stellar radiation (ASR), outgoing longwave radiation (OLR), downwelling longwave radiation (GLR), upwelling longwave radiation from the surface (LWU). These quantities correspond to spatially averaged fluxes over the sphere. High-frequency time oscillations of these quantities may indicate a too large time step for radiative transfer.
To control the convergence of a simulation, it would be interesting to print the time evolution of the usual variables of the energy budget: incident stellar radiation (ISR), absorbed stellar radiation (ASR), outgoing longwave radiation (OLR), downwelling longwave radiation (GLR), upwelling longwave radiation from the surface (LWU). These quantities correspond to spatially averaged fluxes over the sphere. High-frequency time oscillations of these quantities may indicate a too large time step for radiative transfer.