nt">%src = addSource(src,double(c1),rate,'sat',[1,1]);</span>
    <span class="comment">%src = addSource(src,double(c2),-rate,'sat',[1,1]);</span>
</pre><h2>Set up solution structures with only one phase<a name="11"></a></h2><pre class="codeinput">    refSol = initResSol(g_fine, bhp, [0, 1]);
    msSol = refSol;

    gravity <span class="string">off</span>
    verbose = false;
</pre><h2>Set up pressure and transport solvers<a name="12"></a></h2><pre class="codeinput">    <span class="comment">% Reference TP