[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Ër›0Ýç+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 ¹¯Ózn2š“ a8PØ3í^쑵p!!<1õB#A©‹µq)ð2“´ñÂÅ~„Ô^öV"Ä¿
+xÚíXËr›0Ýç+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 ¹¯Ózn2š“ a8PØ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ŠÇ2„ Cç²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²hƒX yßþ¢u.ÉnI/k—ÿáäÒU†»trÖÓDA!‹Y™(içp¥UVÔyÛˆl
-pÄ_^AС‚ #µ:QY¥mÁ,p¬ÐÊ÷6u$ðlx<
-xû¡ežT&›·
-Ob¶þ/'»¯¤qˆç‰zBÝú'nÌTf²¾©ÎµJ3§Œó×ÜZGûyêzÁaúSgÊá;YSÒÖÓSë×qÜ{z(†M½5þo2Jww =3IæšÀIÃU߃¤ª,S͉›•kúÜ,G’¹åò7Q3 ›b,›Ì±ù%Ëû¦:ßþXi]™ò~ gbÀ°ðÑÊ–G>lœ÷§ALüzf½ÖBëöC} KÐ`– Ä"Š¸ÅBôßývëÁñôàH S
+xIËØæ9.&õr©Kc©àxOÅ bKz!$‘°ê ‘¥:"'ÀSiL¦®ƒ5L™¤42 !-B´áeÀ€ŸUu’g¿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 Þ‡ß {œx7vR y7U‰Ë‹¡ eŠPÒùë·ÀCÐx6ûD°qt£ïnò‡Õ³›7YqUúeeW«ÍU ]åqqµZ_ýî½úé寫ۋe”^â/–ix Ià}|ÿþíû qbõÓ-î~ùy÷/³å·»_o_½¹}½^|ZýãæM’Œ%iáq~ ®0‘k®–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Çöô†IŽrÐ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‚:Ž Šm—‚2±ç¹ÃÐ[áä–òõ©!ûå%×Á*á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º#¸Š›Ûߢׯß~p˜gTÆ~bN•ªï›Xžú1àËj])ÞI}@Y)_!ŒjJWso†5/€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 ^=ÑIŠ0ã Hºl)ãJfœ5ÝÒF7¹–íûÚVõ èVÆp AÈŒ9dŽ>È<‹òi¸2ü2™w8 ZÍ=OŸË”OÆD™kh(¤.ªLZã+Ø©®ÙS4‡Á¨8gM€5µ
-M‰Á²â$³þË)¬tÜ`Òþ\À?Á× ÛîÐ
-â°˜EnëP¤†äú¸cÕ†El þù õ½@€*˜ài£¯B¹3î <¤ý2ÑÑbðm󩯇z;:çeKºÙ 3k11A )KG\‰‰U¨ÉÔ/Á!?9kÞ‘ô7ƒrmÝ·Ê 6 ‰åàl¥ÂšãtÀ„4~ßÖ˜|èÞÕ9ƒÄO‹¡{xãÔ€(öË!æû(õ·›Yn5Ô!Cb6o·_ø ÷¢k(g:+ð.ÙÜþ<~Ž ,¶ÓÀÈ~Ä'¶mL>rdÝn:¹#æ(Ð<&{C“ú]/ûØs™´šŽ@ ÚÖôo“iFI|€”t®[gêóÏE ÞirïDµ›uVr²…â[×hs£OóÎÁ ”fšù½a[ÓsŽ´o™ì2ûìÄäÞRK—G¶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@º—²<Ó´7žiYsÓ0mܶŠ3çf|þµà+ØýŒ~–´Ù8½)‰T{yn ÁyËøÜP¹…Cø(êJì ¨—Q/ÉÌØm99Xã¡oÁ´M k‹]É¿fRÓu±öˆ‘ïÌ'¨¾$*q®¤úB¶Æ÷‹ÞZ›†Ü½”zóÊ¥¦± ‡6Ç}ÒºnjË¿Þ¨c •ÔI4í„*ÅXÉ`kà Žežûä\“Ξ|¡èË “‰S p«pVÊ(—a+œ©îe¦³p„®5çœÄ1^þXV P|]‘ "WV–Aêf[†¬q¦w~eÿ‡^J“(ÓI¯nù7hE¨?N.“•ñ#nBå& úßmÀŽ·x´• ª¸ÌÃ&m(Gbm=gƒÖ¨ tæÊÄŽm:c—e&£u˜X™¹)OüdXâèŒÿ €¬ª*³X Yø(K—/ÎOˆü<ËgLÓ,8`žÖ¯žÈ!7PýtLòÂhÕw U†& 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}RežƒRpB黇‡ë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»ŸX„AàýË.ùåîÝí«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€çEšxš™÷骻 I=á°¾kî)YÅì ^2Ž;Í‚7‚Ztk&ñkê…Ëreæv¶ÁÍZÞ,ÂÜÛ!ëðû 7¥Ó—V£ÀÖX27 Ø ßÓa¶Ü
+º©=8Dmµnä±µÔ£cÃp©mxE MÝ`„™'MïF6bÝÖÍ™|@š| KR?+ë>¬ÞþpwãðQâçi?ͧk²‚ùijî™Îï)QGƒÜ;× 62Í\m|U–²#`¢
+вÌ.Ëç~†ÖhXU«hÜ.TÝn'T+!˜Ž2& ýV†zŽ 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õ}Ë?‚óè¢ÆDWšw¨$+í+”UMt5÷VadõØÞ7b¬Ò¬ç¥O_âˆ7^!ö‹ ›*Õj/\Ü €‘3KöÝÏ Þ¿ºuùžÜ’l¸]ÀJë¾ëõ ÜÖݵγ²ÐŠÐnòúåê¥ÓÍ1?Ê
+;kpœdeûº+Í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æÕKOR‰a‚!„+ÌØ·%jêi ®A*i
+£&:Ÿîà’ÛDá}V†‘ŸG½ñéÑLÑο>Oõé Z .¡
+'
+óå“Ùþò? Ìòµ= …KQ? É¢2Àt¤¨,M!<ð²ô‰x¹4s\j[„~öj:NÌý”Ev‚)zß’T!RÓ·Ö`CÍ*ÖuXÂÌ86€»¡6K ·¤; ‹ º'€=‘Ât ñÙYTÊRN’S@»ñïQûašM½Ï,‰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´p—w–¥³, –^¬W2”~>:7ÑX ÜÌ0©6Ú®©ˆq Ê d»7`
+sBíÄ\™O úÀk{²‚‘—[‰éú¥f"Y<Lv¹#?ã ÀxÇ›[Yº|a–ùY_ ¶¹ŸÓ̲têi1éÿ¸íª5:mžùqÜ;ÛWûº¦(©u8ŠCã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‚,O€h¢7å½q*pÈ?MgPa|
+¢ŠÍR«|YJš’%Fб©²â•™¾ð28“ªa—†ŽãD:Ö
+óq‰ FŽµ¬Úe[q®<te_®ƒŽNQÉ 6äÔ3cÎp1…2OY6uŸ–_6¥N‡Aë±çÉsô v1‰0½“èúdž>ÁNM1.°Õ°08tÎY:Æ,+ŠÓ ”8…•Œ+Ǩij:ŸÓçÞFÝe”ƒæ3HÞ;”k
+§½\ï¡™cÐl t ëj£¨CWBà×
+F#•\»3݈N ÍÇDGóÁ·™s08š¦YNˆkÙðv¶Â N)3´¤š ¹HË]¸bbz0ñ‹øéÙð b(j )y×Dw˜,$Œƒ¯E} Ç0ߦ÷õ}Ëeeóœ{WE0 b?ɇW§°È/,ï“Ðßng9ÓP_®yպ¯xS[`)øÒYáæ!ÂåÏ£çT ¤"à #úŸä®â¥%³ÝO÷ÜnŠ'±æ;©Ë™iÿQ‡“Q[é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×8W0ŸrU] ÆÊõüûìà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âœ"ò³¬w—O9ô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’Àë[ üÍŸ–1ƒƒŠhßÇnBcPý¯Ö-Ç¿š;ŽÊŸ¨âˆÃÆ~IC9õ!¶O§Àâú 5ðý…+S{¹me<Ra‚h %!7å± SlÅ›‚§ž=þg Hªìd³ ðJ§/.w`~–f³
+¦0 6˜ÃúÕ90`ýNF /œ£|ÏÕYhìñ£iFæ9>öjžh¼,©ÓÔ¨ 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š&¹·6 T•¶GiŠQE1r30¯°ÁÅô`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ï_‘7 hT[þ>†‡£´ô˜»R.a˜Žê(ÛÊÙr{å¯gW+9vj>&–VÒîjû¥|»9»¸N²EÎò„'‹Ín‘òEf‹Ívñ›wyóönsõq¹â±ïEl¹Šßû¸Œ|ïçÛÛw·ßÓÂææŠë ß/ß÷~qG~]ß]]^_}àß-ÿØüpq E#AQœ1?L>É
-bÜsvµ9 €ä/‚EçŒGñ"æ9K‚hQÔgŸÎ|e!çfËxlÝIK¸xWç‹ïÔÙOðsK+Çu5bû-Za¢\Àc–Eé"öS–ä ©x]>ötÖào–«0ð½KU D[vª¡µ£¯ÞÛÝ^ ätMØ-m)šÇÊ.Šf{rò{Ñw l¡Ù®o
-]:1¢Ð½¨ªšõÜ2käptÏX - at Q–F9]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> ýØøj†e ;¢Á«µ?ãø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>X”yïìñÕÓàÐÊ ä (¡ß~z\)„eŽIì3
-õÞI²jˆ&@Ç;©Ñr›(Œ(ÑŠZbÆyíe‘nnöóúêãýfÎ×¢ˆ¥y0©ùcN˜ÂcðÔÚÐVÍ ÑJÛÏ×›$FaÆM9Ìrø"~¹ï TWêrÉcëXçäÇjVÖ‘`ﮋª7ù
-’çfo™3ÎKûÕ{¡iTã0 é=EõÕÖ u–UôH%\YÂØùk Èž„L[
-÷˜1•½÷ÄÛVÿCÊ2>˜éXéé8'ëÓ ,ù°¦NôÇD4ሞ:ØNo-m¶¼Ö8„€É‡\§ÕVÝWe§ï•ËÝoU£ï[Y«'Éôg=—µòŠúÖw®W°½2›6I6
-á|¶µü¾ìÎizhÕC%k7uÜêŠ'å LB˜7Ç/
-½¸µ—êáB¶[NÓý§©’ÆØlí¨ëÕb „3Q¨¾Óeq!+_ôß,<’a«/qÏÎ^„Ô)Ž™žÏìµJ˜;¬æ’ÛÑ6¸ÂC–[§¹šcFm Xg°îñ~Ô½B É£çe™ÄÞù7‹m£ì@hèûq«æ šË‰èª
-tHÕË %·P„!‹cÛ1F,X®à!â{7
-µy¶Œ}ÛÞŠ¼[™zÏúÆEAìײ~hE#±JõYØôcê¦1ÁyçŒ'VkÌ[NÒìåœ3SM€ÁžªœIxƒu:÷^ØnKmÍÑÛßÌ4Ö>µ¼Ü¾ ÌÅsîÍ= àUCœo–0>±›85ŸÉ˜#UK¢*ŽY—bö —¢õà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‹ÝY‚Mø®¯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ÛQ“d·* 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òcR‰2^ËCsSÎÞWÉâ-?ù ~æÈ5TÝÙ×Â
+á
+PìG‹ G(LB%âEqß7¹±}¹t=‚7¼Ú±¦hyNøF}»Fm·l—OϘ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áú‚Â7”Zß¾.¤˜Â†’ãØÆEŽ]d%wŒÄ©RFj‚"+6›¼±w
+í®!¦¾>8÷Þ¾œñ"a‹ªõ9¸²*òûŽÈ hÿOê ZŪhfû:‘Vs!§%pI€ŸLD¦ÿ]dú¿‰|ýbVÜ„¢$&S o »€î±¨ Ô XiQYµ+‹ÏÓ”
+‰¹èØûyêÊ„¾a잧ôòHè<ç
+W«®¨4L£„¤¶¨úR3û”7MÞ‚ÞY«1ùÁ,’7crÑph
+Й©4R(êðf†JÃEfËû³
+mñ4órÖtÛ¯=û"l’‹ü
+Cg%‰bê´ ‹%1()·³Xc•a_´Ê/„ŒëuBgC¨S“¶†T§ç¯4…IŇÂ1!8ωüØ6¬1|%Ëåå1Š¢Ðv<:#E„‚( á) ¸4ÅI Ζ•`°ºnÖ7ÊJG\cŒ< ˜?Î$€áÛ°îðLä_,Ñl=4Šx¹'ˆZW)_rååJ‚v¬h¦*Eüù؛ƎNøx”ð‰*jè¸&’“€>,IàäÍ^AÓ’·]ˆ @*Ëäq›×j%ã°¨ï†å;ÞO =ç2gKŠ‡¢-²)Zu·7ì˜Ú-#ÀЄE#;ôÈQÒ„c—Ši]*wRGøò»Žu.»làÜ픜m°H0?vÞëë{ޫŮɡ'‰¹MvØh8I™&.ªØ“nÌ[ÃI'¨Ê 3lóN-Ñ…B²b
+«rQqŽ£ÌÇfž
+³·ç Ö«¹Xó}
+#HQ<Ž‰‚Æ"4šO
+.éô´o8a at FiFå –&ðþK°³ãmÑ¢§«À:Uq\óÚÕ¸zìée}ƒâ¹ÚjbãÊ%ö…þš±…€ùeÀ §G)¼/3}ÕX DéF"‰“%¬M¼&Äz†89TÚÒê1c*÷ÄR!B1µfÆKu*ë+ —Æ( æÙÀú¡M(ŠH4Þm®;´öÀ±ÄbjMlëxÆ×eÑvë]™ƒåÖ¯»u“Wü!GÝS7WµECY·ãŸžî±éØ÷¥æº™ß‰„ïïÅd*¶»†ß•ye¶†ZU@óTõ€;YæÍEÅ»"òìÅ{PH„·6hÚoPuÒ@Œ™^µýnÇ1‰ KyßvEz–—L|EüÆÞ T>>¢ŽÞ xjd =ÁÕBH ܹ‰`°8¡ Jô„x>GLÍu` kÝA?5Fû “Èžý2 œÓ±ß´ok®¬ƒ± ÏjGÏ ë
+º
+LHåÈ—¾h”T»ÂóPH}%¿ÈÒ…×/v.¹æQæêÛôšå{ûi=·—& \åÕ]Ãê\t }dºüȾ)Mð¡QH#õfuËpšUÎ³ê °Øª.' žµ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øÍ;Õ‡œ 0™8UÿŸ Œ¦ÆÓÈ å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ä.(QP‡e;‡‚ýr%º*´¾€†ÖøÍ;ÐO3°èk0’Á²Žê€™hð-¥ðg|“££B”$*xlˆñÅgz¯ñ#Ù¸L`:e™ß@@hYwWG½è³ÌÚ *Šî ‡ÅÃvšLáå7X$¯†¦‚žF%mòû¬-}3â N™?¡í®ŽÚÐoÔäýP-ß\&ûÙÜøøûÏMhfÚŠ<Ò] ‚# ÉÐxíê(Ö<SIˆx»ãð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~8‰Pä²¼™dÐp²L'¿N>ìŸ.Ïfs½)C³9÷½éÙŒyÓ‹ããÅñ?ôÆòá œŸ|šaÏ›^ÚOþ}~zxptø™¼Ÿ}]þs #ÆCäÑ`âi^Ø7gè3D™OÔ¡¹=5'|IõY†Èlœ¼é‡
-8Nï4»¦Ò¿)ô .3Y5õ ‡ÓjýàÁºß×tf˜O«<¯f„Oï²òV“9ˆë$«JCT4†xŽôàW …á×]=¹#š•$y,e–ĹžÖ"Ý$"Õ“Ë&Ñ´Ên
-í²jâX D~Ý8ŠØîP03%ž|}§ÔÜ9\î`ÐΛà æ E„O8À ºI±ómÇC”G˜·GúãvÓ~ è VË{‹{“÷ÕÎøoOÌ-ñyú/ʦ0ƒMq 4vµ\m¤9+G(•4ÜŠRÔ-°
-blqµ@¯âñ¦ßÕÍ¿yz• `Êìñr ñMU>–M-¤œÃOl…¬EÞÞŠ\eëwÚšü`hM8B„+ýËlfÅ:0ÓSc“Ê õv©µ´40Š8×4dv[ÄW÷÷úÌÏú'Á¸=Rä‡>|Ð"ßžÿ[»ÁQ„“Þrº¹¿J
-‘¿["¼=À1ŠŸ!ô‰üæq/Ý<Ž µ´ aèÓ– ¦b! ƒ–øc»7§ØC!e.4 GhÐ €ý“Ñ ÷c4ø4˜
-ö<ü)Ì-,ûÿÐàOÑx™ÁQléaÓ󽻕¨ÇáÀUzâZŒ|?+e—‰p¸6f 1Ü%«¾ªøTApÑÓuϬß;„Ƙ" 0K9S áò"€ŸRÄ „šÔGß”úÖ•N{²d“%YóRú Lº_‰ä÷×c($B5·˜7Í6˜;£â04ªžÆõÕM– WA Øsï´áã êss¹ƒ…‰ÇR$m·àõHaĽޟþj3JUˆ¦~гêÆfœTä¶^0ù¸r¥G©IŸŸ\yÿ†R ˜L7®k
-åÅÙñÕåâü`¹<<¾Ø_.NŽ¯NŽ] ZéTÙu0àwð~ñœÖÖŽÔR¤Õõ)-e<²GÎ)(Ä;ž©h³´ ‚4+=úžÉ¤j“<LÒ¸XkëI›Lý`ë)Ø©µQ½Pˆ4ÛÈ%àÈðËš2䓬Lj#bÔGŒ…ÃèÙYø·Mœ+÷Ðæ®|*NdWí5ãr2zWžÜÔâÛF”‰*+ \NzSò®$B!ÔH]TP±Í¥µ5VµRÑÄY.˜fõ3ì[1†&½œQ>«bmjJ(ïHD€qÄüÖË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'9“3ÄÈmd8>?¹8;pY£(d¼ ZÊ•¨Ü×¢ÖãÖ:¶ªÓ‹èôT9V‡Ãð`GŪ'ºÈ€AšÉu ?ˆÔÞF¤Ïzд¡1¤†TH¦ÊZsW9W†·Wö¬¾”£Ð÷û ¨ß²d¥é›¶Ù‚Åk¡µÄŠ@AîZ¨pèŒj” îï÷—û{o‰jÌ ƒB´°e[D çOXõO){WcälD)ÛXc
-tO Ëhˆ¼ "{mŒø6S–fº¿azŒŠÝ®N n¥¥@IaÜÑѹ: ñ!åøl ÍÁçåùɧ ®\ ö3D§p0V0ˆO]袦1Ö:Õ²ëhŸœTþfã}È &#?&Þ·QB$"Sá黎Qz'ÏÊ'ÁæY²ÍiRd²¨ná„;;Dൠ¿RbqÄH4Êû—‡WG'gÿÚ?s½?¬o¼©’®ÈºÜáS5ŠMöq1ºŽøõ ©Ú|Ý{´Hí-sÏUô@ãù]E†šz#Ðs®Þ…6VTõ`AîŒú>ð}@»Pnâ\Š·±Ñ%óÁ’~EÙ4Ý«Ï%‚.ë}D·rÛÄãAç„GÙz¹m ãmoX‹®a’2»Î{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ÿEŽSõÔ(Õ»£j´(T$¦ëÚ‹U£ÖèU C ÌT_£féNn œ®Ëè.ÜB0 »ß°ŠZÛL)ÖÊ–Õ‚zµPóGQWý%†®võ~¬GÙ[í˜æ]íZ=vUà£úU‡ö \AEƒCÔµƒÛ ãgÔfðí”^¼©Ñbº~Çö V¶ò´ëÃö·w ¶Ÿ»ƒ ùÞŠó¬õÒ^奨ëÊ w °”Ä¥^Jª²ÑÏÝJð- D¥x-b;‹¶ZµÒÏM]Á<ˆç£4ÜiÖÅÐ^nlâR —üaæsÛð‚ß\™3?Û^$ÏR›Ji ³Ô 9þ`ê®ë
+xÚµXÛrã6}÷W(Ù‡P›L/S•E#ïxk,;¶<»©É”‹&! ^4¼Ø–¿>EÒôx*›” @`wãt÷é^ ŸxÁ$D¡çx“õfâ;Ÿ“u2ùd-ÞÏ/ÖËËéÌ¡¶å¢éŒz¶u9umëzµ:]ýKý±~¿Tƒ«óSlÛÖGóʯ.–‹“å™ónúyýïã×í(ri€lâOl¥ {bÏdf–gŽ[‰ús½åˆíÕMÉÔD-b«ÈÓ½^Év)ËX^³Dï/J5`iTÕ<V“,ªYÉ£T¿ åzwÅòŠ×|êPëž×Zæïp(‹•9Kµ
+q”«Á-SïÇQ7iÔ*ýͦvÅ´‘ëi`[ÑNÏXžQŠÔ»Ÿ Ûö?˳ۓÆ(¤T Z™0&¬Ž¸T“"WÏŒe·e”35«šr#ÌŒb½ð gSǶZ¬›0’j°‹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þ¥ÑËàjŽC…ë!‘ü0‹äÃv0É!ÄZ¬.¯ÄÐ ³¨¬ ä_®u.ˆÍ>Ž°iWÂÞ{ž° •Ä:Õ_€]!¢ëAo-¤aPzCS\ ¼WKUǡѶê^ÙEe”1j-WE4 JV5i]ÅW±Q¡©Ç
+{©FÇXÅ3ּȑv ÑC\Ïþy„ Ïq•<9Óä¿m½/Dô=(!òdð„ØÒúr^u9ÅUŒÆ…-\ƒ .)÷D˜Tà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–Ú-É6šcx>@)+*Á ËYiÐ^W FWô¶ ³]Y›ç®ŒS×lÏ· ؾ\ÕÐÕ ‘1²d*D«-ß½ ¡%õ0qÊï¦3Ïî
+5Õ1)âQý«§ª\ÃÌ©ø]Ý<>ª=?©GŒ±Ü@J ¼€|W#ûOùE¡ L:ËIóx“h!?!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}Ô‹¨¼Ùð”µé>
+ ³ï
+| 8ˆx´ïÜSÓ4°X2¸¯#
+C±iáýþ ¦¢Ð¹”{ÝÞnLÅIXjú…Ä´¼ÕÖtĪ#=ß92„NÔé[= `æ`Ñ~ꮯ/W7 O¯óõz¹ºž¯OÏW7ç«12E:ësáÞ_ìÑhëvË¢Yð yÖç²D8ÐÐl=¹\Ž
+„ZOi«3a²J{ÄH½U£{^‰ÞVM’(Û©è‰,¦ž/{A˜Bœ S Ð'¼ÉИm€£‹¿~RWfŠÞÂ+}Ñ# rÝ ï‚6¿4Qªn`"ÜÕÕ¦6”PµÝ^=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¤õÔ ÆÇŠðu Qwá2
+kíÄÛÎáí(ù_ÁóÚð܈“û_d3¬>Æ}Àšð ÿ?7VE#+ˆ¹æÝ…¡pH‡ä
+vø¾¸õŠëzhz!ÐÆЊ'z”¦ À¤jvÀ µš ~5×ë`OÌ*4Zœ#Ì°º:¿¾\ŒEƒKPàÒ=Hɵ±(o²[Vª± „V¸é…D èÜ;[)ÚP5QM ^íÒh¯>gOaäÈþ¬¤Æ€hQc‰hMÇÚp>¸ìÅóŠÏë?HØòx«äg¼lÁ¢ü¾Oe‚cL `wÉ Ž²C8-îïæëùñ·°šk: @A ”WØ\6 Yâ ë‘zdP²§àmgÍòJr!ŒÐ“(.#²ý°ì-{ñ˃¹ÜPö¦íû%ûœ¡8Óé8rs ƃ’søâ$¡Yœ¯Î?\ ºË [!Ú³QÛDAŸZê"úb¬ÎTVíöÙN‘o†ïz2ç¯éª$K°˜qAO÷Š£Ô?)ÏŸ‘‡)³Î¡¦UŒWYq;Æ«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¸f fH8`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¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢œËµˆãƒŽ9 6Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢Š(¢¿ÿÙ
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 @@
xœc`Ù
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Éý Êf•fˆsr±< ÌÂYÈR?dÉ,‰"?ŽÓY¹=ùvøy˜3fP÷æÐQàôõ6 g/ÕÉ{ø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
-š÷R‘M{GAT¡ÝΟ°4
-2?Ì#gC=é9¼„ѹÄpo!^Õ"_D at utcÆl(@Ð/| žÖÁʘM fAì YPd!µ
-ÇÀZà’Ö~8^ŽQ¡©¢ìåœ%Þ÷y˜x¢¹›§‰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,l–sèÌBuÊcŸYä_§*ê·`ÃÊûZô«uVßE£JÙK?êb²äæPž‚±ñ¦Qè‡[´5-Ì$+&Ò`&üؽ~6îÜÖZƒ0"‡¶¤V_&õᆘŒ ö•‚
-\o+|X¤`~xì†Z‹oÃè?Sdia(H3¥M.QšÓÛÛð°äM94Ô¼ðèVØuèˆÜÕAk,OBu´ŽÔ1˜f†3"‰ß~†è8ÃÁ <´@ìÑ¡ R&nþœä?HII{§(ŒÚ&…ľÿФõ }ÝÊ~C
-Ø r-ôX%ÒaSXNˆ,´§Jh‡'»^«
-úEÕuAˆ¶¢(N†˜SxÊ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׫™Æ ®²4ŽsÈR á•5ƒíZUw¤’v¡¸)ž…yN•N®IU X©E5™4î*ÂgعéNvî¨åZª³© ß/ò𸪾ãzUËFL½j¡-±q2?£
-›‡¦åA=ë c ¼ !Ô;ÈÇ’noBƒÂ!CÆ„ „Z
- òÓÀôN\ü^ŸÐiA aâñCñæu¯»ˆAÿ’-Áe=†ôc¥òÔö¯¸ÈÕ4â7«16§º ý<- U¥½4ZÕ[ÙÐ ·¤¢w†hvöŒêÀø¦³µÉ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 ]äŒMž—oáÏÎÍÚù†ÞgÖu CÀ˜xŹºéj9›s üx6÷ñôì»È«ŒØeŠ£?mU®ŠÃKËš˜ëZ’\Svu"
+§?]®‰È¢•$mT°Š&?»®W7™5
+H¢V¢¸±&¨ÆH¬$.(䊄Jþ¶‹TNÔƒ¬Klæào Æß˲•xˆh¹h‘â@¡~äÙñò5ÃÒL«Â0od!k‘ÙÕÒÌ®gÁô¾*ÑËaì8DÀ«í ´mÓŠÚî«ÍE&â²¥ÕBÌÏnà梾…œ¢çÈ
+©Ðê@‚ìÜcvkí¸Vi;†E%j‘K< Åáb3Š=‡¡
+âÖ tæ;^ĬˆF‰G «ó|'XôwâýëïNÏFTòÈaq¯Ï# ÕÞ¡÷e—¡ÿtâz7†èœ¸ˆ "ÎÓÖ³72Æ2ß 3pÑHÅÊΒʵÈÒùª«E«Jòi
+¢‚Ö²©œO¹»pXÄé(rx=¸¤pðjî Vm˜vÐìc'ì(u¯´ÛÄò ËÛ`qÃ↓šptm ¹öÒš«Ë* UèªLZ5ó‚é· ¦2»Ÿ…ÁÔ\0߉!
+ ø ×'3/
+RŽ« Dƒ ß /HA®oøåW#o¯·ùÊË•Ìþý2‚a:œ÷ Ѐ ³ ÐàJÔ¤RX‚Ò’UÙ4ê:³÷•EÃ\JkÇñÈYÛû’Û#fN<Äý±Q¨/4*½*ºüZÒ‘0{]âcE#™ãn
+o_ŠªDçݘÈø€7…ü»÷scð€ûb8x’gÞŽg c–à–ˆà–Ã-'þ?çÈè ïF݃@âÞ“Üã½{MOQèTkE^˜áöûwK<ñ@ÄKZü2µ#©—°ØY œ¸Éù&êsoáD1yuŽ÷Tih\TÓÖe¥übp{EÓªäH$e‡„‘Éä7… ¿ºf,‹³0€KYTÆp‹ |ì¼Î¨‘ÙcºpÃ7–éMšaꜙ†y³‘z
+3£ÁÌZ? iqÆ©5A~¡š’ 'œ Ãq£íã'4 WÆö¸
+µ0îú.áQfJ ³Y8ƒ=š_•u¹µýp8½ï N8ÔÈ Doo²c´›Ö ç²YV"
+ÃJÊ¢ºP¢á ‹,Bú+ ¼ÒÖÏ
+²žïB>·ÓaïYŸC©*¨¼ËZQHØÅ89}o®HæÄV§Llm(ÕÞõ1¸.ÝZe„wOXÙßÏ[ êrŸl[š¾FF·W¼’uno;-°ŸI™We
+ãvÁR\« 5š Àr“ÈvKÖþ°ýáðawƒæ_ì°aÿÃööó¿ØOÚo±zí?ìë™ñŸÝöÁÓe]³_O:vì?«£&%)öz?©rô"ð›j’ÇBG™wW6?¶ YKQÿãö}pòhÕ×N ç‘bjI;œÉfÉÀ§Ñ1¾"<Û'¶O wr0l ,ЃÆ0• ëZBG¥“6ü'ŠV?OCYB©À¦2Z¢{m÷í¨m)ë{3ÙH’zUÖ²HP+͘üÄrïÀBu ±OžëF_Æú÷9’ ¤aôèˆnd{µaÎÕ7™•‰j<vÓx´Á v‡B˜e%âp‡¾†±.v8âEêô{?ª‹Ãž²¤ñy´4튄 &𥯠*\“’þÝôª„÷:@OÅ/Á}¶
+CZ˯]ŸÉÄaÈ1”~‰‡ z'ó0¢Ê
+„yIád"²¤Ëè©„SwÒŒX¼´„̓È4
+ÖIb½fƒ-s´¤þÅ@oIú†_ €Ûÿb 4–ÿx.ìÒØÐ íœü©è?¸Š.í}IH›mž.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]o›0}ϯà¤áÚÆ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;Ò @ÈÐ äÁÙÕæû[Øsi–o5„öW¥›²3v] çÉí®â¬•9çõ–ËÞk±¹9j÷s$8Ä Ò“ÛA¥êDSðL´…¬Ä®( ,pC\Æ-ÑÔQ“Ó™»{9yÌÆŸÕu“KÅ;ñ' Z‘ X ‘ä•9\W\Yÿ–¥Ë= ÝQ×z‹¦Vv·²+jè—±kŠ‡™°Ì` bôÿÄ´õ®Éįîþzl,ø°²w¦7²ÃÔ•Ñ3l‹qˆoŒ¿š¿ºúg`¼À’ßW{çZšBB„-Mètž¸ÚTBwÌFi ³á¯‘qõôéiXz‚G]ý|<LŸl¼§Rå2”ÆOÔ‹T™ÌE/ží‘dšÃ‹^C
-¥j;Áó´ªv:
-í!N’tòK&`"
+xÚ•VQo¤6~ϯØGV
+>0Æ@¥{hÓMÓª—‹’½S¥¦Š¼`«¬½5ÐÜö×wŒ
+ÙÝøšë“Íx<óÍçÏc~X_¼»¦ù¢@Åt±®^dI¾XW‹ßƒ«›ïïÖ«ûeˆÓ( h¦4
+î—$
+>ÝÞþ|û“]X߬ìäáã¯Ë8Š‚ÏÓ–ß îVW׫øÇå ë_Þ]r”ˆ¤9Š’lÙ\qî|’ELPB(6Náäâ v&Ö— ² !SÜ(È<Ût½²cÇ{;Ù·Lr;}fKœoÞYs£¢¬µ ¥’|*Ùy G #L©…±^æQ 9“´:±ZÖsûµ×j϶Ë4`&”µ©ÚŽB–j'äÖ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ïFUIO–AVL ¬?À.'(L;¨šÿ5Í+ä+t:Êù {§«FtÓ’dZ¨ï¬6hv,gTä1 vú¸cú©-÷‰?C9Ž'íñ`/èšPÃ{Ä8…{RïžE_6œžA&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Ä …xˆGúºÛý¯›9¹9—fj‡S‹kܤc;~Ú9祆‹mÓÏœí;êWæÛö+7çQ„gò„O×8Ã@Û‘°{®kü+w+OÍçºìu“ ÄM<â~ƒ¹/o3'ÎÙø–JéJÈñ=ðr2Ññò…Ëѳtü’µ\kµ’î œþcÕšO‡1èz PrÇ ¼hEüÿ‰éÔ Mö©?ì}g at s„_Žì½Í ;0cÓ62zUÛ¿Dwñ›û >ÿg€Ä àž Z#`ÿ^GSFQcGSüí<1¹mùø ùh2ó/PÉäéÕƒ² ûpÿ/ó‰O÷T@ -g¥±oÔ ü…ˆŠO’õìÎ$£_nôc”FBv=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