Dear MRST Developers,
I'm working with the code, initially constructed for a two-phase (oil-water) simulation, we are numerically simulating hydraulic fracturing (with one water injecting well) - coupled flow-geomechanical simulation. I am trying to rebuild it to work on a three-phase (blackoil model); Some parts of the code are attached, I added third phase (gas) into the simulation, previous version of code was done for two phase. However, I am encountering some errors, after added gas to the simulation:
Error using cellfun
Non-scalar in Uniform output, at index 1, output 1.
Set 'UniformOutput' to false.
Error in PhysicalModel/stepFunction (line 489)
            if any(cellfun(@(d) ~all(isfinite(d)), dx))
                   ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in NonLinearSolver/solveMinistep (line 356)
                    model.stepFunction(state, state0, dt, drivingForces, ...
                    ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in NonLinearSolver/solveTimestep (line 220)
                    solveMinistep(solver, model, state, state0_inner, dt, drivingForces);
                    ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in MechFluidFixedStressSplitModel/stepFunction (line 107)
            [mstate, mreport] = mechsolver.solveTimestep(mstate0, dt, mechModel, ...
                                ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in NonLinearSolver/solveMinistep (line 356)
                    model.stepFunction(state, state0, dt, drivingForces, ...
                    ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in NonLinearSolver/solveTimestep (line 220)
                    solveMinistep(solver, model, state, state0_inner, dt, drivingForces);
                    ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in simulateScheduleAD (line 250)
            [state, report] = solver.solveTimestep(state0, dt(i), model,...
                              ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in THREE_phase_Base (line 755)
    [w, s, report] = simulateScheduleAD(initState, model, schedule, ...
Thank you in advance!
Timur Sabirov