diff --git a/LDDsimulation/LDDsimulation.py b/LDDsimulation/LDDsimulation.py
index 662d5c41476fa35927182fe924f307bf8553d2a7..4e77d417531fc95e24a4116d364b633b45e33685 100644
--- a/LDDsimulation/LDDsimulation.py
+++ b/LDDsimulation/LDDsimulation.py
@@ -1046,13 +1046,14 @@ class LDDsimulation(object):
         """
         if self.exact_solution is not None:
             for sd_ind, subdomain in self.subdomain.items():
-                mesh = subdomain.mesh
                 subdomain.pressure_exact = dict()
                 for phase in subdomain.has_phases:
                     print(f"initialising exact solutions for subdomain{subdomain.subdomain_index} and phase {phase}")
                     pa_exact = self.exact_solution[sd_ind][phase]
                     subdomain.pressure_exact.update(
-                        {phase: df.Expression(pa_exact, domain = mesh, degree=self.interpolation_degree, t=self.starttime) }
+                        {phase: df.Expression(pa_exact, domain = subdomain.mesh, #
+                                                        degree=self.interpolation_degree, #
+                                                        t=self.starttime) }
                         )
         else:
             pass