[cig-commits] r14485 - cs/cigma/trunk/src

[luis at geodynamics.org]

Author: luis
Date: 2009-03-27 11:53:42 -0700 (Fri, 27 Mar 2009)
New Revision: 14485

Modified:
   cs/cigma/trunk/src/fn_inclusion.cpp
Log:
Adjust constant in benchmark::circular_inclusion::Velocity

This function gives a much better fit to the velocity field.
We use the same formula as in the boundary condition, but the
differences stay nearly constant over the entire domain (probably
because we are not accounting for the variations in the near field).
How do we interpret the formula in the Schmid-Podladchikov paper?

Modified: cs/cigma/trunk/src/fn_inclusion.cpp
===================================================================
--- cs/cigma/trunk/src/fn_inclusion.cpp	2009-03-27 17:09:17 UTC (rev 14484)
+++ cs/cigma/trunk/src/fn_inclusion.cpp	2009-03-27 18:53:42 UTC (rev 14485)
@@ -39,10 +39,11 @@
     const double mu_m = 1;
     const double mu_i = 2;
     const double mag_shear = 1.0;
-    const double C = mag_shear * mu_m / (mu_m + mu_i);
+    //const double C = 2 * mag_shear * mu_m / (mu_m + mu_i);
+    const double C = mag_shear;
 
-    value[0] = -2 * C * x[0];
-    value[1] = +2 * C * x[1];
+    value[0] = -C*x[0];
+    value[1] = +C*x[1];
     value[2] = 0;
 
     return true;