diff --git a/forcemodels.py b/forcemodels.py
index 5b65b271ddee6812c2159f4370be67ddd7215000..1fe5438868a9bcb45a99389fb2e036da8b4cd122 100644
--- a/forcemodels.py
+++ b/forcemodels.py
@@ -90,14 +90,6 @@ import sympy as sp
 from lbmpy.relaxationrates import get_shear_relaxation_rate
 
 
-class ScalarSource(object):
-    def __init__(self, source):
-        self._source = source
-
-    def __call__(self, lb_method, **kwargs):
-        return [w_i * self._source for w_i in lb_method.weights]
-
-
 class Simple(object):
     r"""
     A simple force model which introduces the following additional force term in the
diff --git a/updatekernels.py b/updatekernels.py
index 5ee8f782f49ee91031be79d9b3bd93f03acb2511..3eacdbbe22d4308da06cca04f037ed2f70916f43 100644
--- a/updatekernels.py
+++ b/updatekernels.py
@@ -104,20 +104,3 @@ def create_pdf_array(size, num_directions, ghost_layers=1, layout='fzyx'):
     if isinstance(layout, str):
         layout = layout_string_to_tuple(layout, dim + 1)
     return create_numpy_array_with_layout(size_with_gl + [num_directions], layout)
-
-
-# ------------------------------------------- Add output fields to kernel ----------------------------------------------
-
-
-def add_output_field_for_conserved_quantities(collision_rule, conserved_quantities_to_output_field_dict):
-    method = collision_rule.method
-    cqc = method.conserved_quantity_computation.output_equations_from_pdfs(method.pre_collision_pdf_symbols,
-                                                                           conserved_quantities_to_output_field_dict)
-    return collision_rule.new_merged(cqc)
-
-
-def write_quantities_to_field(collision_rule, symbols, output_field):
-    if not hasattr(symbols, "__len__"):
-        symbols = [symbols]
-    eqs = [Assignment(output_field(i), s) for i, s in enumerate(symbols)]
-    return collision_rule.copy(collision_rule.main_assignments + eqs)