= 1;
    end
    
    % compute the fine-scale transmissibilities
    S = computeTrans_DFM(g_fine, rock,'hybrid',true);
    
    
    % compute the fine-scale hybrid transmissibilities
    if isfield(g_fine,'hybridNeighbors')
        [g_fine,S2] = computeHybridTrans(g_fine,S);
    else
        S2 = [];
    end
    
    %% Define a simple two phase fluid object
    % Linear tracer transport
    fluid = initSimpleFluid('mu' , [   1,  1]*centi*poise     , ...
        'rho', [1014, 859]*kilogram/meter^3, ...
  