      'n'  , [  1,   1]);

    
    %% Extract dual grid from g_fine
    DG = partitionDualDFM(g_fine);
    DG.lineedge = [];
    
    % store the dual grid in the grid structure
    g_fine.cells.level = zeros(g_fine.cells.num,1);
    g_fine.cells.level(DG.nn) = 0;
    g_fine.cells.level(DG.ee) = 1;
    g_fine.cells.level(DG.ii) = 2;
    
    % plot the dual grid
    if doPlotGrid
        clf, plotDual_DFM(g_fine, DG)
    end
    
    %% Generate primal grid
    p  = partitioningByAggregation(g_fine,DG);
  