From b001c580b94104d80d879a1443a6fcb778b939aa Mon Sep 17 00:00:00 2001
From: Kai Bao <kai.bao@sintef.no>
Date: Fri, 16 May 2014 13:32:34 +0200
Subject: [PATCH] Add drsmaxrel for the update of rs.

drsmaxrel is a relative limit.
---
 opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp b/opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp
index 5d3b8e212..3d15f5994 100644
--- a/opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp
+++ b/opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp
@@ -1275,13 +1275,13 @@ namespace {
             }
         }
 
-        const double drsmax = 1e9;
+        const double drsmaxrel = 0.2;
         const double drvmax = 1e9;//% same as in Mrst
         V rs;
         if (disgas) {
             const V rs_old = Eigen::Map<const V>(&state.gasoilratio()[0], nc);
             const V drs = isRs * dxvar;
-            const V drs_limited = sign(drs) * drs.abs().min(drsmax);
+            const V drs_limited = sign(drs) * drs.abs().min(rs_old.abs()*drsmaxrel);
             rs = rs_old - drs_limited;
         }
         V rv;