[cig-commits] [commit] devel: added a test to check that the relaxation times found for viscoelastic materials are positive (f4001a1)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem2d

On branch  : devel
Link       : https://github.com/geodynamics/specfem2d/compare/9ff6f713a2d20b6f88543ef354dbdbb9113f91c4...9bc9f1b010e4097a44a4bbddee81572460ca6e58

>---------------------------------------------------------------

commit f4001a1d6eba252ec03a3c9a4562b9066d2573ba
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Mon Jun 9 00:28:53 2014 +0200

    added a test to check that the relaxation times found for viscoelastic materials are positive


>---------------------------------------------------------------

f4001a1d6eba252ec03a3c9a4562b9066d2573ba
 src/specfem2D/attenuation_model.f90 | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/src/specfem2D/attenuation_model.f90 b/src/specfem2D/attenuation_model.f90
index e473d79..38a3be0 100644
--- a/src/specfem2D/attenuation_model.f90
+++ b/src/specfem2D/attenuation_model.f90
@@ -117,6 +117,9 @@
 
   endif
 
+  if(any(tau_sigma_nu1 < 0.d0) .or. any(tau_sigma_nu2 < 0.d0) .or. any(tau_epsilon_nu1 < 0.d0) .or. any(tau_epsilon_nu2 < 0.d0)) &
+       stop 'error: negative relaxation time found for a viscoelastic material'
+
 ! When implementing viscoelasticity according to the Carcione 1993 paper, attenuation is
 ! non-causal rather than causal. We fixed the problem by using equations in Carcione's
 ! 2004 paper and his 2007 book. See also file doc/problem_attenuation_reference_Specfem2D_fixed_by_Xie_Zhinan.pdf