[cig-commits] r22777 - seismo/3D/SPECFEM3D/trunk/src/generate_databases

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-09-05 06:41:37 -0700 (Thu, 05 Sep 2013)
New Revision: 22777

Modified:
   seismo/3D/SPECFEM3D/trunk/src/generate_databases/calc_jacobian.f90
Log:
clarified the error message about negative Jacobians


Modified: seismo/3D/SPECFEM3D/trunk/src/generate_databases/calc_jacobian.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/generate_databases/calc_jacobian.f90	2013-09-04 19:51:44 UTC (rev 22776)
+++ seismo/3D/SPECFEM3D/trunk/src/generate_databases/calc_jacobian.f90	2013-09-05 13:41:37 UTC (rev 22777)
@@ -91,8 +91,8 @@
              xeta*(yxi*zgamma-ygamma*zxi) + &
              xgamma*(yxi*zeta-yeta*zxi)
 
-! can ignore negative jacobian in mesher if needed when debugging code
-      if(jacobian <= ZERO) call exit_MPI(myrank,'3D Jacobian undefined')
+! check that the Jacobian transform is invertible, i.e. that the Jacobian never becomes negative or null
+      if(jacobian <= ZERO) call exit_MPI(myrank,'error: negative or null 3D Jacobian found')
 
 !     invert the relation (Fletcher p. 50 vol. 2)
       xix = (yeta*zgamma-ygamma*zeta) / jacobian