[cig-commits] commit 2484 by ian.rose to /var/svn/dealii/aspect

[dealii.demon at gmail.com]

Revision 2484

Hack in rebound box

U   branches/freesurface/source/postprocess/topography.cc


http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=2484&peg=2484

Diff:
Modified: branches/freesurface/source/postprocess/topography.cc
===================================================================
--- branches/freesurface/source/postprocess/topography.cc	2014-04-11 17:56:39 UTC (rev 2483)
+++ branches/freesurface/source/postprocess/topography.cc	2014-04-11 18:04:12 UTC (rev 2484)
@@ -21,6 +21,7 @@
 
 #include <aspect/postprocess/topography.h>
 
+#include <aspect/geometry_model/rebound.h>
 #include <aspect/geometry_model/box.h>
 #include <aspect/geometry_model/sphere.h>
 #include <aspect/geometry_model/spherical_shell.h>
@@ -51,6 +52,14 @@
         vertical_gravity = true;
         relevant_boundary = (dim == 2 ? 3 : 5);
       }
+      else if(GeometryModel::ReboundBox<dim> *gm = dynamic_cast<GeometryModel::ReboundBox<dim> *>
+                                             (const_cast<GeometryModel::Interface<dim> *>(&this->get_geometry_model())))
+      {
+        Point<dim> extents = gm->get_extents();
+        reference_height = extents[dim-1];
+        vertical_gravity = true;
+        relevant_boundary = (dim == 2 ? 3 : 5);
+      }
       else if(GeometryModel::Sphere<dim> *gm = dynamic_cast<GeometryModel::Sphere<dim> *>
                                              (const_cast<GeometryModel::Interface<dim> *>(&this->get_geometry_model())))
       {