[cig-commits] r22777 - seismo/3D/SPECFEM3D/trunk/src/generate_databases
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Thu Sep 5 06:41:38 PDT 2013
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
More information about the CIG-COMMITS
mailing list