diff --git a/marlpde/parameters.py b/marlpde/parameters.py
index b8ce451..5c071a4 100644
--- a/marlpde/parameters.py
+++ b/marlpde/parameters.py
@@ -133,7 +133,7 @@ def post_init(self):
         self.rhos = self.rhos0
         self.Xstar = self.D0Ca / self.sedimentationrate
         self.Tstar = self.Xstar / self.sedimentationrate
-        self.b = (self.b/1e4) * 0.8**3 / (0.8*3)
+        self.b = (self.b/1e4)
         self.m2 = self.m1
         self.n2 = self.n1
         self.DCa = self.D0Ca
@@ -205,8 +205,8 @@ class Solver():
     '''
     dt: float = 1.e-6
     # t_range is the integration time in units of T*.
-    t_range: int = 30
-    solver: str = "explicit"
+    t_range: int = 1
+    solver: str = "scipy"
     # Beware that "scheme" and "adaptive" will only be propagated if you have 
     # chosen py-pde's native "explicit" solver above.
     scheme: str = "euler"