Fix adiabatic density gradients

geodynamics · Jul 14, 2023 · bdb1b22 · bdb1b22
1 parent 0620896
commit bdb1b22
Showing 1 changed file with 10 additions and 3 deletions.
diff --git a/cookbooks/vankeken_subduction/plugin/van_Keken_mesh.cc b/cookbooks/vankeken_subduction/plugin/van_Keken_mesh.cc
@@ -274,9 +274,16 @@ namespace aspect
 
     template <int dim>
     double
-    vanKeken<dim>::depth(const Point<dim> &) const
+    vanKeken<dim>::depth(const Point<dim> &position) const
     {
-      return 0.0;
+      // Get the surface x (,y) point
+      Point<dim-1> surface_point;
+      for (unsigned int d=0; d<dim-1; ++d)
+        surface_point[d] = position[d];
+
+      std::vector<double> extents = {660e3, 600e3};
+      const double d = extents[dim-1] - (position(dim-1));
+      return std::min (std::max (d, 0.), maximal_depth());
     }
 
 
@@ -305,7 +312,7 @@ namespace aspect
     double
     vanKeken<dim>::maximal_depth() const
     {
-      return 0.0;
+      return 600e3;
     }