fine.cells.num, 1]);

    <span class="comment">% use the cubic assumption for the fracture permeability</span>
    <span class="comment">% with unit porosity</span>
    kf = g_fine.faces.apertures.^2/12;
    <span class="keyword">for</span> i = 2:-1:1 <span class="comment">% note that hybrid cells with tag 3</span>
        <span class="comment">% i.e the Cartesian grid is omitted from the loop)</span>
        [fcells,ncells] = getFractureCells(g_fine,i);
        rock.perm(fcells) = kf(g_fine.cells.tags(fce