_Incomp_DFM(state, g_fine, CG, S, fluid, ...
        'Reconstruct', true, 'Dual', DG, 'src',src,'Wells', W,...
        'Update', true, 'Iterations',iterations,'Tolerance', tolerance, 'Iterator', 'gmres','c2cTrans',S2);
    
    
    
    %% Implicit transport solver
    tsolve   = @(state, dT) implicitTransport_DFM(state, g_fine, dT, rock, ...
        fluid, 'wells', W, 'src',src, ...
        'verbose', false,'nltol',1e-1);
    %% Solve initial pressure in reservoir
    
    % referance solution 
    refSol