[cig-commits] r8962 - mc/3D/CitcomS/trunk/lib

[tan2 at geodynamics.org]

Author: tan2
Date: 2007-12-20 12:27:41 -0800 (Thu, 20 Dec 2007)
New Revision: 8962

Modified:
   mc/3D/CitcomS/trunk/lib/Full_geometry_cartesian.c
Log:
Fine tune the mesh so that each cap has the same surface area

Modified: mc/3D/CitcomS/trunk/lib/Full_geometry_cartesian.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Full_geometry_cartesian.c	2007-12-20 20:26:04 UTC (rev 8961)
+++ mc/3D/CitcomS/trunk/lib/Full_geometry_cartesian.c	2007-12-20 20:27:41 UTC (rev 8962)
@@ -85,7 +85,9 @@
   E->sphere.caps = 12;
   E->sphere.max_connections = 6;
 
-  offset = 10.0/180.0*M_PI;
+  /* adjust the corner coordinates so that the size (surface area) of
+     each cap is about the same. */
+  offset = 9.736/180.0*M_PI;
 
   for (i=1;i<=4;i++)  {
     E->sphere.cap[(i-1)*3+1].theta[1] = 0.0;