[cig-commits] [commit] devel,master: fixed hdur and f0_for_PML in all the CPML examples; switched to a source that creates S waves as well for the elastic example (f2f9f14)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 16:51:31 PDT 2014
Repository : https://github.com/geodynamics/specfem3d
On branches: devel,master
Link : https://github.com/geodynamics/specfem3d/compare/6026e367984905ab133865f62fa6293b343759b9...47f703851338234f96397e7da9fbff63d8178b8a
>---------------------------------------------------------------
commit f2f9f14a6cca757565eed40b23afd568ce336b34
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date: Tue Mar 26 14:50:24 2013 +0000
fixed hdur and f0_for_PML in all the CPML examples; switched to a source that creates S waves as well for the elastic example
>---------------------------------------------------------------
f2f9f14a6cca757565eed40b23afd568ce336b34
.../HOMO8_FREE_SURFACE/DATA/CMTSOLUTION | 2 +-
.../DATA/MESH/print_min_and_max_mesh_size.f90 | 28 +++++++++++++++++++++
.../HOMO8_FREE_SURFACE/DATA/Par_file | 2 +-
.../HOMO8_NOFREESURFACE/DATA/CMTSOLUTION | 2 +-
.../DATA/MESH/print_min_and_max_mesh_size.f90 | 28 +++++++++++++++++++++
.../HOMO8_NOFREESURFACE/DATA/Par_file | 2 +-
.../geometrical_properties_of_this_PML_mesh.txt | 23 +++++++++++++++++
CPML/purely_acoustic/image_of_this_CPML_mesh.jpg | Bin 0 -> 699041 bytes
.../HOMO8_FREE_SURFACE/DATA/CMTSOLUTION | 18 ++++++-------
.../DATA/MESH/print_min_and_max_mesh_size.f90 | 28 +++++++++++++++++++++
.../HOMO8_FREE_SURFACE/DATA/Par_file | 2 +-
.../HOMO8_NOFREESURFACE/DATA/CMTSOLUTION | 18 ++++++-------
.../DATA/MESH/print_min_and_max_mesh_size.f90 | 28 +++++++++++++++++++++
.../HOMO8_NOFREESURFACE/DATA/Par_file | 2 +-
CPML/purely_elastic/image_of_this_CPML_mesh.jpg | Bin 0 -> 699041 bytes
15 files changed, 159 insertions(+), 24 deletions(-)
diff --git a/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/CMTSOLUTION b/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/CMTSOLUTION
index 563743f..4731cff 100644
--- a/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/CMTSOLUTION
+++ b/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/CMTSOLUTION
@@ -1,7 +1,7 @@
PDE 1999 01 01 00 00 00.00 67000 67000 -25000 4.2 4.2 hom_explosion
event name: hom_explosion
time shift: 0.0000
-half duration: 5.0
+half duration: 15.0
latitude: 67000.0
longitude: 67000.0
depth: 25.0
diff --git a/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/MESH/print_min_and_max_mesh_size.f90 b/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/MESH/print_min_and_max_mesh_size.f90
new file mode 100644
index 0000000..5cd9b35
--- /dev/null
+++ b/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/MESH/print_min_and_max_mesh_size.f90
@@ -0,0 +1,28 @@
+
+ program print_min_and_max_mesh_size
+
+ implicit none
+
+ real, dimension(:), allocatable :: x,y,z
+
+ integer :: ipoin,npoin,idummy
+
+ open(unit=23,file="nodes_coords_file",status="old")
+ read(23,*) npoin
+
+ allocate(x(npoin))
+ allocate(y(npoin))
+ allocate(z(npoin))
+
+ do ipoin = 1,npoin
+ read(23,*) idummy,x(ipoin),y(ipoin),z(ipoin)
+ enddo
+
+ close(23)
+
+ print *, 'X min max = ',minval(x),maxval(x),' deltaX = ',maxval(x) - minval(x)
+ print *, 'Y min max = ',minval(y),maxval(y),' deltaY = ',maxval(y) - minval(y)
+ print *, 'Z min max = ',minval(z),maxval(z),' deltaZ = ',maxval(z) - minval(z)
+
+ end program print_min_and_max_mesh_size
+
diff --git a/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/Par_file b/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/Par_file
index 407b552..1ec6133 100644
--- a/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/Par_file
+++ b/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/Par_file
@@ -65,7 +65,7 @@ PML_WIDTH_MIN = 10000.0d0
PML_WIDTH_MAX = 25000.0d0
# C-PML dominant frequency
-f0_FOR_PML = 12.7
+f0_FOR_PML = 0.09160
# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
# ROTATE_PML_ACTIVATE = .false.
diff --git a/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/CMTSOLUTION b/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/CMTSOLUTION
index 563743f..4731cff 100644
--- a/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/CMTSOLUTION
+++ b/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/CMTSOLUTION
@@ -1,7 +1,7 @@
PDE 1999 01 01 00 00 00.00 67000 67000 -25000 4.2 4.2 hom_explosion
event name: hom_explosion
time shift: 0.0000
-half duration: 5.0
+half duration: 15.0
latitude: 67000.0
longitude: 67000.0
depth: 25.0
diff --git a/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/MESH/print_min_and_max_mesh_size.f90 b/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/MESH/print_min_and_max_mesh_size.f90
new file mode 100644
index 0000000..5cd9b35
--- /dev/null
+++ b/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/MESH/print_min_and_max_mesh_size.f90
@@ -0,0 +1,28 @@
+
+ program print_min_and_max_mesh_size
+
+ implicit none
+
+ real, dimension(:), allocatable :: x,y,z
+
+ integer :: ipoin,npoin,idummy
+
+ open(unit=23,file="nodes_coords_file",status="old")
+ read(23,*) npoin
+
+ allocate(x(npoin))
+ allocate(y(npoin))
+ allocate(z(npoin))
+
+ do ipoin = 1,npoin
+ read(23,*) idummy,x(ipoin),y(ipoin),z(ipoin)
+ enddo
+
+ close(23)
+
+ print *, 'X min max = ',minval(x),maxval(x),' deltaX = ',maxval(x) - minval(x)
+ print *, 'Y min max = ',minval(y),maxval(y),' deltaY = ',maxval(y) - minval(y)
+ print *, 'Z min max = ',minval(z),maxval(z),' deltaZ = ',maxval(z) - minval(z)
+
+ end program print_min_and_max_mesh_size
+
diff --git a/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/Par_file b/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/Par_file
index 0f7eeed..7b8711e 100644
--- a/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/Par_file
+++ b/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/Par_file
@@ -65,7 +65,7 @@ PML_WIDTH_MIN = 10000.0d0
PML_WIDTH_MAX = 25000.0d0
# C-PML dominant frequency
-f0_FOR_PML = 12.7
+f0_FOR_PML = 0.09160
# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
# ROTATE_PML_ACTIVATE = .false.
diff --git a/CPML/purely_acoustic/geometrical_properties_of_this_PML_mesh.txt b/CPML/purely_acoustic/geometrical_properties_of_this_PML_mesh.txt
new file mode 100644
index 0000000..460580a
--- /dev/null
+++ b/CPML/purely_acoustic/geometrical_properties_of_this_PML_mesh.txt
@@ -0,0 +1,23 @@
+
+- regular mesh
+
+- flat model edges aligned with the coordinate axes X, Y and Z
+
+- 12 spectral elements along X, 12 spectral elements along X, 13 spectral elements along Z
+
+- thus total number of spectral elements = 12 * 12 * 13 = 1872
+
+- size of the model:
+
+ X min max = 273500.0 m 352100.0 m deltaX = 78600.00 m (size of the whole model)
+ Y min max = 3614900. m 3691600. m deltaY = 76700.00 m (size of the whole model)
+ Z min max = -79785.46 m -685.4575 m deltaZ = 79100.00 m (size of the whole model)
+
+- size of a spectral element:
+
+ X size = 78600.00 / 12 = 6550 m
+ Y size = 76700.00 / 12 = 6391.67 m
+ Z size = 79100.00 / 13 = 6084.615 m
+
+- acoustic material of constant properties in the whole mesh: rho = 2300 kg/m3, Vp = 2800 m/s
+
diff --git a/CPML/purely_acoustic/image_of_this_CPML_mesh.jpg b/CPML/purely_acoustic/image_of_this_CPML_mesh.jpg
new file mode 100644
index 0000000..b64d6ec
Binary files /dev/null and b/CPML/purely_acoustic/image_of_this_CPML_mesh.jpg differ
diff --git a/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/CMTSOLUTION b/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/CMTSOLUTION
index 563743f..55fd3a7 100644
--- a/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/CMTSOLUTION
+++ b/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/CMTSOLUTION
@@ -1,13 +1,13 @@
-PDE 1999 01 01 00 00 00.00 67000 67000 -25000 4.2 4.2 hom_explosion
-event name: hom_explosion
+PDE 1999 01 01 00 00 00.00 67000 67000 -25000 4.2 4.2 BOLIVIA
+event name: BOLIVIA
time shift: 0.0000
-half duration: 5.0
+half duration: 15.0
latitude: 67000.0
longitude: 67000.0
depth: 25.0
-Mrr: 1.000000e+23
-Mtt: 1.000000e+23
-Mpp: 1.000000e+23
-Mrt: 0.000000
-Mrp: 0.000000
-Mtp: 0.000000
+Mrr: -0.759000E+28
+Mtt: 0.775000E+28
+Mpp: -0.160000E+27
+Mrt: -0.250300E+29
+Mrp: 0.420000E+27
+Mtp: -0.248000E+28
diff --git a/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/MESH/print_min_and_max_mesh_size.f90 b/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/MESH/print_min_and_max_mesh_size.f90
new file mode 100644
index 0000000..5cd9b35
--- /dev/null
+++ b/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/MESH/print_min_and_max_mesh_size.f90
@@ -0,0 +1,28 @@
+
+ program print_min_and_max_mesh_size
+
+ implicit none
+
+ real, dimension(:), allocatable :: x,y,z
+
+ integer :: ipoin,npoin,idummy
+
+ open(unit=23,file="nodes_coords_file",status="old")
+ read(23,*) npoin
+
+ allocate(x(npoin))
+ allocate(y(npoin))
+ allocate(z(npoin))
+
+ do ipoin = 1,npoin
+ read(23,*) idummy,x(ipoin),y(ipoin),z(ipoin)
+ enddo
+
+ close(23)
+
+ print *, 'X min max = ',minval(x),maxval(x),' deltaX = ',maxval(x) - minval(x)
+ print *, 'Y min max = ',minval(y),maxval(y),' deltaY = ',maxval(y) - minval(y)
+ print *, 'Z min max = ',minval(z),maxval(z),' deltaZ = ',maxval(z) - minval(z)
+
+ end program print_min_and_max_mesh_size
+
diff --git a/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/Par_file b/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/Par_file
index 407b552..1ec6133 100644
--- a/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/Par_file
+++ b/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/Par_file
@@ -65,7 +65,7 @@ PML_WIDTH_MIN = 10000.0d0
PML_WIDTH_MAX = 25000.0d0
# C-PML dominant frequency
-f0_FOR_PML = 12.7
+f0_FOR_PML = 0.09160
# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
# ROTATE_PML_ACTIVATE = .false.
diff --git a/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/CMTSOLUTION b/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/CMTSOLUTION
index 563743f..7d91cc6 100644
--- a/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/CMTSOLUTION
+++ b/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/CMTSOLUTION
@@ -1,13 +1,13 @@
-PDE 1999 01 01 00 00 00.00 67000 67000 -25000 4.2 4.2 hom_explosion
-event name: hom_explosion
+PDE 1999 01 01 00 00 00.00 67000 67000 -25000 4.2 4.2 BOLIVIA
+event name: BOLIVIA
time shift: 0.0000
-half duration: 5.0
+half duration: 15.0
latitude: 67000.0
longitude: 67000.0
depth: 25.0
-Mrr: 1.000000e+23
-Mtt: 1.000000e+23
-Mpp: 1.000000e+23
-Mrt: 0.000000
-Mrp: 0.000000
-Mtp: 0.000000
+Mrr: -0.759000E+28
+Mtt: 0.775000E+28
+Mpp: -0.160000E+27
+Mrt: -0.250300E+29
+Mrp: 0.420000E+27
+Mtp: -0.248000E+28
diff --git a/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/MESH/print_min_and_max_mesh_size.f90 b/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/MESH/print_min_and_max_mesh_size.f90
new file mode 100644
index 0000000..5cd9b35
--- /dev/null
+++ b/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/MESH/print_min_and_max_mesh_size.f90
@@ -0,0 +1,28 @@
+
+ program print_min_and_max_mesh_size
+
+ implicit none
+
+ real, dimension(:), allocatable :: x,y,z
+
+ integer :: ipoin,npoin,idummy
+
+ open(unit=23,file="nodes_coords_file",status="old")
+ read(23,*) npoin
+
+ allocate(x(npoin))
+ allocate(y(npoin))
+ allocate(z(npoin))
+
+ do ipoin = 1,npoin
+ read(23,*) idummy,x(ipoin),y(ipoin),z(ipoin)
+ enddo
+
+ close(23)
+
+ print *, 'X min max = ',minval(x),maxval(x),' deltaX = ',maxval(x) - minval(x)
+ print *, 'Y min max = ',minval(y),maxval(y),' deltaY = ',maxval(y) - minval(y)
+ print *, 'Z min max = ',minval(z),maxval(z),' deltaZ = ',maxval(z) - minval(z)
+
+ end program print_min_and_max_mesh_size
+
diff --git a/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/Par_file b/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/Par_file
index 0f7eeed..7b8711e 100644
--- a/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/Par_file
+++ b/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/Par_file
@@ -65,7 +65,7 @@ PML_WIDTH_MIN = 10000.0d0
PML_WIDTH_MAX = 25000.0d0
# C-PML dominant frequency
-f0_FOR_PML = 12.7
+f0_FOR_PML = 0.09160
# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
# ROTATE_PML_ACTIVATE = .false.
diff --git a/CPML/purely_elastic/image_of_this_CPML_mesh.jpg b/CPML/purely_elastic/image_of_this_CPML_mesh.jpg
new file mode 100644
index 0000000..b64d6ec
Binary files /dev/null and b/CPML/purely_elastic/image_of_this_CPML_mesh.jpg differ
More information about the CIG-COMMITS
mailing list