Change default sheath_boundary_simple extrapolation to exponential, free

This was previousle exponential free if quantity was decreasing and Neumann if it was increasing
bendudson · Sep 11, 2024 · 5433abd · 5433abd
1 parent ac9210b
commit 5433abd
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/sheath_boundary_simple.cxx b/src/sheath_boundary_simple.cxx
@@ -128,15 +128,15 @@ SheathBoundarySimple::SheathBoundarySimple(std::string name, Options& alloptions
 
   density_boundary_mode = options["density_boundary_mode"]
     .doc("BC mode: 0=LimitFree, 1=ExponentialFree, 2=LinearFree")
-    .withDefault<BoutReal>(0);
+    .withDefault<BoutReal>(1);
 
   pressure_boundary_mode = options["pressure_boundary_mode"]
     .doc("BC mode: 0=LimitFree, 1=ExponentialFree, 2=LinearFree")
-    .withDefault<BoutReal>(0);
+    .withDefault<BoutReal>(1);
 
   temperature_boundary_mode = options["temperature_boundary_mode"]
     .doc("BC mode: 0=LimitFree, 1=ExponentialFree, 2=LinearFree")
-    .withDefault<BoutReal>(0);
+    .withDefault<BoutReal>(1);
 
 
 }