--- a/+multiblock/stitchSchemes.m	Thu Sep 14 16:00:36 2017 +0200
+++ b/+multiblock/stitchSchemes.m	Fri Sep 15 09:21:29 2017 +0200
@@ -11,7 +11,7 @@
 % Output parameters are cell arrays and cell matrices.
 %
 % Ex: [schms, D, H] = stitchSchemes(schmHand, order, schmParam, blocks, ms, conn, bound)
-function [schms, D, H,penalty] = stitchSchemes(schmHand, order, schmParam, grids, conn, bound,timeDep)
+function [schms, D, H,penalty,J] = stitchSchemes(schmHand, order, schmParam, grids, conn, bound,timeDep)
     default_arg('schmParam',[]);
     default_arg('timeDep','N');
     n_blocks = numel(grids);
@@ -48,6 +48,11 @@
         H{i,i} = schms{i}.H;
     end
     
+    J = cell(n_blocks,n_blocks);
+    for i = 1:n_blocks
+        J{i,i} = schms{i}.J;
+    end
+    
     %% Total system matrix
     
     % Differentiation terms