From ae37eda5ac7caface3b71e41eb87e0c276280894 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Thu, 20 Jan 2011 18:50:32 +0000 Subject: [PATCH] Remove two coefficients that aren't used right now. git-svn-id: https://svn.dealii.org/trunk@23229 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-32/step-32.cc | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) diff --git a/deal.II/examples/step-32/step-32.cc b/deal.II/examples/step-32/step-32.cc index 8ecf15b014..1fb847f0b4 100644 --- a/deal.II/examples/step-32/step-32.cc +++ b/deal.II/examples/step-32/step-32.cc @@ -111,9 +111,6 @@ namespace EquationData const double T0 = 4000+273; /* K */ const double T1 = 700+273; /* K */ - const double g0 = 10.7; /* m / s^2 */ - const double g1 = 9.81; /* m / s^2 */ - const double year_in_seconds = 60*60*24*365.2425; const double end_time = 1e8 * year_in_seconds; @@ -131,6 +128,10 @@ namespace EquationData template Tensor<1,dim> gravity_vector (const Point &p) { +// interpolate the following values with a physically realistic model: +// const double g0 = 10.7; /* m / s^2 */ +// const double g1 = 9.81; /* m / s^2 */ + const double r = p.norm(); return -(1.245e-6 * r + 7.714e13/r/r) * p / p.norm(); } @@ -3413,7 +3414,7 @@ int main (int argc, char *argv[]) Utilities::System::MPI_InitFinalize mpi_initialization(argc, argv); - const int dim = 2; + const int dim = 3; BoussinesqFlowProblem flow_problem; unsigned int ref = (dim == 2 ? 5 : 2); -- 2.39.5