[cig-commits] commit 1879 by dannberg to /var/svn/dealii/aspect

[dealii.demon at gmail.com]

Revision 1879

fix a bug in the adiabatic initial conditions for 3D models

U   branches/j-dannberg/source/initial_conditions/adiabatic.cc


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

Diff:
Modified: branches/j-dannberg/source/initial_conditions/adiabatic.cc
===================================================================
--- branches/j-dannberg/source/initial_conditions/adiabatic.cc	2013-09-11 14:44:03 UTC (rev 1878)
+++ branches/j-dannberg/source/initial_conditions/adiabatic.cc	2013-09-11 14:48:13 UTC (rev 1879)
@@ -93,12 +93,18 @@
                 }
               else if (dim==3)
                 {
-                  mid_point(0) = inner_radius * sin(angle) * cos(angle),
-                  mid_point(1) = inner_radius * sin(angle) * sin(angle);
-                  if (dynamic_cast<const GeometryModel::SphericalShell<dim>&> (*this->geometry_model).opening_angle() == 90)
-                    mid_point(2) = inner_radius * cos(4.0/3.0 * angle);
+            	  if (dynamic_cast<const GeometryModel::SphericalShell<dim>&> (*this->geometry_model).opening_angle() == 90)
+            	  {
+            		mid_point(0) = std::sqrt(inner_radius*inner_radius/3),
+            		mid_point(1) = std::sqrt(inner_radius*inner_radius/3),
+            		mid_point(2) = std::sqrt(inner_radius*inner_radius/3);
+            	  }
                   else
+                  {
+                    mid_point(0) = inner_radius * sin(angle) * cos(angle),
+                    mid_point(1) = inner_radius * sin(angle) * sin(angle),
                     mid_point(2) = inner_radius * cos(angle);
+                  }
                 }
             }
           else if (dynamic_cast <const GeometryModel::Box<dim>*> (this->geometry_model) != 0)