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

[pieyre at geodynamics.org]

Author: pieyre
Date: 2010-05-10 02:59:41 -0700 (Mon, 10 May 2010)
New Revision: 16680

Modified:
   seismo/2D/SPECFEM2D/trunk/specfem2D.F90
Log:
added back the correct formula for the Ricker source time funtion (which is second derivative of a Gaussian)


Modified: seismo/2D/SPECFEM2D/trunk/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2010-05-10 07:04:04 UTC (rev 16679)
+++ seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2010-05-10 09:59:41 UTC (rev 16680)
@@ -3755,10 +3755,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