docs: def _compute_heating_rates(self):
          """Computes energy flux convergences to get heating rates in :math:`W/m^2`.
  
          """
          self._compute_emission()
          for varname, value in self.state.items():
              self.heating_rate[varname] = -self.OLR