--- a/+scheme/Wave2d.m	Sun Dec 23 14:06:26 2018 +0100
+++ b/+scheme/Wave2d.m	Sun Dec 23 14:39:31 2018 +0100
@@ -158,7 +158,7 @@
             end
         end
 
-        function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)
+        function [closure, penalty] = interface(obj, boundary, neighbour_scheme, neighbour_boundary, type)
             % u denotes the solution in the own domain
             % v denotes the solution in the neighbour domain
             [e_u,d_u,s_u,gamm_u, halfnorm_inv] = obj.get_boundary_ops(boundary);