[cig-commits] r21639 - in seismo/3D/SPECFEM3D/trunk/examples/CPML: purely_acoustic purely_acoustic/HOMO8_FREE_SURFACE/DATA purely_acoustic/HOMO8_FREE_SURFACE/DATA/MESH purely_acoustic/HOMO8_NOFREESURFACE/DATA purely_acoustic/HOMO8_NOFREESURFACE/DATA/MESH purely_elastic purely_elastic/HOMO8_FREE_SURFACE/DATA purely_elastic/HOMO8_FREE_SURFACE/DATA/MESH purely_elastic/HOMO8_NOFREESURFACE/DATA purely_elastic/HOMO8_NOFREESURFACE/DATA/MESH
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Tue Mar 26 07:50:25 PDT 2013
Author: dkomati1
Date: 2013-03-26 07:50:24 -0700 (Tue, 26 Mar 2013)
New Revision: 21639
Added:
seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/MESH/print_min_and_max_mesh_size.f90
seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/MESH/print_min_and_max_mesh_size.f90
seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/geometrical_properties_of_this_PML_mesh.txt
seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/image_of_this_CPML_mesh.jpg
seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/MESH/print_min_and_max_mesh_size.f90
seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/MESH/print_min_and_max_mesh_size.f90
seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/image_of_this_CPML_mesh.jpg
Modified:
seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/CMTSOLUTION
seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/CMTSOLUTION
seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/CMTSOLUTION
seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/CMTSOLUTION
seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/Par_file
Log:
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
Modified: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/CMTSOLUTION
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/CMTSOLUTION 2013-03-26 13:35:01 UTC (rev 21638)
+++ seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/CMTSOLUTION 2013-03-26 14:50:24 UTC (rev 21639)
@@ -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
Added: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/MESH/print_min_and_max_mesh_size.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/MESH/print_min_and_max_mesh_size.f90 (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/MESH/print_min_and_max_mesh_size.f90 2013-03-26 14:50:24 UTC (rev 21639)
@@ -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
+
Modified: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/Par_file 2013-03-26 13:35:01 UTC (rev 21638)
+++ seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_FREE_SURFACE/DATA/Par_file 2013-03-26 14:50:24 UTC (rev 21639)
@@ -65,7 +65,7 @@
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.
Modified: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/CMTSOLUTION
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/CMTSOLUTION 2013-03-26 13:35:01 UTC (rev 21638)
+++ seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/CMTSOLUTION 2013-03-26 14:50:24 UTC (rev 21639)
@@ -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
Added: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/MESH/print_min_and_max_mesh_size.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/MESH/print_min_and_max_mesh_size.f90 (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/MESH/print_min_and_max_mesh_size.f90 2013-03-26 14:50:24 UTC (rev 21639)
@@ -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
+
Modified: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/Par_file 2013-03-26 13:35:01 UTC (rev 21638)
+++ seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/HOMO8_NOFREESURFACE/DATA/Par_file 2013-03-26 14:50:24 UTC (rev 21639)
@@ -65,7 +65,7 @@
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.
Added: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/geometrical_properties_of_this_PML_mesh.txt
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/geometrical_properties_of_this_PML_mesh.txt (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/geometrical_properties_of_this_PML_mesh.txt 2013-03-26 14:50:24 UTC (rev 21639)
@@ -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
+
Added: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/image_of_this_CPML_mesh.jpg
===================================================================
(Binary files differ)
Property changes on: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_acoustic/image_of_this_CPML_mesh.jpg
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Modified: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/CMTSOLUTION
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/CMTSOLUTION 2013-03-26 13:35:01 UTC (rev 21638)
+++ seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/CMTSOLUTION 2013-03-26 14:50:24 UTC (rev 21639)
@@ -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
Added: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/MESH/print_min_and_max_mesh_size.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/MESH/print_min_and_max_mesh_size.f90 (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/MESH/print_min_and_max_mesh_size.f90 2013-03-26 14:50:24 UTC (rev 21639)
@@ -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
+
Modified: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/Par_file 2013-03-26 13:35:01 UTC (rev 21638)
+++ seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_FREE_SURFACE/DATA/Par_file 2013-03-26 14:50:24 UTC (rev 21639)
@@ -65,7 +65,7 @@
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.
Modified: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/CMTSOLUTION
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/CMTSOLUTION 2013-03-26 13:35:01 UTC (rev 21638)
+++ seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/CMTSOLUTION 2013-03-26 14:50:24 UTC (rev 21639)
@@ -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
Added: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/MESH/print_min_and_max_mesh_size.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/MESH/print_min_and_max_mesh_size.f90 (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/MESH/print_min_and_max_mesh_size.f90 2013-03-26 14:50:24 UTC (rev 21639)
@@ -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
+
Modified: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/Par_file 2013-03-26 13:35:01 UTC (rev 21638)
+++ seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/HOMO8_NOFREESURFACE/DATA/Par_file 2013-03-26 14:50:24 UTC (rev 21639)
@@ -65,7 +65,7 @@
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.
Added: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/image_of_this_CPML_mesh.jpg
===================================================================
(Binary files differ)
Property changes on: seismo/3D/SPECFEM3D/trunk/examples/CPML/purely_elastic/image_of_this_CPML_mesh.jpg
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
More information about the CIG-COMMITS
mailing list