[cig-commits] [commit] master: Update solidus.cc (c28ec77)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Tue May 20 12:19:59 PDT 2014
Repository : https://github.com/geodynamics/aspect
On branch : master
Link : https://github.com/geodynamics/aspect/compare/a21aaa79d75b89a3d50d5c865e7dc51d309da9b9...b2eaff9e459f8351633e8b4b43c1284c90373873
>---------------------------------------------------------------
commit c28ec773a71fae9cee3d157272b7883198b0cbe2
Author: SiqiZhang <SiqiZhang at users.noreply.github.com>
Date: Fri May 16 11:42:07 2014 -0500
Update solidus.cc
>---------------------------------------------------------------
c28ec773a71fae9cee3d157272b7883198b0cbe2
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 81f613d..6d5063e 100644
--- a/source/initial_conditions/solidus.cc
+++ b/source/initial_conditions/solidus.cc
@@ -116,7 +116,7 @@ namespace aspect
R1=(dynamic_cast< const GeometryModel::SphericalShell<dim> &>(this->get_geometry_model())).R1;
T_min=(this->get_boundary_temperature()).minimal_temperature();
- litho_thick_theta=litho_thick-Magnitude_lith*sin(n*Theta);
+ litho_thick_theta=litho_thick-Magnitude_lith*cos(n*Theta)*cos(n*Phi);
T_litho=Solidus_curve.T(0,R1-litho_thick_theta)+deltaT;
if(litho_thick_theta>0 && Depth<litho_thick_theta)
@@ -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)*cos(n*Phi);
+ T_perturbation=Depth/( this->geometry_model->maximal_depth() )*Magnitude_T*cos(n*Theta)*cos(n*Phi);
return T_solidus+T_perturbation;
}
More information about the CIG-COMMITS
mailing list