def _compute_heating_rates(self):
self._compute_flux()<EOL>self.heating_rate['<STR_LIT>'] = -self._flux<EOL>self.heating_rate['<STR_LIT>'][..., -<NUM_LIT:1>, np.newaxis] = self._flux<EOL>
Compute energy flux convergences to get heating rates in :math:`W/m^2`.