[cig-commits] [commit] devel: Updates SEP example coarser mesh and stations (bb11664)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Thu Oct 2 16:01:19 PDT 2014
Repository : https://github.com/geodynamics/specfem3d
On branch : devel
Link : https://github.com/geodynamics/specfem3d/compare/9473eea83e6cad4b6bd8adbaee3e9c074ce3b975...95b31791a026eedc69a0970b245c691394cfbc71
>---------------------------------------------------------------
commit bb11664d5a102596083a44bcad98541b5d89d05c
Author: Matthieu Lefebvre <ml15 at princeton.edu>
Date: Thu Oct 2 15:49:05 2014 -0400
Updates SEP example coarser mesh and stations
>---------------------------------------------------------------
bb11664d5a102596083a44bcad98541b5d89d05c
EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/CMTSOLUTION | 8 ++++----
.../meshfem3D_examples/sep_bathymetry/DATA/FORCESOLUTION | 12 ++++++------
EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/Par_file | 4 ++--
.../meshfem3D_examples/sep_bathymetry/DATA/STATIONS_FILTERED | 12 +++++++++---
.../sep_bathymetry/DATA/meshfem3D_files/Mesh_Par_file | 12 ++++++------
.../sep_bathymetry/DATA/meshfem3D_files/interfaces.dat | 4 ++--
.../meshfem3D_examples/sep_bathymetry/run_this_example.sh | 6 ++++--
7 files changed, 33 insertions(+), 25 deletions(-)
diff --git a/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/CMTSOLUTION b/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/CMTSOLUTION
index e252010..4538ace 100644
--- a/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/CMTSOLUTION
+++ b/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/CMTSOLUTION
@@ -1,10 +1,10 @@
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
-latitude: 33.6
-longitude: -118.4
-depth: 10.0
+half duration: 0.5
+latitude: 2000.0
+longitude: 5000.0
+depth: 0.01
Mrr: 1.000000e+23
Mtt: 1.000000e+23
Mpp: 1.000000e+23
diff --git a/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/FORCESOLUTION b/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/FORCESOLUTION
index 8c9b49d..c88d980 100644
--- a/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/FORCESOLUTION
+++ b/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/FORCESOLUTION
@@ -1,10 +1,10 @@
FORCE 001
-time shift: 0.0000
+time shift: 1.0000
hdur: 0.5
-latitude: 33.6
-longitude: -118.4
-depth: 10.0
+latitude: 2000.0
+longitude: 5000.0
+depth: 0.01
factor force source: 1.d15
-component dir vect source E: 1.d0
-component dir vect source N: -2.d0
+component dir vect source E: 0.d0
+component dir vect source N: 0.d0
component dir vect source Z_UP: -1.d0
diff --git a/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/Par_file b/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/Par_file
index fdf0f11..086c2e0 100644
--- a/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/Par_file
+++ b/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/Par_file
@@ -15,8 +15,8 @@ SUPPRESS_UTM_PROJECTION = .true.
NPROC = 5
# time step parameters
-NSTEP = 200
-DT = 0.03
+NSTEP = 10000
+DT = 0.0018
# number of nodes for 2D and 3D shape functions for hexahedra
# we use either 8-node mesh elements (bricks) or 27-node elements.
diff --git a/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/STATIONS_FILTERED b/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/STATIONS_FILTERED
index 57ff5c9..0aacd10 100644
--- a/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/STATIONS_FILTERED
+++ b/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/STATIONS_FILTERED
@@ -1,3 +1,9 @@
- K400 CE 0.340370E+02 -0.118178E+03 0.000000E+00 0.100000E+03
- K851 CE 0.340700E+02 -0.118346E+03 0.000000E+00 0.100000E+03
- K853 CE 0.340770E+02 -0.118286E+03 0.000000E+00 0.100000E+03
+ STATION_3 NETWORK 0.200000E+04 0.100000E+04 0.000000E+00 0.500000E+03
+STATION_58 NETWORK 0.200000E+04 0.200000E+04 0.000000E+00 0.500000E+03
+STATION_11 NETWORK 0.200000E+04 0.300000E+04 0.000000E+00 0.500000E+03
+STATION_16 NETWORK 0.200000E+04 0.400000E+04 0.000000E+00 0.500000E+03
+STATION_22 NETWORK 0.200000E+04 0.500000E+04 0.000000E+00 0.500000E+03
+STATION_28 NETWORK 0.200000E+04 0.600000E+04 0.000000E+00 0.500000E+03
+STATION_33 NETWORK 0.200000E+04 0.700000E+04 0.000000E+00 0.500000E+03
+STATION_39 NETWORK 0.200000E+04 0.800000E+04 0.000000E+00 0.500000E+03
+STATION_44 NETWORK 0.200000E+04 0.900000E+04 0.000000E+00 0.500000E+03
diff --git a/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/meshfem3D_files/Mesh_Par_file b/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/meshfem3D_files/Mesh_Par_file
index 9042290..0982f84 100644
--- a/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/meshfem3D_files/Mesh_Par_file
+++ b/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/meshfem3D_files/Mesh_Par_file
@@ -14,8 +14,8 @@ 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 = 200
-NEX_ETA = 8
+NEX_XI = 100
+NEX_ETA = 4
# number of MPI processors along xi and eta (can be different)
NPROC_XI = 5
@@ -44,14 +44,14 @@ NMATERIALS = 2
# Q : quality factor
# anisotropy_flag : 0=no anisotropy/ 1,2,.. check with implementation in aniso_model.f90
# domain_id : 1=acoustic / 2=elastic
-1 10 10000 10000 50.0 0 1
-2 20 20000 20000 40.0 0 2
+1 1000 2500 0 50.0 0 1
+2 2000 4500 2500 40.0 0 2
#1 1100 1600 0 50.0 0 1
#2 1200 1500 750 40.0 0 2
# number of regions
NREGIONS = 2
# 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 200 1 8 1 15 1
-1 200 1 8 16 100 2
+1 100 1 4 1 2 1
+1 100 1 4 3 50 2
diff --git a/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/meshfem3D_files/interfaces.dat b/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/meshfem3D_files/interfaces.dat
index 30ff471..ceb2d94 100644
--- a/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/meshfem3D_files/interfaces.dat
+++ b/EXAMPLES/meshfem3D_examples/sep_bathymetry/DATA/meshfem3D_files/interfaces.dat
@@ -17,6 +17,6 @@
#
# for each layer, we give the number of spectral elements in the vertical direction
# layer number 1 (bottom layer)
- 85
+ 48
# layer number 2 (top layer)
- 15
+ 2
diff --git a/EXAMPLES/meshfem3D_examples/sep_bathymetry/run_this_example.sh b/EXAMPLES/meshfem3D_examples/sep_bathymetry/run_this_example.sh
index 6135546..0843394 100755
--- a/EXAMPLES/meshfem3D_examples/sep_bathymetry/run_this_example.sh
+++ b/EXAMPLES/meshfem3D_examples/sep_bathymetry/run_this_example.sh
@@ -31,7 +31,7 @@ rm -rf OUTPUT_FILES/DATABASES_MPI/*
# compiles executables in root directory
cd ../../../
-#make clean
+make clean
make -j 4 #> $currentdir/tmp.log
cd $currentdir
@@ -51,6 +51,7 @@ echo
echo "running mesher..."
echo
mpirun -np $NPROC ./bin/xmeshfem3D
+#./bin/xmeshfem3D
#mv OUTPUT_FILES/output_mesher.txt OUTPUT_FILES/output_meshfem3D.txt
# runs database generation
@@ -58,12 +59,13 @@ echo
echo "running database generation..."
echo
mpirun -np $NPROC ./bin/xgenerate_databases
+#./bin/xgenerate_databases
# runs simulation
echo
echo " running solver..."
echo
-#mpirun -np $NPROC ./bin/xspecfem3D
+mpirun -np $NPROC ./bin/xspecfem3D
echo
echo "see results in directory: OUTPUT_FILES/"
More information about the CIG-COMMITS
mailing list