def _compute_heating_rates(self):
for varname in list(self.state.keys()):<EOL><INDENT>self.heating_rate[varname] = self.state[varname] * <NUM_LIT:0.><EOL><DEDENT>
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.