[cig-commits] r16680 - seismo/2D/SPECFEM2D/trunk
pieyre at geodynamics.org
pieyre at geodynamics.org
Mon May 10 02:59:43 PDT 2010
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
More information about the CIG-COMMITS
mailing list