[cig-commits] [commit] master: Update solidus.cc (f1b6231)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/aspect

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/a21aaa79d75b89a3d50d5c865e7dc51d309da9b9...b2eaff9e459f8351633e8b4b43c1284c90373873

>---------------------------------------------------------------

commit f1b623146818b33c06a2679746ff332c06ba6c74
Author: SiqiZhang <SiqiZhang at users.noreply.github.com>
Date:   Fri May 16 11:30:52 2014 -0500

    Update solidus.cc


>---------------------------------------------------------------

f1b623146818b33c06a2679746ff332c06ba6c74
 source/initial_conditions/solidus.cc | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/source/initial_conditions/solidus.cc b/source/initial_conditions/solidus.cc
index 89c6ef0..81f613d 100644
--- a/source/initial_conditions/solidus.cc
+++ b/source/initial_conditions/solidus.cc
@@ -103,7 +103,7 @@ namespace aspect
         if(dim==3)
             Phi=atan2(position(0),position(2));
         else
-            Phi=M_PI;
+            Phi=0.;
         double T_solidus,T_perturbation;
         double litho_thick_theta;
         double Depth=this->geometry_model->depth(position);
@@ -124,7 +124,7 @@ namespace aspect
         else
             T_solidus=Solidus_curve.T(0,sqrt(position.square()))+deltaT;
 
-        T_perturbation=Depth/( this->geometry_model->maximal_depth() )*Magnitude_T*sin(n*Theta)*sin(n*Phi);
+        T_perturbation=Depth/( this->geometry_model->maximal_depth() )*Magnitude_T*sin(n*Theta)*cos(n*Phi);
         return T_solidus+T_perturbation;
     }