[cig-commits] r22843 - seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/meshfem3D

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-09-22 17:52:46 -0700 (Sun, 22 Sep 2013)
New Revision: 22843

Modified:
   seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/meshfem3D/add_topography_410_650.f90
Log:
added a stop statement to prevent the add_topography_410_650() routine from being used, since it has a bug that makes the mesher crash


Modified: seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/meshfem3D/add_topography_410_650.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/meshfem3D/add_topography_410_650.f90	2013-09-23 00:41:18 UTC (rev 22842)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/SPECFEM3D_GLOBE_SUNFLOWER/src/meshfem3D/add_topography_410_650.f90	2013-09-23 00:52:46 UTC (rev 22843)
@@ -47,6 +47,9 @@
   double precision r,theta,phi
   double precision gamma
 
+!! DK DK added this safety test for now
+  stop 'there is a currently a bug in this routine, it makes the mesher crash'
+
 ! we loop on all the points of the element
   do ia = 1,NGNOD
 
@@ -135,6 +138,9 @@
   double precision r,theta,phi
   double precision gamma
 
+!! DK DK added this safety test for now
+  stop 'there is a currently a bug in this routine, it makes the mesher crash'
+
   ! we loop on all GLL points of the element
   do k = 1,NGLLZ
      do j = 1,NGLLY