From 6e7980ba47893b70a90a304b03147167ab29998e Mon Sep 17 00:00:00 2001 From: Xavier Raynaud Date: Thu, 14 Jun 2012 14:01:53 +0200 Subject: [PATCH] fixed bug in relperm computation. --- opm/polymer/IncompPropertiesDefaultPolymer.hpp | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/opm/polymer/IncompPropertiesDefaultPolymer.hpp b/opm/polymer/IncompPropertiesDefaultPolymer.hpp index 76198926c..3376780ce 100644 --- a/opm/polymer/IncompPropertiesDefaultPolymer.hpp +++ b/opm/polymer/IncompPropertiesDefaultPolymer.hpp @@ -82,10 +82,10 @@ namespace Opm kr[2*i] = krw(s[2*i]); kr[2*i+1] = kro(s[2*i+1]); if (dkrds != 0) { - dkrds[2*i] = krw_dsw(s[2*i]); - dkrds[2*i+3] = kro_dso(s[2*i+1]); - dkrds[2*i+1] = 0.0; - dkrds[2*i+2] = 0.0; + dkrds[4*i + 0] = krw_dsw(s[2*i]); + dkrds[4*i + 3] = kro_dso(s[2*i+1]); + dkrds[4*i + 1] = 0.0; + dkrds[4*i + 2] = 0.0; } } }