diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-model-change-interface12.py b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-model-change-interface12.py
index 1a3486ee337f365d8ef9941b1a316356052be020..d4016c0b2adf5044a95de81163b2f76cb468f926 100755
--- a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-model-change-interface12.py
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-model-change-interface12.py
@@ -1,7 +1,26 @@
 #!/usr/bin/python3
-"""TP-TP layered soil with inner patch simulation.
-
-This program sets up an LDD simulation
+"""TP-TP layered soil with inner patch simulation, Copyright 2020, David Seus
+
+This program runs an LDD simulation for a layered soil substructuring with
+inner patch.
+This examples uses a model change at interface 12
+
+# LICENCE #####################################################################
+Copyright 2020, David Seus
+david.seus[at]ians.uni-stuttgart.de
+This program is free software: you can redistribute it and/or modify
+it under the terms of the GNU General Public License as published by
+the Free Software Foundation, either version 3 of the License, or
+(at your option) any later version.
+
+This program is distributed in the hope that it will be useful,
+but WITHOUT ANY WARRANTY; without even the implied warranty of
+MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+GNU General Public License for more details.
+
+You should have received a copy of the GNU General Public License
+along with this program.  If not, see <http://www.gnu.org/licenses/>.
+###############################################################################
 """
 import dolfin as df
 import sympy as sym
diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-realistic-pure-dd.py b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-realistic-pure-dd.py
index 0dbf65683f3f7eaec841fc953efd8dde99576e4e..ee8945984490f20669ad6da03ef253946a046655 100755
--- a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-realistic-pure-dd.py
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-realistic-pure-dd.py
@@ -1,7 +1,26 @@
 #!/usr/bin/python3
-"""TP-TP layered soil with inner patch simulation.
-
-This program sets up an LDD simulation
+"""TP-TP layered soil with inner patch simulation, Copyright 2020, David Seus
+
+This program runs an LDD simulation for a layered soil substructuring with
+inner patch.
+This example assumes realistic parameters but remains pure dd
+
+# LICENCE #####################################################################
+Copyright 2020, David Seus
+david.seus[at]ians.uni-stuttgart.de
+This program is free software: you can redistribute it and/or modify
+it under the terms of the GNU General Public License as published by
+the Free Software Foundation, either version 3 of the License, or
+(at your option) any later version.
+
+This program is distributed in the hope that it will be useful,
+but WITHOUT ANY WARRANTY; without even the implied warranty of
+MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+GNU General Public License for more details.
+
+You should have received a copy of the GNU General Public License
+along with this program.  If not, see <http://www.gnu.org/licenses/>.
+###############################################################################
 """
 import dolfin as df
 import sympy as sym
diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py
index 3a45bb7e45fe95ca6535e968259803aa34e0c08a..2e53df413c2f3124aec3eff1eaa3cc551249810f 100755
--- a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py
@@ -1,7 +1,25 @@
 #!/usr/bin/python3
-"""TP-TP layered soil with inner patch simulation.
-
-This program sets up an LDD simulation
+"""TP-TP layered soil with inner patch simulation, Copyright 2020, David Seus
+
+This program runs an LDD simulation for a layered soil substructuring with
+inner patch.
+
+# LICENCE #####################################################################
+Copyright 2020, David Seus
+david.seus[at]ians.uni-stuttgart.de
+This program is free software: you can redistribute it and/or modify
+it under the terms of the GNU General Public License as published by
+the Free Software Foundation, either version 3 of the License, or
+(at your option) any later version.
+
+This program is distributed in the hope that it will be useful,
+but WITHOUT ANY WARRANTY; without even the implied warranty of
+MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+GNU General Public License for more details.
+
+You should have received a copy of the GNU General Public License
+along with this program.  If not, see <http://www.gnu.org/licenses/>.
+###############################################################################
 """
 import dolfin as df
 import sympy as sym
diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/mesh_study/TP-TP-layered_soil_with_inner_patch-realistic-pure-dd.py b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/mesh_study/TP-TP-layered_soil_with_inner_patch-realistic-pure-dd.py
index f0ba034a78caba42ed2cd9bd54ebd4f55821631e..51a27ab65f0a203e503e7a0bbd961404cd2e1322 100755
--- a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/mesh_study/TP-TP-layered_soil_with_inner_patch-realistic-pure-dd.py
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/mesh_study/TP-TP-layered_soil_with_inner_patch-realistic-pure-dd.py
@@ -1,7 +1,27 @@
 #!/usr/bin/python3
-"""TP-TP layered soil with inner patch simulation.
-
-This program sets up an LDD simulation
+"""TP-TP layered soil with inner patch simulation, Copyright 2020, David Seus
+
+This program runs a mesh study for an LDD simulation for a layered soil
+substructuring with inner patch.
+This example assumes realistic parameters but remains pure dd and performs a
+mesh study.
+
+# LICENCE #####################################################################
+Copyright 2020, David Seus
+david.seus[at]ians.uni-stuttgart.de
+This program is free software: you can redistribute it and/or modify
+it under the terms of the GNU General Public License as published by
+the Free Software Foundation, either version 3 of the License, or
+(at your option) any later version.
+
+This program is distributed in the hope that it will be useful,
+but WITHOUT ANY WARRANTY; without even the implied warranty of
+MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+GNU General Public License for more details.
+
+You should have received a copy of the GNU General Public License
+along with this program.  If not, see <http://www.gnu.org/licenses/>.
+###############################################################################
 """
 import dolfin as df
 import sympy as sym
diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/mesh_study/TP-TP-layered_soil_with_inner_patch_mesh_study.py b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/mesh_study/TP-TP-layered_soil_with_inner_patch_mesh_study.py
index a32082c7d8d9e778141292fdaa46b428ef416222..782c6fea05874fa59b0420d3f039de45327f6bc6 100755
--- a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/mesh_study/TP-TP-layered_soil_with_inner_patch_mesh_study.py
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/mesh_study/TP-TP-layered_soil_with_inner_patch_mesh_study.py
@@ -1,7 +1,26 @@
 #!/usr/bin/python3
-"""TP-TP layered soil with inner patch simulation.
-
-This program sets up an LDD simulation
+"""TP-TP layered soil with inner patch simulation, Copyright 2020, David Seus
+
+This program runs a mesh study LDD simulation for a layered soil substructuring
+with inner patch.
+This example assumes realistic parameters and performs a mesh study
+
+# LICENCE #####################################################################
+Copyright 2020, David Seus
+david.seus[at]ians.uni-stuttgart.de
+This program is free software: you can redistribute it and/or modify
+it under the terms of the GNU General Public License as published by
+the Free Software Foundation, either version 3 of the License, or
+(at your option) any later version.
+
+This program is distributed in the hope that it will be useful,
+but WITHOUT ANY WARRANTY; without even the implied warranty of
+MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+GNU General Public License for more details.
+
+You should have received a copy of the GNU General Public License
+along with this program.  If not, see <http://www.gnu.org/licenses/>.
+###############################################################################
 """
 import dolfin as df
 import sympy as sym
@@ -240,23 +259,6 @@ L = {
          'nonwetting': Lnw6}
 }
 
-
-# interface_num : lambda parameter for the L-scheme on that interface.
-# Note that interfaces are numbered starting from 0, because
-# adjacent_subdomains is a list and not a dict. Historic fuckup, I know
-# We have defined above as interfaces
-# # interface_vertices introduces a global numbering of interfaces.
-# interface_def_points = [interface12_vertices,
-#                         interface23_vertices,
-#                         interface24_vertices,
-#                         interface25_vertices,
-#                         interface34_vertices,
-#                         interface36_vertices,
-#                         interface45_vertices_a,
-#                         interface45_vertices_b,
-#                         interface46_vertices,
-#                         interface56_vertices,
-#                         ]
 lambda_param = {
     0: {'wetting': lambda12_w,
          'nonwetting': lambda12_nw},