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

[walter at geodynamics.org]

Author: walter
Date: 2007-12-07 15:49:14 -0800 (Fri, 07 Dec 2007)
New Revision: 8476

Modified:
   seismo/2D/SPECFEM2D/trunk/constants.h
   seismo/2D/SPECFEM2D/trunk/convolve_source_timefunction.f90
Log:
changed 2.628 to 1.628 for decay rate in 2D code, following Qinya's same change in 3D code


Modified: seismo/2D/SPECFEM2D/trunk/constants.h
===================================================================
--- seismo/2D/SPECFEM2D/trunk/constants.h	2005-04-16 22:07:53 UTC (rev 8475)
+++ seismo/2D/SPECFEM2D/trunk/constants.h	2007-12-07 23:49:14 UTC (rev 8476)
@@ -56,7 +56,7 @@
   integer, parameter :: NUM_ITER = 4
 
 ! error function source decay rate for Heaviside
-  double precision, parameter :: SOURCE_DECAY_RATE = 2.628d0
+  double precision, parameter :: SOURCE_DECAY_RATE = 1.628d0
 
 ! display non lineaire pour rehausser les faibles amplitudes sur les images PNM
   double precision, parameter :: POWER_DISPLAY_PNM = 0.30d0

Modified: seismo/2D/SPECFEM2D/trunk/convolve_source_timefunction.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/convolve_source_timefunction.f90	2005-04-16 22:07:53 UTC (rev 8475)
+++ seismo/2D/SPECFEM2D/trunk/convolve_source_timefunction.f90	2007-12-07 23:49:14 UTC (rev 8476)
@@ -38,9 +38,6 @@
   read(33,*) triangle
   close(33)
 
-! for Gaussian use 1.66667*hdur to get roughly a triangle with half-duration hdur
-  if(.not.triangle) hdur = hdur * 5. / 3.
-
 ! allocate arrays
   allocate(time(nlines),sem(nlines),sem_fil(nlines))