mment">% inexact solution</span>
    msSol = psolve(msSol);

    <span class="comment">% store data for plotting</span>
    caseData.msSol{casenr} = msSol;
    caseData.g_fine{casenr} = g_fine;
    caseData.rock{casenr} = rock;
</pre><h2>Transport loop<a name="15"></a></h2><p>Solve a linear transport problem with fixed fluxes computed from the pressure equation</p><pre class="codeinput">    <span class="comment">% Simulation time is 1 PVI</span>
    <span class="keyword">if</span> casenr == 1
    T      = s