From 2d796b3f1fdb77d25720a5987de7c7df19b524dc Mon Sep 17 00:00:00 2001
From: Ethan Coon <coonet@ornl.gov>
Date: Mon, 3 Jul 2023 16:49:32 -0400
Subject: [PATCH] adds perm rescaling, testing mfd in transient works

---
 src/pks/flow/richards.hh    | 2 +-
 src/pks/flow/richards_pk.cc | 5 +++--
 2 files changed, 4 insertions(+), 3 deletions(-)

diff --git a/src/pks/flow/richards.hh b/src/pks/flow/richards.hh
index 88087efc2..0bdf5507c 100644
--- a/src/pks/flow/richards.hh
+++ b/src/pks/flow/richards.hh
@@ -383,7 +383,7 @@ class Richards : public PK_PhysicalBDF_Default {
   Teuchos::RCP<Operators::PDE_DiffusionWithGravity> preconditioner_diff_;
   Teuchos::RCP<Operators::PDE_DiffusionWithGravity> face_matrix_diff_;
   Teuchos::RCP<Operators::PDE_Accumulation> preconditioner_acc_;
-  // double perm_scale_;
+  double perm_scale_;
 
   // flag to do jacobian and therefore coef derivs
   bool precon_used_;
diff --git a/src/pks/flow/richards_pk.cc b/src/pks/flow/richards_pk.cc
index c52d196b9..8bdf46452 100644
--- a/src/pks/flow/richards_pk.cc
+++ b/src/pks/flow/richards_pk.cc
@@ -93,8 +93,9 @@ Richards::ParseParameterList_()
     deform_key_ = Keys::readKey(*plist_, domain_, "deformation indicator", "base_porosity");
 
   // scaling for permeability for better "nondimensionalization"
-  // perm_scale_ = plist_->get<double>("permeability rescaling", 1.e7);
-  // S_->GetEvaluatorList(coef_key_).set<double>("permeability rescaling", perm_scale_);
+  perm_scale_ = plist_->get<double>("permeability rescaling", 1.e7);
+  S_->GetEvaluatorList(coef_key_).set<double>("permeability rescaling", perm_scale_);
+  S_->GetEvaluatorList(perm_key_).set<double>("rescaling factor", perm_scale_);
 }
 
 // -------------------------------------------------------------