]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Remove two coefficients that aren't used right now.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Thu, 20 Jan 2011 18:50:32 +0000 (18:50 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Thu, 20 Jan 2011 18:50:32 +0000 (18:50 +0000)
git-svn-id: https://svn.dealii.org/trunk@23229 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-32/step-32.cc

index 8ecf15b0143207f787c0421a54bad69b07ed389e..1fb847f0b42780387476498c54e405b616a9276a 100644 (file)
@@ -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 <int dim>
   Tensor<1,dim> gravity_vector (const Point<dim> &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<dim> flow_problem;
 
       unsigned int ref = (dim == 2 ? 5 : 2);

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.