[cig-commits] r18217 - in seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D: . example_deep example_utm
carltape at geodynamics.org
carltape at geodynamics.org
Sun Apr 10 11:22:53 PDT 2011
Author: carltape
Date: 2011-04-10 11:22:52 -0700 (Sun, 10 Apr 2011)
New Revision: 18217
Added:
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/Mesh_Par_file_deep
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/Mesh_Par_file_deep_honor60km
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/README
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/honor60km_or_not.png
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/interface0.dat
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/interfaces_deep.dat
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/interfaces_deep_honor60km.dat
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_utm/
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_utm/Mesh_Par_file_utm
Removed:
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_utm/Mesh_Par_file_001_utm
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/
Modified:
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_utm/README
Log:
added socal1D sub-example to illustrate how to extend the identical model downward in depth; because layers are ordered from bottom up, this requires several minor changes to ensure that the meshes are identical in the overlapping region.
Added: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/Mesh_Par_file_deep
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/Mesh_Par_file_deep (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/Mesh_Par_file_deep 2011-04-10 18:22:52 UTC (rev 18217)
@@ -0,0 +1,57 @@
+
+# coordinates of mesh block in latitude/longitude and depth in km
+LATITUDE_MIN = 33.790d0
+LATITUDE_MAX = 34.157d0
+LONGITUDE_MIN = -118.728d0
+LONGITUDE_MAX = -118.237d0
+DEPTH_BLOCK_KM = 200.d0
+UTM_PROJECTION_ZONE = 11
+SUPPRESS_UTM_PROJECTION = .false.
+
+# file that contains the interfaces of the model / mesh
+INTERFACES_FILE = interfaces.dat
+
+# number of elements at the surface along edges of the mesh at the surface
+# (must be 8 * multiple of NPROC below if mesh is not regular and contains mesh doublings)
+# (must be multiple of NPROC below if mesh is regular)
+NEX_XI = 24
+NEX_ETA = 24
+
+# number of MPI processors along xi and eta (can be different)
+NPROC_XI = 1
+NPROC_ETA = 1
+
+# Regular/irregular mesh
+USE_REGULAR_MESH = .false.
+# Only for irregular meshes, number of doubling layers (1 or 2) and their position
+NDOUBLINGS = 2
+# NZ_DOUGLING_1 is the parameter to set up if there is only one doubling layer
+NZ_DOUGLING_1 = 20
+NZ_DOUGLING_2 = 23
+
+# create mesh files for visualisation or further checking
+CREATE_ABAQUS_FILES = .true.
+CREATE_DX_FILES = .true.
+
+# path to store the databases files
+LOCAL_PATH = ../in_out_files/DATABASES_MPI
+
+# number of materials
+NMATERIALS = 4
+# define the different materials in the model as :
+# #material_id #rho #vp #vs #Q #anisotropy_flag #domain_id
+# Q : quality factor
+# anisotropy_flag : 0=no anisotropy/ 1,2,.. check with implementation in aniso_model.f90
+# domain_id : 1=acoustic / 2=elastic
+1 3000 7800 4500 0.0 0 2
+2 2800 6700 3870 0.0 0 2
+3 2670 6300 3640 0.0 0 2
+4 2400 5500 3180 0.0 0 2
+# number of regions
+NREGIONS = 4
+# define the different regions of the model as :
+#NEX_XI_BEGIN #NEX_XI_END #NEX_ETA_BEGIN #NEX_ETA_END #NZ_BEGIN #NZ_END #material_id
+1 24 1 24 1 17 1
+1 24 1 24 18 20 2
+1 24 1 24 21 23 3
+1 24 1 24 24 26 4
\ No newline at end of file
Added: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/Mesh_Par_file_deep_honor60km
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/Mesh_Par_file_deep_honor60km (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/Mesh_Par_file_deep_honor60km 2011-04-10 18:22:52 UTC (rev 18217)
@@ -0,0 +1,58 @@
+
+# coordinates of mesh block in latitude/longitude and depth in km
+LATITUDE_MIN = 33.790d0
+LATITUDE_MAX = 34.157d0
+LONGITUDE_MIN = -118.728d0
+LONGITUDE_MAX = -118.237d0
+DEPTH_BLOCK_KM = 200.d0
+UTM_PROJECTION_ZONE = 11
+SUPPRESS_UTM_PROJECTION = .false.
+
+# file that contains the interfaces of the model / mesh
+INTERFACES_FILE = interfaces.dat
+
+# number of elements at the surface along edges of the mesh at the surface
+# (must be 8 * multiple of NPROC below if mesh is not regular and contains mesh doublings)
+# (must be multiple of NPROC below if mesh is regular)
+NEX_XI = 24
+NEX_ETA = 24
+
+# number of MPI processors along xi and eta (can be different)
+NPROC_XI = 1
+NPROC_ETA = 1
+
+# Regular/irregular mesh
+USE_REGULAR_MESH = .false.
+# Only for irregular meshes, number of doubling layers (1 or 2) and their position
+NDOUBLINGS = 2
+# NZ_DOUGLING_1 is the parameter to set up if there is only one doubling layer
+NZ_DOUGLING_1 = 20
+NZ_DOUGLING_2 = 23
+
+# create mesh files for visualisation or further checking
+CREATE_ABAQUS_FILES = .true.
+CREATE_DX_FILES = .true.
+
+# path to store the databases files
+LOCAL_PATH = ../in_out_files/DATABASES_MPI
+
+# number of materials
+NMATERIALS = 4
+# define the different materials in the model as :
+# #material_id #rho #vp #vs #Q #anisotropy_flag #domain_id
+# Q : quality factor
+# anisotropy_flag : 0=no anisotropy/ 1,2,.. check with implementation in aniso_model.f90
+# domain_id : 1=acoustic / 2=elastic
+1 3000 7800 4500 0.0 0 2
+2 2800 6700 3870 0.0 0 2
+3 2670 6300 3640 0.0 0 2
+4 2400 5500 3180 0.0 0 2
+# number of regions
+NREGIONS = 5
+# define the different regions of the model as :
+#NEX_XI_BEGIN #NEX_XI_END #NEX_ETA_BEGIN #NEX_ETA_END #NZ_BEGIN #NZ_END #material_id
+1 24 1 24 1 14 1
+1 24 1 24 15 17 1
+1 24 1 24 18 20 2
+1 24 1 24 21 23 3
+1 24 1 24 24 26 4
\ No newline at end of file
Added: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/README
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/README (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/README 2011-04-10 18:22:52 UTC (rev 18217)
@@ -0,0 +1,29 @@
+README
+(SPECFEM3D/examples/meshfem3D_examples/socal1D/example_deep/README)
+
+This example illustrates two options for extending the bottom layer of the socal1D model. In some scenarios, it is helpful to extend the boundaries to eliminate spurious reflections from boundaries.
+
+First run the default socal1D example with 1 core.
+Then try examples A and B.
+
+EXAMPLE A: Change bottom interface from -60 km to -200 km; use uniform element
+ spacing in bottom layer.
+
+ > cp Mesh_Par_file_deep ../Mesh_Par_file
+ > cd interfaces_deep.dat ../interfaces.dat
+ > cd ..
+ > ./process.sh
+
+EXAMPLE B: Change bottom interface from -60 km to -200 km; add -60 km as a
+ new interface, such that the mesh above -60 km remains identical to
+ the default model that extends to 60 km depth only.
+
+ > cp Mesh_Par_file_deep_honor60km ../Mesh_Par_file
+ > cp interfaces_deep_honor60km.dat ../interfaces.dat
+ > cp interface0.dat ..
+ > cd ..
+ > ./process.sh
+
+The models for A and B are identical, but the meshes are subtley different: as shown in honor60km_or_not.png, mesh A does not honor the 60 km interface, whereas mesh B does. Since there is no material contrast at 60 km depth, the synthetics are expected to be nearly identical.
+
+=======================
\ No newline at end of file
Added: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/honor60km_or_not.png
===================================================================
(Binary files differ)
Property changes on: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/honor60km_or_not.png
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/interface0.dat
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/interface0.dat (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/interface0.dat 2011-04-10 18:22:52 UTC (rev 18217)
@@ -0,0 +1,4 @@
+-60000
+-60000
+-60000
+-60000
Added: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/interfaces_deep.dat
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/interfaces_deep.dat (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/interfaces_deep.dat 2011-04-10 18:22:52 UTC (rev 18217)
@@ -0,0 +1,30 @@
+# number of interfaces
+ 4
+#
+# We describe each interface below, structured as a 2D-grid, with several parameters :
+# number of points along XI and ETA, minimal XI ETA coordinates
+# and spacing between points which must be constant.
+# Then the records contain the Z coordinates of the NXI x NETA points.
+#
+# interface number 1 (moho)
+ .false. 2 2 -121.d0 32.d0 6.0d0 6.0d0
+ interface1.dat
+# interface number 2
+ .false. 2 2 -121.d0 32.d0 6.0d0 6.0d0
+ interface2.dat
+# interface number 3
+ .false. 2 2 -121.d0 32.d0 6.0d0 6.0d0
+ interface3.dat
+# interface number 4 (topography, top of the mesh)
+ .false. 2 2 -121.d0 32.d0 6.0d0 6.0d0
+ interface4.dat
+#
+# for each layer, we give the number of spectral elements in the vertical direction
+# layer number 1 (bottom layer)
+ 17
+# layer number 2
+ 3
+# layer number 3
+ 3
+# layer number 4 (top layer)
+ 3
\ No newline at end of file
Added: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/interfaces_deep_honor60km.dat
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/interfaces_deep_honor60km.dat (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_deep/interfaces_deep_honor60km.dat 2011-04-10 18:22:52 UTC (rev 18217)
@@ -0,0 +1,35 @@
+# number of interfaces
+ 5
+#
+# We describe each interface below, structured as a 2D-grid, with several parameters :
+# number of points along XI and ETA, minimal XI ETA coordinates
+# and spacing between points which must be constant.
+# Then the records contain the Z coordinates of the NXI x NETA points.
+#
+# interface number 0
+ .false. 2 2 -121.d0 32.d0 6.0d0 6.0d0
+ interface0.dat
+# interface number 1 (moho)
+ .false. 2 2 -121.d0 32.d0 6.0d0 6.0d0
+ interface1.dat
+# interface number 2
+ .false. 2 2 -121.d0 32.d0 6.0d0 6.0d0
+ interface2.dat
+# interface number 3
+ .false. 2 2 -121.d0 32.d0 6.0d0 6.0d0
+ interface3.dat
+# interface number 4 (topography, top of the mesh)
+ .false. 2 2 -121.d0 32.d0 6.0d0 6.0d0
+ interface4.dat
+#
+# for each layer, we give the number of spectral elements in the vertical direction
+# layer number 0 (bottom layer)
+ 14
+# layer number 1
+ 3
+# layer number 2
+ 3
+# layer number 4
+ 3
+# layer number 5 (top layer)
+ 3
\ No newline at end of file
Copied: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_utm (from rev 18216, seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example)
Deleted: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_utm/Mesh_Par_file_001_utm
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/Mesh_Par_file_001_utm 2011-04-10 05:14:33 UTC (rev 18216)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_utm/Mesh_Par_file_001_utm 2011-04-10 18:22:52 UTC (rev 18217)
@@ -1,58 +0,0 @@
-
-# coordinates of mesh block in latitude/longitude and depth in km
-LATITUDE_MIN = 3739440.d0
-LATITUDE_MAX = 3781380.d0
-LONGITUDE_MIN = 340080.d0
-LONGITUDE_MAX = 385720.d0
-DEPTH_BLOCK_KM = 60.d0
-UTM_PROJECTION_ZONE = 11
-SUPPRESS_UTM_PROJECTION = .true.
-
-# file that contains the interfaces of the model / mesh
-INTERFACES_FILE = interfaces.dat
-
-# number of elements at the surface along edges of the mesh at the surface
-# (must be 8 * multiple of NPROC below if mesh is not regular and contains mesh doublings)
-# (must be multiple of NPROC below if mesh is regular)
-NEX_XI = 24
-NEX_ETA = 24
-
-# number of MPI processors along xi and eta (can be different)
-NPROC_XI = 1
-NPROC_ETA = 1
-
-# Regular/irregular mesh
-USE_REGULAR_MESH = .false.
-# Only for irregular meshes, number of doubling layers (1 or 2) and their position
-NDOUBLINGS = 2
-# NZ_DOUGLING_1 is the parameter to set up if there is only one doubling layer
-NZ_DOUGLING_1 = 6
-NZ_DOUGLING_2 = 9
-
-# create mesh files for visualisation or further checking
-CREATE_ABAQUS_FILES = .true.
-CREATE_DX_FILES = .true.
-
-# path to store the databases files
-LOCAL_PATH = ../in_out_files/DATABASES_MPI
-
-# number of materials
-NMATERIALS = 4
-# define the different materials in the model as :
-# #material_id #rho #vp #vs #Q #anisotropy_flag #domain_id
-# Q : quality factor
-# anisotropy_flag : 0=no anisotropy/ 1,2,.. check with implementation in aniso_model.f90
-# domain_id : 1=acoustic / 2=elastic
-1 3000 7800 4500 0.0 0 2
-2 2800 6700 3870 0.0 0 2
-3 2670 6300 3640 0.0 0 2
-4 2400 5500 3180 0.0 0 2
-# number of regions
-NREGIONS = 4
-# define the different regions of the model as :
-#NEX_XI_BEGIN #NEX_XI_END #NEX_ETA_BEGIN #NEX_ETA_END #NZ_BEGIN #NZ_END #material_id
-1 24 1 24 1 3 1
-1 24 1 24 4 6 2
-1 24 1 24 7 9 3
-1 24 1 24 10 12 4
-
Copied: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_utm/Mesh_Par_file_utm (from rev 18216, seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/Mesh_Par_file_001_utm)
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_utm/Mesh_Par_file_utm (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_utm/Mesh_Par_file_utm 2011-04-10 18:22:52 UTC (rev 18217)
@@ -0,0 +1,58 @@
+
+# coordinates of mesh block in latitude/longitude and depth in km
+LATITUDE_MIN = 3739440.d0
+LATITUDE_MAX = 3781380.d0
+LONGITUDE_MIN = 340080.d0
+LONGITUDE_MAX = 385720.d0
+DEPTH_BLOCK_KM = 60.d0
+UTM_PROJECTION_ZONE = 11
+SUPPRESS_UTM_PROJECTION = .true.
+
+# file that contains the interfaces of the model / mesh
+INTERFACES_FILE = interfaces.dat
+
+# number of elements at the surface along edges of the mesh at the surface
+# (must be 8 * multiple of NPROC below if mesh is not regular and contains mesh doublings)
+# (must be multiple of NPROC below if mesh is regular)
+NEX_XI = 24
+NEX_ETA = 24
+
+# number of MPI processors along xi and eta (can be different)
+NPROC_XI = 1
+NPROC_ETA = 1
+
+# Regular/irregular mesh
+USE_REGULAR_MESH = .false.
+# Only for irregular meshes, number of doubling layers (1 or 2) and their position
+NDOUBLINGS = 2
+# NZ_DOUGLING_1 is the parameter to set up if there is only one doubling layer
+NZ_DOUGLING_1 = 6
+NZ_DOUGLING_2 = 9
+
+# create mesh files for visualisation or further checking
+CREATE_ABAQUS_FILES = .true.
+CREATE_DX_FILES = .true.
+
+# path to store the databases files
+LOCAL_PATH = ../in_out_files/DATABASES_MPI
+
+# number of materials
+NMATERIALS = 4
+# define the different materials in the model as :
+# #material_id #rho #vp #vs #Q #anisotropy_flag #domain_id
+# Q : quality factor
+# anisotropy_flag : 0=no anisotropy/ 1,2,.. check with implementation in aniso_model.f90
+# domain_id : 1=acoustic / 2=elastic
+1 3000 7800 4500 0.0 0 2
+2 2800 6700 3870 0.0 0 2
+3 2670 6300 3640 0.0 0 2
+4 2400 5500 3180 0.0 0 2
+# number of regions
+NREGIONS = 4
+# define the different regions of the model as :
+#NEX_XI_BEGIN #NEX_XI_END #NEX_ETA_BEGIN #NEX_ETA_END #NZ_BEGIN #NZ_END #material_id
+1 24 1 24 1 3 1
+1 24 1 24 4 6 2
+1 24 1 24 7 9 3
+1 24 1 24 10 12 4
+
Modified: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_utm/README
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/README 2011-04-10 05:14:33 UTC (rev 18216)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/example_utm/README 2011-04-10 18:22:52 UTC (rev 18217)
@@ -1,5 +1,5 @@
README
-(SPECFEM3D/examples/meshfem3D_examples/socal1D/utm/README)
+(SPECFEM3D/examples/meshfem3D_examples/socal1D/example_utm/README)
This example is mainly to show the ease of using either UTM or lon-lat coordinates.
The utm example is better suited for direct comparison with external meshes for
@@ -10,7 +10,7 @@
EXAMPLE A: Input geographic bounds in UTM, everything else in lon-lat.
- > cp Mesh_Par_file_001_utm ../Mesh_Par_file
+ > cp Mesh_Par_file_utm ../Mesh_Par_file
> cd ..
> ./process.sh
@@ -19,7 +19,7 @@
> cp CMTSOLUTION_utm ../CMTSOLUTION
> cp STATIONS_utm ../STATIONS
> cp interfaces_utm.dat ../interfaces.dat
- > cp Mesh_Par_file_001_utm ../Mesh_Par_file
+ > cp Mesh_Par_file_utm ../Mesh_Par_file
> cp Par_file_utm ../Par_file
> cd ..
> ./process.sh
More information about the CIG-COMMITS
mailing list