docs: def _compute_heating_rates(self):
          """Computes energy flux convergences to get heating rates in unit
          :math:`\\textrm{W}/ \\textrm{m}^2`.
  
          This method should be over-ridden by daughter classes.
  
          """
          for varname in list(self.state.keys()):
              self.heating_rate[varname] = self.state[varname] * 0.