[cig-commits] [commit] devel: updates stations and source in CPML acoustic examples; out-comments unused array potential_dot_dot_acoustic_old (eeec5ae)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Fri Nov 28 14:14:16 PST 2014
Repository : https://github.com/geodynamics/specfem3d
On branch : devel
Link : https://github.com/geodynamics/specfem3d/compare/9a84d06c76e869f5276019e4f84affce23830a4d...8dce71b713c1fd0b510e7538fea0f2307c7b29e8
>---------------------------------------------------------------
commit eeec5aef347463ba75854408a1970e3b6e85f024
Author: daniel peter <peterda at ethz.ch>
Date: Fri Nov 28 13:50:33 2014 +0100
updates stations and source in CPML acoustic examples; out-comments unused array potential_dot_dot_acoustic_old
>---------------------------------------------------------------
eeec5aef347463ba75854408a1970e3b6e85f024
.../DATA/CMTSOLUTION | 12 ++++++------
.../DATA/STATIONS | 8 ++++----
.../DATA/CMTSOLUTION | 12 ++++++------
.../DATA/STATIONS | 8 ++++----
.../DATA/STATIONS | 8 ++++----
.../DATA/STATIONS | 8 ++++----
src/specfem3D/compute_forces_acoustic_calling_routine.f90 | 4 ++--
src/specfem3D/pml_allocate_arrays.f90 | 8 ++++----
src/specfem3D/pml_par.f90 | 2 +-
9 files changed, 35 insertions(+), 35 deletions(-)
diff --git a/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_5sides/DATA/CMTSOLUTION b/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_5sides/DATA/CMTSOLUTION
index fca749f..6c3e90c 100644
--- a/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_5sides/DATA/CMTSOLUTION
+++ b/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_5sides/DATA/CMTSOLUTION
@@ -5,9 +5,9 @@ half duration: 5.0
latorUTM: 67000.0
longorUTM: 67000.0
depth: 30.0
-Mrr: 1.000000e+24
-Mtt: 1.000000e+24
-Mpp: 1.000000e+24
-Mrt: 0.
-Mrp: 0.
-Mtp: 0.
+Mrr: 7.700000e+27
+Mtt: 7.700000e+27
+Mpp: 7.700000e+27
+Mrt: 0.0
+Mrp: 0.0
+Mtp: 0.0
diff --git a/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_5sides/DATA/STATIONS b/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_5sides/DATA/STATIONS
index 1ee47ff..0f18aea 100644
--- a/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_5sides/DATA/STATIONS
+++ b/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_5sides/DATA/STATIONS
@@ -1,4 +1,4 @@
-X20 DB 67000.00 22732.14 0.0 0.0
-X30 DB 67000.00 34696.43 0.0 0.0
-X40 DB 67000.00 46660.71 0.0 0.0
-X50 DB 67000.00 58625.00 0.0 0.0
+X20 DB 67000.00 22732.14 0.0 50.0
+X30 DB 67000.00 34696.43 0.0 50.0
+X40 DB 67000.00 46660.71 0.0 50.0
+X50 DB 67000.00 58625.00 0.0 50.0
diff --git a/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_6sides/DATA/CMTSOLUTION b/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_6sides/DATA/CMTSOLUTION
index fca749f..6c3e90c 100644
--- a/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_6sides/DATA/CMTSOLUTION
+++ b/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_6sides/DATA/CMTSOLUTION
@@ -5,9 +5,9 @@ half duration: 5.0
latorUTM: 67000.0
longorUTM: 67000.0
depth: 30.0
-Mrr: 1.000000e+24
-Mtt: 1.000000e+24
-Mpp: 1.000000e+24
-Mrt: 0.
-Mrp: 0.
-Mtp: 0.
+Mrr: 7.700000e+27
+Mtt: 7.700000e+27
+Mpp: 7.700000e+27
+Mrt: 0.0
+Mrp: 0.0
+Mtp: 0.0
diff --git a/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_6sides/DATA/STATIONS b/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_6sides/DATA/STATIONS
index 1ee47ff..0f18aea 100644
--- a/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_6sides/DATA/STATIONS
+++ b/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_6sides/DATA/STATIONS
@@ -1,4 +1,4 @@
-X20 DB 67000.00 22732.14 0.0 0.0
-X30 DB 67000.00 34696.43 0.0 0.0
-X40 DB 67000.00 46660.71 0.0 0.0
-X50 DB 67000.00 58625.00 0.0 0.0
+X20 DB 67000.00 22732.14 0.0 50.0
+X30 DB 67000.00 34696.43 0.0 50.0
+X40 DB 67000.00 46660.71 0.0 50.0
+X50 DB 67000.00 58625.00 0.0 50.0
diff --git a/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_5sides/DATA/STATIONS b/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_5sides/DATA/STATIONS
index 1ee47ff..0f18aea 100644
--- a/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_5sides/DATA/STATIONS
+++ b/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_5sides/DATA/STATIONS
@@ -1,4 +1,4 @@
-X20 DB 67000.00 22732.14 0.0 0.0
-X30 DB 67000.00 34696.43 0.0 0.0
-X40 DB 67000.00 46660.71 0.0 0.0
-X50 DB 67000.00 58625.00 0.0 0.0
+X20 DB 67000.00 22732.14 0.0 50.0
+X30 DB 67000.00 34696.43 0.0 50.0
+X40 DB 67000.00 46660.71 0.0 50.0
+X50 DB 67000.00 58625.00 0.0 50.0
diff --git a/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_6sides/DATA/STATIONS b/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_6sides/DATA/STATIONS
index 1ee47ff..0f18aea 100644
--- a/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_6sides/DATA/STATIONS
+++ b/EXAMPLES/CPML_examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_6sides/DATA/STATIONS
@@ -1,4 +1,4 @@
-X20 DB 67000.00 22732.14 0.0 0.0
-X30 DB 67000.00 34696.43 0.0 0.0
-X40 DB 67000.00 46660.71 0.0 0.0
-X50 DB 67000.00 58625.00 0.0 0.0
+X20 DB 67000.00 22732.14 0.0 50.0
+X30 DB 67000.00 34696.43 0.0 50.0
+X40 DB 67000.00 46660.71 0.0 50.0
+X50 DB 67000.00 58625.00 0.0 50.0
diff --git a/src/specfem3D/compute_forces_acoustic_calling_routine.f90 b/src/specfem3D/compute_forces_acoustic_calling_routine.f90
index a2d1d4e..0bc9b6d 100644
--- a/src/specfem3D/compute_forces_acoustic_calling_routine.f90
+++ b/src/specfem3D/compute_forces_acoustic_calling_routine.f90
@@ -59,7 +59,7 @@ subroutine compute_forces_acoustic()
use specfem_par_elastic
use specfem_par_poroelastic
use pml_par,only: spec_to_CPML,is_CPML,rmemory_coupling_ac_el_displ,nglob_interface_PML_acoustic,&
- b_PML_potential,b_reclen_PML_potential,potential_dot_dot_acoustic_old,potential_acoustic_old
+ b_PML_potential,b_reclen_PML_potential,potential_acoustic_old ! potential_dot_dot_acoustic_old
implicit none
! local parameters
@@ -235,7 +235,7 @@ subroutine compute_forces_acoustic()
potential_dot_acoustic(iglob) = 0.0
potential_acoustic(iglob) = 0.0
if (ELASTIC_SIMULATION) then
- potential_dot_dot_acoustic_old(iglob) = 0.0
+ !potential_dot_dot_acoustic_old(iglob) = 0.0
potential_acoustic_old(iglob) = 0.0
endif
enddo
diff --git a/src/specfem3D/pml_allocate_arrays.f90 b/src/specfem3D/pml_allocate_arrays.f90
index 40a7e01..727bfde 100644
--- a/src/specfem3D/pml_allocate_arrays.f90
+++ b/src/specfem3D/pml_allocate_arrays.f90
@@ -176,8 +176,8 @@ subroutine pml_allocate_arrays()
if (ier /= 0) stop 'error allocating potential_acoustic_new array'
! store the potential acoustic field at n-1 time step
- allocate(potential_dot_dot_acoustic_old(NGLOB_AB),stat=ier)
- if (ier /= 0) stop 'error allocating potential_dot_dot_acoustic_old array'
+! allocate(potential_dot_dot_acoustic_old(NGLOB_AB),stat=ier)
+! if (ier /= 0) stop 'error allocating potential_dot_dot_acoustic_old array'
! stores derivatives of potential with respect to x, y and z
allocate(PML_dpotential_dxl(NGLLX,NGLLY,NGLLZ),stat=ier)
@@ -302,7 +302,7 @@ subroutine pml_allocate_arrays()
if (ACOUSTIC_SIMULATION) then
potential_acoustic_old(:) = 0._CUSTOM_REAL
potential_acoustic_new(:) = 0._CUSTOM_REAL
- potential_dot_dot_acoustic_old(:) = 0._CUSTOM_REAL
+! potential_dot_dot_acoustic_old(:) = 0._CUSTOM_REAL
PML_dpotential_dxl(:,:,:) = 0._CUSTOM_REAL
PML_dpotential_dyl(:,:,:) = 0._CUSTOM_REAL
PML_dpotential_dzl(:,:,:) = 0._CUSTOM_REAL
@@ -529,7 +529,7 @@ subroutine pml_allocate_arrays_dummy()
if (.not. allocated(rmemory_dpotential_dzl)) allocate(rmemory_dpotential_dzl(1,1,1,1,3))
if (.not. allocated(potential_acoustic_old)) allocate(potential_acoustic_old(1))
if (.not. allocated(potential_acoustic_new)) allocate(potential_acoustic_new(1))
- if (.not. allocated(potential_dot_dot_acoustic_old)) allocate(potential_dot_dot_acoustic_old(1))
+! if (.not. allocated(potential_dot_dot_acoustic_old)) allocate(potential_dot_dot_acoustic_old(1))
if (.not. allocated(rmemory_potential_acoustic)) allocate(rmemory_potential_acoustic(1,1,1,1,3))
if (.not. allocated(potential_dot_dot_acoustic_CPML)) allocate(potential_dot_dot_acoustic_CPML(1,1,1))
if (.not. allocated(rmemory_coupling_ac_el_displ)) allocate(rmemory_coupling_ac_el_displ(3,1,1,1,1,2))
diff --git a/src/specfem3D/pml_par.f90 b/src/specfem3D/pml_par.f90
index 69e2444..e6e7b96 100644
--- a/src/specfem3D/pml_par.f90
+++ b/src/specfem3D/pml_par.f90
@@ -126,7 +126,7 @@ module pml_par
real(kind=CUSTOM_REAL), dimension(:,:,:,:,:), allocatable :: rmemory_duz_dxl_z,rmemory_duz_dyl_z,rmemory_duz_dzl_z
!store the field accel at time step n-1 for second order convolution scheme
- real(kind=CUSTOM_REAL), dimension(:), allocatable :: potential_dot_dot_acoustic_old
+! real(kind=CUSTOM_REAL), dimension(:), allocatable :: potential_dot_dot_acoustic_old
real(kind=CUSTOM_REAL), dimension(:,:,:,:,:), allocatable :: rmemory_dpotential_dxl
real(kind=CUSTOM_REAL), dimension(:,:,:,:,:), allocatable :: rmemory_dpotential_dyl
More information about the CIG-COMMITS
mailing list