From af3b78460b0450260f5fcef23ad2690e1ac1c93b Mon Sep 17 00:00:00 2001
From: Kai Bao <kai.bao@sintef.no>
Date: Fri, 5 Jun 2015 10:41:54 +0200
Subject: [PATCH] adding SHRATE to PolymerPropsAd

---
 opm/polymer/fullyimplicit/PolymerPropsAd.cpp | 10 ++++++++++
 opm/polymer/fullyimplicit/PolymerPropsAd.hpp |  3 +++
 2 files changed, 13 insertions(+)

diff --git a/opm/polymer/fullyimplicit/PolymerPropsAd.cpp b/opm/polymer/fullyimplicit/PolymerPropsAd.cpp
index 0fb23f6e7..e215562a5 100644
--- a/opm/polymer/fullyimplicit/PolymerPropsAd.cpp
+++ b/opm/polymer/fullyimplicit/PolymerPropsAd.cpp
@@ -102,6 +102,16 @@ namespace Opm {
         return polymer_props_.plyshlogRefTemp();
     }
 
+    double PolymerPropsAd::shrate() const
+    {
+        if (polymer_props_.hasShrate()) {
+            return polymer_props_.shrate();
+        } else {
+            OPM_THROW(std::logic_error, "the SHRATE keyword is not specified while requested \n");
+        }
+    }
+
+
     double
     PolymerPropsAd::viscMult(double c) const
     {
diff --git a/opm/polymer/fullyimplicit/PolymerPropsAd.hpp b/opm/polymer/fullyimplicit/PolymerPropsAd.hpp
index 07ddc7a11..4630a9d7b 100644
--- a/opm/polymer/fullyimplicit/PolymerPropsAd.hpp
+++ b/opm/polymer/fullyimplicit/PolymerPropsAd.hpp
@@ -62,6 +62,9 @@ namespace Opm {
         /// \ return    The reference temperature in PLYSHLOG keyword
         double plyshlogRefTemp() const;
 
+        /// \ return   the value of SHRATE
+        double shrate() const;
+
         double viscMult(double c) const; // multipler interpolated from PLYVISC table
 
 		typedef AutoDiffBlock<double> ADB;