code: def flux_components(self, kwargs_light, n_grid=400, delta_grid=0.01, deltaPix=0.05, type="lens"):
          
          flux_list = []
          R_h_list = []
          x_grid, y_grid = util.make_grid(numPix=n_grid, deltapix=delta_grid)
          kwargs_copy = copy.deepcopy(kwargs_light)
          for k, kwargs in enumerate(kwargs_light):
              if 'center_x' in kwargs_copy[k]:
                  kwargs_copy[k]['center_x'] = 0
                  kwargs_copy[k]['center_y'] = 0
              if type == 'lens':
                  light = self.LensLightModel.surface_brightness(x_grid, y_grid, kwargs_copy, k=k)
              elif type == 'source':
                  light = self.SourceModel.surface_brightness(x_grid, y_grid, kwargs_copy, k=k)
              else:
                  raise ValueError("type %s not supported!" % type)
              flux = np.sum(light)*delta_grid**2 / deltaPix**2
              R_h = analysis_util.half_light_radius(light, x_grid, y_grid)
              flux_list.append(flux)
              R_h_list.append(R_h)
          return flux_list, R_h_list
docs: computes the total flux in each component of the model
  
          :param kwargs_light:
          :param n_grid:
          :param delta_grid:
          :return: