[cig-commits] r16015 - seismo/2D/SPECFEM2D/trunk

[cmorency at geodynamics.org]

Author: cmorency
Date: 2009-11-20 07:51:39 -0800 (Fri, 20 Nov 2009)
New Revision: 16015

Modified:
   seismo/2D/SPECFEM2D/trunk/specfem2D.F90
Log:
Fixed flag for attenuation and anisotropic calculation with p_sv = true, following Dylan Mikesell email.


Modified: seismo/2D/SPECFEM2D/trunk/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2009-11-19 23:11:05 UTC (rev 16014)
+++ seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2009-11-20 15:51:39 UTC (rev 16015)
@@ -1039,7 +1039,7 @@
   print*, '*************** WARNING ***************'
   stop
   endif
-  if(p_sv .and. (TURN_ATTENUATION_ON .or. TURN_ANISOTROPY_ON)) then
+  if(.not. p_sv .and. (TURN_ATTENUATION_ON .or. TURN_ANISOTROPY_ON)) then
   print*, '*************** WARNING ***************'
   print*, 'Attenuation and anisotropy are not implemented for surface (membrane) waves calculation'
   print*, '*************** WARNING ***************'
@@ -3796,10 +3796,10 @@
 
 ! Ricker (second derivative of a Gaussian) source time function
       if(time_function_type(i_source) == 1) then
-!        source_time_function(i_source,it) = - factor(i_source) * (ONE-TWO*aval(i_source)*(time-t0(i_source))**2) * &
-!                                           exp(-aval(i_source)*(time-t0(i_source))**2)
-        source_time_function(i_source,it) = - factor(i_source) * TWO*aval(i_source)*sqrt(aval(i_source))*&
-                                            (time-t0(i_source))/pi * exp(-aval(i_source)*(time-t0(i_source))**2)
+        source_time_function(i_source,it) = - factor(i_source) * (ONE-TWO*aval(i_source)*(time-t0(i_source))**2) * &
+                                           exp(-aval(i_source)*(time-t0(i_source))**2)
+!        source_time_function(i_source,it) = - factor(i_source) * TWO*aval(i_source)*sqrt(aval(i_source))*&
+!                                            (time-t0(i_source))/pi * exp(-aval(i_source)*(time-t0(i_source))**2)
 
 ! first derivative of a Gaussian source time function
       else if(time_function_type(i_source) == 2) then