adds perm rescaling, testing mfd in transient works

amanzi · Jul 3, 2023 · 2d796b3 · 2d796b3
1 parent db8df78
commit 2d796b3
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 3 deletions.
diff --git 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
@@ -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_);
 }
 
 // -------------------------------------------------------------