[cig-commits] r14895 - in seismo/2D/SPECFEM2D/branches/BIOT: . DATA

[cmorency at geodynamics.org]

Author: cmorency
Date: 2009-05-06 12:12:06 -0700 (Wed, 06 May 2009)
New Revision: 14895

Modified:
   seismo/2D/SPECFEM2D/branches/BIOT/DATA/CMTSOLUTION
   seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90
Log:
BIOT2D: Heaviside source time function error fixed


Modified: seismo/2D/SPECFEM2D/branches/BIOT/DATA/CMTSOLUTION
===================================================================
--- seismo/2D/SPECFEM2D/branches/BIOT/DATA/CMTSOLUTION	2009-05-06 19:10:14 UTC (rev 14894)
+++ seismo/2D/SPECFEM2D/branches/BIOT/DATA/CMTSOLUTION	2009-05-06 19:12:06 UTC (rev 14895)
@@ -3,7 +3,7 @@
 xs                              = 1600.          # source location x in meters
 zs                              = 2900.          # source location z in meters
 source_type                     = 1              # elastic force or acoustic pressure = 1 or moment tensor = 2
-time_function_type              = 1              # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
+time_function_type              = 5              # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
 f0                              = 15.0           # dominant source frequency (Hz) if not Dirac or Heaviside
 t0                              = 0.0            # time shift when multi sources (if one source, must be zero)
 angleforce                      = 0.             # angle of the source (for a force only)

Modified: seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90	2009-05-06 19:10:14 UTC (rev 14894)
+++ seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90	2009-05-06 19:12:06 UTC (rev 14895)
@@ -612,6 +612,8 @@
   allocate( t0(NSOURCE) )
   allocate( factor(NSOURCE) )
   allocate( angleforce(NSOURCE) )
+  allocate( hdur(NSOURCE) )
+  allocate( hdur_gauss(NSOURCE) )
   allocate( Mxx(NSOURCE) )
   allocate( Mxz(NSOURCE) )
   allocate( Mzz(NSOURCE) )