diff --git a/LDDsimulation.py b/LDDsimulation.py
index 19d4a14651356d01def2a0a7b0385ca316f2a4a7..b86e73a3c9bc3270654bdc10dbebcfba493d4572 100644
--- a/LDDsimulation.py
+++ b/LDDsimulation.py
@@ -135,6 +135,7 @@ class LDDsimulation(object):
         self._init_interfaces()
         self._init_subdomains()
         self._init_initial_values()
+        self._eval_sources()
 
     ## Private methods
     def _init_meshes_and_markers(self, subdomain_def_points: tp.List[tp.List[df.Point]] = None,#