= @(state, dT) implicitTransport_DFM(state, g_fine, dT, rock, <span class="keyword">...</span>
        fluid, <span class="string">'wells'</span>, W, <span class="string">'src'</span>,src, <span class="keyword">...</span>
        <span class="string">'verbose'</span>, false,<span class="string">'nltol'</span>,1e-1);
</pre><h2>Solve initial pressure in reservoir<a name="14"></a></h2><pre class="codeinput">    <span class="comment">% referance solution</span>
    refSol= r_psolve(refSol);

    <span class="co