--- a/+scheme/Scheme.m	Tue Feb 07 15:47:51 2017 +0100
+++ b/+scheme/Scheme.m	Tue Feb 07 16:09:02 2017 +0100
@@ -20,13 +20,12 @@
         %                           'l','r' or 'e','w','n','s'.
         %       type                is a string specifying the type of
         %                           boundary condition if there are several.
-        %       data                is a function returning the data that
-        %                           should be applied at the boundary.
         %       neighbour_scheme    is an instance of Scheme that should be
         %                           interfaced to.
         %       neighbour_boundary  is a string specifying which boundary to
         %                           interface to.
-        [closure, penalty] = boundary_condition(obj,boundary,type,data)
+        %       penalty  may be a cell array if there are several penalties with different weights
+        [closure, penalty] = boundary_condition(obj,boundary,type)
         [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)
 
         % Returns the number of degrees of freedom.
@@ -42,4 +41,4 @@
             [vv,vu] = schm_v.interface(bound_v,schm_u,bound_u);
         end
     end
-end
\ No newline at end of file
+end