rewrite calc_gli_term

LDDsimulation/boundary_and_interface.py

@@ -241,6 +241,8 @@ class Interface(BoundaryPart):
         # global index ob the subdomain
         self.global_index = global_index
         self.marker_value = self.global_index+1
+        # dictionary holding the facet objects
+        self.facets = dict()
         # dictionary containing for each facet belonging to the interface (with
         # respect to a subdomain mesh) the coordinates of the the vertices
         # belonging to that facet. This dictionary is created by
@@ -420,6 +422,7 @@ class Interface(BoundaryPart):
         """
         self.facet_to_vertex_coordinates.update({subdomain_index: dict()})
         self.outer_normals.update({subdomain_index: dict()})
+        self.facets.update({subdomain_index: dict()})
         # Get a subset iterator for mesh entities along the facets marked by
         # interface_marker and value interface_marker_value.
         mesh_entity_iterator = df.SubsetIterator(interface_marker, interface_marker_value)
@@ -442,6 +445,7 @@ class Interface(BoundaryPart):
             self.facet_to_vertex_coordinates[subdomain_index].update(
                 {facet.index(): facet_vertices}
             )
+            self.facets[subdomain_index].update({facet.index(): facet})
             if subdomain_index == 0:
                 self.outer_normals[subdomain_index].update(
                     {facet.index(): None}