def _compute_fluxes(self):
self.emission = self._compute_emission()<EOL>self.emission_sfc = self._compute_emission_sfc()<EOL>fromspace = self._from_space()<EOL>self.flux_down = self.trans.flux_down(fromspace, self.emission)<EOL>self.flux_reflected_up = self.trans.flux_reflected_up(self.flux_down, self.albedo_sfc)<EOL>self.flux_to_sfc = self.flux_down[..., -<NUM_LIT:1>, np.newaxis]<EOL>self.flux_from_sfc = (self.emission_sfc +<EOL>self.flux_reflected_up[..., -<NUM_LIT:1>, np.newaxis])<EOL>self.flux_up = self.trans.flux_up(self.flux_from_sfc,<EOL>self.emission + self.flux_reflected_up[...,<NUM_LIT:0>:-<NUM_LIT:1>])<EOL>self.flux_net = self.flux_up - self.flux_down<EOL>self.absorbed = np.diff(self.flux_net, axis=-<NUM_LIT:1>)<EOL>self.absorbed_total = np.sum(self.absorbed, axis=-<NUM_LIT:1>)<EOL>self.flux_to_space = self._compute_flux_top()<EOL>
All fluxes are band by band