code: def _compute_heating_rates(self):
          '''Compute energy flux convergences to get heating rates in :math:`W/m^2`.'''
          self._compute_flux()
          self.heating_rate['Ts'] = -self._flux
          # Modify only the lowest model level
          self.heating_rate['Tatm'][..., -1, np.newaxis] = self._fluf _compute_heating_rates(self):
          '''Compute energy flux convergences to get heating rates in :math:`W/m^2`.'''
          self._compute_flux()
          self.heating_rate['Ts'] = -self._flux
          # Modify only the lowest model level
          self.heating_rate['Tatm'][..., -1, np.newaxis] = self._flux
docs: Compute energy flux convergences to get heating rates in :math:`W/m^2`.