t(0,'DefaultFigureWindowStyle','docked')

%% set up the case

% bounding box of the domain 
box = [0,0; 100,100];

% the two set of fractures are drawn in open-office 
fset{1} = struct('fractures','coarseFractures.odp',...
    'box',box, ...
    'precision',1e-3, ...
    'aperture',0.001, ...
    'tag', 1 ...
    );

fset{2} = struct('fractures','fineFractures.odp',...
    'box',box, ...
    'precision',1e-3, ...
    'aperture',0.0005, ...
    'tag', 2 ...
    );


% make an nonconforming coarse grid for co