[cig-commits] r16339 - in seismo/2D/SPECFEM2D/trunk: . DATA
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Thu Feb 25 06:17:07 PST 2010
Author: dkomati1
Date: 2010-02-25 06:17:06 -0800 (Thu, 25 Feb 2010)
New Revision: 16339
Added:
seismo/2D/SPECFEM2D/trunk/DATA/SOURCE
seismo/2D/SPECFEM2D/trunk/DATA/SOURCE.adj
seismo/2D/SPECFEM2D/trunk/DATA/SOURCE.regular
Removed:
seismo/2D/SPECFEM2D/trunk/DATA/CMTSOLUTION
seismo/2D/SPECFEM2D/trunk/DATA/CMTSOLUTION.adj
seismo/2D/SPECFEM2D/trunk/DATA/CMTSOLUTION.regular
Modified:
seismo/2D/SPECFEM2D/trunk/README_MANUAL.txt
seismo/2D/SPECFEM2D/trunk/constants.h
seismo/2D/SPECFEM2D/trunk/meshfem2D.F90
Log:
renamed DATA/CMTSOLUTION to DATA/SOURCE (because that file is not at all in the classical CMTSOLUTION format,
and the units for the moment tensor source are N.m, not dyne.cm)
Deleted: seismo/2D/SPECFEM2D/trunk/DATA/CMTSOLUTION
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/CMTSOLUTION 2010-02-25 05:02:27 UTC (rev 16338)
+++ seismo/2D/SPECFEM2D/trunk/DATA/CMTSOLUTION 2010-02-25 14:17:06 UTC (rev 16339)
@@ -1,13 +0,0 @@
-#source 1
-source_surf = .false. # source inside the medium or at the surface
-xs = 0. # source location x in meters
-zs = -100. # 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
-f0 = 5.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.0 # angle of the source (for a force only)
-Mxx = 1. # Mxx component (for a moment tensor source only)
-Mzz = 1. # Mzz component (for a moment tensor source only)
-Mxz = 0. # Mxz component (for a moment tensor source only)
-factor = 1.d10 # amplification factor
Deleted: seismo/2D/SPECFEM2D/trunk/DATA/CMTSOLUTION.adj
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/CMTSOLUTION.adj 2010-02-25 05:02:27 UTC (rev 16338)
+++ seismo/2D/SPECFEM2D/trunk/DATA/CMTSOLUTION.adj 2010-02-25 14:17:06 UTC (rev 16339)
@@ -1,13 +0,0 @@
-# source 1
-source_surf = .false. # source inside the medium or at the surface
-xs = 2500. # source location x in meters
-zs = 2500. # 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
-f0 = 3 # dominant source frequency (Hz) if not Dirac or Heaviside
-t0 = 0. # offset of the source, irrelevant if NSOURCE=1
-angleforce = 00. # angle of the source (for a force only)
-Mxx = 1. # Mxx component (for a moment tensor source only)
-Mzz = 1. # Mzz component (for a moment tensor source only)
-Mxz = 0. # Mxz component (for a moment tensor source only)
-factor = 1.d10 # amplification factor
Deleted: seismo/2D/SPECFEM2D/trunk/DATA/CMTSOLUTION.regular
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/CMTSOLUTION.regular 2010-02-25 05:02:27 UTC (rev 16338)
+++ seismo/2D/SPECFEM2D/trunk/DATA/CMTSOLUTION.regular 2010-02-25 14:17:06 UTC (rev 16339)
@@ -1,13 +0,0 @@
-# source 1
-source_surf = .false. # source inside the medium or at the surface
-xs = 1600. # source location x in meters
-zs = 2900. # source location z in meters
-source_type = 2 # 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
-f0 = 15 # dominant source frequency (Hz) if not Dirac or Heaviside
-t0 = 0.
-angleforce = 00. # angle of the source (for a force only)
-Mxx = 1. # Mxx component (for a moment tensor source only)
-Mzz = 1. # Mzz component (for a moment tensor source only)
-Mxz = 0. # Mxz component (for a moment tensor source only)
-factor = 1.d10 # amplification factor
Copied: seismo/2D/SPECFEM2D/trunk/DATA/SOURCE (from rev 16331, seismo/2D/SPECFEM2D/trunk/DATA/CMTSOLUTION)
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/SOURCE (rev 0)
+++ seismo/2D/SPECFEM2D/trunk/DATA/SOURCE 2010-02-25 14:17:06 UTC (rev 16339)
@@ -0,0 +1,17 @@
+#source 1
+source_surf = .false. # source inside the medium or at the surface
+xs = 0. # source location x in meters
+zs = -100. # 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
+f0 = 5.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.0 # angle of the source (for a force only)
+Mxx = 1. # Mxx component (for a moment tensor source only)
+Mzz = 1. # Mzz component (for a moment tensor source only)
+Mxz = 0. # Mxz component (for a moment tensor source only)
+factor = 1.d10 # amplification factor
+
+# The components of a moment tensor source must be given in N.m,
+# not in dyne.cm as in the DATA/CMTSOLUTION source file of the 3D version of the code.
+
Copied: seismo/2D/SPECFEM2D/trunk/DATA/SOURCE.adj (from rev 16331, seismo/2D/SPECFEM2D/trunk/DATA/CMTSOLUTION.adj)
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/SOURCE.adj (rev 0)
+++ seismo/2D/SPECFEM2D/trunk/DATA/SOURCE.adj 2010-02-25 14:17:06 UTC (rev 16339)
@@ -0,0 +1,17 @@
+# source 1
+source_surf = .false. # source inside the medium or at the surface
+xs = 2500. # source location x in meters
+zs = 2500. # 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
+f0 = 3 # dominant source frequency (Hz) if not Dirac or Heaviside
+t0 = 0. # offset of the source, irrelevant if NSOURCE=1
+angleforce = 00. # angle of the source (for a force only)
+Mxx = 1. # Mxx component (for a moment tensor source only)
+Mzz = 1. # Mzz component (for a moment tensor source only)
+Mxz = 0. # Mxz component (for a moment tensor source only)
+factor = 1.d10 # amplification factor
+
+# The components of a moment tensor source must be given in N.m,
+# not in dyne.cm as in the DATA/CMTSOLUTION source file of the 3D version of the code.
+
Copied: seismo/2D/SPECFEM2D/trunk/DATA/SOURCE.regular (from rev 16331, seismo/2D/SPECFEM2D/trunk/DATA/CMTSOLUTION.regular)
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/SOURCE.regular (rev 0)
+++ seismo/2D/SPECFEM2D/trunk/DATA/SOURCE.regular 2010-02-25 14:17:06 UTC (rev 16339)
@@ -0,0 +1,17 @@
+# source 1
+source_surf = .false. # source inside the medium or at the surface
+xs = 1600. # source location x in meters
+zs = 2900. # source location z in meters
+source_type = 2 # 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
+f0 = 15 # dominant source frequency (Hz) if not Dirac or Heaviside
+t0 = 0.
+angleforce = 00. # angle of the source (for a force only)
+Mxx = 1. # Mxx component (for a moment tensor source only)
+Mzz = 1. # Mzz component (for a moment tensor source only)
+Mxz = 0. # Mxz component (for a moment tensor source only)
+factor = 1.d10 # amplification factor
+
+# The components of a moment tensor source must be given in N.m,
+# not in dyne.cm as in the DATA/CMTSOLUTION source file of the 3D version of the code.
+
Modified: seismo/2D/SPECFEM2D/trunk/README_MANUAL.txt
===================================================================
--- seismo/2D/SPECFEM2D/trunk/README_MANUAL.txt 2010-02-25 05:02:27 UTC (rev 16338)
+++ seismo/2D/SPECFEM2D/trunk/README_MANUAL.txt 2010-02-25 14:17:06 UTC (rev 16339)
@@ -240,8 +240,9 @@
In section "# source parameters":
The code now support multi sources.
NSOURCE is the number of source.
-Parameters of the sources are displayed in the file CMTSOLUTION, which must be
-in the directory DATA.
+Parameters of the sources are displayed in the file SOURCE, which must be
+in the directory DATA. The components of a moment tensor source must be given in N.m,
+not in dyne.cm as in the DATA/CMTSOLUTION source file of the 3D version of the code.
In section "# receiver line parameters for seismograms":
SAVE_FORWARD determines if the last frame of a forward simulation is saved (.true.) or not (.false)
Modified: seismo/2D/SPECFEM2D/trunk/constants.h
===================================================================
--- seismo/2D/SPECFEM2D/trunk/constants.h 2010-02-25 05:02:27 UTC (rev 16338)
+++ seismo/2D/SPECFEM2D/trunk/constants.h 2010-02-25 14:17:06 UTC (rev 16339)
@@ -85,8 +85,8 @@
! uncomment this to write to file instead
! integer, parameter :: IOUT = 41
-! number of lines per source in CMTSOLUTION file
- integer, parameter :: NLINES_PER_CMTSOLUTION_SOURCE = 13
+! number of lines per source in SOURCE file
+ integer, parameter :: NLINES_PER_SOURCE = 13
! flags for absorbing boundaries
integer, parameter :: IBOTTOM = 1
Modified: seismo/2D/SPECFEM2D/trunk/meshfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/meshfem2D.F90 2010-02-25 05:02:27 UTC (rev 16338)
+++ seismo/2D/SPECFEM2D/trunk/meshfem2D.F90 2010-02-25 14:17:06 UTC (rev 16339)
@@ -604,22 +604,22 @@
allocate(Mzz(NSOURCE))
allocate(factor(NSOURCE))
- open(unit=IIN_SOURCE,file='DATA/CMTSOLUTION',iostat=ios,status='old',action='read')
- if(ios /= 0) stop 'error opening CMTSOLUTION file'
+ open(unit=IIN_SOURCE,file='DATA/SOURCE',iostat=ios,status='old',action='read')
+ if(ios /= 0) stop 'error opening SOURCE file'
icounter = 0
do while(ios == 0)
read(IIN_SOURCE,"(a)",iostat=ios) dummystring
if(ios == 0) icounter = icounter + 1
enddo
close(IIN_SOURCE)
- if(mod(icounter,NLINES_PER_CMTSOLUTION_SOURCE) /= 0) &
- stop 'total number of lines in CMTSOLUTION file should be a multiple of NLINES_PER_CMTSOLUTION_SOURCE'
- NSOURCES = icounter / NLINES_PER_CMTSOLUTION_SOURCE
- if(NSOURCES < 1) stop 'need at least one source in CMTSOLUTION file'
+ if(mod(icounter,NLINES_PER_SOURCE) /= 0) &
+ stop 'total number of lines in SOURCE file should be a multiple of NLINES_PER_SOURCE'
+ NSOURCES = icounter / NLINES_PER_SOURCE
+ if(NSOURCES < 1) stop 'need at least one source in SOURCE file'
if(NSOURCES /= NSOURCE) &
stop 'total number of sources read is different than declared in Par_file'
- open(unit=IIN_SOURCE,file='DATA/CMTSOLUTION',status='old',action='read')
+ open(unit=IIN_SOURCE,file='DATA/SOURCE',status='old',action='read')
do i_source=1,NSOURCE
call read_value_logical(IIN_SOURCE,IGNORE_JUNK,source_surf(i_source))
call read_value_double_precision(IIN_SOURCE,IGNORE_JUNK,xs(i_source))
More information about the CIG-COMMITS
mailing list