def METHOD_NAME(self, sd: pp.Grid, data: dict):
    """Return the rhs for a discretization of the integrated source term. Also
    discretize the necessary operators if the data dictionary does not contain a
    source term.
    Parameters:
        sd (Grid): Computational grid, with geometry fields computed.
        data (dictionary): With data stored.
    Returns:
        np.array (self.ndof): Right hand side vector representing the
            source.
    """
    matrix_dictionary = data[pp.DISCRETIZATION_MATRICES][self.keyword]
    parameter_dictionary = data[pp.PARAMETERS][self.keyword]
    sources = parameter_dictionary["source"]
    assert sources.size == self.ndof(
        sd
    ), "There should be one source value for each cell"
    return matrix_dictionary["bound_source"] * sources