def _compute_heating_rates(self):
self.LW_to_atm = self._compute_emission()<EOL>self.SW_to_atm = self._compute_reflected_flux()<EOL>self.heating_rate['<STR_LIT>'] = ( self.LW_from_atm - self.LW_to_atm<EOL>+ self.SW_from_atm - self.SW_to_atm )<EOL>
Compute energy flux convergences to get heating rates in :math:`W/m^2`.