OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Use S_o=1 as initial condition.

Browse Source 
This is the case in 'sim_simple.cpp', too.
This commit is contained in:
Bård Skaflestad 2013-05-03 20:28:36 +02:00
parent 9b65338461
commit fc5a3f181e
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
mrst/sim_simple.m
View File

@ -14,8 +14,9 @@ function x = sim_simple(cartDims, physDims, tf, verb)
src = addSource([], [1, g.cells.num], [1, -1], 'sat', [ 1, 0 ; 0, 1]);
state = initState(g, [], 0, repmat([ 0.5, 0.5 ], [g.cells.num, 1]));
state = incompTPFA(state, g, T, fluid, 'src', src);
s0 = [0, 1];
state = incompTPFA(initState(g, [], 0, s0), ...
g, T, fluid, 'src', src);
if nargin < 4, verb = false; end
state = implicitTransport(state, g, tf, rock, fluid, ...