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

[tan2 at geodynamics.org]

Author: tan2
Date: 2007-12-20 12:26:04 -0800 (Thu, 20 Dec 2007)
New Revision: 8961

Modified:
   mc/3D/CitcomS/trunk/lib/Regional_version_dependent.c
Log:
Removed obsolete code

Modified: mc/3D/CitcomS/trunk/lib/Regional_version_dependent.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Regional_version_dependent.c	2007-12-20 20:25:33 UTC (rev 8960)
+++ mc/3D/CitcomS/trunk/lib/Regional_version_dependent.c	2007-12-20 20:26:04 UTC (rev 8961)
@@ -229,22 +229,6 @@
     }          /* lev   */
 
 
-/*    do not need to rotate for the mesh grid for one regional problem   */
-
-
-  ro = -0.5*(M_PI/4.0)/E->lmesh.elx;
-  fo = 0.0;
-
-  E->sphere.dircos[1][1] = cos(ro)*cos(fo);
-  E->sphere.dircos[1][2] = cos(ro)*sin(fo);
-  E->sphere.dircos[1][3] = -sin(ro);
-  E->sphere.dircos[2][1] = -sin(fo);
-  E->sphere.dircos[2][2] = cos(fo);
-  E->sphere.dircos[2][3] = 0.0;
-  E->sphere.dircos[3][1] = sin(ro)*cos(fo);
-  E->sphere.dircos[3][2] = sin(ro)*sin(fo);
-  E->sphere.dircos[3][3] = cos(ro);
-
   for (j=1;j<=E->sphere.caps_per_proc;j++)   {
      regional_coord_of_cap(E,j,0);
      }
@@ -262,12 +246,6 @@
     }
     fflush(E->fp_out);
   }
-                   /* rotate the mesh to avoid two poles on mesh points */
-/*
-  for (j=1;j<=E->sphere.caps_per_proc;j++)   {
-     rotate_mesh(E,j,0);
-     }
-*/
 
   compute_angle_surf_area (E);   /* used for interpolation */