[cig-commits] r19346 - 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/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 src/meshfem2D src/specfem2D
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Mon Jan 9 15:34:10 PST 2012
Author: dkomati1
Date: 2012-01-09 15:34:09 -0800 (Mon, 09 Jan 2012)
New Revision: 19346
Modified:
seismo/2D/SPECFEM2D/trunk/DATA/Par_file
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/Par_file
seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/Par_file_Tape2007_onerec
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/Par_file_Tromp2005
seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/Par_file_acoustic_poroelastic
seismo/2D/SPECFEM2D/trunk/EXAMPLES/attenuation/Par_file_attenuation_2D
seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon
seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave
seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_for
seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_kernel
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_fair
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_good
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_best
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_fair
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_good
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_fair
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_good
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_1
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_2
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_3
seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90
seismo/2D/SPECFEM2D/trunk/src/meshfem2D/save_databases.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/checkgrid.F90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_fluid.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_solid.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/get_poroelastic_velocities.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/gmat01.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
Log:
changed ATTENUATION_POROELASTIC_SOLID to ATTENUATION_PORO_FLUID_PART to further clarify
Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external model from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in 2012-01-09 23:34:09 UTC (rev 19346)
@@ -18,7 +18,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in 2012-01-09 23:34:09 UTC (rev 19346)
@@ -18,7 +18,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/Par_file
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/Par_file 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/Par_file 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external model from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .true. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .true. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .false. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .false. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/Par_file_Tape2007_onerec
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/Par_file_Tape2007_onerec 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/Par_file_Tape2007_onerec 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .false. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/Par_file_Tromp2005
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/Par_file_Tromp2005 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/Par_file_Tromp2005 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/Par_file_acoustic_poroelastic
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/Par_file_acoustic_poroelastic 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/Par_file_acoustic_poroelastic 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/attenuation/Par_file_attenuation_2D
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/attenuation/Par_file_attenuation_2D 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/attenuation/Par_file_attenuation_2D 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .true. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external model from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external model from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_for
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_for 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_for 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_kernel
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_kernel 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_kernel 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_fair
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_fair 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_fair 2012-01-09 23:34:09 UTC (rev 19346)
@@ -16,7 +16,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .false. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_good
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_good 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_good 2012-01-09 23:34:09 UTC (rev 19346)
@@ -16,7 +16,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .false. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_best
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_best 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_best 2012-01-09 23:34:09 UTC (rev 19346)
@@ -16,7 +16,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .false. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_fair
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_fair 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_fair 2012-01-09 23:34:09 UTC (rev 19346)
@@ -16,7 +16,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .false. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_good
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_good 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_good 2012-01-09 23:34:09 UTC (rev 19346)
@@ -16,7 +16,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .false. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_fair
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_fair 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_fair 2012-01-09 23:34:09 UTC (rev 19346)
@@ -16,7 +16,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .false. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_good
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_good 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_good 2012-01-09 23:34:09 UTC (rev 19346)
@@ -16,7 +16,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .false. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_1
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_1 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_1 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .false. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_2
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_2 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_2 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .false. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_3
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_3 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_3 2012-01-09 23:34:09 UTC (rev 19346)
@@ -17,7 +17,7 @@
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
-ATTENUATION_POROELASTIC_SOLID = .false. # turn viscous attenuation on or off for poroelastic solid parts of the model
+ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
freq0 = 10 # frequency for viscous attenuation
p_sv = .false. # set the type of calculation (P-SV or SH/membrane waves)
Modified: seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90 2012-01-09 23:34:09 UTC (rev 19346)
@@ -68,7 +68,7 @@
integer :: nx,ngnod
logical :: initialfield,add_Bielak_conditions,assign_external_model, &
- READ_EXTERNAL_SEP_FILE,ATTENUATION_VISCOELASTIC_SOLID,ATTENUATION_POROELASTIC_SOLID
+ READ_EXTERNAL_SEP_FILE,ATTENUATION_VISCOELASTIC_SOLID,ATTENUATION_PORO_FLUID_PART
double precision :: Q0,freq0
@@ -211,7 +211,7 @@
if(err_occurred() /= 0) stop 'error reading parameter 11 in Par_file'
! read viscous attenuation parameters (poroelastic media)
- call read_value_logical_p(ATTENUATION_POROELASTIC_SOLID, 'solver.ATTENUATION_POROELASTIC_SOLID')
+ call read_value_logical_p(ATTENUATION_PORO_FLUID_PART, 'solver.ATTENUATION_PORO_FLUID_PART')
if(err_occurred() /= 0) stop 'error reading parameter 12 in Par_file'
call read_value_double_precision_p(Q0, 'solver.Q0')
Modified: seismo/2D/SPECFEM2D/trunk/src/meshfem2D/save_databases.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/meshfem2D/save_databases.f90 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/src/meshfem2D/save_databases.f90 2012-01-09 23:34:09 UTC (rev 19346)
@@ -130,8 +130,8 @@
write(15,*) 'output_grid output_energy output_wavefield_snapshot'
write(15,*) output_grid,output_energy,output_wavefield_snapshot
- write(15,*) 'ATTENUATION_VISCOELASTIC_SOLID ATTENUATION_POROELASTIC_SOLID'
- write(15,*) ATTENUATION_VISCOELASTIC_SOLID,ATTENUATION_POROELASTIC_SOLID
+ write(15,*) 'ATTENUATION_VISCOELASTIC_SOLID ATTENUATION_PORO_FLUID_PART'
+ write(15,*) ATTENUATION_VISCOELASTIC_SOLID,ATTENUATION_PORO_FLUID_PART
write(15,*) 'Q0 freq0'
write(15,*) Q0,freq0
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/checkgrid.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/checkgrid.F90 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/checkgrid.F90 2012-01-09 23:34:09 UTC (rev 19346)
@@ -50,7 +50,7 @@
coorg,xinterp,zinterp,shapeint,knods,simulation_title, &
npgeo,pointsdisp,ngnod,any_elastic,any_poroelastic,all_anisotropic, &
myrank,nproc,NSOURCES,poroelastic, &
- freq0,Q0,ATTENUATION_POROELASTIC_SOLID,US_LETTER,output_postscript_snapshot)
+ freq0,Q0,ATTENUATION_PORO_FLUID_PART,US_LETTER,output_postscript_snapshot)
! check the mesh, stability and number of points per wavelength
@@ -102,7 +102,7 @@
double precision :: deltat
logical :: assign_external_model,initialfield,any_elastic,any_poroelastic,all_anisotropic, &
- ATTENUATION_POROELASTIC_SOLID,output_postscript_snapshot
+ ATTENUATION_PORO_FLUID_PART,output_postscript_snapshot
integer :: myrank,nproc
@@ -276,7 +276,7 @@
M_biot = kappa_s*kappa_s/(D_biot - kappa_fr)
call get_poroelastic_velocities(cpIsquare,cpIIsquare,cssquare,H_biot,C_biot,M_biot,mu_fr,phi, &
- tort,denst_s,denst_f,eta_f,perm,f0(1),freq0,Q0,w_c,ATTENUATION_POROELASTIC_SOLID)
+ tort,denst_s,denst_f,eta_f,perm,f0(1),freq0,Q0,w_c,ATTENUATION_PORO_FLUID_PART)
cpIloc = sqrt(cpIsquare)
cpIIloc = sqrt(cpIIsquare)
@@ -592,7 +592,7 @@
M_biot = kappa_s*kappa_s/(D_biot - kappa_fr)
call get_poroelastic_velocities(cpIsquare,cpIIsquare,cssquare,H_biot,C_biot,M_biot,mu_fr,phi, &
- tort,denst_s,denst_f,eta_f,perm,f0(1),freq0,Q0,w_c,ATTENUATION_POROELASTIC_SOLID)
+ tort,denst_s,denst_f,eta_f,perm,f0(1),freq0,Q0,w_c,ATTENUATION_PORO_FLUID_PART)
cpIloc = sqrt(cpIsquare)
cpIIloc = sqrt(cpIIsquare)
@@ -1076,7 +1076,7 @@
M_biot = kappa_s*kappa_s/(D_biot - kappa_fr)
call get_poroelastic_velocities(cpIsquare,cpIIsquare,cssquare,H_biot,C_biot,M_biot,mu_fr,phi, &
- tort,denst_s,denst_f,eta_f,perm,f0(1),freq0,Q0,w_c,ATTENUATION_POROELASTIC_SOLID)
+ tort,denst_s,denst_f,eta_f,perm,f0(1),freq0,Q0,w_c,ATTENUATION_PORO_FLUID_PART)
cpIloc = sqrt(cpIsquare)
else
@@ -1430,7 +1430,7 @@
M_biot = kappa_s*kappa_s/(D_biot - kappa_fr)
call get_poroelastic_velocities(cpIsquare,cpIIsquare,cssquare,H_biot,C_biot,M_biot,mu_fr,phi, &
- tort,denst_s,denst_f,eta_f,perm,f0(1),freq0,Q0,w_c,ATTENUATION_POROELASTIC_SOLID)
+ tort,denst_s,denst_f,eta_f,perm,f0(1),freq0,Q0,w_c,ATTENUATION_PORO_FLUID_PART)
cpIloc = sqrt(cpIsquare)
csloc = sqrt(cssquare)
@@ -1848,7 +1848,7 @@
M_biot = kappa_s*kappa_s/(D_biot - kappa_fr)
call get_poroelastic_velocities(cpIsquare,cpIIsquare,cssquare,H_biot,C_biot,M_biot,mu_fr,phi, &
- tort,denst_s,denst_f,eta_f,perm,f0(1),freq0,Q0,w_c,ATTENUATION_POROELASTIC_SOLID)
+ tort,denst_s,denst_f,eta_f,perm,f0(1),freq0,Q0,w_c,ATTENUATION_PORO_FLUID_PART)
cpIloc = sqrt(cpIsquare)
else
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_fluid.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_fluid.f90 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_fluid.f90 2012-01-09 23:34:09 UTC (rev 19346)
@@ -44,7 +44,7 @@
subroutine compute_forces_poro_fluid(nglob,nspec,myrank,nelemabs,numat, &
ispec_selected_source,ispec_selected_rec,is_proc_source,which_proc_receiver,&
source_type,it,NSTEP,anyabs, &
- initialfield,ATTENUATION_VISCOELASTIC_SOLID,ATTENUATION_POROELASTIC_SOLID,deltatcube, &
+ initialfield,ATTENUATION_VISCOELASTIC_SOLID,ATTENUATION_PORO_FLUID_PART,deltatcube, &
deltatfourth,twelvedeltat,fourdeltatsquare,ibool,kmato,numabs,poroelastic,codeabs, &
accelw_poroelastic,velocw_poroelastic,displw_poroelastic,velocs_poroelastic,displs_poroelastic,&
b_accelw_poroelastic,b_displw_poroelastic,b_displs_poroelastic,&
@@ -123,7 +123,7 @@
double precision, dimension(NGLLX,NGLLZ,nspec) :: rx_viscous
double precision, dimension(NGLLX,NGLLZ,nspec) :: rz_viscous
double precision :: theta_e,theta_s
- logical ATTENUATION_POROELASTIC_SOLID
+ logical ATTENUATION_PORO_FLUID_PART
double precision, dimension(3):: bl_unrelaxed,bl_relaxed
! derivatives of Lagrange polynomials
@@ -475,7 +475,7 @@
bl_relaxed(2) = etal_f*invpermlxz
bl_relaxed(3) = etal_f*invpermlzz
- if(ATTENUATION_POROELASTIC_SOLID) then
+ if(ATTENUATION_PORO_FLUID_PART) then
bl_unrelaxed(1) = etal_f*invpermlxx*theta_e/theta_s
bl_unrelaxed(2) = etal_f*invpermlxz*theta_e/theta_s
bl_unrelaxed(3) = etal_f*invpermlzz*theta_e/theta_s
@@ -486,7 +486,7 @@
iglob = ibool(i,j,ispec)
- if(ATTENUATION_POROELASTIC_SOLID) then
+ if(ATTENUATION_PORO_FLUID_PART) then
! compute the viscous damping term with the unrelaxed viscous coef and add memory variable
viscodampx = velocw_poroelastic(1,iglob)*bl_unrelaxed(1) + velocw_poroelastic(2,iglob)*bl_unrelaxed(2)&
- rx_viscous(i,j,ispec)
@@ -552,7 +552,7 @@
M_biot = kappal_s*kappal_s/(D_biot - kappal_fr)
call get_poroelastic_velocities(cpIsquare,cpIIsquare,cssquare,H_biot,C_biot,M_biot,mul_fr,phil, &
- tortl,rhol_s,rhol_f,etal_f,permlxx,f0,freq0,Q0,w_c,ATTENUATION_POROELASTIC_SOLID)
+ tortl,rhol_s,rhol_f,etal_f,permlxx,f0,freq0,Q0,w_c,ATTENUATION_PORO_FLUID_PART)
cpIl = sqrt(cpIsquare)
cpIIl = sqrt(cpIIsquare)
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_solid.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_solid.f90 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_solid.f90 2012-01-09 23:34:09 UTC (rev 19346)
@@ -44,7 +44,7 @@
subroutine compute_forces_poro_solid(nglob,nspec,myrank,nelemabs,numat, &
ispec_selected_source,ispec_selected_rec,is_proc_source,which_proc_receiver,&
source_type,it,NSTEP,anyabs, &
- initialfield,ATTENUATION_VISCOELASTIC_SOLID,ATTENUATION_POROELASTIC_SOLID,deltatcube, &
+ initialfield,ATTENUATION_VISCOELASTIC_SOLID,ATTENUATION_PORO_FLUID_PART,deltatcube, &
deltatfourth,twelvedeltat,fourdeltatsquare,ibool,kmato,numabs,poroelastic,codeabs, &
accels_poroelastic,velocs_poroelastic,velocw_poroelastic,displs_poroelastic,displw_poroelastic,&
b_accels_poroelastic,b_displs_poroelastic,b_displw_poroelastic,&
@@ -123,7 +123,7 @@
double precision, dimension(NGLLX,NGLLZ,nspec) :: rx_viscous
double precision, dimension(NGLLX,NGLLZ,nspec) :: rz_viscous
double precision :: theta_e,theta_s
- logical ATTENUATION_POROELASTIC_SOLID
+ logical ATTENUATION_PORO_FLUID_PART
double precision, dimension(3):: bl_unrelaxed,bl_relaxed
! derivatives of Lagrange polynomials
@@ -495,7 +495,7 @@
bl_relaxed(2) = etal_f*invpermlxz
bl_relaxed(3) = etal_f*invpermlzz
- if(ATTENUATION_POROELASTIC_SOLID) then
+ if(ATTENUATION_PORO_FLUID_PART) then
bl_unrelaxed(1) = etal_f*invpermlxx*theta_e/theta_s
bl_unrelaxed(2) = etal_f*invpermlxz*theta_e/theta_s
bl_unrelaxed(3) = etal_f*invpermlzz*theta_e/theta_s
@@ -506,7 +506,7 @@
iglob = ibool(i,j,ispec)
- if(ATTENUATION_POROELASTIC_SOLID) then
+ if(ATTENUATION_PORO_FLUID_PART) then
! compute the viscous damping term with the unrelaxed viscous coef and add memory variable
viscodampx = velocw_poroelastic(1,iglob)*bl_unrelaxed(1) + velocw_poroelastic(2,iglob)*bl_unrelaxed(2)&
- rx_viscous(i,j,ispec)
@@ -571,7 +571,7 @@
M_biot = kappal_s*kappal_s/(D_biot - kappal_fr)
call get_poroelastic_velocities(cpIsquare,cpIIsquare,cssquare,H_biot,C_biot,M_biot,mul_fr,phil, &
- tortl,rhol_s,rhol_f,etal_f,permlxx,f0,freq0,Q0,w_c,ATTENUATION_POROELASTIC_SOLID)
+ tortl,rhol_s,rhol_f,etal_f,permlxx,f0,freq0,Q0,w_c,ATTENUATION_PORO_FLUID_PART)
cpIl = sqrt(cpIsquare)
cpIIl = sqrt(cpIIsquare)
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/get_poroelastic_velocities.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/get_poroelastic_velocities.f90 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/get_poroelastic_velocities.f90 2012-01-09 23:34:09 UTC (rev 19346)
@@ -46,7 +46,7 @@
!----
subroutine get_poroelastic_velocities(cpIsquare,cpIIsquare,cssquare,H_biot,C_biot,M_biot,mul_fr,phil, &
- tortl,rhol_s,rhol_f,etal_f,perm,fi,f0,Q0,w_c,ATTENUATION_POROELASTIC_SOLID)
+ tortl,rhol_s,rhol_f,etal_f,perm,fi,f0,Q0,w_c,ATTENUATION_PORO_FLUID_PART)
implicit none
@@ -64,7 +64,7 @@
double precision :: wi,fi,taus,taue,Q0,bbr,bbi
double precision :: gA,gB,sa,sb,xxs,yys
- logical :: ATTENUATION_POROELASTIC_SOLID
+ logical :: ATTENUATION_PORO_FLUID_PART
rhol_bar = (1.d0 - phil)*rhol_s + phil*rhol_f
@@ -78,7 +78,7 @@
taue = (sqrt(Q0*Q0+1) +1)/(w0il*Q0)
taus = (sqrt(Q0*Q0+1) -1)/(w0il*Q0)
- if(ATTENUATION_POROELASTIC_SOLID) then
+ if(ATTENUATION_PORO_FLUID_PART) then
! high frequency, with memory variables
bbr = etal_f/perm*(1.d0+alpha*alpha*taus*taue)/(1.d0 + alpha*alpha*taus*taus)
bbi = etal_f/perm*alpha*(taue-taus)/(1.d0 + alpha*alpha*taus*taus)
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/gmat01.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/gmat01.f90 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/gmat01.f90 2012-01-09 23:34:09 UTC (rev 19346)
@@ -45,7 +45,7 @@
subroutine gmat01(density_array,porosity_array,tortuosity_array, &
aniso_array,permeability,poroelastcoef, &
numat,myrank,ipass,QKappa_array,Qmu_array, &
- freq0,Q0,f0,ATTENUATION_POROELASTIC_SOLID)
+ freq0,Q0,f0,ATTENUATION_PORO_FLUID_PART)
! reads properties of a 2D isotropic or anisotropic linear elastic element
@@ -57,7 +57,7 @@
double precision :: aniso_array(6,numat),tortuosity_array(numat),permeability(3,numat)
double precision :: QKappa_array(numat),Qmu_array(numat)
double precision :: f0,Q0,freq0
- logical :: ATTENUATION_POROELASTIC_SOLID
+ logical :: ATTENUATION_PORO_FLUID_PART
! local parameters
double precision :: lambdaplus2mu,kappa
@@ -199,7 +199,7 @@
call get_poroelastic_velocities(cpIsquare,cpIIsquare,cssquare, &
H_biot,C_biot,M_biot,mu_fr,phi, &
tortuosity,density(1),density(2),eta_f, &
- val4,f0,freq0,Q0,w_c,ATTENUATION_POROELASTIC_SOLID)
+ val4,f0,freq0,Q0,w_c,ATTENUATION_PORO_FLUID_PART)
porosity_array(n) = val2
tortuosity_array(n) = val3
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90 2012-01-09 23:34:09 UTC (rev 19346)
@@ -52,7 +52,7 @@
anglerec,initialfield,add_Bielak_conditions, &
seismotype,imagetype,assign_external_model,READ_EXTERNAL_SEP_FILE, &
output_grid,output_energy,output_wavefield_snapshot,ATTENUATION_VISCOELASTIC_SOLID, &
- ATTENUATION_POROELASTIC_SOLID,Q0,freq0,p_sv, &
+ ATTENUATION_PORO_FLUID_PART,Q0,freq0,p_sv, &
NSTEP,deltat,NTSTEP_BETWEEN_OUTPUT_SEISMO,NSOURCES, &
factor_subsample_image,USE_SNAPSHOT_NUMBER_IN_FILENAME,DRAW_WATER_CONSTANT_BLUE_IN_JPG,US_LETTER, &
POWER_DISPLAY_COLOR,PERFORM_CUTHILL_MCKEE,SU_FORMAT,USER_T0, &
@@ -72,7 +72,7 @@
logical :: meshvect,modelvect,boundvect,initialfield,add_Bielak_conditions, &
assign_external_model,READ_EXTERNAL_SEP_FILE, &
output_grid,output_energy,output_wavefield_snapshot,p_sv
- logical :: ATTENUATION_VISCOELASTIC_SOLID,ATTENUATION_POROELASTIC_SOLID
+ logical :: ATTENUATION_VISCOELASTIC_SOLID,ATTENUATION_PORO_FLUID_PART
double precision :: cutsnaps,sizemax_arrows,anglerec
double precision :: Q0,freq0
@@ -197,7 +197,7 @@
read(IIN,*) output_grid,output_energy,output_wavefield_snapshot
read(IIN,"(a80)") datlin
- read(IIN,*) ATTENUATION_VISCOELASTIC_SOLID,ATTENUATION_POROELASTIC_SOLID
+ read(IIN,*) ATTENUATION_VISCOELASTIC_SOLID,ATTENUATION_PORO_FLUID_PART
read(IIN,"(a80)") datlin
read(IIN,*) Q0,freq0
@@ -255,7 +255,7 @@
if (myrank == 0 .and. ipass == 1) write(IOUT,703) NSTEP,deltat,NSTEP*deltat
if( SIMULATION_TYPE == 1 .and. SAVE_FORWARD .and. &
- (ATTENUATION_VISCOELASTIC_SOLID .or. ATTENUATION_POROELASTIC_SOLID) ) then
+ (ATTENUATION_VISCOELASTIC_SOLID .or. ATTENUATION_PORO_FLUID_PART) ) then
print*, '*************** WARNING ***************'
print*, 'Anisotropy & Attenuation & Viscous damping are not presently implemented for adjoint calculations'
stop
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90 2012-01-09 22:20:27 UTC (rev 19345)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90 2012-01-09 23:34:09 UTC (rev 19346)
@@ -535,7 +535,7 @@
double precision :: theta_e,theta_s
double precision :: Q0,freq0
double precision :: alphaval,betaval,gammaval,thetainv
- logical :: ATTENUATION_POROELASTIC_SOLID
+ logical :: ATTENUATION_PORO_FLUID_PART
double precision, dimension(NGLLX,NGLLZ) :: viscox_loc,viscoz_loc
double precision :: Sn,Snp1,etal_f
double precision, dimension(3):: bl_relaxed
@@ -908,7 +908,7 @@
anglerec,initialfield,add_Bielak_conditions, &
seismotype,imagetype,assign_external_model,READ_EXTERNAL_SEP_FILE, &
output_grid,output_energy,output_wavefield_snapshot,ATTENUATION_VISCOELASTIC_SOLID, &
- ATTENUATION_POROELASTIC_SOLID,Q0,freq0,p_sv, &
+ ATTENUATION_PORO_FLUID_PART,Q0,freq0,p_sv, &
NSTEP,deltat,NTSTEP_BETWEEN_OUTPUT_SEISMO,NSOURCES, &
factor_subsample_image,USE_SNAPSHOT_NUMBER_IN_FILENAME,DRAW_WATER_CONSTANT_BLUE_IN_JPG,US_LETTER, &
POWER_DISPLAY_COLOR,PERFORM_CUTHILL_MCKEE,SU_FORMAT,USER_T0, &
@@ -940,7 +940,7 @@
anglerec,initialfield,add_Bielak_conditions, &
seismotype,imagetype,assign_external_model,READ_EXTERNAL_SEP_FILE, &
output_grid,output_energy,output_wavefield_snapshot,ATTENUATION_VISCOELASTIC_SOLID, &
- ATTENUATION_POROELASTIC_SOLID,Q0,freq0,p_sv, &
+ ATTENUATION_PORO_FLUID_PART,Q0,freq0,p_sv, &
NSTEP,deltat,NTSTEP_BETWEEN_OUTPUT_SEISMO,NSOURCES, &
factor_subsample_image,USE_SNAPSHOT_NUMBER_IN_FILENAME,DRAW_WATER_CONSTANT_BLUE_IN_JPG,US_LETTER, &
POWER_DISPLAY_COLOR,PERFORM_CUTHILL_MCKEE,SU_FORMAT,USER_T0, &
@@ -1054,7 +1054,7 @@
!---- read the material properties
!
call gmat01(density,porosity,tortuosity,anisotropy,permeability,poroelastcoef,numat,&
- myrank,ipass,QKappa_attenuation,Qmu_attenuation,freq0,Q0,f0(1),ATTENUATION_POROELASTIC_SOLID)
+ myrank,ipass,QKappa_attenuation,Qmu_attenuation,freq0,Q0,f0(1),ATTENUATION_PORO_FLUID_PART)
!
!---- read spectral macrobloc data
!
@@ -1198,7 +1198,7 @@
! allocate memory variables for viscous attenuation (poroelastic media)
if(ipass == 1) then
- if(ATTENUATION_POROELASTIC_SOLID) then
+ if(ATTENUATION_PORO_FLUID_PART) then
allocate(rx_viscous(NGLLX,NGLLZ,nspec))
allocate(rz_viscous(NGLLX,NGLLZ,nspec))
allocate(viscox(NGLLX,NGLLZ,nspec))
@@ -2805,7 +2805,7 @@
coorg,xinterp,zinterp,shape2D_display,knods,simulation_title, &
npgeo,pointsdisp,ngnod,any_elastic,any_poroelastic,all_anisotropic, &
myrank,nproc,NSOURCES,poroelastic, &
- freq0,Q0,ATTENUATION_POROELASTIC_SOLID,US_LETTER,output_postscript_snapshot)
+ freq0,Q0,ATTENUATION_PORO_FLUID_PART,US_LETTER,output_postscript_snapshot)
! convert receiver angle to radians
anglerec = anglerec * pi / 180.d0
@@ -3553,7 +3553,7 @@
! the common nodes forming the edge are computed here
if(coupled_elastic_poro) then
- if(ATTENUATION_VISCOELASTIC_SOLID .or. ATTENUATION_POROELASTIC_SOLID) &
+ if(ATTENUATION_VISCOELASTIC_SOLID .or. ATTENUATION_PORO_FLUID_PART) &
stop 'Attenuation not supported for mixed elastic/poroelastic simulations'
if ( myrank == 0 ) then
@@ -3912,7 +3912,7 @@
seismo_current = 0
! Precompute Runge Kutta coefficients if viscous attenuation
- if(ATTENUATION_POROELASTIC_SOLID) then
+ if(ATTENUATION_PORO_FLUID_PART) then
theta_e = (sqrt(Q0**2+1.d0) +1.d0)/(2.d0*pi*freq0*Q0)
theta_s = (sqrt(Q0**2+1.d0) -1.d0)/(2.d0*pi*freq0*Q0)
@@ -4166,7 +4166,7 @@
!--------------------------------------------------------------------------------------------
! implement viscous attenuation for poroelastic media
!
- if(ATTENUATION_POROELASTIC_SOLID .and. any_poroelastic) then
+ if(ATTENUATION_PORO_FLUID_PART .and. any_poroelastic) then
! update memory variables with fourth-order Runge-Kutta time scheme for attenuation
! loop over spectral elements
do ispec = 1,nspec
@@ -5320,7 +5320,7 @@
call compute_forces_poro_solid(nglob,nspec,myrank,nelemabs,numat, &
ispec_selected_source,ispec_selected_rec,is_proc_source,which_proc_receiver,&
source_type,it,NSTEP,anyabs, &
- initialfield,ATTENUATION_VISCOELASTIC_SOLID,ATTENUATION_POROELASTIC_SOLID,deltatcube, &
+ initialfield,ATTENUATION_VISCOELASTIC_SOLID,ATTENUATION_PORO_FLUID_PART,deltatcube, &
deltatfourth,twelvedeltat,fourdeltatsquare,ibool,kmato,numabs,poroelastic,codeabs, &
accels_poroelastic,velocs_poroelastic,velocw_poroelastic,displs_poroelastic,displw_poroelastic,&
b_accels_poroelastic,b_displs_poroelastic,b_displw_poroelastic,&
@@ -5343,7 +5343,7 @@
call compute_forces_poro_fluid(nglob,nspec,myrank,nelemabs,numat, &
ispec_selected_source,ispec_selected_rec,is_proc_source,which_proc_receiver,&
source_type,it,NSTEP,anyabs, &
- initialfield,ATTENUATION_VISCOELASTIC_SOLID,ATTENUATION_POROELASTIC_SOLID,deltatcube, &
+ initialfield,ATTENUATION_VISCOELASTIC_SOLID,ATTENUATION_PORO_FLUID_PART,deltatcube, &
deltatfourth,twelvedeltat,fourdeltatsquare,ibool,kmato,numabs,poroelastic,codeabs, &
accelw_poroelastic,velocw_poroelastic,displw_poroelastic,velocs_poroelastic,displs_poroelastic,&
b_accelw_poroelastic,b_displw_poroelastic,b_displs_poroelastic,&
More information about the CIG-COMMITS
mailing list