docs: def _compute_fluxes(self):
          ''' All fluxes are band by band'''
          self.emission = self._compute_emission()
          self.emission_sfc = self._compute_emission_sfc()
          fromspace = self._from_space()
          self.flux_down = self.trans.flux_down(fromspace, self.emission)
          self.flux_reflected_up = self.trans.flux_reflected_up(self.flux_down, self.albedo_sfc)
          # this ensure same dimensions as other fields
          self.flux_to_sfc = self.flux_down[..., -1, np.newaxis]
          self.flux_from_sfc = (self.emission_sfc +
                                self.flux_reflected_up[..., -1, np.newaxis])
          self.flux_up = self.trans.flux_up(self.flux_from_sfc,
                              self.emission + self.flux_reflected_up[...,0:-1])
          self.flux_net = self.flux_up - self.flux_down
          # absorbed radiation (flux convergence) in W / m**2 (per band)
          self.absorbed = np.diff(self.flux_net, axis=-1)
          self.absorbed_total = np.sum(self.absorbed, axis=-1)
          self.flux_to_space = self._compute_flux_top()