  CG = generateCoarseGrid(g_fine, p);
    CG = coarsenGeometry(CG);
    
    %% Setup a producer / injector pair of wells
    
    % Place the producer / injector in the large-scale fractures
    [fcells,ncells] = getFractureCells(g_fine,1);
    
    % Injector in lower left corner
    c1 = findClosestCell(g_fine,box(1,1),box(1,2),fcells);
    % Producer in upper right corner
    c2 = findClosestCell(g_fine,box(2,1),box(2,2),fcells);
    
    % set up the well
    rate = 10*meter^3/day;
    bhp  = 1*barsa;
