[cig-commits] r20673 - in seismo/2D/SPECFEM2D/trunk: DATA EXAMPLES/Abel_Brest EXAMPLES/Gmsh_example_MPI EXAMPLES/Gmsh_example_serial EXAMPLES/INDUSTRIAL_FORMAT EXAMPLES/M2_UPPA EXAMPLES/Rayleigh_wave_no_crack EXAMPLES/Rayleigh_wave_with_crack EXAMPLES/Tape2007 EXAMPLES/Tape2007_kernel EXAMPLES/Tromp2005 EXAMPLES/Tromp2005_kernel EXAMPLES/acoustic_poroelastic EXAMPLES/attenuation EXAMPLES/canyon EXAMPLES/canyon/mesh EXAMPLES/fluid_solid/fluid_solid_external_mesh EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom EXAMPLES/init_plane EXAMPLES/noise_layered/model_0 EXAMPLES/noise_layered/model_1 EXAMPLES/noise_layered/model_2 EXAMPLES/noise_uniform EXAMPLES/semi_infinite_homo doc/USER_MANUAL setup src/meshfem2D src/specfem2D
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Sun Sep 2 19:09:57 PDT 2012
Author: dkomati1
Date: 2012-09-02 19:09:56 -0700 (Sun, 02 Sep 2012)
New Revision: 20673
Modified:
seismo/2D/SPECFEM2D/trunk/DATA/SOURCE
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/SOURCE_Abel_Balanche_bathy_source_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/SOURCE.SqrCirc
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/SOURCE.SqrCirc
seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/SOURCE
seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/SOURCE_M2_UPPA
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_no_crack/SOURCE_Rayleigh_2D
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_with_crack/SOURCE_Rayleigh_2D
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_001
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_002
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_003
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_004
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_005
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_006
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_007
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_008
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_009
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_010
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_011
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_012
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_013
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_014
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_015
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_016
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_017
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_018
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_019
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_020
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_021
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_022
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_023
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_024
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_025
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/SOURCE_001
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/SOURCE_Tromp2005
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/SOURCE_Tromp2005_s100
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/SOURCE_Tromp2005
seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/SOURCE_acoustic_poroelastic
seismo/2D/SPECFEM2D/trunk/EXAMPLES/attenuation/SOURCE_attenuation_2D
seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/SOURCE_canyon
seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/mesh/canyon_absorbing_surface_file
seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/fluid_solid_external_mesh/SOURCE_fluid_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/SOURCE_fluid_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/SOURCE_fluid_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/SOURCE_Slave
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/SOURCE_noise
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/SOURCE_noise
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/SOURCE_noise
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/SOURCE_noise
seismo/2D/SPECFEM2D/trunk/EXAMPLES/semi_infinite_homo/SOURCE_elastic_2D
seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.pdf
seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.tex
seismo/2D/SPECFEM2D/trunk/setup/constants.h.in
seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_source_file.f90
seismo/2D/SPECFEM2D/trunk/src/meshfem2D/save_databases.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_Bielak_conditions.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/paco_beyond_critical.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/prepare_initialfield.F90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/set_sources.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
Log:
switched from angleforce to anglesource; added flexible reading of DATA/SOURCE; added possibility to use the incident field only for a plane wave
Modified: seismo/2D/SPECFEM2D/trunk/DATA/SOURCE
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/SOURCE 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/DATA/SOURCE 2012-09-03 02:09:56 UTC (rev 20673)
@@ -1,12 +1,16 @@
-#source 1. 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.
+# 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.
+# 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
+# source type: elastic force or acoustic pressure = 1 or moment tensor = 2;
+# for a plane wave including converted and reflected waves at the free surface, P wave = 1, S wave = 2, Rayleigh wave = 3
+# for a plane wave without converted nor reflected waves at the free surface, i.e. with the incident wave only, P wave = 4, S wave = 5
+source_type = 1
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 10.0 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.0 # angle of the source (for a force only)
+anglesource = 0.0 # angle of the source (for a force only); for a plane wave, this is the incidence angle; for moment tensor sources this is unused
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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/SOURCE_Abel_Balanche_bathy_source_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/SOURCE_Abel_Balanche_bathy_source_solid 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/SOURCE_Abel_Balanche_bathy_source_solid 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.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)
+anglesource = 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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/SOURCE.SqrCirc
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/SOURCE.SqrCirc 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/SOURCE.SqrCirc 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 6.d3 # 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)
+anglesource = 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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/SOURCE.SqrCirc
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/SOURCE.SqrCirc 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/SOURCE.SqrCirc 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 6.d3 # 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)
+anglesource = 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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/SOURCE
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/SOURCE 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/SOURCE 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 10.0 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.0 # angle of the source (for a force only)
+anglesource = 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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/SOURCE_M2_UPPA
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/SOURCE_M2_UPPA 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/SOURCE_M2_UPPA 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 10.0 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.0 # angle of the source (for a force only)
+anglesource = 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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_no_crack/SOURCE_Rayleigh_2D
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_no_crack/SOURCE_Rayleigh_2D 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_no_crack/SOURCE_Rayleigh_2D 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 1.0 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 0. # angle of the source (for a force only)
+anglesource = 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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_with_crack/SOURCE_Rayleigh_2D
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_with_crack/SOURCE_Rayleigh_2D 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_with_crack/SOURCE_Rayleigh_2D 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 1.0 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 0. # angle of the source (for a force only)
+anglesource = 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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_001
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_001 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_001 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_002
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_002 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_002 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_003
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_003 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_003 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_004
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_004 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_004 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_005
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_005 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_005 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_006
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_006 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_006 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_007
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_007 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_007 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_008
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_008 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_008 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_009
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_009 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_009 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_010
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_010 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_010 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_011
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_011 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_011 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_012
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_012 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_012 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_013
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_013 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_013 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_014
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_014 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_014 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_015
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_015 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_015 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_016
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_016 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_016 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_017
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_017 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_017 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_018
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_018 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_018 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_019
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_019 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_019 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_020
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_020 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_020 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_021
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_021 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_021 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_022
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_022 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_022 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_023
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_023 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_023 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_024
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_024 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_024 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_025
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_025 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/SOURCE_025 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/SOURCE_001
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/SOURCE_001 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/SOURCE_001 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.400e-02 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.000e+00 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.00 # angle of the source (for a force only)
+anglesource = 0.00 # angle of the source (for a force only)
Mxx = 1.000000 # Mxx component (for a moment tensor source only)
Mzz = -1.000000 # Mzz component (for a moment tensor source only)
Mxz = 0.000000 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/SOURCE_Tromp2005
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/SOURCE_Tromp2005 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/SOURCE_Tromp2005 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 0.42 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 270.0 # angle of the source (for a force only)
+anglesource = 270.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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/SOURCE_Tromp2005_s100
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/SOURCE_Tromp2005_s100 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/SOURCE_Tromp2005_s100 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 42. # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 270.0 # angle of the source (for a force only)
+anglesource = 270.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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/SOURCE_Tromp2005
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/SOURCE_Tromp2005 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/SOURCE_Tromp2005 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 0.42 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 270.0 # angle of the source (for a force only)
+anglesource = 270.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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/SOURCE_acoustic_poroelastic
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/SOURCE_acoustic_poroelastic 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/SOURCE_acoustic_poroelastic 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 1 # 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)
+anglesource = 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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/attenuation/SOURCE_attenuation_2D
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/attenuation/SOURCE_attenuation_2D 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/attenuation/SOURCE_attenuation_2D 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 18.0 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 0. # angle of the source (for a force only)
+anglesource = 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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/SOURCE_canyon
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/SOURCE_canyon 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/SOURCE_canyon 2012-09-03 02:09:56 UTC (rev 20673)
@@ -1,12 +1,16 @@
-# source parameters
+# 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.
+# source 1
source_surf = .false. # source inside the medium or at the surface
-xs = 8. # source location x in meters
-zs = 9. # source location z in meters
-source_type = 2 # 1 for plane P waves, 2 for plane SV waves, 3 for Rayleigh wave
+xs = 8. # source location x in meters
+zs = 9. # source location z in meters
+# source type: elastic force or acoustic pressure = 1 or moment tensor = 2;
+# for a plane wave including converted and reflected waves at the free surface, P wave = 1, S wave = 2, Rayleigh wave = 3
+# for a plane wave without converted nor reflected waves at the free surface, i.e. with the incident wave only, P wave = 4, S wave = 5
+source_type = 2
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
-f0 = 1.0 # dominant source frequency (Hz) if not Dirac or Heaviside
+f0 = 1.0 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 30. # angle of the source (for a force only)
+anglesource = 30. # angle of the source (for a force only); for a plane wave, this is the incidence angle; for moment tensor sources this is unused
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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/mesh/canyon_absorbing_surface_file
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/mesh/canyon_absorbing_surface_file 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/mesh/canyon_absorbing_surface_file 2012-09-03 02:09:56 UTC (rev 20673)
@@ -18,7 +18,7 @@
417 2 468 469 1
418 2 469 470 1
419 2 470 471 1
- 420 2 471 472 1 **********************
+ 420 2 471 472 1
420 2 472 493 2
440 2 493 514 2
460 2 514 535 2
@@ -79,7 +79,7 @@
1218 2 1309 1310 1
1219 2 1310 1311 1
1220 2 1311 452 1
- 2001 2 2132 2133 1 **********************
+ 2001 2 2132 2133 1
2001 2 2173 2132 4
2002 2 2133 2134 1
2003 2 2134 2135 1
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/fluid_solid_external_mesh/SOURCE_fluid_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/fluid_solid_external_mesh/SOURCE_fluid_solid 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/fluid_solid_external_mesh/SOURCE_fluid_solid 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
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
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.0 # angle of the source (for a force only)
+anglesource = 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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/SOURCE_fluid_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/SOURCE_fluid_solid 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/SOURCE_fluid_solid 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 10.0 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.0 # angle of the source (for a force only)
+anglesource = 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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/SOURCE_fluid_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/SOURCE_fluid_solid 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/SOURCE_fluid_solid 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 10.0 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 0.0 # angle of the source (for a force only)
+anglesource = 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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/SOURCE_Slave
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/SOURCE_Slave 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/SOURCE_Slave 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 0.4 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 23
+anglesource = 23
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)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/SOURCE_noise
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/SOURCE_noise 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/SOURCE_noise 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 3 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 10.0 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 0. # angle of the source (for a force only)
+anglesource = 0. # angle of the source (for a force only)
Mxx = 0.d0 # Mxx component (for a moment tensor source only)
Mzz = 0.d0 # Mzz component (for a moment tensor source only)
Mxz = 0.d0 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/SOURCE_noise
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/SOURCE_noise 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/SOURCE_noise 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 3 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 10.0 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 0. # angle of the source (for a force only)
+anglesource = 0. # angle of the source (for a force only)
Mxx = 0.d0 # Mxx component (for a moment tensor source only)
Mzz = 0.d0 # Mzz component (for a moment tensor source only)
Mxz = 0.d0 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/SOURCE_noise
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/SOURCE_noise 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/SOURCE_noise 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 3 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 10.0 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 0. # angle of the source (for a force only)
+anglesource = 0. # angle of the source (for a force only)
Mxx = 0.d0 # Mxx component (for a moment tensor source only)
Mzz = 0.d0 # Mzz component (for a moment tensor source only)
Mxz = 0.d0 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/SOURCE_noise
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/SOURCE_noise 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/SOURCE_noise 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 3 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 10.0 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 0. # angle of the source (for a force only)
+anglesource = 0. # angle of the source (for a force only)
Mxx = 0.d0 # Mxx component (for a moment tensor source only)
Mzz = 0.d0 # Mzz component (for a moment tensor source only)
Mxz = 0.d0 # Mxz component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/semi_infinite_homo/SOURCE_elastic_2D
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/semi_infinite_homo/SOURCE_elastic_2D 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/semi_infinite_homo/SOURCE_elastic_2D 2012-09-03 02:09:56 UTC (rev 20673)
@@ -6,7 +6,7 @@
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 8.0 # dominant source frequency (Hz) if not Dirac or Heaviside
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
-angleforce = 0. # angle of the source (for a force only)
+anglesource = 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)
Modified: seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.pdf
===================================================================
--- seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.pdf 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.pdf 2012-09-03 02:09:56 UTC (rev 20673)
@@ -4788,18 +4788,14 @@
/ProcSet [ /PDF /Text ]
>> endobj
225 0 obj <<
-/Length 791
+/Length 792
/Filter /FlateDecode
>>
stream
-xÚíXËr0Ýç+Xª7Ò2Móh;fjO6mØ(1\yôë#@8¶K¦3íL¯$t5>çÞs~?=xwB°HózÓ+/Fà ÂyÓÔûæie¤2Uðcú©Ý+"Ák¶/d0ïv ÿ£2e ¹¯Ózn2a8PØ3í^ìµp!!<1õB#A©µq)ð2´ñÂÅ~Ô^öV"Ä¿
+xÚíXËr0Ýç+Xª7Ò2Móh;fjO6mØ(1\yôë#@8¶K¦3íL¯$t5>çs/~?=xwB°HózÓ+/Fà ÂyÓÔûæie¤2Uðcú©Ý+"Ák¶/d0ïv ÿ£2e ¹¯Ózn2a8PØ3í^ìµp!!<1õB#A©µq)ð2´ñÂÅ~Ô^öV"Ä¿
-xIËØæ9.&õr©Kc©àxOÅ b«ôBH"a«'Dè
-¹<Ædê:QÃIJ#Ó¡"d£ð2`ÀϪ:ɳ߫;ÌBvªº¶»õ<ÉóÇ]¥tõÕö9]óþNíqø"HlÇ2Cç²Zt¿öT*Y&OÉ$Í^1t¦Dý[°îƺrèN.NÏÑfö%ø³|CY5
õiѳܿºîØvïÕ2Ô¿3²TIÞ=)Úüà¾!²`B0ÔTpûoF©×£×_ÞûMõx>!ÈëL×L8&y
-ÌÅF®öÈ+%ªbËjëým -ÒþM³{ú $óÔ¶ÿ¾xT×à£áÄiø×Z©?3u¢ó2Yvë»j)çW²hX yßþ¢u.ÉnI/kÿáäÒU»trÖÓDA!Y(içp¥UVÔyÛl
-pÄ_^AС #µ:QY¥mÁ,p¬ÐÊ÷6u$ðlx<
-xû¡eT&·
-Ob¶þ/'»¯¤qçzBÝú'nÌTf²¾©ÎµJ3§ó×ÜZGûyêzÁaúSgÊá;YSÒÖÓSë×qÜ{z(M½5þo2Jww =3IæÀIÃUߤª,SÍkúÜ,G¹åò7Q3b,̱ù%Ëû¦:ßþXi]ò~ gbÀ°ðÑÊG>l÷§ALüzf½ÖBëöC} KÐ` Ä"¸ÅBôßývëÁñôàH S
+xIËØæ9.&õr©Kc©àxOÅ bKz!$°ê ¥:"'ÀSiL¦®5L¤42!-B´áeÀUug¿Ww
ì&U=+um/wëyç»Jéë«ísº.çýÚãðEØ2e;Îeµè~í©T²L²I>!¼bèLú·aÝuåÐ\£Íí%ø³|ÃY5
õiѳܿºîØvïÕ2Ô¿3²TIÞ=)Úüà¾!²`B0TTpûo[R¯;Fï¿¿1÷îñ|B×®pLò\í-VNT%Å2ÕÖûÛ Z¤üf÷ô; Hæ©-ÿ1|ñ¨®Á9ÆÃóð¯µRfêDç7
+d²ìÖwÕRίdÑ.±@>ê}ûÖ¹$»-½¬]þ_K§wé䬧B³2QÒÎá.J«¬¨ó¶Ùà¿<AÐ!A^¨¬ÒV0Ëk4ò}:x6<¼ýÐ2O*Í['1Ûÿ]WÒ8ÄóD=áîÎý7f*3Y_TçZ¥óÆùk.|l- ®¦?u¦¾5'm{zjûu÷==ÃM½müÞdî® zfÌ«¾IUY¦7«®és³´=Ì-¿ Ítl±l2Çæ,ïBêúöGE¤ueÊûÃÂG«¶l=òaÓ)qÄÄïþ«g¶×ZhÝ~¨d ÌXD·Â¢Xþ»ßn=8< KH T
endstream
endobj
224 0 obj <<
@@ -14008,34 +14004,34 @@
/ProcSet [ /PDF /Text /ImageC ]
>> endobj
651 0 obj <<
-/Length 3145
+/Length 3358
/Filter /FlateDecode
>>
stream
-xÚ¥ÙÛ6òÝ_1奶FÞß{x7vRy7UË¡ ePÒùë·ÀCÐx6ûD°qt£ïnòÕ³7YqUúeeW«ÍU]åqqµZ_ýî½úé寫Ûe^â/ixIà}|ÿþíûqbõÓ-î~ùy÷/³å·»_o_½¹}½^|ZýãæM%iáq~ ®0k®IúiY^-£Æ8ùG¥Î>ÒÅ2*Ö
-ïp°#(ð¾.ÂÔÓK~}÷3Ƚàí=k¶øZñfÍ:Æïl¨é6]k@Ïn§Oú#âo¯ÙÚ²Ö'Úg2`Ê¿E`÷ÂîM^¯i«iiô5Bݸð®~¦x]ºYD©ÄUè#BØÞ´äZ¾åÞ«:Þ,oí=®á->ïø~OÛ
-×øNÿì¹5¼²Ft¬¥PØ[í<7õ !÷úYêF±ÁgÇöôIrÐDòzÁôjD.+=`ÜÕü ¸àã©(½ªo[à R¢Ba
,©xû¡(9¹Å$·êQa`æ¯Í±-0â´ÈROq#÷ÞHºxk6HRrZݵðk*:Cd%ã¸óXá¤Yâ@éÍ$~P¸,zmoDï%í",¼d|ßÃMûÒhØZÎ
-(6è-úÈ
-Å|À!ªeÎP/5X4º4hy¡zÓÐ5kiÕñö> ËG> J3?/Sã>®Þþ|wãðqê]æã5£2ò³Lß3ßI
-ïÄ{9ÈsñFY¶Â̶8Ûr³ûª¤a5Ày³QôÛ-[`:B,ÀÕädÐ=a5: m2±ç¹ÃÐ[áäòõ©!ûå%×Á*ásÔ¸ª*ÌmOìÈj-ï
OÉ»{Ì~v2^K×pNkÑaM_/²` OPT+
-UN$Hð1ÑA¹F<LY~òý¡ïPÙ+¬#
-©O«/¢áu/Mß6VªÊãà=]3b!z²72&bB`K7`oè/`k. îíö¤ýûÃÍG>÷p3¸ã ߤçUϾÅiá5Õ$Uò)5Úana
-ÂØÑ®ë/nnÇ£/¬¢KÉx¯íÐ íMÍî[Ò èmº#¸Ûߢׯß~pgTÆ~bNªïXú1àÃj])ÞI}@Y)_!jJWso5/maàý@Ç*Mõ~Rû8ø=,Ëäö
-_ùT©V;êâ$EdȾûåãW·.ßSøq·¢Ú¸o^
-ÊmÜ]çÄ~\æ×/W/n.òã¼4«ÇV¶ã}Ýk
-@# ;çu±KEKeâyäôBãG¸aøa03G!
-8_ÑÏ¢o7ÓÄÏëï¤ldì0¡HßÉVC¹7Yìg¹U¿k{ê»òÅÐ/rËBeݲTe[m[Æg Úõº{=ge&'%V¯Vw,ªôréÏVaVæyäqKLÐá¥!«;±oCÔÔÓA ¹)1+Lé|ú½Kn
×6ö c¿ñéQLQÎO¾=LõiO;"\NæÛÌÑûGÅEEíNêRÔÔÏÃp¢¨Èaã|¾J[!uolòåÚÁsµ-²^¥··VE[È*c¢À
-Þ, ïrľ¥×&Ôí$.,¥QÙÛFÈ+ÏÈÛsÛèi#¥&|8bÈâÒtæ9dmðyÓ7ï%¹õ!âÕs4¦H
Ï$Ôhò0!0$eü¿Æb¯NöDA ÌG¤Wú tlß?0¨À]@tòïÚåÛr¿HÇÄ ñ^m§©a¦únP<2LÙäÍ:6¤o.yéC
-cN¿èP
-a.˧úGÒÁHãÂ"«ôÀªò=õ×ÆkÖʪ»ãT%-}â©?QÌÑÙà7'z£,MëSÃ(
JbæQ6Á÷<»<3Ê'GRSòT k¾g
-iôÂM^=ÑI0ãHºl)ãJf5ÝÒF7¹íûÚVõ èVÆp AÈ9d>È<òi¸2ü2w8ZÍ=OËOÆDkh(¤.ªLZã+Ø©®ÙS4Á¨8gM5µ
-MÁ²â$³þË)¬tÜ`Òþ\À?Á×ÛîÐ
-â°EnëP¤äú¸cÕElþù õ½@*ài£¯B¹3î <¤ý2ÑÑbðmó©¯z;:çeKºÙ3k11A)KG\U¨ÉÔ/Á!?9kÞô7rmÝ·Ê6åàl¥ÂãtÀ4~ßÖ|èÞÕ9ÄO¡{xãÔ(öË!æû(õ·Yn5Ô!Cb6o·_ø÷¢k(g:+ð.ÙÜþ<~,¶ÓÀÈ~Ä'¶mL>rdÝn:¹#æ(Ð<&{Cú]/ûØs´@ÚÖôoiFI|t®[gêóÏEÞirïDµuVr²
â[×hs£OóÎÁfù½a[Ós´oì2ûìÄäÞRKG¶6\ãG)³¡Uz.A'×~DÕ÷ÑLßL[X§s:NÛ7g÷«u¯d$o
-ÄúC$ûë¶8¤o´å3UpN¢²c`\mazÁöµrv¨þY$ñÍêSëm^®-8øY»ðÀE·<´MaëV²é,×,Ò¾¹6ETE¬½±.võñÝÔÅP"f=®¹Å£ÛoíNcÛÀE]ZüÕHÖ2òíìlßIJ á̽¥
-mÉåfÀ)ª½ÎøúËõ¾Îõ%ËÛvºV
-OC±;yߤÒjÏ'Â,ëætâRÈ~À>¸ESÏdµ¦YôO¹ÈéÆT@º²<Ó´7iYsÓ0mܶ3çf|þµà+Øý~´Ù8½)T{ynÁyËøÜP¹
Cø(êJì¨Q/ÉÌØm99Xã¡oÁ´Mk]É¿fRÓu±öïÌ'¨¾$*q®¤úB¶Æ÷ÞZܽzóÊ¥¦±6Ç}ÒºnjË¿Þ¨cÔI4í*ÅXÉ`kà eûä\Î|¡èËSp«pVÊ(a+©îe¦³p®5çÄ1^þXVP|]"WVAêf[¬q¦w~eÿ^J(ÓI¯nù7hE¨?N.ñ#nBå& úßmÀ·x´ª¸ÌÃ&m(Gbm=gÖ¨tæÊÄm:ce&£uX¹)OüdXâèÿ ¬ª*³XYø(K/ÎOü<ËgLÓ,8`Ö¯È!7PýtLòÂhÕwU&9èa¬?Û%ÂÓ)©kê&4^¹¿B%æÃS2äAߥç!k];ö8OVñª¡g$*-2
y ®¿39ÚÌýéÜ"ñc['Á°æ
-EÊ) ªYª&å¤úWµ[²Õ¦ \øÞ,u
-Ð#:"Ô{lªÏrêìG:.ÓT%úÜÎuÇØ¥`
_ÿå`|ã¬EÊÍ>2jñê½M}ReRpBé»ëwß¾]¿{øæ 1ËÆÞ¾m}ÛR}¬ü ?WËwI¢nhÈ_Y$ÿô£Ùe»Àe69mðMª©£½1±0ý<z®c^Þ{Þi¸V¯2¾$¡³ºªoùÉ`@Ø"}t ÈSúÖ³ÿÖLÉ_}í¦ óðÇÒ4QÍǯ]ßÐi?[;_áÃW¯³ãç.\oÎ{Eû÷Å©¡¯Ô£g·«gÿl
+xÚ¥ZYÛ6~÷¯òËRU# ÞÞÚÇ'ÎæpyäÝT%)/F$)R!ÈÈò¯ßn4ÀKñìæI ÎF_7ôÍêÙÍ4¿*ü"eéÕj{±«,ʯV«_½Wß½|·º}¿X²$ðb±LÒÀ{¿ïÃO?½ýé[X}wK»XAàýË.ùåîÝí«7·?²×ßWßß¼ãÑAqûA]tVã«e~RWKÁÔc,ÓaïX,# bSÓoU·ÔØó¼?aâ)ï~üü^Õͽ¬vô¹®«le])úc)¢jÅÆtTôÛîÍN¿AÔÉͪK¹±lÌàY´ÀòïEØÃ;Õ¯MºÜÆÐR
+J¡¾qîÇ)\7ý"Ièºb»` ÊÖè-BXÄrïÈþYò¿2ïUÉÏm]-o«ÝT-Í©ú½«ÑlE¹¡oñGÇíáSVª|B#oµ¶_ÑæuU©çÞü|ýɱ¥ßVÄ54ãä]1àõéÕ8N+<`ܦõQsºOûzLEỦ)Å *&ôJ<·k©$ëIL¸ÔìH
+¶]i·mçEx÷骻 I=á°¾kî)YÅì ^2;Í7Ztk&ñkê
Ëreæv¶ÁÍZÞ,ÂÜÛ!ëðû 7¥ÓV£ÀÖX27 ØßÓa¶Ü
+º©=8Dmµnä±µÔ£cÃp©mxEMÝ`'MïF6bÝÖÍ|@| KR?+ë>¬ÞþpwãðQâçi?ͧk²ùijîÎï)QGÜ;×62Í\m|U²#`¢
+вÌ.Ëç~ÖhXU«hÜ.TÝn'T+!2&ýVz 2Y=âÀeI]mM]k¹[*!Õ¡¦þ0½
+îD½9Wü°|ÈuȵòkÒ¸´ªjÌmÏì$K#ïE~wûm³éZØÒ\£1£EÇ
+×~?èÈF%<AQ{ ¬YÅD~&:sÁÃPÏÏúpìZÛ©íæñçV®µ/EÑÔe&J_Û^ªÚãÐØHn·á-µNf°³2æjB`#¶`oä/`¬Ü=LÚxó ÖþàfpÇ3}¡çÕ¿]C
+ÞÀ%Ja!Uøí0·0óÜÚѾm/nnN§¯r-Èx¯i+Ð ÝM)ïÞ!èmÛ¸Û_Øë×oß;ÌÛ]Qõ}Ë?óè¢ÆDWw¨$+í+UMt5÷VadõØÞ7b¬Ò¬ç¥O_â7^!ö *Õj/\Ü 3KöÝÏÞ¿ºuùÜl¸]ÀJë¾ëõ ÜÖݵγ²ÐÐnòúåê¥ÓÍ1?Ê
+;kpdeûº+Ía÷ÐH'äÎ벤إ£¥6qÀg?¿0gG£³aøa03G
_ªk¶ØOÞ5ß¡l0F¶(bíw#¾3½µ7iä§Y¯>~ÛtÂwáÅÐϳ
Úºcm©Ú¶à ²mÏ@tßiæݱ^f8È{bÍl}GdÑÚL×Hg¼·³2æÕKORa!+ÌØ·%jêi ®A*i
+£&:îàÛDá}VG½ñéÑLÑο>Oõé Z.¡
+'
+óåÙþò?Ìòµ=
KQ?â2Àt¤¨,M!<ð²ôx¹4s\j[~öj:NÌýEv)zßT!RÓ·Ö`CÍ*ÖuXÂÌ86»¡6K·¤; º'=ÂtñÙYTÊRNS@»ñïQûaM½Ï,YF¨Þ &¹ 1¦Â^@Ôe·1~*BÄY¦ÆÅEq0ªèUÀ;Úù®Ê,hé×QämaZ#S×h;öÞ=Fà?è'¼¦ß»'Ìefß÷ü\
+¹Û»RË
+ÂáayÑßaÀ~Z0ØÏ©ßp{ÜëO²
`ÑÒ¸Ie7ðÒµ§¢Üì§y
+â%¶f¼bOúcsÓAÄàÖªvs®NJ´pw¥³, ^¬W2~>:7ÑXÜÌ0©6Ú®©q Êd»7`
+sBíÄ\OúÀk{²[éú¥f"Y<Lv¹#?ã ÀxÇ[Yº|aùY_ ¶¹Ó̲têi1éÿ¸íª5:îmùqÜ;ÛWûº¦(©u8CãK¸õÛñÅ DØ 6öÏøéôy,tõè¯\øC$tÀzÓ8}¿ëO(xε´d~çùV<g+{C5âÇF»mD#)õrÈË\Gf
Ä=¾v$7ðëÜEË;¥$¯\§`ù-è)v8ÍSt¬~-¾¦fíR÷-ëIK¸ëw,Oh¢7å½q*pÈ?MgPa|
+¢ÍR«|YJ%Fб©²â¾ð28ªaãD:Ö
+óq Fµ¬Úe[q®<te_®NQÉ6äÔ3cÎp1
2OY6u_6¥NAë±çÉsôv10½èúÇ>ÁNM1.°Õ°08tÎY:Æ,+Ó8
+Ǩij:ÓçÞFÝeæ3HÞ;k
+§½\ï¡cÐl tëj£¨CWBà×
+F#\»3ÝNÍÇDGóÁ·s08¦YNkÙðv¶ÂN)3´¤¹HË]¸bbz0ñøéÙðb(j )y×Dw,$¯E}Ç0ߦ÷õ}Ëeeó{WE0b?ÉW§°È/,ïÐßng9ÓP_®yպ¯xS[`)øÒYáæ!ÂåÏ£çT ¤"à#úä®â¥%³ÝO÷Ün'±æ;©ËiÿQQ[ékàB»RüÍdÌhÜ20õ¨öüss:Ü[µÞϪ}+lêÃÑÔ^ ä²ßçÁ=ªiF÷Fîì[3®e"°Ùoe§&÷F-]äÆrÞºJÅ)×ät÷Ä5>¢ê1{6Ó7ûÜ£ów]ç
+§éªûäí°H«8ú¶¨lÌÒòE4õLF}^ AÄÒ5£\`Wj_Gê2<oV×8W0rU]ÆÊõüûìàX«vyljH U_âÛ¶çZhW]ÛâÈ÷öjÝ¡e¶o§.Fp5«]Ï-¼~#E{kÍÜÔEeÄÿÿÉßa0Á¶¯ö,¨jiÞJ4üá"ß ÆNu~nÄr³®ÇKGýå,»m½Jñ(ÏÃã@ w²4ê̳.I«±>+;;ÚݹK!»áôÁ-ÚÛìz¡düSîrº©[Ò®jyäö1Î^xfdýAÍMCÁ°uÛ:Î\ñå+à`gâ#øiëÃô¶pâ"ò³¬wO9ô0ç8;Iý0íÝÁ£G¯ÕÞqôÅð³x¬¹5»L[\÷àßæÕ$¥±êÒÛ#E¾ ßHi¬såëO|g}¿êzk3=w/Qo^¹4ãÈÂcàûGYV¥¿«ÿ¼QÆ8+H#b6͵b¬0ô)ðä%$ÙäÀÜãÀ? òåÉÄ)¸U8Ëd´KéKêWÔpêÝ.R«hÐGÓåa¹çPØC r¡² [òh"ü¥A/Ñ$dÀë[ üÍ1hßÇnBcPý¯Ö-Ç¿;ʨâÃÆ~IC9õ!¶O§Àâú 5ðý
+S{¹me<Rah%!7å±Slŧ=þg Hªìd³ ðJ§/.w`~f³
+¦06ÃúÕ90`ýNF /£|ÏÕYhìñ£iFæ9>öjh¼,©ÓÔ¨ eÎÅôºÛåxÀA¾Ñóð5ä;Ö8wÖuñ¬¡dµÙ¼<Ì ®¿·mæþ¶"zDgQL:6µ®ÈCÆ-CcÕPAbþ"¤' ¢l½h$×yU]-
+
Þ9"è·Ðú¹]ïMöòdÓÔ)Ïðä~
+ÓQ4%ÓÌéUeÖæ#Qâ^NPúÌ«£^[gÖêØÏ´/ùe)=AÙP&¹·6T¶GiQE1r30¯°ÁÅô`Cêø,I§NÐ>âÌÃøÔlwWA¸Þh<»]=û/óiá
endstream
endobj
650 0 obj <<
@@ -14069,42 +14065,42 @@
642 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [95.91 336.274 189.628 347.178]
+/Rect [95.91 300.409 189.628 311.312]
/Subtype /Link
/A << /S /GoTo /D (cite.KoTr02a) >>
>> endobj
643 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [193.444 336.274 215.361 347.178]
+/Rect [193.444 300.409 215.361 311.312]
/Subtype /Link
/A << /S /GoTo /D (cite.KoTr02a) >>
>> endobj
644 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[1 0 0]
-/Rect [290.33 336.274 304.776 347.178]
+/Rect [290.33 300.409 304.776 311.312]
/Subtype /Link
/A << /S /GoTo /D (figure.4.2) >>
>> endobj
645 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [188.094 300.409 281.811 311.312]
+/Rect [188.094 264.543 281.811 275.447]
/Subtype /Link
/A << /S /GoTo /D (cite.KoTr02a) >>
>> endobj
646 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [285.627 300.409 307.545 311.312]
+/Rect [285.627 264.543 307.545 275.447]
/Subtype /Link
/A << /S /GoTo /D (cite.KoTr02a) >>
>> endobj
647 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 1]
-/Rect [99.228 240.633 196.862 251.537]
+/Rect [99.228 204.767 196.862 215.671]
/Subtype/Link/A<</Type/Action/S/URI/URI(www.llnl.gov/sac)>>
>> endobj
652 0 obj <<
@@ -14114,39 +14110,38 @@
/Font << /F68 529 0 R /F44 199 0 R /F67 351 0 R /F43 198 0 R /F8 653 0 R /F11 357 0 R >>
/ProcSet [ /PDF /Text ]
>> endobj
-666 0 obj <<
-/Length 1968
+660 0 obj <<
+/Length 1750
/Filter /FlateDecode
>>
stream
-xÚXÝsã4ï_7hT[þ>£´ô»R.aê(ÛÊÙr{å¯gW+9vj>&VÒîjû¥|»9»¸N²EÎò'ÍnòEfÍvñwyóönsõq¹â±ïEl¹ßû¸|ïçÛÛw·ßÓÂææëß/ß÷~qG~]ß]]^_}àß-ÿØüpqE#AQ1?L>É
-bÜsvµ9ä/EçGñ"æ9KhQÔgÎ|e!çfËxlÝIK¸xWçïÔÙOðsK+Çu5bû-Za¢\ÀcEé"öSä ©x]>ötÖào«0ð½KUD[vª¡µ£¯ÞÛÝ^ätMØ-m)ÇÊ.f{rò{Ñwl¡Ù®o
-]:1¢Ð½¨ªõÜ2käptÏX- at QF9]c åÍ* FÒ
-;â8U°Vµl´%ʦS-°j`©³Ê:³9·_ÖsÚæå¥iJ:ìà¦À÷µª!à#«ª¨Uù»ïóB´ôó`:üÎÝ*-ÏÉ0¤k
Aâý%[E#]ÖÐÛÔu_YY8/TÛÊ®½íìNur¢ ÈvÌ.=.:_¸p ¥¸§Ú.*·öü¬xY®°<éR´zÿ©¢I$Är$ÞÆ0õ¹×iXÆaà8Ñp¯v¯ó²#X`:yÀÂÁ4ídFAkÖ¿±&Á°èè©?Ç$²!w̨'
-RÞ¼fÀYÉÏh²8ÍFJ´M÷änÏ^T`°à¶oÉJ¯¤f>ýØøje ;¢Á«µ?ãøAãj±? 4rx3 ",¡\[G°Àb^£ùáÔwl¼g£xÏ=¸¢¥Sn$!õiÄl_ZTª³ÛQ$Q!åy/?,GÚ1È=:Þ×@HBïF%÷½§²+·'Ztt¶wâÍÇA2g>X.8G~Á?
-AÊóä©;ÂW=hQ6óv{/D
-.v>XyïìñÕÓàÐÊä (¡ß~z\)eIì3
-õÞI²j&@Ç;©Ñr((ÑZbÆyíennöóúêãýfÎ×¢¥y0©ùcNÂcðÔÚÐVÍ ÑJÛÏ×$FaÆM9Ìrø"~¹ïTWêrÉcëXçäÇjVÖ`ﮪ7ù
-çfo3ÎKûÕ{¡iTã0 é=EõÕÖuUôH%\YÂØùkÈL[
-÷1½÷ÄÛVÿCÊ2>éXéé8'ëÓ,ù°¦NôÇD4á:ØNo-m¶¼Ö8É\§ÕVÝWe§ïËÝoU£ï[Y«'Éôg=µòúÖw®W°½26I6
-á|¶µü¾ìÎizhÕC%k7uÜê'å LB7Ç/
-½¸µêáB¶[NÓý§©ÆØlí¨ëÕb3Q¨¾Óeq!+_ôß,<a«/qÏÎ^Ô)ÏìµJ;¬æÛÑ6¸ÂC[§¹cFmXg°îñ~Ô½BÉ£çeÄÞù7m£ì@hèûq«æËèª
-tHÕË%·P!cÛ1F,X®à!â{7
-µy¶}ÛÞ¼[zÏúÆEAìײ~hE#±JõYØôcê¦1Áyç'VkÌ[NÒìå3SMÁªIxu:÷^ØnKmÍÑÛßÌ4Ö>µ¼Ü¾ ÌÅsîÍ= àUCo0>±85É#UK¢*Ybö ¢õàEîE0ÀþyõÝ]õE ÂÍä$À{!:¡Ã}÷D
-Ø7ôaºÅRðúÐ¥yzúʯjº-¨w¢½ß.Cþ¯}R ó,¶{ofqX½Æ¢pÀ4@ã´ðqZ Îaä)6YïÑ&ä°¯´²<00ðX¥ÞOTzY
-/@"ßZ¡DÜoÔb&ö ¨¦5ÀhAvÒ:¤w=³P¯~é4 6ÔW%v¤Ðô¨ O5æ]xUJìyMMH\JÆÎÓK>A·m®¥Ô'5Þm#+«u~0xtÞã"²´JW P
-M_ÌæìoÜ÷Ó?&ù>9©¸[ ]keÔŦóÛ[Û0Ëqô0*øO`yÝYMø®¯h|÷¢÷ôpL½®hË&:\+îÌKû¡uõÕûN<A«Q¶ìð2Õp24X+AÝç¶ÄRC²°aFùwxá4»l¡÷(f»VÕ4z/týªl%Leã$@V~ê-Ü:89qÐío©HÃvôû0ô.o?®qÈsïh;YÌRä]
ÛgÜ6AÚD}ੳe(2ôÞiZÛQd·* L±Ã
-½Vàeaÿ+R7-*Ä%º$"ÝÛÁò%A+»¾ÒÝ\?qòçÖã¨YoNË%söý
-tÞ^C
+xÚÛnÛ6ô=_aìIî"¤îú°¶IÓaɲÆ-¬ÁHr,T]]8_¿Ã«$[[;`ðÈÃÃs¿Ñ¯W'ga¼HPÒp±Ú,"º¼x±Ê:o.¾YXº4ÀnbçÃÒÇÎÇëë÷×ïÔÁêò\-nûuI0v>+ÜÞ¿¹8¿¢o~9»ðý#?ö¢V¼H pNÎW'@xA$IõE@V'_O0òcR2^ËCsSÎÞWÉâ-?ù~æÈ5TÝÙ×Â
+á
+PìG G(LB%âEqß7¹±}¹t=7¼Ú±¦hyNøF}»Fm·lOÏFi
+VßúÕÙÁÍw¬o[@Q»M_§]aØ°´ëYYîÕ®oói#{#=ÜÀz¸E´z\==^=?^==Ïø%ã+^åuLÈéòºåZ§ 5¯á¨UûÏ8ÀK8¬,2âµVs«²Ô µ÷M·§`QèݺÛK³Áñô O ý8×¼ÓðnËh`2PeáúÂ7Zß¾.¤ÂãØÆE]d%wÄ©RFj"+6¼±w
+í®!¦¾>8÷Þ¾ñ"aªõ9¸²*òûÈhÿOêZŪhfû:Vs!§%pILD¦ÿ]dú¿|ýbVÜ¢$&So »î±¨ÔXiQYµ+ÏÓ
+¹èØûyêʾaì§ôòHè<ç
+W«®¨4L£¤¶¨úR3û7MÞÞY«1ùÁ,7crÑph
+Щ4R(êðfJÃEfËû³
+mñ4órÖtÛ¯=û"lü
+Cg%bê´%1()·³Xca_´Ê/ëuBgC¨S¶T§ç¯4
IÅÂ1!8ÏüØ6¬1|%Ëåå1¢Ðv<:#E( á)¸4ÅIÎ`°ºnÖ7ÊJG\c<?Î$áÛ°îðLä_,Ñl=4x¹'ZW)_rååJv¬h¦*EüùØÆNøxð*jè¸&>,IàäÍ^AÓ·] @*Ëäq×j%ã°¨ïå;ÞO =ç2gK¢-²)Zu·7ìÚ-#ÀÐE#;ôÈQÒci]*wRGøò»u.»làÜím°H0?vÞëë{ޫŮɡ'¹MvØh8I&.ªØnÌ[ÃI'¨Ê3lóN-Ñ
B²b
+«rQq£ÌÇf
+³·çÖ«¹Xó}
+#HQ<Æ"4O
+.éô´o8a at FiFå&ðþK°³ãmÑ¢§«À:Uq\óÚÕ¸zìée}â¹ÚjbãÊ%ö
þ±
ùeÀ §G)¼/3}ÕXDéF"%¬M¼&Äz89TÚÒê1c*÷ÄR!B1µfÆKu*ë+Æ( æÙÀú¡M(H4Þm®;´öÀ±ÄbjMlëxÆ×eÑvë]åÖ¯»uWü!GÝS7WµECY·ãî±éØ÷¥æºßïïÅd*¶»ßye¶ZU@óTõ;YæÍEÅ»"òìÅ{PH·6hÚoPuÒ@^µýnÇ1KyßvEzL|EüÆÞ T>>¢Þxjd=ÁÕBHܹ`°8¡Jôx>GLÍu`kÝA?5FûÈý2ӱߴok®¬±ÏjGÏ ë
+º
+LHåȾhT»ÂóPH}%¿ÈÒ
×/v.¹æQæêÛôå{ûi=·&\åÕ]Ãê\t }dºüȾ)Mð¡QH#õfuËpUÎ³ê °Øª.'µNkÞ:YVÖ²ÐÑþ¯©xB¹ <Åbó|µðÝØ¡ùdŪC"ee:T]³;v¿ÝÙ(2Oëö'÷Yc`}mèp35ð`ðí$´[·j {¥>¨kD+8VÆ°(9ó]´
m§¼aÍzS
+ù]5ú A&qpð¼õAA|ìÈ~åQÔG-|g s¾ðÔ7±ïNrÀ+4/åX(÷ÀBºGó(ºíD¤½k
+ k:|+.8
+§BøÍ;Õ 08Uÿ ¦ÆÓÈåjñ¿Õß|ܬ
endstream
endobj
-665 0 obj <<
+659 0 obj <<
/Type /Page
-/Contents 666 0 R
-/Resources 664 0 R
+/Contents 660 0 R
+/Resources 658 0 R
/MediaBox [0 0 612 792]
/Parent 616 0 R
-/Annots [ 656 0 R 657 0 R ]
>> endobj
648 0 obj <<
/Type /XObject
@@ -14154,15 +14149,15 @@
/FormType 1
/PTEX.FileName (./figures/gauss_vs_triangle_mod.pdf)
/PTEX.PageNumber 1
-/PTEX.InfoDict 671 0 R
+/PTEX.InfoDict 665 0 R
/BBox [0 0 447 394]
/Resources <<
/ProcSet [ /PDF /ImageB /Text ]
/ExtGState <<
-/R7 672 0 R
->>/Font << /R8 673 0 R>>
+/R7 666 0 R
+>>/Font << /R8 667 0 R>>
>>
-/Length 674 0 R
+/Length 668 0 R
/Filter /FlateDecode
>>
stream
@@ -14197,7 +14192,7 @@
|÷¿£YÕZ*§Ho©°jä.(QPe;ýr%º*´¾ÖøÍ;ÐO3°èk0Á²êhð-¥ðg|££B$*xlñÅgz¯ñ#Ù¸L`:eß@@hYwWG½è³ÌÚ *îÅÃvLáå7X$¯¦F%mòû¬-}3â N?¡í®ÚÐoÔäýP-ß\&ûÙÜøøûÏMhfÚ<Ò]# ÉÐxíê(Ö<SIx»ãðløM±yzmáèáIÆ#cFóäéñÚíâds< kUl^L;eÙ:¨Ç¦L¶t àÕÈ;÷lòôu%óø¦LOXiÊâëéìðPõ_})Ë©
endstream
endobj
-671 0 obj
+665 0 obj
<<
/Producer (GPL Ghostscript 8.57)
/CreationDate (D:20110307214926Z)
@@ -14207,17 +14202,17 @@
/Author (Seismo Lab)
>>
endobj
-672 0 obj
+666 0 obj
<<
/Type /ExtGState
/OPM 1
/SA true
>>
endobj
-673 0 obj
+667 0 obj
<<
/BaseFont /HSZFNI#2BMyriad-Roman
-/FontDescriptor 675 0 R
+/FontDescriptor 669 0 R
/Type /Font
/FirstChar 32
/LastChar 121
@@ -14226,10 +14221,10 @@
/Subtype /Type1
>>
endobj
-674 0 obj
+668 0 obj
3794
endobj
-675 0 obj
+669 0 obj
<<
/Type /FontDescriptor
/FontName /HSZFNI#2BMyriad-Roman
@@ -14242,10 +14237,10 @@
/StemV 111
/MissingWidth 212
/CharSet (/C/M/T/a/c/d/e/f/h/hyphen/i/l/n/o/r/s/space/t/u/y)
-/FontFile3 676 0 R
+/FontFile3 670 0 R
>>
endobj
-676 0 obj
+670 0 obj
<<
/Filter /FlateDecode
/Subtype /Type1C
@@ -14266,84 +14261,79 @@
f~Ýåhþ®·<Üåy¯
ÇG³Êã% þ:`'
endstream
endobj
-656 0 obj <<
-/Type /Annot
-/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [468.671 116.668 511.843 127.572]
-/Subtype /Link
-/A << /S /GoTo /D (cite.TaLiTr07) >>
+661 0 obj <<
+/D [659 0 R /XYZ 72 738 null]
>> endobj
-657 0 obj <<
-/Type /Annot
-/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [515.761 116.668 537.679 127.572]
-/Subtype /Link
-/A << /S /GoTo /D (cite.TaLiTr07) >>
->> endobj
-667 0 obj <<
-/D [665 0 R /XYZ 72 738 null]
->> endobj
655 0 obj <<
-/D [665 0 R /XYZ 172.661 519.651 null]
+/D [659 0 R /XYZ 172.661 519.651 null]
>> endobj
-669 0 obj <<
-/D [665 0 R /XYZ 72 341.445 null]
+663 0 obj <<
+/D [659 0 R /XYZ 72 292.818 null]
>> endobj
-670 0 obj <<
-/D [665 0 R /XYZ 72 292.03 null]
+664 0 obj <<
+/D [659 0 R /XYZ 72 243.336 null]
>> endobj
50 0 obj <<
-/D [665 0 R /XYZ 72 228.855 null]
+/D [659 0 R /XYZ 72 180.095 null]
>> endobj
-664 0 obj <<
-/Font << /F68 529 0 R /F44 199 0 R /F43 198 0 R /F11 357 0 R /F8 653 0 R /F14 668 0 R /F67 351 0 R >>
+658 0 obj <<
+/Font << /F68 529 0 R /F44 199 0 R /F43 198 0 R /F11 357 0 R /F8 653 0 R /F14 662 0 R /F67 351 0 R >>
/XObject << /Im9 648 0 R >>
/ProcSet [ /PDF /Text ]
>> endobj
-684 0 obj <<
-/Length 2149
+683 0 obj <<
+/Length 2083
/Filter /FlateDecode
>>
stream
-xÚµYYwÛ¸~÷¯P§/RkÁÄÆ%çÌDZ5íزÛ-v¸(/¿¾@4¼N&y6ÞåÃ]á_;{G~8Pä²¼dÐp²L'¿N>ì.Ïfs½)C³9÷½éÙyÓããÅñ?ôÆòá|aÏ^ÚOþ}~zxptø¼}]þsï±#ÆCäÑ`âi^Ø7gè3DOÔ¡¹=5'|IõYÈl¼é
-8Nï4»¦Ò¿)ô .3Y5õÓjýàÁºß×tfO«<¯fOï²òV9ë$«JCT4xôàW
á×]=¹#$y,eĹÖ"Ý$"ÕË&Ñ´Ên
-í²jâXD~Ý8ØîP03%|}§ÔÜ9\î`ÐÎà æEO8À ºI±ómÇCG·GúãvÓ~ èVË{{÷ÕÎøoOÌ-ñyú/ʦ0Mq 4vµ\m¤9+G(4ÜRÔ-°
-blqµ@¯âñ¦ßÕÍ¿yz `ÊìñrñMU>M-¤ÃOl
¬EÞÞ\eëwÚü`hM8B+ýËlîfÅ:0ÓScÊõv©µ´408×4dv[ÄW÷÷úÌÏú'Á¸=Rä>|Ð"ßÿ[»ÁQÞrº¹¿J
-¿["¼=À1!ôüæq/Ý< µ´aèÓ¦b!øc»7§ØC!e.4GhÐýÑ ÷c4ø4
-ö<ü)Ì-,ûÿÐàOÑxÁQléaÓ󽻨ÇáÀUzâZ|?+ep¸6f1Ü%«¾ªøTApÑÓuϬß;Æ"0K9Sáò"RÄ ÔGßúÖN{²d%YóRúLº_ä÷×c($B5·7Í6;£â04ªÆõÕMWA Øsï´áã êss¹
ÇR$m·àõHaĽÞþj3JU¦~гêÆfTä¶^0ù¸r¥G©I\yÿR L7®k
-åÅÙñÕåâü`¹<<¾Ø_.N¯N]ZéTÙu0àwð~ñÖÖÔR¤Õõ)-e<²GÎ)(Ä;©h³´4+=úɤj<LÒ¸XkëILý`ë)Ø©µQ½P4ÛÈ%àÈðË2ä¬Lj#bÔG
ÃèÙYø·M+÷Ðæ®|*NdWí5ãr2zWÜÔâÛF*+\NzSò®$B!ÔH]TP±Í¥µ5VµRÑÄY.fõ3ì[1&½Q>«bmjJ(ïHDqÄüÖËxý1/k²¢+¡á¶ë\OÔÒéD^&L_t"¨!àá/>_|ÒγüÏé¡tD¨[KýÍóH)äè.µSìi¡{ñ±ø&ܵæS§c"®ÒN}M¶_Ç髬llüiVU:ºÁßqQÂÞ:ôIº¾¨kè»FYmÚ·\a ñqð9 $lí¯zJñÔÊÊ&ÈÍbB£'_u'93ÄÈmd8>?¹8;pY£(d¼ZʨÜ×¢ÖãÖ:¶ªÓèôT9VÃð`GŪ'ºÈAÉu?ÔÞF¤Ïzд¡1¤TH¦ÊZsW9W·Wö¬¾£Ð÷û¨ß²d¥é¶ÙÅk¡µÄ@AîZ¨pèj îï÷û{ojÌ B´°e[D çOXõO){WcälD)ÛXc
-tOËh¼ "{mø6Sfº¿azÝ®Nn¥¥@IaÜÑѹ: ñ!åøl ÍÁçåùɧ®\ö3D§p0V0O]袦1Ö:Õ²ëhTþfã}È &#?&Þ·QB$"Sáé»Qz'ÏÊ'ÁæY²ÍiRd²¨ná;;DൿRbqÄH4ÊûWG'gÿÚ?s½?¬o¼©®ÈºÜáS5Möq1ºøõ ©Ú|Ý{´Hí-sÏUô@ãù]Ez#Ðs®Þ
6VTõ`Aîú>ð}@»Pnâ\·±Ñ%óÁ~EÙ4Ý«Ï%.ë}D·rÛÄãAçGÙz¹m ãmoX®a2»Î{tô ò¹w£
CLÝÏvRv7ѽêzÞ3^µ2Ý¥)¡.8ÞàçËÇx½6Ö÷¥Ø¸¸¶
Xã´ÀC½¯ö/jFF©g,3 vÀÍÚjW8ÿqÊ¡-ÙÂ%oRvêÉndLg½Ê¦SzzðhfðÑ|üñÆþÖöÙ9¶¸±kpãJUÕ¬ê¥< AÿESõÔ(Õ»£j´(T$¦ëÚU£ÖèU CÌT_£féNn®Ëè.ÜB0»ß°ZÛL)ÖÊÕzµPóGQWý%®võ~¬GÙ[íæ]íZ=vUà£úUö\AECÔµÛãgÔfðí^¼©Ñbº~Çö V¶ò´ëÃö·w ¶» ùÞó¬õÒ^奨ëÊw°Ä¥^Jª²ÑÏÝJð- D¥x-b;¶ZµÒÏM]Á<ç£4ÜiÖÅÐ^nlâR üaæsÛðß\3?Û^$ÏRJi³Ô9þ`ê®ë
+xÚµXÛrã6}÷W(ÙPL/SE#ïxk,;¶<»©É&!^4¼Ø¿>EÒôx*@`wãt÷é^xÁ$D¡çxõfâ;u2ùd-ÞÏ/ÖËËéÌ¡¶å¢éz¶u9umëzµ:]ýKý±~¿T«óSlÛÖGóʯ.åónúyýïã×í(rilâOl¥{bÏdfg[ús½åíÕMÉÔD-b«ÈÓ½^Év)ËX^³Dï/J5`iTÕ<V,ªYÉ£T¿åzwÅò×|êPë×Zæïp(9Kµ
+q«Á-SïÇQ7iÔ*ýͦvÅ´ëi`[ÑNÏXQÔ»Ûö?˳ÛÆ(¤TZ0&¬¸T"WÏe·e35«r#Ìb½ð gSǶî§Z¬0j°BðDO×|:#¶g5Û4©_ìëPFlߪâïjµMá Ò¡ßq(vl)P;ôjyvSE÷ì&á%ÚíG" ç¹fÿæ>¼®Y®tÝîþqÈX¶(9¸2çj¶)L>D·EÕñ!¦^?cà2þ¥ÑËàjC
ë!ü0äÃv0É!ÄZ¬.¯ÄÐ ³¨¬ ä_®u.Í>°iWÂÞ{° Ä:Õ_]!¢ëAo-¤aPzCS\ ¼WKUǡѶê^ÙEe1j-WE4JV5i]ÅW±Q¡©Ç
+{©FÇXÅ3Ö¼ÈvÑC\Ïþy Ïq<9Óä¿m½/Dô=(!òdðØÒúr^u9ÅUÆ
-\.)÷DTàx~§Ä,¢2æ
I WÔä>á0>¿ÑÙÀ6(à`-Tm5-YÒÄ2aòq
øµC+xÏÒ[~§5à0tßô
+ÓSßoÅ1ë#§³'x)d¡C'.qPÞ$ξÙÐS¹¥;7 ßÆbùø4Ãöä]qôüÌ>ëHÿY0}f`zRÚ-É6cx>@)+*ÁËYiÐ^W FWô¶³]Yç®S×lϷؾ\ÕÐÕ1²d*D«-ß½¡%õ0qÊï¦3Ïî
+5Õ1)âQý«§ª\ÃÌ©ø]Ý<>ª=?©G±Ü@J¼|W#ûOùE¡L:ËIóxh!?!TîõaÁ/ñ»BDI§¡ )ûIþ@>Q5øÈ
+è·(Âcl£¸ch<=
+Ð äoFü]h<Ñ ÏÐpÇÑp_F>GÃ7ÖýshÐçh|]Á·t° %Éí[VéÀª¦¨dÜçyU";Öq¸.rqÛBvªO׳ëõÇ£1&úmC2ZóÀ$èKC<,}äJß®PeO7£¢Å|µü-¶,þýu
B(æ³]SÈ|9"
+C}Ô¨¼Ùðµé>
+³ï
+|8x´ïÜSÓ4°X2¸¯#
+C±iáýþ¦¢Ð¹{ÝÞnLÅIXjú
Ä´¼ÕÖtĪ#=ß92NÔé[= `æ`Ñ~ꮯ/W7O¯óõz¹º¯OÏW7ç«12E:ësáÞ_ìÑhëvË¢YðyÖç²D8ÐÐl=¹\
+ZOi«3a²J{ÄH½U£{^ÞVM(Û©è,¦/{ABS Ð'¼ÉÐm£¿~RWfÞÂ+}Ñ#rÝ ï6¿4Qªn`"ÜÕÕ¦6PµÝ^=lç ¢·íɦdÐÄç±h+öcIØç5\Ap`Ámº)-M°Öê]ÕÔgð¿1£Òë)¡\`vº§f1øÌõlp½~ûòg³¬æYÛB1ípiM";&_M"è¡à¾¯NÏ®?¨äYÿz±:ZÔC¤þfÛNΪ/ëýõIápÑ1¤õÔÆÇðuQwá2
+kíÄÛÎáí(ù_ÁóÚðÜû_d3¬>Æ}Àðÿ?7VE#+¹æÝ
¡pHä
+vø¾¸õëzhz!ÐÆÐ'z¦ À¤jvÀ µ ~5×ë`OÌ*4Z#Ì°º:¿¾\EKPàÒ=Hɵ±(o²[Vª±V¸é
DèÜ;[)ÚP5QM^íÒh¯>gOaäÈþ¬¤ÆhQchMÇÚp>¸ìÅóÏë?HØòx«äg¼lÁ¢ü¾OecL `wɲC8-îïæëùñ·°k:@AWØ\6 Yâ ëzdP²§àmgÍòJr!Ð(.#²ý°ì-{ñ˹ÜPö¦íû%û¡8Óé8rsÆsøâ$¡Y¯Î?\ºË [!Ú³QÛDAZê"úb¬ÎTVíöÙNoïz2ç¯éª$K°qAO÷£Ô?)Ï)³Î¡¦UWYq;Æ«CYKñ+-E®ªÃüãòæäüò?ó˱¯Â¾Û
+ÞDXñ¶6ð!|¢G1Å>Êî^¯ Jª©×ñ2µÇ¸x^Û¡ºlz)Õ[jSeEtä6¨_Ðïû¤¥
+´Ò}ý¹é³í
+
endstream
endobj
-683 0 obj <<
+682 0 obj <<
/Type /Page
-/Contents 684 0 R
-/Resources 682 0 R
+/Contents 683 0 R
+/Resources 681 0 R
/MediaBox [0 0 612 792]
/Parent 616 0 R
-/Annots [ 658 0 R 659 0 R 660 0 R 661 0 R 662 0 R 663 0 R 679 0 R 680 0 R ]
+/Annots [ 656 0 R 657 0 R 671 0 R 672 0 R 673 0 R 674 0 R 675 0 R 676 0 R 678 0 R 679 0 R ]
>> endobj
-678 0 obj <<
+677 0 obj <<
/Type /XObject
/Subtype /Form
/FormType 1
/PTEX.FileName (./figures/Carcione_et_al_aniso_1988_image.pdf)
/PTEX.PageNumber 1
-/PTEX.InfoDict 686 0 R
+/PTEX.InfoDict 685 0 R
/Matrix [0 -1 1 0 0 595]
/BBox [0 0 595 842]
/Resources <<
/ColorSpace <<
-/R13 687 0 R
+/R13 686 0 R
>>/ExtGState <<
-/R16 688 0 R
-/R15 689 0 R
-/R11 690 0 R
+/R16 687 0 R
+/R15 688 0 R
+/R11 689 0 R
>>/ProcSet [ /PDF /ImageC ]
/XObject <<
-/R14 691 0 R
+/R14 690 0 R
>>>>
-/Length 692 0 R
+/Length 691 0 R
/Filter /FlateDecode
>>
stream
@@ -14351,7 +14341,7 @@
%ÇæªÿÍ3H¸ffH8`a#Ojäs®î-%éðúßMQZ±âmñØõáÑïöoxô÷{þ YpKb
endstream
endobj
-686 0 obj
+685 0 obj
<<
/Author (Dimitri)
/CreationDate (D:20120520013856+02'00')
@@ -14361,34 +14351,34 @@
/Title (Untitled)
>>
endobj
+686 0 obj
+[/Separation/Black/DeviceCMYK 692 0 R]
+endobj
687 0 obj
-[/Separation/Black/DeviceCMYK 693 0 R]
-endobj
-688 0 obj
<<
/Type /ExtGState
/SA true
>>
endobj
-689 0 obj
+688 0 obj
<<
/Type /ExtGState
-/BG 694 0 R
-/UCR 695 0 R
+/BG 693 0 R
+/UCR 694 0 R
>>
endobj
-690 0 obj
+689 0 obj
<<
/Type /ExtGState
-/BG 694 0 R
-/UCR 695 0 R
+/BG 693 0 R
+/UCR 694 0 R
/OPM 1
>>
endobj
-691 0 obj
+690 0 obj
<<
/Subtype /Image
-/ColorSpace 687 0 R
+/ColorSpace 686 0 R
/Width 791
/Height 1263
/BitsPerComponent 8
@@ -15534,10 +15524,10 @@
*X¥xdI"vGF¬§ÐØÓîe3ÌÓX"ásÜã Éçx qUéÊv°lÅMuq%ÝÃÍ3ïºcØ à ªôQV$I ;ÈPNHÏ\g''Ôæ½bIä{x¢fýÜYØ '$ñÔ9<àÐW¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢Ëµã96(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢(¢¿ÿÙ
endstream
endobj
-692 0 obj
+691 0 obj
200
endobj
-693 0 obj
+692 0 obj
<<
/Filter /FlateDecode
/FunctionType 4
@@ -15549,7 +15539,7 @@
x«N)-P0PÈ-ÍQHHÎPÀÏ53kß
endstream
endobj
-694 0 obj
+693 0 obj
<<
/Filter /FlateDecode
/FunctionType 0
@@ -15563,7 +15553,7 @@
xc`Ù
endstream
endobj
-695 0 obj
+694 0 obj
<<
/Filter /FlateDecode
/FunctionType 0
@@ -15578,95 +15568,114 @@
x«¯Ù Ä1
endstream
endobj
-658 0 obj <<
+656 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [114.062 684.618 173.877 695.522]
+/Rect [468.671 706.884 511.843 717.788]
/Subtype /Link
+/A << /S /GoTo /D (cite.TaLiTr07) >>
+>> endobj
+657 0 obj <<
+/Type /Annot
+/Border[0 0 0]/H/I/C[0 1 0]
+/Rect [515.761 706.884 537.679 717.788]
+/Subtype /Link
+/A << /S /GoTo /D (cite.TaLiTr07) >>
+>> endobj
+671 0 obj <<
+/Type /Annot
+/Border[0 0 0]/H/I/C[0 1 0]
+/Rect [114.062 603.173 173.877 614.077]
+/Subtype /Link
/A << /S /GoTo /D (cite.CaKoKo88) >>
>> endobj
-659 0 obj <<
+672 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [177.693 684.618 199.61 695.522]
+/Rect [177.693 603.173 199.61 614.077]
/Subtype /Link
/A << /S /GoTo /D (cite.CaKoKo88) >>
>> endobj
-660 0 obj <<
+673 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [383.09 684.618 412.201 695.522]
+/Rect [383.09 603.173 412.201 614.077]
/Subtype /Link
/A << /S /GoTo /D (cite.Hel94) >>
>> endobj
-661 0 obj <<
+674 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [415.19 684.618 437.107 695.522]
+/Rect [415.19 603.173 437.107 614.077]
/Subtype /Link
/A << /S /GoTo /D (cite.Hel94) >>
>> endobj
-662 0 obj <<
+675 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [440.096 684.618 478.053 695.522]
+/Rect [440.096 603.173 478.053 614.077]
/Subtype /Link
/A << /S /GoTo /D (cite.Car07) >>
>> endobj
-663 0 obj <<
+676 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [481.042 684.618 502.96 695.522]
+/Rect [481.042 603.173 502.96 614.077]
/Subtype /Link
/A << /S /GoTo /D (cite.Car07) >>
>> endobj
-679 0 obj <<
+678 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [394.368 326.362 478.313 337.266]
+/Rect [394.368 241.908 478.313 252.812]
/Subtype /Link
/A << /S /GoTo /D (cite.MoTr08) >>
>> endobj
-680 0 obj <<
+679 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [482.129 326.362 504.046 337.266]
+/Rect [482.129 241.908 504.046 252.812]
/Subtype /Link
/A << /S /GoTo /D (cite.MoTr08) >>
>> endobj
-685 0 obj <<
-/D [683 0 R /XYZ 72 738 null]
+684 0 obj <<
+/D [682 0 R /XYZ 72 738 null]
>> endobj
54 0 obj <<
-/D [683 0 R /XYZ 72 720 null]
+/D [682 0 R /XYZ 72 644.369 null]
>> endobj
58 0 obj <<
-/D [683 0 R /XYZ 72 411.567 null]
+/D [682 0 R /XYZ 72 327.786 null]
>> endobj
-682 0 obj <<
-/Font << /F68 529 0 R /F44 199 0 R /F43 198 0 R /F67 351 0 R >>
-/XObject << /Im10 678 0 R >>
+681 0 obj <<
+/Font << /F68 529 0 R /F44 199 0 R /F67 351 0 R /F43 198 0 R >>
+/XObject << /Im10 677 0 R >>
/ProcSet [ /PDF /Text ]
>> endobj
703 0 obj <<
-/Length 1677
+/Length 1665
/Filter /FlateDecode
>>
stream
-xÚkoÛ6ð{~
?Ê@HÔ;@?¬i²vXÓ;h¨,ºÔ4ùõ½#²(+
-BêxoÞ~±<9½LóYá)KgËz±Yå³e5ûì¿úíÝòâÃ|ÁÀýù"IïÃ<¼OWW¯¯~·ËWvóñíó0¼¿Éßß]_^¼a/ç_^Æñ 8Éý Êffsr±< ÌÂYÈR?dÉ,"?ÓY¹=ùvøy3fP÷æÐQàôõ6g/ÕÉ{øsgÇvqÀ÷ºáH;pC~ÂÈòzÐb¾¢<ÍqVx?øv×XÕ¸Æ^/·²½¶°Zin´ ¼N
-ºqgì-7.EÛ¶¬ ¢Ã/AÀt×ÛSË0´äíµSAv£HЯÜÒî^hå£ãÁ7°·Hkïê^"jÎ{ÜE°CþsëÕ[¥=^VhÞ8jaO7óÄ»Û)´rz<Ã-8O;$¶ë¹vrºD¿qu®@à ¶\ß@zh9ÂxG,;ÀÁíÛ$q#ë~Ê;®ùVàå(N³Ã(f~¤.û`"ÐÃØgyèP"kècv,ölÌo?}8¿`å~X,ñ®1í-Z¯í§7ÒBôÎó =2¯gnY«Yì® %ï!ñ¶r§ÄrÃzQ
-÷RM{GAT¡Ýΰ4
-2?Ì#gC=é9¼Ñ¹Äpo!^Õ"_D at utcÆl(@Ð/|ÖÁÊM fAì YPd!µ
-ÇÀZàÖ~8^Q¡©¢ìå%Þ÷yx¢¹§g®¶²DĤ§Þ¨UmÛî2de?*Ñv²¿É<¢Û)á¦ÊfÒA=ÉZèíGî³TÛj!
í§"ÝBðµl£ +ý$÷ûcî÷Ïî
-èþûc îóÐÿïv|Õ±ïTã®l=g747ök«ª¡º§ùÔS×þkéMLj̾©èù?EàwÙOÊv ÙG´tÁõ/Å¿ßÇÜS}øþ>jLz^öû;aé¾Ya]g.s¹©ÉÑ𮥸ªdhAgÒZ¨¥&áo{Ó02ïÚ2O!1{."1ÕÖ§$]¾³ ¬7J¶D®tbK,lsèÌBuÊcYä_§*ê·`ÃÊûZô«uVßE£JÙK?êb²äæP±ñ¦Qè[´5-Ì$+&Ò`&üؽ~6îÜÖZ0"¶¤V_&õáö
-\o+|X¤`~xìZoÃè?Sdia(H3¥M.QÓÛÛð°äM94Ô¼ðèVØuèÜÕAk,OBu´Ô1f3"ß~è8ÃÁ<´@ìÑ¡R&nþä?HII{§(Ú&
ľÿФõ }ÝÊ~C
-Ø r-ôX%ÒaSXN,´§Jh'»^«
-úEÕuA¶¢(NSxÊbkbìÇó<!ï¿uýÊFJîÞ¹·Ì
¡ÏiL«l!Ò@ÙJ¢ÿº¨ÇzÂ|édÔXÎshp+ Ä¥Übûµù5².pʸʶT4ù#¦U1J;±p¬ÈöV:oW {Q#c1.[<Ý'Æ k(iiCļåÍäLÄ$ÊB©MÆ"Ý
- lfDX7j«`ÎWCg)¡ÄƵÂß°3×«Æ ®²4sÈRá5íZUw¤v¡¸)
yNN®IUX©E54î*ÂgعéNvî¨åZª³© ß/ò𸪾ãzUËFL½j¡-±q2?£
-¦åA=ëc ¼ !Ô;ÈÇnoBÂ!CÆ Z
-òÓÀôN\ü^ÂÐiAaâñCñæu¯»Aÿ-Áe=ôc¥òÔö¯¸ÈÕ4â7«16§ºý<-U¥½4ZÕ[Ùз¤¢whvöêÀø¦³µÉ5< å4¨h×>
-?ËÙ£ø4¼ËÞôlÂñK?èVTã°¥Òw?SüÍæö
+xÚXmoÛ8þÞ_Ü}IƵ,Û±ôÃko½Ýº·l8`
+Õ¡~_Öµ¿þH²ÔÙuw·MQù"©<[Ñ$vâÐ'Ët²ð&M«É§ééßÞ,ÏÞÍæ^àN}g6Bwúnæ»Ó¿å3C¼ýç¹îô£]ò×û7g§çg¯¼ç³/Ë?Î}c#?/&®Ù-PæàlyÀåNØy¡Ã¼`pîø~8Iò¯®±Èó´È&'íJb]äMoáÏÎÍÚùÞgÃu CÀxŹºéj9s üx6÷ñôì»È«Øe£?mU®ÃKËëZ\Svu"
+§?]®È¢$mT°&?»®W75
+H¢V¢¸±&¨ÆH¬$.(äJþ¶TNÔ¬KlæàoÆß˲xh¹hâ@¡~äÙñò5ÃÒL«Â0od!kÙÕÒÌ®gÁô¾*ÑËaì8DÀ«í´mÓÚî«ÍE&â²¥ÕBÌÏnà梾
¢çÈ
+©Ðê@ìÜcvkí¸Vi;E%jK<Åáb3=¡
+âÖ tæ;^ĬFG«ó|'XôwâýëïNÏFTòÈaq¯Ï# ÕÞ¡÷e¡ÿtâz7è¸ "ÎÓÖ³72Æ2ß3pÑHÅÊÎʵÈÒùª«E«Jòi
+¢Ö²©O¹»pXÄé(rx=¸¤pðjîVmvÐìc'ì(u¯´ÛÄòËÛ`qÃâptm ¹öÒ«Ë* UèªLZ5óé·¦2»
ÁÔ\0ß!
+ ø×'3/
+R« D ß/HA®oøåW#o¯·ùÊËÌþý2a:÷гÐàJÔ¤RXÒUÙ4ê:³÷EÃ\JkÇñÈYÛûÛ#fN<Äý±Q¨/4*½*ºüZÒ0{]âcE#ãn
+o_ªDçÝÈø7
ü»÷scðûb8xgÞg cààÃ-'þ?çÈèïFÝ@âÞÜã½{MOQèTkE^áöûwK<ñ@ÄKZü2µ#©°ØY¸Éù&êsoáD1yu÷Tih\TÓÖe¥übp{EÓªäH$eÉä7
¿ºf,³0KYTÆp |ì¼Î¨ÙcºpÃ7éMaêy³z
+3£ÁÌZ? iqÆ©5A~¡' Ãq£íã'4 WÆö¸
+µ0îú.áQfJ³Y8=_u¹µýp8½ïN8ÔÈDoo²c´Ö ç²YV"
+ÃJÊ¢ºP¢á ,Bú+¼ÒÖÏ
+²ïB>·ÓaïYC©*¨¼ËZQHØÅ89}o®HæÄV§Llm(ÕÞõ1¸.ÝZewOXÙßÏ[êrl[¾FF·W¼uno;-°IWe
+ãvÁR\«5ÀrÈvKÖþ°ýáðawæ_ì°aÿÃööó¿ØOÚo±zí?ìëñÝöÁÓe]³_O:vì?«£&%)öz?©rô"ðjÇBGwW6?¶ YKQÿãö}pòhÕ×NçbjI;ÉfÉÀ§Ñ1¾"<Û'¶Owr0l,ÐÆ0ëZBG¥6ü'V?OCYB©À¦2Z¢{m÷í¨m)ë{3ÙHzUÖ²HP+ÍüÄrïÀBu±OëF_Æú÷9¤aôènd{µaÎÕ7j<vÓx´Á vBe%âp¾±.v8âEêô{?ªÃ²¤ñy´4í&ð¥¯*\þÝôª÷:@OÅ/Á}¶
+CZ˯]ÉÄaÈ1~z'ó0¢Ê
+yIád"²¤Ëè©SwÒX¼´ÍÈ4
+ÖIb½f-s´¤þÅ@oIú_Ûÿb 4ÿx.ìÒØÐíü©è?¸.í}IHm.CËAïúÇÒj×T`@
+§ðÌ°9Ñz¢¬4
åÈPzµ²Í$6-}°TíMúN.èµ
endstream
endobj
702 0 obj <<
@@ -15677,7 +15686,7 @@
/Parent 706 0 R
/Annots [ 699 0 R 700 0 R ]
>> endobj
-681 0 obj <<
+680 0 obj <<
/Type /XObject
/Subtype /Form
/FormType 1
@@ -15771,14 +15780,14 @@
699 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [456.935 281.062 540.996 291.966]
+/Rect [456.935 169.481 540.996 180.385]
/Subtype /Link
/A << /S /GoTo /D (cite.MoTr08) >>
>> endobj
700 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [99.228 269.665 121.146 280.011]
+/Rect [99.228 158.084 121.146 168.43]
/Subtype /Link
/A << /S /GoTo /D (cite.MoTr08) >>
>> endobj
@@ -15788,22 +15797,22 @@
705 0 obj <<
/D [702 0 R /XYZ 121.344 523.483 null]
>> endobj
-62 0 obj <<
-/D [702 0 R /XYZ 72 210.577 null]
->> endobj
701 0 obj <<
/Font << /F68 529 0 R /F44 199 0 R /F67 351 0 R /F43 198 0 R >>
-/XObject << /Im11 681 0 R >>
+/XObject << /Im11 680 0 R >>
/ProcSet [ /PDF /Text ]
>> endobj
712 0 obj <<
-/Length 523
+/Length 980
/Filter /FlateDecode
>>
stream
-xÚT]o0}ϯà¤áÚÆ60i[FÖM[%´´L&X"¦Òýú9µ!mJ×ìëëûuÎ=æC:9±ÈAÌ0sÒ ±æÎwzþ~&KÏǺx>eÐ]zºóùçù's'ÆX]|õîUò}µH¦³äþèýL¿ÍyÔÐÀ t é
"ÃÂG1>¡9>µ30¡«Ëå4)P1ÖW|»ÐÄÈñÃÔܬ1¦#Ý"@46ûO;Ò @ÈÐ äÁÙÕæû[Øsio5öW¥²3v]çÉí®â¬9çõËÞk±¹9j÷s$8Ä ÒÛA¥êDSðL´
¬Ä®( ,pC\Æ-ÑÔQÓ»{9yÌÆÕuKÅ;ñ' ZX ä9\W\Yÿ¥Ë=ÝQ×z¦Vv·²+jè±k°Ì`bôÿÄ´õ®Éįîþzl,ø°²w¦7²ÃÔÑ3lqo¿¿ºúg`¼ÀßW{çZBB-Mèt¸ÚTBwÌFi ³á¯qõôéiXzG]ý|<Ll¼§Rå2ÆOÔTÌE/ídÃ^C
-¥j;Áó´ªv:
-í!NtòK&`"
+xÚVQo¤6~ϯØGV
+>0Æ@¥{hÓMÓª½S¥¦¼`«¬½5ÐÜö×w
+ÙÝøëÍx<óÍçÏc~X_¼»¦ù¢@Åt±®^dI¾XWß«ïïÖ«ûeÓ( h¦4
+î$
+>ÝÞþ|û]X߬ìäáã¯Ë8ÏÓßîVW׫øÇåë_Þ]r¤9lÙ\qî|ELPB(6Náäâv&Ö ²!SÜ(È<Ût½²cÇ{;Ù·Lr;}fKoÞYs£¢¬µ¥|*Ùy G #L©
±^æQ 9´:±ZÖsûµ×j϶Ë4`&µ©ÚBj'äÖyZ 5U 8éÓúÛ±oøIª9p¥vLÈKøÈéK1få¨1Ãu¼rhÜf&Y{èE [ ^Sa£"Mm
µÒ&á° c¸HgZÁäÁl̳µS[.ÕÐÙ;UñÖ.In Ùxl0n¸uã©·ÊZXeË`nïFUIOAVL¬?À.'(L;¨ÿ5Í+ä+t:Êù{§«FtÓdZ¨ï¬6hv,gTä1vú¸cú©-÷?C9'íñ`/èPÃ{Ä8
{RïE_6A&SL'Zðx¦à½P¯<B STóôÌ$ öJHgÊz#é× r|.òiØO³6=Hb£Ó®J y¤êmîóãìØðvoקÆÚئSz3Þ5³àD"xwi¹Rü¢@Y_é_º.áDÑ5j0ùÆ5~Ð;\h45òp(úÛH>~º¿ZyHHRàI#RÄ
xGúºÛý¯9¹9fjSkܤc;~Ú9ç¥mÓÏí;êWæÛö+7çQgòO×8Ã@Û°{®kü+w+OÍçºìuÄM<â~¹/o3'ÎÙøJéJÈñ=ðr2Ññò
Ëѳtüµ\kµî þcÕO1èzPrǼhEüÿéÔ Mö©?ì}g at s_ì½Í;0cÓ62zUÛ¿Dwñû>ÿgÄàZ#`ÿ^GSFQcGSüí<1¹mùø ùh2ó/PÉäéÕ² ûpÿ/óO÷T@-g¥±oÔü
OõìÎ$£_nôcFBv=gÕW´
+öà/VëÇ©
endstream
endobj
711 0 obj <<
@@ -15816,8 +15825,11 @@
713 0 obj <<
/D [711 0 R /XYZ 72 738 null]
>> endobj
+62 0 obj <<
+/D [711 0 R /XYZ 72 720 null]
+>> endobj
710 0 obj <<
-/Font << /F68 529 0 R /F44 199 0 R /F67 351 0 R >>
+/Font << /F68 529 0 R /F44 199 0 R /F43 198 0 R /F67 351 0 R >>
/ProcSet [ /PDF /Text ]
>> endobj
716 0 obj <<
@@ -15898,7 +15910,7 @@
/D [719 0 R /XYZ 72 392.335 null]
>> endobj
718 0 obj <<
-/Font << /F68 529 0 R /F44 199 0 R /F43 198 0 R /F90 723 0 R /F8 653 0 R /F11 357 0 R /F13 724 0 R /F14 668 0 R /F67 351 0 R >>
+/Font << /F68 529 0 R /F44 199 0 R /F43 198 0 R /F90 723 0 R /F8 653 0 R /F11 357 0 R /F13 724 0 R /F14 662 0 R /F67 351 0 R >>
/ProcSet [ /PDF /Text ]
>> endobj
728 0 obj <<
@@ -15932,7 +15944,7 @@
/D [727 0 R /XYZ 72 720 null]
>> endobj
726 0 obj <<
-/Font << /F43 198 0 R /F44 199 0 R /F90 723 0 R /F14 668 0 R /F8 653 0 R /F11 357 0 R /F7 730 0 R /F10 731 0 R /F13 724 0 R >>
+/Font << /F43 198 0 R /F44 199 0 R /F90 723 0 R /F14 662 0 R /F8 653 0 R /F11 357 0 R /F7 730 0 R /F10 731 0 R /F13 724 0 R >>
/ProcSet [ /PDF /Text ]
>> endobj
736 0 obj <<
@@ -16470,7 +16482,7 @@
535 0 obj <<
/D [765 0 R /XYZ 47.093 412.073 null]
>> endobj
-677 0 obj <<
+695 0 obj <<
/D [765 0 R /XYZ 47.093 356.384 null]
>> endobj
408 0 obj <<
@@ -17934,7 +17946,7 @@
/LastChar 50
/Widths 777 0 R
>> endobj
-668 0 obj <<
+662 0 obj <<
/Type /Font
/Subtype /Type1
/BaseFont /RZHJNT+CMSY10
@@ -18048,7 +18060,7 @@
/Type /Pages
/Count 6
/Parent 820 0 R
-/Kids [612 0 R 628 0 R 635 0 R 650 0 R 665 0 R 683 0 R]
+/Kids [612 0 R 628 0 R 635 0 R 650 0 R 659 0 R 682 0 R]
>> endobj
706 0 obj <<
/Type /Pages
@@ -18276,7 +18288,7 @@
/Limits [(cite.PuAmKa09) (cite.SePr94)]
>> endobj
837 0 obj <<
-/Names [(cite.SePrPr95) 553 0 R (cite.TaCrEtViBeSa10) 535 0 R (cite.TaLiTr07) 677 0 R (cite.TaWi00) 408 0 R (cite.TrKoLi08) 365 0 R (cite.TsKoChTr03) 399 0 R]
+/Names [(cite.SePrPr95) 553 0 R (cite.TaCrEtViBeSa10) 535 0 R (cite.TaLiTr07) 695 0 R (cite.TaWi00) 408 0 R (cite.TrKoLi08) 365 0 R (cite.TsKoChTr03) 399 0 R]
/Limits [(cite.SePrPr95) (cite.TsKoChTr03)]
>> endobj
838 0 obj <<
@@ -18288,7 +18300,7 @@
/Limits [(figure.3.3) (page.10)]
>> endobj
840 0 obj <<
-/Names [(page.11) 614 0 R (page.12) 630 0 R (page.13) 637 0 R (page.14) 652 0 R (page.15) 667 0 R (page.16) 685 0 R]
+/Names [(page.11) 614 0 R (page.12) 630 0 R (page.13) 637 0 R (page.14) 652 0 R (page.15) 661 0 R (page.16) 684 0 R]
/Limits [(page.11) (page.16)]
>> endobj
841 0 obj <<
@@ -18304,11 +18316,11 @@
/Limits [(page.28) (page.6)]
>> endobj
844 0 obj <<
-/Names [(page.7) 581 0 R (page.8) 588 0 R (page.9) 594 0 R (section*.1) 203 0 R (section*.10) 772 0 R (section*.3) 669 0 R]
+/Names [(page.7) 581 0 R (page.8) 588 0 R (page.9) 594 0 R (section*.1) 203 0 R (section*.10) 772 0 R (section*.3) 663 0 R]
/Limits [(page.7) (section*.3)]
>> endobj
845 0 obj <<
-/Names [(section*.4) 670 0 R (section*.5) 722 0 R (section*.6) 725 0 R (section.1.1) 10 0 R (section.1.2) 14 0 R (section.2.1) 22 0 R]
+/Names [(section*.4) 664 0 R (section*.5) 722 0 R (section*.6) 725 0 R (section.1.1) 10 0 R (section.1.2) 14 0 R (section.2.1) 22 0 R]
/Limits [(section*.4) (section.2.1)]
>> endobj
846 0 obj <<
@@ -18356,8 +18368,8 @@
>> endobj
856 0 obj <<
/Author()/Title()/Subject()/Creator(LaTeX with hyperref package)/Producer(pdfTeX-1.40.3)/Keywords()
-/CreationDate (D:20120903015511+02'00')
-/ModDate (D:20120903015511+02'00')
+/CreationDate (D:20120903035538+02'00')
+/ModDate (D:20120903035538+02'00')
/Trapped /False
/PTEX.Fullbanner (This is pdfTeX using libpoppler, Version 3.141592-1.40.3-2.2 (Web2C 7.5.6) kpathsea version 3.5.6)
>> endobj
@@ -18369,80 +18381,80 @@
0000000004 00000 f
0000000000 00000 f
0000000015 00000 n
-0000504205 00000 n
-0001792473 00000 n
+0000504206 00000 n
+0001792849 00000 n
0000000060 00000 n
0000000092 00000 n
-0000526912 00000 n
-0001792401 00000 n
+0000526913 00000 n
+0001792777 00000 n
0000000139 00000 n
0000000170 00000 n
-0000526968 00000 n
-0001792328 00000 n
+0000526969 00000 n
+0001792704 00000 n
0000000218 00000 n
0000000248 00000 n
-0000530491 00000 n
-0001792204 00000 n
+0000530492 00000 n
+0001792580 00000 n
0000000294 00000 n
0000000330 00000 n
-0000533089 00000 n
-0001792143 00000 n
+0000533090 00000 n
+0001792519 00000 n
0000000378 00000 n
0000000459 00000 n
-0000895008 00000 n
-0001792018 00000 n
+0000895009 00000 n
+0001792394 00000 n
0000000505 00000 n
0000000541 00000 n
-0000895060 00000 n
-0001791944 00000 n
+0000895061 00000 n
+0001792320 00000 n
0000000589 00000 n
0000000632 00000 n
-0000967711 00000 n
-0001791857 00000 n
+0000967712 00000 n
+0001792233 00000 n
0000000680 00000 n
0000000747 00000 n
-0001038072 00000 n
-0001791770 00000 n
+0001038073 00000 n
+0001792146 00000 n
0000000795 00000 n
0000000861 00000 n
-0001038128 00000 n
-0001791696 00000 n
+0001038129 00000 n
+0001792072 00000 n
0000000909 00000 n
0000000979 00000 n
-0001391645 00000 n
-0001791571 00000 n
+0001391646 00000 n
+0001791947 00000 n
0000001025 00000 n
0000001075 00000 n
-0001406603 00000 n
-0001791497 00000 n
+0001406260 00000 n
+0001791873 00000 n
0000001123 00000 n
0000001197 00000 n
-0001609235 00000 n
-0001791410 00000 n
+0001609154 00000 n
+0001791786 00000 n
0000001245 00000 n
0000001289 00000 n
-0001609287 00000 n
-0001791323 00000 n
+0001609210 00000 n
+0001791699 00000 n
0000001337 00000 n
0000001385 00000 n
-0001617194 00000 n
-0001791249 00000 n
+0001618470 00000 n
+0001791625 00000 n
0000001433 00000 n
0000001500 00000 n
-0001620222 00000 n
-0001791124 00000 n
+0001620598 00000 n
+0001791500 00000 n
0000001546 00000 n
0000001586 00000 n
-0001620274 00000 n
-0001791063 00000 n
+0001620650 00000 n
+0001791439 00000 n
0000001634 00000 n
0000001697 00000 n
-0001624254 00000 n
-0001790975 00000 n
+0001624630 00000 n
+0001791351 00000 n
0000001743 00000 n
0000001796 00000 n
-0001652053 00000 n
-0001790900 00000 n
+0001652429 00000 n
+0001791276 00000 n
0000001843 00000 n
0000001879 00000 n
0000002071 00000 n
@@ -18451,7 +18463,7 @@
0000001930 00000 n
0000474188 00000 n
0000474239 00000 n
-0001790065 00000 n
+0001790441 00000 n
0000006189 00000 n
0000006380 00000 n
0000006428 00000 n
@@ -18561,670 +18573,670 @@
0000474810 00000 n
0000474699 00000 n
0000474357 00000 n
-0001789359 00000 n
-0001789532 00000 n
+0001789735 00000 n
+0001789908 00000 n
0000475767 00000 n
0000475599 00000 n
0000474895 00000 n
0000475710 00000 n
-0000476999 00000 n
-0000477150 00000 n
-0000477303 00000 n
-0000477456 00000 n
-0000477605 00000 n
-0000477758 00000 n
-0000477909 00000 n
-0000478062 00000 n
-0000478215 00000 n
-0000478367 00000 n
-0000478520 00000 n
-0000478670 00000 n
-0000478823 00000 n
-0000478976 00000 n
-0000479128 00000 n
-0000479281 00000 n
-0000479431 00000 n
-0000479583 00000 n
-0000479734 00000 n
-0000479995 00000 n
-0000476724 00000 n
+0000477000 00000 n
+0000477151 00000 n
+0000477304 00000 n
+0000477457 00000 n
+0000477606 00000 n
+0000477759 00000 n
+0000477910 00000 n
+0000478063 00000 n
+0000478216 00000 n
+0000478368 00000 n
+0000478521 00000 n
+0000478671 00000 n
+0000478824 00000 n
+0000478977 00000 n
+0000479129 00000 n
+0000479282 00000 n
+0000479432 00000 n
+0000479584 00000 n
+0000479735 00000 n
+0000479996 00000 n
+0000476725 00000 n
0000475852 00000 n
-0000479885 00000 n
-0000479938 00000 n
-0000484324 00000 n
-0000484664 00000 n
-0000484831 00000 n
-0000485006 00000 n
-0000485160 00000 n
-0000485314 00000 n
-0000485624 00000 n
-0000485779 00000 n
-0000485935 00000 n
-0000486089 00000 n
-0000486242 00000 n
-0000486400 00000 n
-0000486558 00000 n
-0000486714 00000 n
-0000486870 00000 n
-0000487220 00000 n
-0000487395 00000 n
-0000487820 00000 n
-0000487976 00000 n
-0000488132 00000 n
-0000488307 00000 n
-0000488483 00000 n
-0000488639 00000 n
-0000488795 00000 n
-0000488947 00000 n
-0000489099 00000 n
-0000489255 00000 n
-0000489411 00000 n
-0000489567 00000 n
-0000489723 00000 n
-0000489879 00000 n
-0000490034 00000 n
-0000490187 00000 n
-0000490338 00000 n
-0000490494 00000 n
-0000490652 00000 n
-0000490806 00000 n
-0000490960 00000 n
-0000491114 00000 n
-0000491267 00000 n
-0000491420 00000 n
-0000491573 00000 n
-0000491727 00000 n
-0000491881 00000 n
-0000492037 00000 n
-0000492193 00000 n
-0000492349 00000 n
-0000492505 00000 n
-0000492659 00000 n
-0000492815 00000 n
-0000492969 00000 n
-0000493122 00000 n
-0000493276 00000 n
-0000493430 00000 n
-0000493586 00000 n
-0000493742 00000 n
-0000493898 00000 n
-0000494054 00000 n
-0000494209 00000 n
-0000494365 00000 n
-0000494519 00000 n
-0000494672 00000 n
-0000494828 00000 n
-0000494984 00000 n
-0000495140 00000 n
-0000495295 00000 n
-0000495453 00000 n
-0000495611 00000 n
-0000495764 00000 n
-0000495918 00000 n
-0000496076 00000 n
-0000496234 00000 n
-0000496390 00000 n
-0000496547 00000 n
-0000496700 00000 n
-0000496853 00000 n
-0000497007 00000 n
-0000497161 00000 n
-0000497315 00000 n
-0000497469 00000 n
-0000497623 00000 n
-0000497777 00000 n
-0000498086 00000 n
-0000498240 00000 n
-0000498388 00000 n
-0000498540 00000 n
-0000498693 00000 n
-0000498845 00000 n
-0000498998 00000 n
-0000499151 00000 n
-0000499309 00000 n
-0000499467 00000 n
-0000499625 00000 n
-0000499783 00000 n
-0000499940 00000 n
-0000500097 00000 n
-0000500261 00000 n
-0000500426 00000 n
-0000500581 00000 n
-0000500736 00000 n
-0000501045 00000 n
-0000501199 00000 n
-0000501356 00000 n
-0000501514 00000 n
-0000501667 00000 n
-0000501820 00000 n
-0000501975 00000 n
-0000502130 00000 n
-0000502287 00000 n
-0000502443 00000 n
-0000502600 00000 n
-0000502757 00000 n
-0000502911 00000 n
-0000503065 00000 n
-0000503225 00000 n
-0000503385 00000 n
-0000503538 00000 n
-0000503691 00000 n
-0000503999 00000 n
-0000504256 00000 n
-0000483193 00000 n
-0000480080 00000 n
-0000504152 00000 n
-0001788838 00000 n
-0000484494 00000 n
-0000485470 00000 n
-0000487046 00000 n
-0000487608 00000 n
-0000497933 00000 n
-0001787809 00000 n
-0000500891 00000 n
-0000503847 00000 n
-0001639857 00000 n
-0001644837 00000 n
-0001644898 00000 n
-0001644715 00000 n
-0001644959 00000 n
-0001650100 00000 n
-0001649430 00000 n
-0001630472 00000 n
-0001630777 00000 n
-0001640647 00000 n
-0001634658 00000 n
-0001644472 00000 n
-0001640587 00000 n
-0001650344 00000 n
-0001630046 00000 n
-0001649551 00000 n
-0001630168 00000 n
-0001635451 00000 n
-0001639675 00000 n
-0001634780 00000 n
-0001635512 00000 n
-0001649612 00000 n
-0001645142 00000 n
-0001635694 00000 n
-0001630107 00000 n
-0001630411 00000 n
-0001635268 00000 n
-0001635085 00000 n
-0001635146 00000 n
-0001635024 00000 n
-0001635329 00000 n
-0001645324 00000 n
-0001644655 00000 n
-0001649795 00000 n
-0001635207 00000 n
-0001634841 00000 n
-0001634963 00000 n
-0001649734 00000 n
-0001640161 00000 n
-0001650161 00000 n
-0001640283 00000 n
-0001630655 00000 n
-0001640526 00000 n
-0001640344 00000 n
-0001640466 00000 n
-0001640405 00000 n
-0001639736 00000 n
-0001650405 00000 n
-0001650039 00000 n
-0001640100 00000 n
-0001645020 00000 n
-0000508361 00000 n
-0000508521 00000 n
-0000508681 00000 n
-0000508839 00000 n
-0000508995 00000 n
-0000509150 00000 n
-0000509312 00000 n
-0000509473 00000 n
-0000509628 00000 n
-0000509783 00000 n
-0000509945 00000 n
-0000510106 00000 n
-0000510268 00000 n
-0000510430 00000 n
-0000510590 00000 n
-0000510748 00000 n
-0000510924 00000 n
-0000511098 00000 n
-0000511254 00000 n
-0000511410 00000 n
-0000511564 00000 n
-0000511718 00000 n
-0000511872 00000 n
-0000512026 00000 n
-0000512188 00000 n
-0000512350 00000 n
-0000512508 00000 n
-0000512666 00000 n
-0000512826 00000 n
-0000512988 00000 n
-0000513154 00000 n
-0000513320 00000 n
-0000513477 00000 n
-0000513635 00000 n
-0000513797 00000 n
-0000513959 00000 n
-0000514112 00000 n
-0000514265 00000 n
-0000514572 00000 n
-0000514726 00000 n
-0000514880 00000 n
-0000515034 00000 n
-0000515188 00000 n
-0000515342 00000 n
-0000515495 00000 n
-0000515648 00000 n
-0000515802 00000 n
-0000515955 00000 n
-0000516111 00000 n
-0000516265 00000 n
-0000516417 00000 n
-0000516569 00000 n
-0000516727 00000 n
-0000516885 00000 n
-0000517039 00000 n
-0000517192 00000 n
-0000517346 00000 n
-0000517500 00000 n
-0000517652 00000 n
-0000517804 00000 n
-0000517958 00000 n
-0000518112 00000 n
-0000518268 00000 n
-0000518423 00000 n
-0000518579 00000 n
-0000518734 00000 n
-0000518889 00000 n
-0000519045 00000 n
-0000519198 00000 n
-0000519351 00000 n
-0000519505 00000 n
-0000519660 00000 n
-0000519816 00000 n
-0000519973 00000 n
-0000520127 00000 n
-0000520280 00000 n
-0000520433 00000 n
-0000520586 00000 n
-0000520741 00000 n
-0000520897 00000 n
-0000521073 00000 n
-0000521249 00000 n
-0000521408 00000 n
-0000521568 00000 n
-0000521728 00000 n
-0000521888 00000 n
-0000522049 00000 n
-0000522207 00000 n
-0000522364 00000 n
-0000522522 00000 n
-0000522679 00000 n
-0000522833 00000 n
-0000522991 00000 n
-0000523148 00000 n
-0000523310 00000 n
-0000523472 00000 n
-0000523626 00000 n
-0000523780 00000 n
-0000523934 00000 n
-0000524088 00000 n
-0000524411 00000 n
-0000524573 00000 n
-0000524727 00000 n
-0000524881 00000 n
-0000525262 00000 n
-0000525417 00000 n
-0000525573 00000 n
-0000525749 00000 n
-0000525924 00000 n
-0000526078 00000 n
-0000526232 00000 n
-0000526388 00000 n
-0000526544 00000 n
-0000526701 00000 n
-0000527024 00000 n
-0000507302 00000 n
-0000504367 00000 n
-0000526859 00000 n
-0001789882 00000 n
-0000514419 00000 n
-0000524250 00000 n
-0000525072 00000 n
-0001630228 00000 n
-0001630289 00000 n
-0001649917 00000 n
-0001639917 00000 n
-0001640222 00000 n
-0001640708 00000 n
-0001644350 00000 n
-0001644411 00000 n
-0001640039 00000 n
-0001639978 00000 n
-0001634719 00000 n
-0001635573 00000 n
-0001650283 00000 n
-0001649673 00000 n
-0001645446 00000 n
-0001634902 00000 n
-0001639796 00000 n
-0001635390 00000 n
-0001630716 00000 n
-0001649491 00000 n
-0001649856 00000 n
-0001635634 00000 n
-0001645385 00000 n
-0001645081 00000 n
-0001650222 00000 n
-0001644594 00000 n
-0001645203 00000 n
-0001644533 00000 n
-0001645264 00000 n
-0001644776 00000 n
-0000527647 00000 n
-0000527482 00000 n
-0000527135 00000 n
-0000527594 00000 n
-0001790180 00000 n
-0000529943 00000 n
-0000530096 00000 n
-0000530250 00000 n
-0000532629 00000 n
-0000530543 00000 n
-0000529795 00000 n
-0000527732 00000 n
-0000530438 00000 n
-0001788667 00000 n
-0001649369 00000 n
-0000533145 00000 n
-0000532489 00000 n
-0000530654 00000 n
-0000533036 00000 n
-0000532833 00000 n
-0000535038 00000 n
-0000894766 00000 n
-0000895178 00000 n
-0000534906 00000 n
-0000533256 00000 n
-0000894955 00000 n
-0000895116 00000 n
-0000898042 00000 n
-0000964389 00000 n
-0000897930 00000 n
-0000895312 00000 n
-0000964274 00000 n
-0000964327 00000 n
-0000962258 00000 n
-0000962520 00000 n
-0000963698 00000 n
-0000963721 00000 n
-0000963909 00000 n
-0000967505 00000 n
-0000969830 00000 n
-0000978341 00000 n
-0000967827 00000 n
-0000967373 00000 n
-0000964528 00000 n
-0000967658 00000 n
-0001789009 00000 n
-0001789184 00000 n
-0000967767 00000 n
-0001038184 00000 n
-0000969718 00000 n
-0000967990 00000 n
-0001037957 00000 n
-0001038010 00000 n
-0001790297 00000 n
-0000970137 00000 n
-0000970317 00000 n
-0000978271 00000 n
-0000978291 00000 n
-0000978319 00000 n
-0000978648 00000 n
-0000978828 00000 n
-0001037886 00000 n
-0001037906 00000 n
-0001037934 00000 n
-0001039091 00000 n
-0001038926 00000 n
-0001038349 00000 n
-0001039038 00000 n
-0001041456 00000 n
-0001120215 00000 n
-0001233238 00000 n
-0001391758 00000 n
-0001041344 00000 n
-0001039176 00000 n
-0001391592 00000 n
-0001391697 00000 n
-0001395340 00000 n
-0001395537 00000 n
-0001395695 00000 n
-0001395853 00000 n
-0001396006 00000 n
-0001396161 00000 n
-0001396313 00000 n
-0001396468 00000 n
-0001396623 00000 n
-0001399164 00000 n
-0001396839 00000 n
-0001395144 00000 n
-0001391918 00000 n
-0001396786 00000 n
-0001788097 00000 n
-0001630350 00000 n
-0001406428 00000 n
-0001406063 00000 n
-0001406219 00000 n
-0001607954 00000 n
-0001608110 00000 n
-0001608265 00000 n
-0001608417 00000 n
-0001608569 00000 n
-0001608722 00000 n
-0001406659 00000 n
-0001399024 00000 n
-0001396975 00000 n
-0001406375 00000 n
-0001788380 00000 n
-0001406490 00000 n
-0001406547 00000 n
-0001403293 00000 n
-0001403504 00000 n
-0001403560 00000 n
-0001403954 00000 n
-0001403976 00000 n
-0001404257 00000 n
-0001649978 00000 n
-0001409254 00000 n
-0001608874 00000 n
-0001609028 00000 n
-0001611381 00000 n
-0001609343 00000 n
-0001409066 00000 n
-0001406836 00000 n
-0001609182 00000 n
-0001409879 00000 n
-0001410112 00000 n
-0001410168 00000 n
-0001410217 00000 n
-0001410282 00000 n
-0001410354 00000 n
-0001607431 00000 n
-0001607452 00000 n
-0001607611 00000 n
-0001607771 00000 n
-0001630594 00000 n
-0001635755 00000 n
-0001630533 00000 n
-0001616772 00000 n
-0001616926 00000 n
-0001617250 00000 n
-0001611241 00000 n
-0001609483 00000 n
-0001617079 00000 n
-0001617132 00000 n
-0001790414 00000 n
-0001616503 00000 n
-0001616703 00000 n
-0001616750 00000 n
-0001618159 00000 n
-0001617994 00000 n
-0001617390 00000 n
-0001618106 00000 n
-0001620330 00000 n
-0001620057 00000 n
-0001618257 00000 n
-0001620169 00000 n
-0001622573 00000 n
-0001622294 00000 n
-0001620428 00000 n
-0001622406 00000 n
-0001622459 00000 n
-0001787664 00000 n
-0001788523 00000 n
-0001622516 00000 n
-0001624306 00000 n
-0001624089 00000 n
-0001622748 00000 n
-0001624201 00000 n
-0001788239 00000 n
-0001787953 00000 n
-0001625407 00000 n
-0001625746 00000 n
-0001626028 00000 n
-0001625259 00000 n
-0001624480 00000 n
-0001625918 00000 n
-0001625971 00000 n
-0001625577 00000 n
-0001626699 00000 n
-0001626477 00000 n
-0001626126 00000 n
-0001626589 00000 n
-0001626642 00000 n
-0001790531 00000 n
-0001630838 00000 n
-0001629824 00000 n
-0001626784 00000 n
-0001629936 00000 n
-0001629989 00000 n
-0001789705 00000 n
-0001635816 00000 n
-0001634493 00000 n
-0001630936 00000 n
-0001634605 00000 n
-0001640769 00000 n
-0001639510 00000 n
-0001635914 00000 n
-0001639622 00000 n
-0001645505 00000 n
-0001644185 00000 n
-0001640867 00000 n
-0001644297 00000 n
-0001650466 00000 n
-0001649204 00000 n
-0001645603 00000 n
-0001649316 00000 n
-0001652162 00000 n
-0001651888 00000 n
-0001650564 00000 n
-0001652000 00000 n
-0001652105 00000 n
-0001790648 00000 n
-0001787053 00000 n
-0001652260 00000 n
-0001653156 00000 n
-0001653181 00000 n
-0001653206 00000 n
-0001653231 00000 n
-0001653256 00000 n
-0001653914 00000 n
-0001654041 00000 n
-0001654156 00000 n
-0001654543 00000 n
-0001654774 00000 n
-0001654969 00000 n
-0001655562 00000 n
-0001656027 00000 n
-0001657001 00000 n
-0001657478 00000 n
-0001660005 00000 n
-0001660227 00000 n
-0001664535 00000 n
-0001664811 00000 n
-0001666541 00000 n
-0001666760 00000 n
-0001669297 00000 n
-0001669560 00000 n
-0001671293 00000 n
-0001671513 00000 n
-0001673317 00000 n
-0001673589 00000 n
-0001674971 00000 n
-0001675197 00000 n
-0001685142 00000 n
-0001685428 00000 n
-0001703425 00000 n
-0001703938 00000 n
-0001715123 00000 n
-0001715500 00000 n
-0001720776 00000 n
-0001721050 00000 n
-0001737222 00000 n
-0001737669 00000 n
-0001758267 00000 n
-0001758942 00000 n
-0001770091 00000 n
-0001770408 00000 n
-0001786603 00000 n
-0001790725 00000 n
-0001790826 00000 n
-0001792582 00000 n
-0001792773 00000 n
-0001792955 00000 n
-0001793169 00000 n
-0001793409 00000 n
-0001793623 00000 n
-0001793827 00000 n
-0001794059 00000 n
-0001794280 00000 n
-0001794507 00000 n
-0001794718 00000 n
-0001794946 00000 n
-0001795164 00000 n
-0001795370 00000 n
-0001795596 00000 n
-0001795802 00000 n
-0001796028 00000 n
-0001796250 00000 n
-0001796433 00000 n
-0001796603 00000 n
-0001796772 00000 n
-0001796942 00000 n
-0001797107 00000 n
-0001797286 00000 n
-0001797481 00000 n
-0001797677 00000 n
-0001797810 00000 n
-0001797925 00000 n
-0001798050 00000 n
-0001798163 00000 n
-0001798276 00000 n
-0001798361 00000 n
-0001798468 00000 n
-0001798506 00000 n
-0001798716 00000 n
+0000479886 00000 n
+0000479939 00000 n
+0000484325 00000 n
+0000484665 00000 n
+0000484832 00000 n
+0000485007 00000 n
+0000485161 00000 n
+0000485315 00000 n
+0000485625 00000 n
+0000485780 00000 n
+0000485936 00000 n
+0000486090 00000 n
+0000486243 00000 n
+0000486401 00000 n
+0000486559 00000 n
+0000486715 00000 n
+0000486871 00000 n
+0000487221 00000 n
+0000487396 00000 n
+0000487821 00000 n
+0000487977 00000 n
+0000488133 00000 n
+0000488308 00000 n
+0000488484 00000 n
+0000488640 00000 n
+0000488796 00000 n
+0000488948 00000 n
+0000489100 00000 n
+0000489256 00000 n
+0000489412 00000 n
+0000489568 00000 n
+0000489724 00000 n
+0000489880 00000 n
+0000490035 00000 n
+0000490188 00000 n
+0000490339 00000 n
+0000490495 00000 n
+0000490653 00000 n
+0000490807 00000 n
+0000490961 00000 n
+0000491115 00000 n
+0000491268 00000 n
+0000491421 00000 n
+0000491574 00000 n
+0000491728 00000 n
+0000491882 00000 n
+0000492038 00000 n
+0000492194 00000 n
+0000492350 00000 n
+0000492506 00000 n
+0000492660 00000 n
+0000492816 00000 n
+0000492970 00000 n
+0000493123 00000 n
+0000493277 00000 n
+0000493431 00000 n
+0000493587 00000 n
+0000493743 00000 n
+0000493899 00000 n
+0000494055 00000 n
+0000494210 00000 n
+0000494366 00000 n
+0000494520 00000 n
+0000494673 00000 n
+0000494829 00000 n
+0000494985 00000 n
+0000495141 00000 n
+0000495296 00000 n
+0000495454 00000 n
+0000495612 00000 n
+0000495765 00000 n
+0000495919 00000 n
+0000496077 00000 n
+0000496235 00000 n
+0000496391 00000 n
+0000496548 00000 n
+0000496701 00000 n
+0000496854 00000 n
+0000497008 00000 n
+0000497162 00000 n
+0000497316 00000 n
+0000497470 00000 n
+0000497624 00000 n
+0000497778 00000 n
+0000498087 00000 n
+0000498241 00000 n
+0000498389 00000 n
+0000498541 00000 n
+0000498694 00000 n
+0000498846 00000 n
+0000498999 00000 n
+0000499152 00000 n
+0000499310 00000 n
+0000499468 00000 n
+0000499626 00000 n
+0000499784 00000 n
+0000499941 00000 n
+0000500098 00000 n
+0000500262 00000 n
+0000500427 00000 n
+0000500582 00000 n
+0000500737 00000 n
+0000501046 00000 n
+0000501200 00000 n
+0000501357 00000 n
+0000501515 00000 n
+0000501668 00000 n
+0000501821 00000 n
+0000501976 00000 n
+0000502131 00000 n
+0000502288 00000 n
+0000502444 00000 n
+0000502601 00000 n
+0000502758 00000 n
+0000502912 00000 n
+0000503066 00000 n
+0000503226 00000 n
+0000503386 00000 n
+0000503539 00000 n
+0000503692 00000 n
+0000504000 00000 n
+0000504257 00000 n
+0000483194 00000 n
+0000480081 00000 n
+0000504153 00000 n
+0001789214 00000 n
+0000484495 00000 n
+0000485471 00000 n
+0000487047 00000 n
+0000487609 00000 n
+0000497934 00000 n
+0001788185 00000 n
+0000500892 00000 n
+0000503848 00000 n
+0001640233 00000 n
+0001645213 00000 n
+0001645274 00000 n
+0001645091 00000 n
+0001645335 00000 n
+0001650476 00000 n
+0001649806 00000 n
+0001630848 00000 n
+0001631153 00000 n
+0001641023 00000 n
+0001635034 00000 n
+0001644848 00000 n
+0001640963 00000 n
+0001650720 00000 n
+0001630422 00000 n
+0001649927 00000 n
+0001630544 00000 n
+0001635827 00000 n
+0001640051 00000 n
+0001635156 00000 n
+0001635888 00000 n
+0001649988 00000 n
+0001645518 00000 n
+0001636070 00000 n
+0001630483 00000 n
+0001630787 00000 n
+0001635644 00000 n
+0001635461 00000 n
+0001635522 00000 n
+0001635400 00000 n
+0001635705 00000 n
+0001645700 00000 n
+0001645031 00000 n
+0001650171 00000 n
+0001635583 00000 n
+0001635217 00000 n
+0001635339 00000 n
+0001650110 00000 n
+0001640537 00000 n
+0001650537 00000 n
+0001640659 00000 n
+0001631031 00000 n
+0001640902 00000 n
+0001640720 00000 n
+0001640842 00000 n
+0001640781 00000 n
+0001640112 00000 n
+0001650781 00000 n
+0001650415 00000 n
+0001640476 00000 n
+0001645396 00000 n
+0000508362 00000 n
+0000508522 00000 n
+0000508682 00000 n
+0000508840 00000 n
+0000508996 00000 n
+0000509151 00000 n
+0000509313 00000 n
+0000509474 00000 n
+0000509629 00000 n
+0000509784 00000 n
+0000509946 00000 n
+0000510107 00000 n
+0000510269 00000 n
+0000510431 00000 n
+0000510591 00000 n
+0000510749 00000 n
+0000510925 00000 n
+0000511099 00000 n
+0000511255 00000 n
+0000511411 00000 n
+0000511565 00000 n
+0000511719 00000 n
+0000511873 00000 n
+0000512027 00000 n
+0000512189 00000 n
+0000512351 00000 n
+0000512509 00000 n
+0000512667 00000 n
+0000512827 00000 n
+0000512989 00000 n
+0000513155 00000 n
+0000513321 00000 n
+0000513478 00000 n
+0000513636 00000 n
+0000513798 00000 n
+0000513960 00000 n
+0000514113 00000 n
+0000514266 00000 n
+0000514573 00000 n
+0000514727 00000 n
+0000514881 00000 n
+0000515035 00000 n
+0000515189 00000 n
+0000515343 00000 n
+0000515496 00000 n
+0000515649 00000 n
+0000515803 00000 n
+0000515956 00000 n
+0000516112 00000 n
+0000516266 00000 n
+0000516418 00000 n
+0000516570 00000 n
+0000516728 00000 n
+0000516886 00000 n
+0000517040 00000 n
+0000517193 00000 n
+0000517347 00000 n
+0000517501 00000 n
+0000517653 00000 n
+0000517805 00000 n
+0000517959 00000 n
+0000518113 00000 n
+0000518269 00000 n
+0000518424 00000 n
+0000518580 00000 n
+0000518735 00000 n
+0000518890 00000 n
+0000519046 00000 n
+0000519199 00000 n
+0000519352 00000 n
+0000519506 00000 n
+0000519661 00000 n
+0000519817 00000 n
+0000519974 00000 n
+0000520128 00000 n
+0000520281 00000 n
+0000520434 00000 n
+0000520587 00000 n
+0000520742 00000 n
+0000520898 00000 n
+0000521074 00000 n
+0000521250 00000 n
+0000521409 00000 n
+0000521569 00000 n
+0000521729 00000 n
+0000521889 00000 n
+0000522050 00000 n
+0000522208 00000 n
+0000522365 00000 n
+0000522523 00000 n
+0000522680 00000 n
+0000522834 00000 n
+0000522992 00000 n
+0000523149 00000 n
+0000523311 00000 n
+0000523473 00000 n
+0000523627 00000 n
+0000523781 00000 n
+0000523935 00000 n
+0000524089 00000 n
+0000524412 00000 n
+0000524574 00000 n
+0000524728 00000 n
+0000524882 00000 n
+0000525263 00000 n
+0000525418 00000 n
+0000525574 00000 n
+0000525750 00000 n
+0000525925 00000 n
+0000526079 00000 n
+0000526233 00000 n
+0000526389 00000 n
+0000526545 00000 n
+0000526702 00000 n
+0000527025 00000 n
+0000507303 00000 n
+0000504368 00000 n
+0000526860 00000 n
+0001790258 00000 n
+0000514420 00000 n
+0000524251 00000 n
+0000525073 00000 n
+0001630604 00000 n
+0001630665 00000 n
+0001650293 00000 n
+0001640293 00000 n
+0001640598 00000 n
+0001641084 00000 n
+0001644726 00000 n
+0001644787 00000 n
+0001640415 00000 n
+0001640354 00000 n
+0001635095 00000 n
+0001635949 00000 n
+0001650659 00000 n
+0001650049 00000 n
+0001645822 00000 n
+0001635278 00000 n
+0001640172 00000 n
+0001635766 00000 n
+0001631092 00000 n
+0001649867 00000 n
+0001650232 00000 n
+0001636010 00000 n
+0001645761 00000 n
+0001645457 00000 n
+0001650598 00000 n
+0001644970 00000 n
+0001645579 00000 n
+0001644909 00000 n
+0001645640 00000 n
+0001645152 00000 n
+0000527648 00000 n
+0000527483 00000 n
+0000527136 00000 n
+0000527595 00000 n
+0001790556 00000 n
+0000529944 00000 n
+0000530097 00000 n
+0000530251 00000 n
+0000532630 00000 n
+0000530544 00000 n
+0000529796 00000 n
+0000527733 00000 n
+0000530439 00000 n
+0001789043 00000 n
+0001649745 00000 n
+0000533146 00000 n
+0000532490 00000 n
+0000530655 00000 n
+0000533037 00000 n
+0000532834 00000 n
+0000535039 00000 n
+0000894767 00000 n
+0000895179 00000 n
+0000534907 00000 n
+0000533257 00000 n
+0000894956 00000 n
+0000895117 00000 n
+0000898043 00000 n
+0000964390 00000 n
+0000897931 00000 n
+0000895313 00000 n
+0000964275 00000 n
+0000964328 00000 n
+0000962259 00000 n
+0000962521 00000 n
+0000963699 00000 n
+0000963722 00000 n
+0000963910 00000 n
+0000967506 00000 n
+0000969831 00000 n
+0000978342 00000 n
+0000967828 00000 n
+0000967374 00000 n
+0000964529 00000 n
+0000967659 00000 n
+0001789385 00000 n
+0001789560 00000 n
+0000967768 00000 n
+0001038185 00000 n
+0000969719 00000 n
+0000967991 00000 n
+0001037958 00000 n
+0001038011 00000 n
+0001790673 00000 n
+0000970138 00000 n
+0000970318 00000 n
+0000978272 00000 n
+0000978292 00000 n
+0000978320 00000 n
+0000978649 00000 n
+0000978829 00000 n
+0001037887 00000 n
+0001037907 00000 n
+0001037935 00000 n
+0001039092 00000 n
+0001038927 00000 n
+0001038350 00000 n
+0001039039 00000 n
+0001041457 00000 n
+0001120216 00000 n
+0001233239 00000 n
+0001391759 00000 n
+0001041345 00000 n
+0001039177 00000 n
+0001391593 00000 n
+0001391698 00000 n
+0001395554 00000 n
+0001395751 00000 n
+0001395909 00000 n
+0001396067 00000 n
+0001396220 00000 n
+0001396375 00000 n
+0001396527 00000 n
+0001396682 00000 n
+0001396837 00000 n
+0001399132 00000 n
+0001397053 00000 n
+0001395358 00000 n
+0001391919 00000 n
+0001397000 00000 n
+0001788473 00000 n
+0001630726 00000 n
+0001406084 00000 n
+0001607561 00000 n
+0001607717 00000 n
+0001406316 00000 n
+0001399020 00000 n
+0001397189 00000 n
+0001406031 00000 n
+0001788756 00000 n
+0001406146 00000 n
+0001406203 00000 n
+0001403261 00000 n
+0001403472 00000 n
+0001403528 00000 n
+0001403922 00000 n
+0001403944 00000 n
+0001404225 00000 n
+0001607873 00000 n
+0001608029 00000 n
+0001608184 00000 n
+0001608336 00000 n
+0001608488 00000 n
+0001608641 00000 n
+0001408861 00000 n
+0001608793 00000 n
+0001608947 00000 n
+0001611292 00000 n
+0001609266 00000 n
+0001408657 00000 n
+0001406493 00000 n
+0001609101 00000 n
+0001409486 00000 n
+0001409719 00000 n
+0001409775 00000 n
+0001409824 00000 n
+0001409889 00000 n
+0001409961 00000 n
+0001607038 00000 n
+0001607059 00000 n
+0001607218 00000 n
+0001607378 00000 n
+0001650354 00000 n
+0001630970 00000 n
+0001636131 00000 n
+0001630909 00000 n
+0001616683 00000 n
+0001616837 00000 n
+0001617104 00000 n
+0001611152 00000 n
+0001609406 00000 n
+0001616989 00000 n
+0001617042 00000 n
+0001790790 00000 n
+0001616414 00000 n
+0001616614 00000 n
+0001616661 00000 n
+0001618522 00000 n
+0001618305 00000 n
+0001617244 00000 n
+0001618417 00000 n
+0001620706 00000 n
+0001620433 00000 n
+0001618633 00000 n
+0001620545 00000 n
+0001622949 00000 n
+0001622670 00000 n
+0001620804 00000 n
+0001622782 00000 n
+0001622835 00000 n
+0001788040 00000 n
+0001788899 00000 n
+0001622892 00000 n
+0001624682 00000 n
+0001624465 00000 n
+0001623124 00000 n
+0001624577 00000 n
+0001788615 00000 n
+0001788329 00000 n
+0001625783 00000 n
+0001626122 00000 n
+0001626404 00000 n
+0001625635 00000 n
+0001624856 00000 n
+0001626294 00000 n
+0001626347 00000 n
+0001625953 00000 n
+0001627075 00000 n
+0001626853 00000 n
+0001626502 00000 n
+0001626965 00000 n
+0001627018 00000 n
+0001790907 00000 n
+0001631214 00000 n
+0001630200 00000 n
+0001627160 00000 n
+0001630312 00000 n
+0001630365 00000 n
+0001790081 00000 n
+0001636192 00000 n
+0001634869 00000 n
+0001631312 00000 n
+0001634981 00000 n
+0001641145 00000 n
+0001639886 00000 n
+0001636290 00000 n
+0001639998 00000 n
+0001645881 00000 n
+0001644561 00000 n
+0001641243 00000 n
+0001644673 00000 n
+0001650842 00000 n
+0001649580 00000 n
+0001645979 00000 n
+0001649692 00000 n
+0001652538 00000 n
+0001652264 00000 n
+0001650940 00000 n
+0001652376 00000 n
+0001652481 00000 n
+0001791024 00000 n
+0001787429 00000 n
+0001652636 00000 n
+0001653532 00000 n
+0001653557 00000 n
+0001653582 00000 n
+0001653607 00000 n
+0001653632 00000 n
+0001654290 00000 n
+0001654417 00000 n
+0001654532 00000 n
+0001654919 00000 n
+0001655150 00000 n
+0001655345 00000 n
+0001655938 00000 n
+0001656403 00000 n
+0001657377 00000 n
+0001657854 00000 n
+0001660381 00000 n
+0001660603 00000 n
+0001664911 00000 n
+0001665187 00000 n
+0001666917 00000 n
+0001667136 00000 n
+0001669673 00000 n
+0001669936 00000 n
+0001671669 00000 n
+0001671889 00000 n
+0001673693 00000 n
+0001673965 00000 n
+0001675347 00000 n
+0001675573 00000 n
+0001685518 00000 n
+0001685804 00000 n
+0001703801 00000 n
+0001704314 00000 n
+0001715499 00000 n
+0001715876 00000 n
+0001721152 00000 n
+0001721426 00000 n
+0001737598 00000 n
+0001738045 00000 n
+0001758643 00000 n
+0001759318 00000 n
+0001770467 00000 n
+0001770784 00000 n
+0001786979 00000 n
+0001791101 00000 n
+0001791202 00000 n
+0001792958 00000 n
+0001793149 00000 n
+0001793331 00000 n
+0001793545 00000 n
+0001793785 00000 n
+0001793999 00000 n
+0001794203 00000 n
+0001794435 00000 n
+0001794656 00000 n
+0001794883 00000 n
+0001795094 00000 n
+0001795322 00000 n
+0001795540 00000 n
+0001795746 00000 n
+0001795972 00000 n
+0001796178 00000 n
+0001796404 00000 n
+0001796626 00000 n
+0001796809 00000 n
+0001796979 00000 n
+0001797148 00000 n
+0001797318 00000 n
+0001797483 00000 n
+0001797662 00000 n
+0001797857 00000 n
+0001798053 00000 n
+0001798186 00000 n
+0001798301 00000 n
+0001798426 00000 n
+0001798539 00000 n
+0001798652 00000 n
+0001798737 00000 n
+0001798844 00000 n
+0001798882 00000 n
+0001799092 00000 n
trailer
<< /Size 857
/Root 855 0 R
/Info 856 0 R
-/ID [<1D684A5DEC0D8F7BD54A2644EE6D88D8> <1D684A5DEC0D8F7BD54A2644EE6D88D8>] >>
+/ID [<5014EB306AC0929EE8491BF934324FC4> <5014EB306AC0929EE8491BF934324FC4>] >>
startxref
-1799047
+1799423
%%EOF
Modified: seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.tex
===================================================================
--- seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.tex 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.tex 2012-09-03 02:09:56 UTC (rev 20673)
@@ -621,6 +621,9 @@
\item[source\_type] Set this value equal to \texttt{1} for elastic forces or acoustic pressure,
set this to \texttt{2} for moment tensor sources.
+For a plane wave including converted and reflected waves at the free surface, P wave = 1, S wave = 2, Rayleigh wave = 3;
+for a plane wave without converted nor reflected waves at the free surface, i.e. the incident wave only, P wave = 4, S wave = 5.
+(incident plane waves are turned on by parameter \texttt{initialfield} in \texttt{DATA/Par\_file}).
\item[time\_function\_type] Choose a source-time function: set this value to \texttt{1} to use a Ricker,
\texttt{2} the first derivative, \texttt{3} a Gaussian, \texttt{4} a Dirac or \texttt{5} a Heaviside source-time function.
@@ -668,7 +671,7 @@
an overall time shift to the synthetics (according to the time shift of the first source), something that can be done
in the post-processing. This time shift parameter can be non-zero when using multiple sources.
-\item[angleforce] Angle of the source (for a force only)
+\item[anglesource] angle of the source (for a force only); for a plane wave, this is the incidence angle. For moment tensor sources this parameter is unused.
\item[Mxx,Mzz,Mxz] Moment tensor components (valid only for moment tensor sources, source\_type "2").
Note that the units for the components of a moment tensor source are different in SPECFEM2D and in SPECFEM3D:
Modified: seismo/2D/SPECFEM2D/trunk/setup/constants.h.in
===================================================================
--- seismo/2D/SPECFEM2D/trunk/setup/constants.h.in 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/setup/constants.h.in 2012-09-03 02:09:56 UTC (rev 20673)
@@ -99,7 +99,7 @@
! integer, parameter :: IOUT = 41
! number of lines per source in SOURCE file
- integer, parameter :: NLINES_PER_SOURCE = 13
+ integer, parameter :: NLINES_PER_SOURCE = 12
! flags for absorbing boundaries
integer, parameter :: IBOTTOM = 1
Modified: seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_source_file.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_source_file.f90 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_source_file.f90 2012-09-03 02:09:56 UTC (rev 20673)
@@ -48,7 +48,7 @@
! source parameters
integer, dimension(:),pointer :: source_type,time_function_type
- double precision, dimension(:),pointer :: xs,zs,f0,tshift_src,angleforce, &
+ double precision, dimension(:),pointer :: xs,zs,f0,tshift_src,anglesource, &
Mxx,Mzz,Mxz,factor
logical, dimension(:),pointer :: source_surf
@@ -65,7 +65,7 @@
! local parameters
integer :: ios,icounter,i_source,num_sources
- character(len=150) dummystring
+ character(len=150) string_read
integer, parameter :: IIN_SOURCE = 22
! allocates memory arrays
@@ -76,7 +76,7 @@
allocate(time_function_type(NSOURCES))
allocate(f0(NSOURCES))
allocate(tshift_src(NSOURCES))
- allocate(angleforce(NSOURCES))
+ allocate(anglesource(NSOURCES))
allocate(Mxx(NSOURCES))
allocate(Mxz(NSOURCES))
allocate(Mzz(NSOURCES))
@@ -88,14 +88,28 @@
icounter = 0
do while(ios == 0)
- read(IIN_SOURCE,"(a)",iostat=ios) dummystring
- if(ios == 0) icounter = icounter + 1
+ read(IIN_SOURCE,"(a)",iostat=ios) string_read
+
+ if(ios == 0) then
+
+! suppress trailing carriage return (ASCII code 13) if any (e.g. if input text file coming from Windows/DOS)
+ if(index(string_read,achar(13)) > 0) string_read = string_read(1:index(string_read,achar(13))-1)
+
+! suppress leading and trailing white spaces, if any
+ string_read = adjustl(string_read)
+ string_read = string_read(1:len_trim(string_read))
+
+! if the line is not empty and is not a comment, count it
+ if(len_trim(string_read) > 0 .and. (index(string_read,'#') == 0 .or. index(string_read,'#') > 1)) icounter = icounter + 1
+
+ endif
+
enddo
close(IIN_SOURCE)
! checks counter
if(mod(icounter,NLINES_PER_SOURCE) /= 0) &
- stop 'total number of lines in SOURCE file should be a multiple of NLINES_PER_SOURCE'
+ stop 'total number of non blank and non comment lines in SOURCE file should be a multiple of NLINES_PER_SOURCE'
! total number of sources
num_sources = icounter / NLINES_PER_SOURCE
@@ -114,7 +128,7 @@
call read_value_integer(IIN_SOURCE,IGNORE_JUNK,time_function_type(i_source))
call read_value_double_precision(IIN_SOURCE,IGNORE_JUNK,f0(i_source))
call read_value_double_precision(IIN_SOURCE,IGNORE_JUNK,tshift_src(i_source))
- call read_value_double_precision(IIN_SOURCE,IGNORE_JUNK,angleforce(i_source))
+ call read_value_double_precision(IIN_SOURCE,IGNORE_JUNK,anglesource(i_source))
call read_value_double_precision(IIN_SOURCE,IGNORE_JUNK,Mxx(i_source))
call read_value_double_precision(IIN_SOURCE,IGNORE_JUNK,Mzz(i_source))
call read_value_double_precision(IIN_SOURCE,IGNORE_JUNK,Mxz(i_source))
@@ -129,7 +143,7 @@
print *,'Frequency, delay = ',f0(i_source),tshift_src(i_source)
print *,'Source type (1=force, 2=explosion): ',source_type(i_source)
print *,'Time function type (1=Ricker, 2=First derivative, 3=Gaussian, 4=Dirac, 5=Heaviside): ',time_function_type(i_source)
- print *,'Angle of the source if force = ',angleforce(i_source)
+ print *,'Angle of the source if force = ',anglesource(i_source)
print *,'Mxx of the source if moment tensor = ',Mxx(i_source)
print *,'Mzz of the source if moment tensor = ',Mzz(i_source)
print *,'Mxz of the source if moment tensor = ',Mxz(i_source)
Modified: seismo/2D/SPECFEM2D/trunk/src/meshfem2D/save_databases.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/meshfem2D/save_databases.f90 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/src/meshfem2D/save_databases.f90 2012-09-03 02:09:56 UTC (rev 20673)
@@ -225,7 +225,7 @@
write(15,*) 'source', i_source
write(15,*) source_type(i_source),time_function_type(i_source), &
xs(i_source),zs(i_source),f0(i_source),tshift_src(i_source), &
- factor(i_source),angleforce(i_source), &
+ factor(i_source),anglesource(i_source), &
Mxx(i_source),Mzz(i_source),Mxz(i_source)
enddo
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_Bielak_conditions.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_Bielak_conditions.f90 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_Bielak_conditions.f90 2012-09-03 02:09:56 UTC (rev 20673)
@@ -45,7 +45,7 @@
! compute analytical initial plane wave for Bielak's conditions
subroutine compute_Bielak_conditions(coord,iglob,nglob,it,deltat,dxUx,dxUz,dzUx,dzUz,veloc_horiz,veloc_vert, &
- x0_source, z0_source, A_plane, B_plane, C_plane, angleforce, angleforce_refl, &
+ x0_source, z0_source, A_plane, B_plane, C_plane, anglesource, anglesource_refl, &
c_inc, c_refl, time_offset,f0)
implicit none
@@ -64,7 +64,7 @@
double precision, external :: ricker_Bielak_veloc
- double precision x0_source, z0_source, angleforce, angleforce_refl
+ double precision x0_source, z0_source, anglesource, anglesource_refl
double precision c_inc, c_refl, time_offset, f0
double precision, dimension(NDIM) :: A_plane, B_plane, C_plane
@@ -124,25 +124,25 @@
! to compute the derivative of the displacement, we take the velocity ricker expression and we multiply by
! the derivative of the interior argument of ricker_Bielak_veloc
- dxUx = A_plane(1) * ricker_Bielak_veloc(t - sin(angleforce)*x/c_inc + cos(angleforce)*z/c_inc,f0) * (-sin(angleforce)/c_inc)&
- + B_plane(1) * ricker_Bielak_veloc(t - sin(angleforce)*x/c_inc - cos(angleforce)*z/c_inc,f0) * (-sin(angleforce)/c_inc)&
- + C_plane(1) * ricker_Bielak_veloc(t - sin(angleforce_refl)*x/c_refl - cos(angleforce_refl)*z/c_refl,f0)&
- * (-sin(angleforce_refl)/c_refl)
+ dxUx = A_plane(1) * ricker_Bielak_veloc(t - sin(anglesource)*x/c_inc + cos(anglesource)*z/c_inc,f0) * (-sin(anglesource)/c_inc)&
+ + B_plane(1) * ricker_Bielak_veloc(t - sin(anglesource)*x/c_inc - cos(anglesource)*z/c_inc,f0) * (-sin(anglesource)/c_inc)&
+ + C_plane(1) * ricker_Bielak_veloc(t - sin(anglesource_refl)*x/c_refl - cos(anglesource_refl)*z/c_refl,f0)&
+ * (-sin(anglesource_refl)/c_refl)
- dzUx = A_plane(1) * ricker_Bielak_veloc(t - sin(angleforce)*x/c_inc + cos(angleforce)*z/c_inc,f0) * (-cos(angleforce)/c_inc)&
- + B_plane(1) * ricker_Bielak_veloc(t - sin(angleforce)*x/c_inc - cos(angleforce)*z/c_inc,f0) * (cos(angleforce)/c_inc)&
- + C_plane(1) * ricker_Bielak_veloc(t - sin(angleforce_refl)*x/c_refl - cos(angleforce_refl)*z/c_refl,f0)&
- * (cos(angleforce_refl)/c_refl)
+ dzUx = A_plane(1) * ricker_Bielak_veloc(t - sin(anglesource)*x/c_inc + cos(anglesource)*z/c_inc,f0) * (-cos(anglesource)/c_inc)&
+ + B_plane(1) * ricker_Bielak_veloc(t - sin(anglesource)*x/c_inc - cos(anglesource)*z/c_inc,f0) * (cos(anglesource)/c_inc)&
+ + C_plane(1) * ricker_Bielak_veloc(t - sin(anglesource_refl)*x/c_refl - cos(anglesource_refl)*z/c_refl,f0)&
+ * (cos(anglesource_refl)/c_refl)
- dxUz = A_plane(2) * ricker_Bielak_veloc(t - sin(angleforce)*x/c_inc + cos(angleforce)*z/c_inc,f0) * (-sin(angleforce)/c_inc)&
- + B_plane(2) * ricker_Bielak_veloc(t - sin(angleforce)*x/c_inc - cos(angleforce)*z/c_inc,f0) * (-sin(angleforce)/c_inc)&
- + C_plane(2) * ricker_Bielak_veloc(t - sin(angleforce_refl)*x/c_refl - cos(angleforce_refl)*z/c_refl,f0)&
- * (-sin(angleforce_refl)/c_refl)
+ dxUz = A_plane(2) * ricker_Bielak_veloc(t - sin(anglesource)*x/c_inc + cos(anglesource)*z/c_inc,f0) * (-sin(anglesource)/c_inc)&
+ + B_plane(2) * ricker_Bielak_veloc(t - sin(anglesource)*x/c_inc - cos(anglesource)*z/c_inc,f0) * (-sin(anglesource)/c_inc)&
+ + C_plane(2) * ricker_Bielak_veloc(t - sin(anglesource_refl)*x/c_refl - cos(anglesource_refl)*z/c_refl,f0)&
+ * (-sin(anglesource_refl)/c_refl)
- dzUz = A_plane(2) * ricker_Bielak_veloc(t - sin(angleforce)*x/c_inc + cos(angleforce)*z/c_inc,f0) * (-cos(angleforce)/c_inc)&
- + B_plane(2) * ricker_Bielak_veloc(t - sin(angleforce)*x/c_inc - cos(angleforce)*z/c_inc,f0) * (cos(angleforce)/c_inc)&
- + C_plane(2) * ricker_Bielak_veloc(t - sin(angleforce_refl)*x/c_refl - cos(angleforce_refl)*z/c_refl,f0)&
- * (cos(angleforce_refl)/c_refl)
+ dzUz = A_plane(2) * ricker_Bielak_veloc(t - sin(anglesource)*x/c_inc + cos(anglesource)*z/c_inc,f0) * (-cos(anglesource)/c_inc)&
+ + B_plane(2) * ricker_Bielak_veloc(t - sin(anglesource)*x/c_inc - cos(anglesource)*z/c_inc,f0) * (cos(anglesource)/c_inc)&
+ + C_plane(2) * ricker_Bielak_veloc(t - sin(anglesource_refl)*x/c_refl - cos(anglesource_refl)*z/c_refl,f0)&
+ * (cos(anglesource_refl)/c_refl)
t = time_veloc
@@ -154,12 +154,12 @@
!!$ veloc_vert = - HALF * ricker_Bielak_veloc(t - x/2.d0 + (9 - z) * (sqrt(3.d0)/2.d0)) &
!!$ + HALF * ricker_Bielak_veloc(t - x/2.d0 - (9 - z) * (sqrt(3.d0)/2.d0))
- veloc_horiz = A_plane(1) * ricker_Bielak_veloc(t - sin(angleforce)*x/c_inc + cos(angleforce)*z/c_inc,f0) &
- + B_plane(1) * ricker_Bielak_veloc(t - sin(angleforce)*x/c_inc - cos(angleforce)*z/c_inc,f0) &
- + C_plane(1) * ricker_Bielak_veloc(t - sin(angleforce_refl)*x/c_refl - cos(angleforce_refl)*z/c_refl,f0)
- veloc_vert = A_plane(2) * ricker_Bielak_veloc(t - sin(angleforce)*x/c_inc + cos(angleforce)*z/c_inc,f0) &
- + B_plane(2) * ricker_Bielak_veloc(t - sin(angleforce)*x/c_inc - cos(angleforce)*z/c_inc,f0) &
- + C_plane(2) * ricker_Bielak_veloc(t - sin(angleforce_refl)*x/c_refl - cos(angleforce_refl)*z/c_refl,f0)
+ veloc_horiz = A_plane(1) * ricker_Bielak_veloc(t - sin(anglesource)*x/c_inc + cos(anglesource)*z/c_inc,f0) &
+ + B_plane(1) * ricker_Bielak_veloc(t - sin(anglesource)*x/c_inc - cos(anglesource)*z/c_inc,f0) &
+ + C_plane(1) * ricker_Bielak_veloc(t - sin(anglesource_refl)*x/c_refl - cos(anglesource_refl)*z/c_refl,f0)
+ veloc_vert = A_plane(2) * ricker_Bielak_veloc(t - sin(anglesource)*x/c_inc + cos(anglesource)*z/c_inc,f0) &
+ + B_plane(2) * ricker_Bielak_veloc(t - sin(anglesource)*x/c_inc - cos(anglesource)*z/c_inc,f0) &
+ + C_plane(2) * ricker_Bielak_veloc(t - sin(anglesource_refl)*x/c_refl - cos(anglesource_refl)*z/c_refl,f0)
end subroutine compute_Bielak_conditions
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90 2012-09-03 02:09:56 UTC (rev 20673)
@@ -45,7 +45,7 @@
subroutine compute_forces_viscoelastic(p_sv,nglob,nspec,myrank,nelemabs,numat, &
ispec_selected_source,ispec_selected_rec,is_proc_source,which_proc_receiver, &
source_type,it,NSTEP,anyabs,assign_external_model, &
- initialfield,ATTENUATION_VISCOELASTIC_SOLID,angleforce,deltatcube, &
+ initialfield,ATTENUATION_VISCOELASTIC_SOLID,anglesource,deltatcube, &
deltatfourth,twelvedeltat,fourdeltatsquare,ibool,kmato,numabs,elastic,codeabs, &
accel_elastic,veloc_elastic,displ_elastic,b_accel_elastic,b_displ_elastic, &
density,poroelastcoef,xix,xiz,gammax,gammaz, &
@@ -55,7 +55,7 @@
dux_dxl_np1,duz_dzl_np1,duz_dxl_np1,dux_dzl_np1,hprime_xx,hprimewgll_xx, &
hprime_zz,hprimewgll_zz,wxgll,wzgll,inv_tau_sigma_nu1,phi_nu1,inv_tau_sigma_nu2,phi_nu2,Mu_nu1,Mu_nu2,N_SLS, &
deltat,coord,add_Bielak_conditions, &
- x0_source, z0_source, A_plane, B_plane, C_plane, angleforce_refl, c_inc, c_refl, time_offset,f0, &
+ x0_source, z0_source, A_plane, B_plane, C_plane, anglesource_refl, c_inc, c_refl, time_offset,f0, &
v0x_left,v0z_left,v0x_right,v0z_right,v0x_bot,v0z_bot,t0x_left,t0z_left,t0x_right,t0z_right,t0x_bot,t0z_bot,&
nleft,nright,nbot,over_critical_angle,NSOURCES,nrec,SIMULATION_TYPE,SAVE_FORWARD,b_absorb_elastic_left,&
b_absorb_elastic_right,b_absorb_elastic_bottom,b_absorb_elastic_top,nspec_left,nspec_right,&
@@ -95,7 +95,7 @@
logical :: SAVE_FORWARD
double precision :: deltatcube,deltatfourth,twelvedeltat,fourdeltatsquare
- double precision, dimension(NSOURCES) :: angleforce
+ double precision, dimension(NSOURCES) :: anglesource
integer, dimension(NGLLX,NGLLZ,nspec) :: ibool
integer, dimension(nspec) :: kmato
@@ -201,7 +201,7 @@
! for analytical initial plane wave for Bielak's conditions
double precision :: veloc_horiz,veloc_vert,dxUx,dzUx,dxUz,dzUz,traction_x_t0,traction_z_t0,deltat
double precision, dimension(NDIM,nglob), intent(in) :: coord
- double precision x0_source, z0_source, angleforce_refl, c_inc, c_refl, time_offset, f0
+ double precision x0_source, z0_source, anglesource_refl, c_inc, c_refl, time_offset, f0
double precision, dimension(NDIM) :: A_plane, B_plane, C_plane
!over critical angle
logical :: over_critical_angle
@@ -1025,9 +1025,9 @@
cpl = sqrt((kappal + 4._CUSTOM_REAL*mul_unrelaxed_elastic/3._CUSTOM_REAL)/rhol)
csl = sqrt(mul_unrelaxed_elastic/rhol)
-!!! DK DK
+!!! DK DK
c_inc = csl
-!!! DK DK
+!!! DK DK
!--- left absorbing boundary
if(codeabs(IEDGE4,ispecabs)) then
@@ -1043,7 +1043,7 @@
if(add_Bielak_conditions .and. initialfield) then
if (.not.over_critical_angle) then
call compute_Bielak_conditions(coord,iglob,nglob,it,deltat,dxUx,dxUz,dzUx,dzUz,veloc_horiz,veloc_vert, &
- x0_source, z0_source, A_plane, B_plane, C_plane, angleforce(1), angleforce_refl, &
+ x0_source, z0_source, A_plane, B_plane, C_plane, anglesource(1), anglesource_refl, &
c_inc, c_refl, time_offset,f0)
traction_x_t0 = (lambdal_unrelaxed_elastic+2*mul_unrelaxed_elastic)*dxUx + lambdal_unrelaxed_elastic*dzUz
traction_z_t0 = mul_unrelaxed_elastic*(dxUz + dzUx)
@@ -1158,7 +1158,7 @@
if(add_Bielak_conditions .and. initialfield) then
if (.not.over_critical_angle) then
call compute_Bielak_conditions(coord,iglob,nglob,it,deltat,dxUx,dxUz,dzUx,dzUz,veloc_horiz,veloc_vert, &
- x0_source, z0_source, A_plane, B_plane, C_plane, angleforce(1), angleforce_refl, &
+ x0_source, z0_source, A_plane, B_plane, C_plane, anglesource(1), anglesource_refl, &
c_inc, c_refl, time_offset,f0)
traction_x_t0 = (lambdal_unrelaxed_elastic+2*mul_unrelaxed_elastic)*dxUx + lambdal_unrelaxed_elastic*dzUz
traction_z_t0 = mul_unrelaxed_elastic*(dxUz + dzUx)
@@ -1279,7 +1279,7 @@
if(add_Bielak_conditions .and. initialfield) then
if (.not.over_critical_angle) then
call compute_Bielak_conditions(coord,iglob,nglob,it,deltat,dxUx,dxUz,dzUx,dzUz,veloc_horiz,veloc_vert, &
- x0_source, z0_source, A_plane, B_plane, C_plane, angleforce(1), angleforce_refl, &
+ x0_source, z0_source, A_plane, B_plane, C_plane, anglesource(1), anglesource_refl, &
c_inc, c_refl, time_offset,f0)
traction_x_t0 = mul_unrelaxed_elastic*(dxUz + dzUx)
traction_z_t0 = lambdal_unrelaxed_elastic*dxUx + (lambdal_unrelaxed_elastic+2*mul_unrelaxed_elastic)*dzUz
@@ -1416,7 +1416,7 @@
! top or bottom edge, vertical normal vector
if(add_Bielak_conditions .and. initialfield) then
call compute_Bielak_conditions(coord,iglob,nglob,it,deltat,dxUx,dxUz,dzUx,dzUz,veloc_horiz,veloc_vert, &
- x0_source, z0_source, A_plane, B_plane, C_plane, angleforce(1), angleforce_refl, &
+ x0_source, z0_source, A_plane, B_plane, C_plane, anglesource(1), anglesource_refl, &
c_inc, c_refl, time_offset,f0)
traction_x_t0 = mul_unrelaxed_elastic*(dxUz + dzUx)
traction_z_t0 = lambdal_unrelaxed_elastic*dxUx + (lambdal_unrelaxed_elastic+2*mul_unrelaxed_elastic)*dzUz
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/paco_beyond_critical.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/paco_beyond_critical.f90 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/paco_beyond_critical.f90 2012-09-03 02:09:56 UTC (rev 20673)
@@ -14,7 +14,7 @@
! modified by Dimitri Komatitsch and Ronan Madec in March 2008
! in particular, converted to Fortran90 and to double precision
-subroutine paco_beyond_critical(coord,nglob,deltat,NSTEP_global,angleforce,&
+subroutine paco_beyond_critical(coord,nglob,deltat,NSTEP_global,anglesource,&
f0,cp_local,cs_local,INCLUDE_ATTENUATION,QD,source_type,v0x_left,v0z_left,v0x_right,v0z_right,&
v0x_bot,v0z_bot,t0x_left,t0z_left,t0x_right,t0z_right,t0x_bot,t0z_bot,left_bound,right_bound,&
bot_bound,nleft,nright,nbot,displ_elastic,veloc_elastic,accel_elastic,x_source)
@@ -23,7 +23,7 @@
include "constants.h"
- double precision :: f0,cp_local,cs_local,deltat,dt,TP,angleforce,QD,delta_in_period
+ double precision :: f0,cp_local,cs_local,deltat,dt,TP,anglesource,QD,delta_in_period
logical :: INCLUDE_ATTENUATION
integer :: npt,NSTEP_global,source_type,nleft,nright,nbot,nglob
@@ -126,7 +126,7 @@
UR=(1.0d0, 0.0d0)
! convert angle to radians
- GAMR = angleforce
+ GAMR = anglesource
BEALF=SQRT((1.0d0-2.0d0*ANU)/(2.0d0*(1.0d0-ANU)))
ALFBE=1.0d0/BEALF
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/prepare_initialfield.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/prepare_initialfield.F90 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/prepare_initialfield.F90 2012-09-03 02:09:56 UTC (rev 20673)
@@ -45,9 +45,9 @@
subroutine prepare_initialfield(myrank,any_acoustic,any_poroelastic,over_critical_angle, &
- NSOURCES,source_type,angleforce,x_source,z_source,f0,t0, &
+ NSOURCES,source_type,anglesource,x_source,z_source,f0,t0, &
nglob,numat,poroelastcoef,density,coord, &
- angleforce_refl,c_inc,c_refl,cploc,csloc,time_offset, &
+ anglesource_refl,c_inc,c_refl,cploc,csloc,time_offset, &
A_plane, B_plane, C_plane, &
accel_elastic,veloc_elastic,displ_elastic)
@@ -62,7 +62,7 @@
integer :: NSOURCES
integer, dimension(NSOURCES) :: source_type
- double precision, dimension(NSOURCES) :: angleforce,x_source,z_source,f0
+ double precision, dimension(NSOURCES) :: anglesource,x_source,z_source,f0
double precision :: t0
integer :: nglob,numat
@@ -70,7 +70,7 @@
double precision, dimension(2,numat) :: density
double precision, dimension(NDIM,nglob) :: coord
- double precision :: angleforce_abs, angleforce_refl,c_inc,c_refl,cploc,csloc
+ double precision :: anglesource_abs, anglesource_refl,c_inc,c_refl,cploc,csloc
double precision :: time_offset,x0_source,z0_source
double precision, dimension(2) :: A_plane, B_plane, C_plane
@@ -114,26 +114,26 @@
write(IOUT,*) '*** calculation of the initial plane wave ***'
write(IOUT,*)
write(IOUT,*) 'To change the initial plane wave, change source_type in DATA/SOURCE'
- write(IOUT,*) 'and use 1 for a plane P wave, 2 for a plane SV wave, 3 for a Rayleigh wave'
+ write(IOUT,*) 'and use 1 or 4 for a plane P wave, 2 or 5 for a plane SV wave, 3 for a Rayleigh wave'
write(IOUT,*)
! only implemented for one source
if(NSOURCES > 1) call exit_MPI('calculation of the initial wave is only implemented for one source')
- if (source_type(1) == 1) then
- write(IOUT,*) 'initial P wave of', angleforce(1)*180.d0/pi, 'degrees introduced.'
- else if (source_type(1) == 2) then
- write(IOUT,*) 'initial SV wave of', angleforce(1)*180.d0/pi, ' degrees introduced.'
+ if (source_type(1) == 1 .or. source_type(1) == 4) then
+ write(IOUT,*) 'initial P wave of', anglesource(1)*180.d0/pi, 'degrees introduced.'
+ else if (source_type(1) == 2 .or. source_type(1) == 5) then
+ write(IOUT,*) 'initial SV wave of', anglesource(1)*180.d0/pi, ' degrees introduced.'
else if (source_type(1) == 3) then
write(IOUT,*) 'Rayleigh wave introduced.'
else
- call exit_MPI('Unrecognized source_type: should be 1 for plane P waves, 2 for plane SV waves, 3 for Rayleigh wave')
+ call exit_MPI('Unrecognized source_type: should be 1 or 4 for plane P waves, 2 or 5 for plane SV waves, 3 for Rayleigh wave')
endif
endif
- ! allow negative angleforce(1): incidence from the right side of the domain
- angleforce_abs=abs(angleforce(1))
- if (angleforce_abs > pi/2.d0 .and. source_type(1) /= 3) &
- call exit_MPI("incorrect angleforce: must have 0 <= angleforce < 90")
+ ! allow negative anglesource(1): incidence from the right side of the domain
+ anglesource_abs=abs(anglesource(1))
+ if (anglesource_abs > pi/2.d0 .and. source_type(1) /= 3) &
+ call exit_MPI("incorrect anglesource: must have 0 <= anglesource < 90")
! only implemented for homogeneous media therefore only 1 material supported
numat_local = numat
@@ -155,75 +155,75 @@
csloc = sqrt(mu/denst)
! P wave case
- if (source_type(1) == 1) then
+ if (source_type(1) == 1 .or. source_type(1) == 4) then
- p=sin(angleforce_abs)/cploc
+ p=sin(anglesource_abs)/cploc
c_inc = cploc
c_refl = csloc
- angleforce_refl = asin(p*c_refl)
+ anglesource_refl = asin(p*c_refl)
! from formulas (5.27) and (5.28) p 134 in Aki & Richards (2002)
- PP = (- cos(2.d0*angleforce_refl)**2/csloc**3 &
- + 4.d0*p**2*cos(angleforce_abs)*cos(angleforce_refl)/cploc) / &
- ( cos(2.d0*angleforce_refl)**2/csloc**3 &
- + 4.d0*p**2*cos(angleforce_abs)*cos(angleforce_refl)/cploc)
+ PP = (- cos(2.d0*anglesource_refl)**2/csloc**3 &
+ + 4.d0*p**2*cos(anglesource_abs)*cos(anglesource_refl)/cploc) / &
+ ( cos(2.d0*anglesource_refl)**2/csloc**3 &
+ + 4.d0*p**2*cos(anglesource_abs)*cos(anglesource_refl)/cploc)
- PS = 4.d0*p*cos(angleforce_abs)*cos(2.d0*angleforce_refl) / &
- (csloc**2*(cos(2.d0*angleforce_refl)**2/csloc**3 &
- +4.d0*p**2*cos(angleforce_abs)*cos(angleforce_refl)/cploc))
+ PS = 4.d0*p*cos(anglesource_abs)*cos(2.d0*anglesource_refl) / &
+ (csloc**2*(cos(2.d0*anglesource_refl)**2/csloc**3 &
+ +4.d0*p**2*cos(anglesource_abs)*cos(anglesource_refl)/cploc))
if (myrank == 0) then
- write(IOUT,*) 'reflected convert plane wave angle: ', angleforce_refl*180.d0/pi
+ write(IOUT,*) 'reflected convert plane wave angle: ', anglesource_refl*180.d0/pi
endif
! from Table 5.1 p141 in Aki & Richards (1980)
! we put the opposite sign on z coefficients because z axis is oriented from bottom to top
- A_plane(1) = sin(angleforce_abs); A_plane(2) = cos(angleforce_abs)
- B_plane(1) = PP * sin(angleforce_abs); B_plane(2) = - PP * cos(angleforce_abs)
- C_plane(1) = PS * cos(angleforce_refl); C_plane(2) = PS * sin(angleforce_refl)
+ A_plane(1) = sin(anglesource_abs); A_plane(2) = cos(anglesource_abs)
+ B_plane(1) = PP * sin(anglesource_abs); B_plane(2) = - PP * cos(anglesource_abs)
+ C_plane(1) = PS * cos(anglesource_refl); C_plane(2) = PS * sin(anglesource_refl)
! SV wave case
- else if (source_type(1) == 2) then
+ else if (source_type(1) == 2 .or. source_type(1) == 5) then
- p=sin(angleforce_abs)/csloc
+ p=sin(anglesource_abs)/csloc
c_inc = csloc
c_refl = cploc
! if this coefficient is greater than 1, we are beyond the critical SV wave angle and there cannot be a converted P wave
if (p*c_refl<=1.d0) then
- angleforce_refl = asin(p*c_refl)
+ anglesource_refl = asin(p*c_refl)
! from formulas (5.30) and (5.31) p 140 in Aki & Richards (1980)
- SS = (cos(2.d0*angleforce_abs)**2/csloc**3 &
- - 4.d0*p**2*cos(angleforce_abs)*cos(angleforce_refl)/cploc) / &
- (cos(2.d0*angleforce_abs)**2/csloc**3 &
- + 4.d0*p**2*cos(angleforce_abs)*cos(angleforce_refl)/cploc)
- SP = 4.d0*p*cos(angleforce_abs)*cos(2*angleforce_abs) / &
- (cploc*csloc*(cos(2.d0*angleforce_abs)**2/csloc**3&
- +4.d0*p**2*cos(angleforce_refl)*cos(angleforce_abs)/cploc))
+ SS = (cos(2.d0*anglesource_abs)**2/csloc**3 &
+ - 4.d0*p**2*cos(anglesource_abs)*cos(anglesource_refl)/cploc) / &
+ (cos(2.d0*anglesource_abs)**2/csloc**3 &
+ + 4.d0*p**2*cos(anglesource_abs)*cos(anglesource_refl)/cploc)
+ SP = 4.d0*p*cos(anglesource_abs)*cos(2*anglesource_abs) / &
+ (cploc*csloc*(cos(2.d0*anglesource_abs)**2/csloc**3&
+ +4.d0*p**2*cos(anglesource_refl)*cos(anglesource_abs)/cploc))
if (myrank == 0) then
- write(IOUT,*) 'reflected convert plane wave angle: ', angleforce_refl*180.d0/pi
+ write(IOUT,*) 'reflected convert plane wave angle: ', anglesource_refl*180.d0/pi
endif
! SV45 degree incident plane wave is a particular case
- else if (angleforce_abs>pi/4.d0-1.0d-11 .and. angleforce_abs<pi/4.d0+1.0d-11) then
- angleforce_refl = 0.d0
+ else if (anglesource_abs>pi/4.d0-1.0d-11 .and. anglesource_abs<pi/4.d0+1.0d-11) then
+ anglesource_refl = 0.d0
SS = -1.0d0
SP = 0.d0
else
over_critical_angle=.true.
- angleforce_refl = 0.d0
+ anglesource_refl = 0.d0
SS = 0.0d0
SP = 0.d0
endif
! from Table 5.1 p141 in Aki & Richards (1980)
! we put the opposite sign on z coefficients because z axis is oriented from bottom to top
- A_plane(1) = cos(angleforce_abs); A_plane(2) = - sin(angleforce_abs)
- B_plane(1) = SS * cos(angleforce_abs); B_plane(2) = SS * sin(angleforce_abs)
- C_plane(1) = SP * sin(angleforce_refl); C_plane(2) = - SP * cos(angleforce_refl)
+ A_plane(1) = cos(anglesource_abs); A_plane(2) = - sin(anglesource_abs)
+ B_plane(1) = SS * cos(anglesource_abs); B_plane(2) = SS * sin(anglesource_abs)
+ C_plane(1) = SP * sin(anglesource_refl); C_plane(2) = - SP * cos(anglesource_refl)
! Rayleigh case
else if (source_type(1) == 3) then
@@ -234,11 +234,17 @@
endif
! correct A_plane and B_plane according to incident direction
- if (angleforce(1) < 0.) then
+ if (anglesource(1) < 0.) then
A_plane(1)=-A_plane(1); B_plane(1)=-B_plane(1)
C_plane(1)=-C_plane(1)
endif
+! to suppress the reflected and converted plane wave fields
+ if(source_type(1) == 4 .or. source_type(1) == 5) then
+ B_plane(:) = 0
+ C_plane(:) = 0
+ endif
+
! get minimum and maximum values of mesh coordinates
xmin = minval(coord(1,:))
zmin = minval(coord(2,:))
@@ -263,7 +269,7 @@
! initialize the time offset to put the plane wave not too close to the irregularity on the free surface
! add -t0 to match with the actual traveltime of plane waves
- if (abs(angleforce(1))<1.d0*pi/180.d0 .and. source_type(1)/=3) then
+ if (abs(anglesource(1))<1.d0*pi/180.d0 .and. source_type(1)/=3) then
time_offset = -1.d0*(zmax-zmin)/2.d0/c_inc - t0
else
time_offset = 0.d0 - t0
@@ -290,7 +296,7 @@
! z is from bottom to top therefore we take -z to make parallel with Aki & Richards
z = z0_source - z
- if (angleforce(1) >= 0.) then
+ if (anglesource(1) >= 0.) then
x = x - x0_source
else
x = x0_source -x
@@ -300,33 +306,33 @@
! formulas for the initial displacement for a plane wave from Aki & Richards (1980)
displ_elastic(1,i) = &
- A_plane(1) * ricker_Bielak_displ(t - sin(angleforce_abs)*x/c_inc + cos(angleforce_abs)*z/c_inc,f0(1)) &
- + B_plane(1) * ricker_Bielak_displ(t - sin(angleforce_abs)*x/c_inc - cos(angleforce_abs)*z/c_inc,f0(1)) &
- + C_plane(1) * ricker_Bielak_displ(t - sin(angleforce_refl)*x/c_refl - cos(angleforce_refl)*z/c_refl,f0(1))
+ A_plane(1) * ricker_Bielak_displ(t - sin(anglesource_abs)*x/c_inc + cos(anglesource_abs)*z/c_inc,f0(1)) &
+ + B_plane(1) * ricker_Bielak_displ(t - sin(anglesource_abs)*x/c_inc - cos(anglesource_abs)*z/c_inc,f0(1)) &
+ + C_plane(1) * ricker_Bielak_displ(t - sin(anglesource_refl)*x/c_refl - cos(anglesource_refl)*z/c_refl,f0(1))
displ_elastic(3,i) = &
- A_plane(2) * ricker_Bielak_displ(t - sin(angleforce_abs)*x/c_inc + cos(angleforce_abs)*z/c_inc,f0(1)) &
- + B_plane(2) * ricker_Bielak_displ(t - sin(angleforce_abs)*x/c_inc - cos(angleforce_abs)*z/c_inc,f0(1)) &
- + C_plane(2) * ricker_Bielak_displ(t - sin(angleforce_refl)*x/c_refl - cos(angleforce_refl)*z/c_refl,f0(1))
+ A_plane(2) * ricker_Bielak_displ(t - sin(anglesource_abs)*x/c_inc + cos(anglesource_abs)*z/c_inc,f0(1)) &
+ + B_plane(2) * ricker_Bielak_displ(t - sin(anglesource_abs)*x/c_inc - cos(anglesource_abs)*z/c_inc,f0(1)) &
+ + C_plane(2) * ricker_Bielak_displ(t - sin(anglesource_refl)*x/c_refl - cos(anglesource_refl)*z/c_refl,f0(1))
! formulas for the initial velocity for a plane wave (first derivative in time of the displacement)
veloc_elastic(1,i) = &
- A_plane(1) * ricker_Bielak_veloc(t - sin(angleforce_abs)*x/c_inc + cos(angleforce_abs)*z/c_inc,f0(1)) &
- + B_plane(1) * ricker_Bielak_veloc(t - sin(angleforce_abs)*x/c_inc - cos(angleforce_abs)*z/c_inc,f0(1)) &
- + C_plane(1) * ricker_Bielak_veloc(t - sin(angleforce_refl)*x/c_refl - cos(angleforce_refl)*z/c_refl,f0(1))
+ A_plane(1) * ricker_Bielak_veloc(t - sin(anglesource_abs)*x/c_inc + cos(anglesource_abs)*z/c_inc,f0(1)) &
+ + B_plane(1) * ricker_Bielak_veloc(t - sin(anglesource_abs)*x/c_inc - cos(anglesource_abs)*z/c_inc,f0(1)) &
+ + C_plane(1) * ricker_Bielak_veloc(t - sin(anglesource_refl)*x/c_refl - cos(anglesource_refl)*z/c_refl,f0(1))
veloc_elastic(3,i) = &
- A_plane(2) * ricker_Bielak_veloc(t - sin(angleforce_abs)*x/c_inc + cos(angleforce_abs)*z/c_inc,f0(1)) &
- + B_plane(2) * ricker_Bielak_veloc(t - sin(angleforce_abs)*x/c_inc - cos(angleforce_abs)*z/c_inc,f0(1)) &
- + C_plane(2) * ricker_Bielak_veloc(t - sin(angleforce_refl)*x/c_refl - cos(angleforce_refl)*z/c_refl,f0(1))
+ A_plane(2) * ricker_Bielak_veloc(t - sin(anglesource_abs)*x/c_inc + cos(anglesource_abs)*z/c_inc,f0(1)) &
+ + B_plane(2) * ricker_Bielak_veloc(t - sin(anglesource_abs)*x/c_inc - cos(anglesource_abs)*z/c_inc,f0(1)) &
+ + C_plane(2) * ricker_Bielak_veloc(t - sin(anglesource_refl)*x/c_refl - cos(anglesource_refl)*z/c_refl,f0(1))
! formulas for the initial acceleration for a plane wave (second derivative in time of the displacement)
accel_elastic(1,i) = &
- A_plane(1) * ricker_Bielak_accel(t - sin(angleforce_abs)*x/c_inc + cos(angleforce_abs)*z/c_inc,f0(1)) &
- + B_plane(1) * ricker_Bielak_accel(t - sin(angleforce_abs)*x/c_inc - cos(angleforce_abs)*z/c_inc,f0(1)) &
- + C_plane(1) * ricker_Bielak_accel(t - sin(angleforce_refl)*x/c_refl - cos(angleforce_refl)*z/c_refl,f0(1))
+ A_plane(1) * ricker_Bielak_accel(t - sin(anglesource_abs)*x/c_inc + cos(anglesource_abs)*z/c_inc,f0(1)) &
+ + B_plane(1) * ricker_Bielak_accel(t - sin(anglesource_abs)*x/c_inc - cos(anglesource_abs)*z/c_inc,f0(1)) &
+ + C_plane(1) * ricker_Bielak_accel(t - sin(anglesource_refl)*x/c_refl - cos(anglesource_refl)*z/c_refl,f0(1))
accel_elastic(3,i) = &
- A_plane(2) * ricker_Bielak_accel(t - sin(angleforce_abs)*x/c_inc + cos(angleforce_abs)*z/c_inc,f0(1)) &
- + B_plane(2) * ricker_Bielak_accel(t - sin(angleforce_abs)*x/c_inc - cos(angleforce_abs)*z/c_inc,f0(1)) &
- + C_plane(2) * ricker_Bielak_accel(t - sin(angleforce_refl)*x/c_refl - cos(angleforce_refl)*z/c_refl,f0(1))
+ A_plane(2) * ricker_Bielak_accel(t - sin(anglesource_abs)*x/c_inc + cos(anglesource_abs)*z/c_inc,f0(1)) &
+ + B_plane(2) * ricker_Bielak_accel(t - sin(anglesource_abs)*x/c_inc - cos(anglesource_abs)*z/c_inc,f0(1)) &
+ + C_plane(2) * ricker_Bielak_accel(t - sin(anglesource_refl)*x/c_refl - cos(anglesource_refl)*z/c_refl,f0(1))
enddo
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90 2012-09-03 02:09:56 UTC (rev 20673)
@@ -389,7 +389,7 @@
!
subroutine read_databases_sources(NSOURCES,source_type,time_function_type, &
- x_source,z_source,Mxx,Mzz,Mxz,f0,tshift_src,factor,angleforce)
+ x_source,z_source,Mxx,Mzz,Mxz,f0,tshift_src,factor,anglesource)
! reads source parameters
@@ -399,7 +399,7 @@
integer :: NSOURCES
integer, dimension(NSOURCES) :: source_type,time_function_type
double precision, dimension(NSOURCES) :: x_source,z_source, &
- Mxx,Mzz,Mxz,f0,tshift_src,factor,angleforce
+ Mxx,Mzz,Mxz,f0,tshift_src,factor,anglesource
! local parameters
integer :: i_source
@@ -416,14 +416,14 @@
f0(:) = 0.d0
tshift_src(:) = 0.d0
factor(:) = 0.d0
- angleforce(:) = 0.d0
+ anglesource(:) = 0.d0
! reads in source info from Database file
do i_source=1,NSOURCES
read(IIN,"(a80)") datlin
read(IIN,*) source_type(i_source),time_function_type(i_source), &
x_source(i_source),z_source(i_source),f0(i_source),tshift_src(i_source), &
- factor(i_source),angleforce(i_source),Mxx(i_source),Mzz(i_source),Mxz(i_source)
+ factor(i_source),anglesource(i_source),Mxx(i_source),Mzz(i_source),Mxz(i_source)
enddo
end subroutine read_databases_sources
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/set_sources.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/set_sources.f90 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/set_sources.f90 2012-09-03 02:09:56 UTC (rev 20673)
@@ -45,7 +45,7 @@
subroutine set_sources(myrank,NSOURCES,source_type,time_function_type, &
- x_source,z_source,Mxx,Mzz,Mxz,f0,tshift_src,factor,angleforce,aval, &
+ x_source,z_source,Mxx,Mzz,Mxz,f0,tshift_src,factor,anglesource,aval, &
t0,initialfield,ipass,deltat,USER_T0)
! gets source parameters
@@ -57,7 +57,7 @@
integer :: NSOURCES
integer, dimension(NSOURCES) :: source_type,time_function_type
double precision, dimension(NSOURCES) :: x_source,z_source, &
- Mxx,Mzz,Mxz,f0,tshift_src,factor,angleforce
+ Mxx,Mzz,Mxz,f0,tshift_src,factor,anglesource
double precision, dimension(NSOURCES) :: aval
double precision :: t0
double precision :: deltat
@@ -83,7 +83,7 @@
if ( myrank == 0 .and. ipass == 1 ) then
! user output
write(IOUT,212) x_source(i_source),z_source(i_source),f0(i_source),tshift_src(i_source), &
- factor(i_source),angleforce(i_source)
+ factor(i_source),anglesource(i_source)
endif
else if(source_type(i_source) == 2) then
if ( myrank == 0 .and. ipass == 1 ) then
@@ -120,7 +120,7 @@
aval(i_source) = PI*PI*f0(i_source)*f0(i_source)
! convert angle from degrees to radians
- angleforce(i_source) = angleforce(i_source) * PI / 180.d0
+ anglesource(i_source) = anglesource(i_source) * PI / 180.d0
enddo ! do i_source=1,NSOURCES
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90 2012-09-02 23:58:03 UTC (rev 20672)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90 2012-09-03 02:09:56 UTC (rev 20673)
@@ -367,7 +367,7 @@
integer NSOURCES,i_source
integer, dimension(:), allocatable :: source_type,time_function_type
double precision, dimension(:), allocatable :: x_source,z_source,xi_source,gamma_source,&
- Mxx,Mzz,Mxz,f0,tshift_src,factor,angleforce
+ Mxx,Mzz,Mxz,f0,tshift_src,factor,anglesource
integer, dimension(:), allocatable :: ix_image_color_source,iy_image_color_source
real(kind=CUSTOM_REAL), dimension(:,:,:,:),allocatable :: sourcearray
double precision :: t0
@@ -804,7 +804,7 @@
integer, dimension(:), allocatable :: recloc, which_proc_receiver
! to compute analytical initial plane wave field
- double precision :: angleforce_refl, c_inc, c_refl, cploc, csloc
+ double precision :: anglesource_refl, c_inc, c_refl, cploc, csloc
double precision, dimension(2) :: A_plane, B_plane, C_plane
double precision :: time_offset
@@ -866,7 +866,7 @@
integer :: n1_tangential_detection_curve
integer, dimension(4) :: n_tangential_detection_curve
integer, dimension(:), allocatable :: rec_tangential_detection_curve
- double precision :: distmin, dist_current, angleforce_recv
+ double precision :: distmin, dist_current, anglesource_recv
double precision, dimension(:), allocatable :: dist_tangential_detection_curve
double precision :: x_final_receiver_dummy, z_final_receiver_dummy
@@ -1103,7 +1103,7 @@
allocate( f0(NSOURCES) )
allocate( tshift_src(NSOURCES) )
allocate( factor(NSOURCES) )
- allocate( angleforce(NSOURCES) )
+ allocate( anglesource(NSOURCES) )
allocate( Mxx(NSOURCES) )
allocate( Mxz(NSOURCES) )
allocate( Mzz(NSOURCES) )
@@ -1120,11 +1120,11 @@
! reads in source infos
call read_databases_sources(NSOURCES,source_type,time_function_type, &
- x_source,z_source,Mxx,Mzz,Mxz,f0,tshift_src,factor,angleforce)
+ x_source,z_source,Mxx,Mzz,Mxz,f0,tshift_src,factor,anglesource)
! sets source parameters
call set_sources(myrank,NSOURCES,source_type,time_function_type, &
- x_source,z_source,Mxx,Mzz,Mxz,f0,tshift_src,factor,angleforce,aval, &
+ x_source,z_source,Mxx,Mzz,Mxz,f0,tshift_src,factor,anglesource,aval, &
t0,initialfield,ipass,deltat,USER_T0)
@@ -2345,7 +2345,7 @@
! users can give an angle with respect to the normal to the topography surface,
! in which case we must compute the normal to the topography
! and add it the existing rotation angle
- call compute_normal_vector( angleforce(i_source), &
+ call compute_normal_vector( anglesource(i_source), &
nodes_tangential_curve(1,n_tangential_detection_curve(1)), &
nodes_tangential_curve(1,n_tangential_detection_curve(2)), &
nodes_tangential_curve(1,n_tangential_detection_curve(3)), &
@@ -2358,24 +2358,24 @@
source_courbe_eros(i_source) = n1_tangential_detection_curve
if ( myrank == 0 .and. is_proc_source(i_source) == 1 .and. nb_proc_source(i_source) == 1 ) then
source_courbe_eros(i_source) = n1_tangential_detection_curve
- angleforce_recv = angleforce(i_source)
+ anglesource_recv = anglesource(i_source)
#ifdef USE_MPI
else if ( myrank == 0 ) then
do i = 1, nb_proc_source(i_source) - is_proc_source(i_source)
call MPI_recv(source_courbe_eros(i_source),1,MPI_INTEGER, &
MPI_ANY_SOURCE,42,MPI_COMM_WORLD,request_mpi_status,ier)
- call MPI_recv(angleforce_recv,1,MPI_DOUBLE_PRECISION, &
+ call MPI_recv(anglesource_recv,1,MPI_DOUBLE_PRECISION, &
MPI_ANY_SOURCE,43,MPI_COMM_WORLD,request_mpi_status,ier)
enddo
else if ( is_proc_source(i_source) == 1 ) then
call MPI_send(n1_tangential_detection_curve,1,MPI_INTEGER,0,42,MPI_COMM_WORLD,ier)
- call MPI_send(angleforce(i_source),1,MPI_DOUBLE_PRECISION,0,43,MPI_COMM_WORLD,ier)
+ call MPI_send(anglesource(i_source),1,MPI_DOUBLE_PRECISION,0,43,MPI_COMM_WORLD,ier)
#endif
endif
#ifdef USE_MPI
- call MPI_bcast(angleforce_recv,1,MPI_DOUBLE_PRECISION,0,MPI_COMM_WORLD,ier)
- angleforce(i_source) = angleforce_recv
+ call MPI_bcast(anglesource_recv,1,MPI_DOUBLE_PRECISION,0,MPI_COMM_WORLD,ier)
+ anglesource(i_source) = anglesource_recv
#endif
endif ! if (is_proc_source(i_source) == 1)
enddo ! do i_source=1,NSOURCES
@@ -3574,9 +3574,9 @@
! Calculation of the initial field for a plane wave
if( any_elastic ) then
call prepare_initialfield(myrank,any_acoustic,any_poroelastic,over_critical_angle, &
- NSOURCES,source_type,angleforce,x_source,z_source,f0,t0, &
+ NSOURCES,source_type,anglesource,x_source,z_source,f0,t0, &
nglob,numat,poroelastcoef,density,coord, &
- angleforce_refl,c_inc,c_refl,cploc,csloc,time_offset, &
+ anglesource_refl,c_inc,c_refl,cploc,csloc,time_offset, &
A_plane, B_plane, C_plane, &
accel_elastic,veloc_elastic,displ_elastic)
endif
@@ -3608,7 +3608,7 @@
allocate(t0z_bot(count_bottom,NSTEP))
! call Paco's routine to compute in frequency and convert to time by Fourier transform
- call paco_beyond_critical(coord,nglob,deltat,NSTEP,angleforce(1),&
+ call paco_beyond_critical(coord,nglob,deltat,NSTEP,anglesource(1),&
f0(1),cploc,csloc,ATTENUATION_VISCOELASTIC_SOLID,QKappa_attenuation(1),source_type(1),v0x_left,v0z_left, &
v0x_right,v0z_right,v0x_bot,v0z_bot,t0x_left,t0z_left,t0x_right,t0z_right, &
t0x_bot,t0z_bot,left_bound(1:count_left),right_bound(1:count_right),bot_bound(1:count_bottom), &
@@ -5363,7 +5363,7 @@
call compute_forces_viscoelastic(p_sv,nglob,nspec,myrank,nelemabs,numat, &
ispec_selected_source,ispec_selected_rec,is_proc_source,which_proc_receiver, &
source_type,it,NSTEP,anyabs,assign_external_model, &
- initialfield,ATTENUATION_VISCOELASTIC_SOLID,angleforce,deltatcube, &
+ initialfield,ATTENUATION_VISCOELASTIC_SOLID,anglesource,deltatcube, &
deltatfourth,twelvedeltat,fourdeltatsquare,ibool,kmato,numabs,elastic,codeabs, &
accel_elastic,veloc_elastic,displ_elastic,b_accel_elastic,b_displ_elastic, &
density,poroelastcoef,xix,xiz,gammax,gammaz, &
@@ -5374,7 +5374,7 @@
hprime_zz,hprimewgll_zz,wxgll,wzgll,inv_tau_sigma_nu1, &
phi_nu1,inv_tau_sigma_nu2,phi_nu2,Mu_nu1,Mu_nu2,N_SLS, &
deltat,coord,add_Bielak_conditions, x_source(1), z_source(1), &
- A_plane, B_plane, C_plane, angleforce_refl, c_inc, c_refl, time_offset, f0(1),&
+ A_plane, B_plane, C_plane, anglesource_refl, c_inc, c_refl, time_offset, f0(1),&
v0x_left(1,it),v0z_left(1,it),v0x_right(1,it),v0z_right(1,it),v0x_bot(1,it),v0z_bot(1,it), &
t0x_left(1,it),t0z_left(1,it),t0x_right(1,it),t0z_right(1,it),t0x_bot(1,it),t0z_bot(1,it), &
count_left,count_right,count_bottom,over_critical_angle, &
@@ -5871,9 +5871,9 @@
iglob = ibool(i,j,ispec_selected_source(i_source))
hlagrange = hxis_store(i_source,i) * hgammas_store(i_source,j)
accel_elastic(1,iglob) = accel_elastic(1,iglob) &
- - sin(angleforce(i_source))*source_time_function(i_source,it,i_stage)*hlagrange
+ - sin(anglesource(i_source))*source_time_function(i_source,it,i_stage)*hlagrange
accel_elastic(3,iglob) = accel_elastic(3,iglob) &
- + cos(angleforce(i_source))*source_time_function(i_source,it,i_stage)*hlagrange
+ + cos(anglesource(i_source))*source_time_function(i_source,it,i_stage)*hlagrange
enddo
enddo
else ! SH (membrane) calculation
@@ -5893,10 +5893,10 @@
iglob = ibool(i,j,ispec_selected_source(i_source))
hlagrange = hxis_store(i_source,i) * hgammas_store(i_source,j)
b_accel_elastic(1,iglob) = b_accel_elastic(1,iglob) &
- - sin(angleforce(i_source))*source_time_function(i_source,NSTEP-it+1,stage_time_scheme-i_stage+1) &
+ - sin(anglesource(i_source))*source_time_function(i_source,NSTEP-it+1,stage_time_scheme-i_stage+1) &
*hlagrange
b_accel_elastic(3,iglob) = b_accel_elastic(3,iglob) &
- + cos(angleforce(i_source))*source_time_function(i_source,NSTEP-it+1,stage_time_scheme-i_stage+1) &
+ + cos(anglesource(i_source))*source_time_function(i_source,NSTEP-it+1,stage_time_scheme-i_stage+1) &
*hlagrange
enddo
enddo
@@ -6645,14 +6645,14 @@
hlagrange = hxis_store(i_source,i) * hgammas_store(i_source,j)
! s
accels_poroelastic(1,iglob) = accels_poroelastic(1,iglob) - hlagrange * &
- (1._CUSTOM_REAL - phil/tortl)*sin(angleforce(i_source))*source_time_function(i_source,it,i_stage)
+ (1._CUSTOM_REAL - phil/tortl)*sin(anglesource(i_source))*source_time_function(i_source,it,i_stage)
accels_poroelastic(2,iglob) = accels_poroelastic(2,iglob) + hlagrange * &
- (1._CUSTOM_REAL - phil/tortl)*cos(angleforce(i_source))*source_time_function(i_source,it,i_stage)
+ (1._CUSTOM_REAL - phil/tortl)*cos(anglesource(i_source))*source_time_function(i_source,it,i_stage)
! w
accelw_poroelastic(1,iglob) = accelw_poroelastic(1,iglob) - hlagrange * &
- (1._CUSTOM_REAL - rhol_f/rhol_bar)*sin(angleforce(i_source))*source_time_function(i_source,it,i_stage)
+ (1._CUSTOM_REAL - rhol_f/rhol_bar)*sin(anglesource(i_source))*source_time_function(i_source,it,i_stage)
accelw_poroelastic(2,iglob) = accelw_poroelastic(2,iglob) + hlagrange * &
- (1._CUSTOM_REAL - rhol_f/rhol_bar)*cos(angleforce(i_source))*source_time_function(i_source,it,i_stage)
+ (1._CUSTOM_REAL - rhol_f/rhol_bar)*cos(anglesource(i_source))*source_time_function(i_source,it,i_stage)
enddo
enddo
else ! backward wavefield
@@ -6662,17 +6662,17 @@
hlagrange = hxis_store(i_source,i) * hgammas_store(i_source,j)
! b_s
b_accels_poroelastic(1,iglob) = b_accels_poroelastic(1,iglob) - hlagrange * &
- (1._CUSTOM_REAL - phil/tortl)*sin(angleforce(i_source))* &
+ (1._CUSTOM_REAL - phil/tortl)*sin(anglesource(i_source))* &
source_time_function(i_source,NSTEP-it+1,stage_time_scheme-i_stage+1)
b_accels_poroelastic(2,iglob) = b_accels_poroelastic(2,iglob) + hlagrange * &
- (1._CUSTOM_REAL - phil/tortl)*cos(angleforce(i_source))* &
+ (1._CUSTOM_REAL - phil/tortl)*cos(anglesource(i_source))* &
source_time_function(i_source,NSTEP-it+1,stage_time_scheme-i_stage+1)
!b_w
b_accelw_poroelastic(1,iglob) = b_accelw_poroelastic(1,iglob) - hlagrange * &
- (1._CUSTOM_REAL - rhol_f/rhol_bar)*sin(angleforce(i_source))* &
+ (1._CUSTOM_REAL - rhol_f/rhol_bar)*sin(anglesource(i_source))* &
source_time_function(i_source,NSTEP-it+1,stage_time_scheme-i_stage+1)
b_accelw_poroelastic(2,iglob) = b_accelw_poroelastic(2,iglob) + hlagrange * &
- (1._CUSTOM_REAL - rhol_f/rhol_bar)*cos(angleforce(i_source))* &
+ (1._CUSTOM_REAL - rhol_f/rhol_bar)*cos(anglesource(i_source))* &
source_time_function(i_source,NSTEP-it+1,stage_time_scheme-i_stage+1)
enddo
enddo
More information about the CIG-COMMITS
mailing list