    radius = 0.05;
    
    W = addWell([], g_fine, rock, c1 ,      ...
        'Type', 'rate' , 'Val', rate, ...
        'Radius', radius, 'InnerProduct', 'ip_tpf', ...
        'Comp_i', [1, 0]);
    
    W = addWell(W, g_fine, rock, c2,      ...
        'Type', 'bhp' , 'Val', bhp, ...
        'Radius', radius, 'InnerProduct', 'ip_tpf', ...
        'Comp_i', [0, 1]);
    
    %W = [];
    
    % Alternatively a simple source
    src =[];
    %src = addSource(src,double(c1),rate,'sat',[1,1]);
    %src = add