fset{2};
    end
    
    
    t = tic;
    
    % filename for to save the grids
    grid_filename = ['grids-',int2str(casenr),'.mat'];
    
    % create the grid hierarchy
    [g_dual,g_fine] = createGridHierarchy(coarseGrid,fineGrid,'save',grid_filename,...
        'numElements', numberOfFinecells,'precision',1e-3,'box',box);
    toc(t)
    
    
    %% Set up permeability
    
    % uniform matrix permeability
    rock.perm = repmat(100*milli*darcy, [g_fine.cells.num, 1]);
    rock.poro = repmat(0.1    