def _compute_heating_rates(self):
for varname, value in self.state.items():<EOL><INDENT>self.heating_rate[varname] = - self.b * (value - global_mean(value))<EOL><DEDENT>
Computes energy flux convergences to get heating rates in :math:`W/m^2`.