uiae

LDDsimulation/LDDsimulation.py

@@ -1045,7 +1045,7 @@ class LDDsimulation(object):
                         subdomain.function_space["pressure"]['wetting']
                         )
                     saturation_nw = df.Function(
-                        subdomain.function_space["pressure"]['wetting']
+                        subdomain.function_space["pressure"]['nonwetting']
                         )
 
                 for phase in subdomain.has_phases:
@@ -1062,7 +1062,6 @@ class LDDsimulation(object):
                         "source_"+"{}".format(phase)
                         )
                     file.write(source[phase], timestep)
-
                     if self.exact_solution:
                         pa_exact = subdomain.pressure_exact[phase]
                         pa_exact.t = timestep