[cig-commits] [commit] devel, master: added two sub-examples to handle utm suppression for the socal1D internal mesh example (493663d)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 16:48:50 PDT 2014
Repository : https://github.com/geodynamics/specfem3d
On branches: devel,master
Link : https://github.com/geodynamics/specfem3d/compare/6026e367984905ab133865f62fa6293b343759b9...47f703851338234f96397e7da9fbff63d8178b8a
>---------------------------------------------------------------
commit 493663d1b796d0e6dffae588297e3a5a75fdefb4
Author: Carl Tape <carltape at gi.alaska.edu>
Date: Sun Apr 10 05:14:33 2011 +0000
added two sub-examples to handle utm suppression for the socal1D internal mesh example
>---------------------------------------------------------------
493663d1b796d0e6dffae588297e3a5a75fdefb4
meshfem3D_examples/socal1D/process.sh | 1 +
.../{CMTSOLUTION => utm_example/CMTSOLUTION_utm} | 4 +--
.../Mesh_Par_file_001_utm} | 10 ++++----
.../socal1D/{Par_file => utm_example/Par_file_utm} | 2 +-
meshfem3D_examples/socal1D/utm_example/README | 27 +++++++++++++++++++
.../socal1D/utm_example/STATIONS_utm | 30 ++++++++++++++++++++++
.../interfaces_utm.dat} | 8 +++---
7 files changed, 70 insertions(+), 12 deletions(-)
diff --git a/meshfem3D_examples/socal1D/process.sh b/meshfem3D_examples/socal1D/process.sh
index f6d8eea..d518b9f 100755
--- a/meshfem3D_examples/socal1D/process.sh
+++ b/meshfem3D_examples/socal1D/process.sh
@@ -64,6 +64,7 @@ cd ../
mv in_out_files/OUTPUT_FILES/output_mesher.txt in_out_files/OUTPUT_FILES/output_meshfem3D.txt
# stores setup
+cp -r in_data_files/meshfem3D_files/ in_out_files/OUTPUT_FILES/
cp in_data_files/Par_file in_out_files/OUTPUT_FILES/
cp in_data_files/CMTSOLUTION in_out_files/OUTPUT_FILES/
cp in_data_files/STATIONS in_out_files/OUTPUT_FILES/
diff --git a/meshfem3D_examples/socal1D/CMTSOLUTION b/meshfem3D_examples/socal1D/utm_example/CMTSOLUTION_utm
similarity index 86%
copy from meshfem3D_examples/socal1D/CMTSOLUTION
copy to meshfem3D_examples/socal1D/utm_example/CMTSOLUTION_utm
index 10609f9..8c03466 100644
--- a/meshfem3D_examples/socal1D/CMTSOLUTION
+++ b/meshfem3D_examples/socal1D/utm_example/CMTSOLUTION_utm
@@ -2,8 +2,8 @@ PDE 2001 9 9 23 59 17.78 34.0745 -118.3792 6.4 4.2 4.2 HOLLYWOOD POLET SEM
event name: 9703873
time shift: 0.0000
half duration: 0.3114
-latitude: 34.0745
-longitude: -118.3792
+latitude: 3771080.0
+longitude: 372738.0
depth: 5.4000
Mrr: -0.002000e+23
Mtt: -0.064000e+23
diff --git a/meshfem3D_examples/socal1D/Mesh_Par_file b/meshfem3D_examples/socal1D/utm_example/Mesh_Par_file_001_utm
similarity index 91%
copy from meshfem3D_examples/socal1D/Mesh_Par_file
copy to meshfem3D_examples/socal1D/utm_example/Mesh_Par_file_001_utm
index acb120f..d0eca99 100644
--- a/meshfem3D_examples/socal1D/Mesh_Par_file
+++ b/meshfem3D_examples/socal1D/utm_example/Mesh_Par_file_001_utm
@@ -1,12 +1,12 @@
# 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
+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 = .false.
+SUPPRESS_UTM_PROJECTION = .true.
# file that contains the interfaces of the model / mesh
INTERFACES_FILE = interfaces.dat
diff --git a/meshfem3D_examples/socal1D/Par_file b/meshfem3D_examples/socal1D/utm_example/Par_file_utm
similarity index 97%
copy from meshfem3D_examples/socal1D/Par_file
copy to meshfem3D_examples/socal1D/utm_example/Par_file_utm
index 1b89c1f..6fbdfc9 100644
--- a/meshfem3D_examples/socal1D/Par_file
+++ b/meshfem3D_examples/socal1D/utm_example/Par_file_utm
@@ -6,7 +6,7 @@ SAVE_FORWARD = .false.
# UTM projection parameters
UTM_PROJECTION_ZONE = 11
-SUPPRESS_UTM_PROJECTION = .false.
+SUPPRESS_UTM_PROJECTION = .true.
# number of MPI processors
NPROC = 1
diff --git a/meshfem3D_examples/socal1D/utm_example/README b/meshfem3D_examples/socal1D/utm_example/README
new file mode 100644
index 0000000..3e865a7
--- /dev/null
+++ b/meshfem3D_examples/socal1D/utm_example/README
@@ -0,0 +1,27 @@
+README
+(SPECFEM3D/examples/meshfem3D_examples/socal1D/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
+"flat models" that do not need to take into account the sphericity of Earth
+(for example, the 3D socal model used in Komatitsch et al. 2004).
+
+First, run the default socal1D example with 1 core.
+
+EXAMPLE A: Input geographic bounds in UTM, everything else in lon-lat.
+
+ > cp Mesh_Par_file_001_utm ../Mesh_Par_file
+ > cd ..
+ > ./process.sh
+
+EXAMPLE B: All input values in UTM.
+
+ > cp CMTSOLUTION_utm ../CMTSOLUTION
+ > cp STATIONS_utm ../STATIONS
+ > cp interfaces_utm.dat ../interfaces.dat
+ > cp Mesh_Par_file_001_utm ../Mesh_Par_file
+ > cp Par_file_utm ../Par_file
+ > cd ..
+ > ./process.sh
+
+=======================
\ No newline at end of file
diff --git a/meshfem3D_examples/socal1D/utm_example/STATIONS_utm b/meshfem3D_examples/socal1D/utm_example/STATIONS_utm
new file mode 100644
index 0000000..e86c6c9
--- /dev/null
+++ b/meshfem3D_examples/socal1D/utm_example/STATIONS_utm
@@ -0,0 +1,30 @@
+ 24851 CE 3770540.194093 375794.778375 0.00 0.00
+ 24853 CE 3771245.142806 381341.534900 0.00 0.00
+ K851 CE 3770540.194093 375794.778375 0.00 0.00
+ K853 CE 3771245.142806 381341.534900 0.00 0.00
+ BVH CI 3771287.223592 371199.161741 0.00 0.00
+ CAB CI 3780444.660252 348729.637807 0.00 0.00
+ DJJ CI 3774690.021889 365788.275263 0.00 0.00
+ DJJB CI 3774690.021889 365788.275263 0.00 0.00
+ GR2 CI 3775840.420798 380163.280642 0.00 0.00
+ GRF CI 3775941.393569 380072.320783 0.00 0.00
+ LAF CI 3748222.391683 376849.747924 0.00 0.00
+ LAX CI 3756632.948247 369335.533371 0.00 0.00
+ LCG CI 3762853.919174 372743.220601 0.00 0.00
+ PDR CI 3758759.514975 367228.127192 0.00 0.00
+ SMF CI 3765299.026424 366421.222454 0.00 0.00
+ SMS CI 3764511.850697 365540.030956 0.00 0.00
+ SPF CI 3769762.226999 348077.027970 0.00 0.00
+ TCF CI 3772403.408916 352470.414102 0.00 0.00
+ USC CI 3764835.171833 381232.340025 0.00 0.00
+ WNS CI 3775572.000924 372724.507168 0.00 0.00
+ WTT CI 3756982.464007 383984.272469 0.00 0.00
+ FABS FA 3770467.255447 355014.474196 0.00 0.00
+ 5081 NP 3772430.358145 352378.540706 0.00 0.00
+ 5399 NP 3756538.331182 369602.304617 0.00 0.00
+ 5402 NP 3771266.130124 371281.939298 0.00 0.00
+ 5404 NP 3780363.811004 348721.882056 0.00 0.00
+ BVH NP 3771287.223592 371199.161741 0.00 0.00
+ GRF NP 3775941.393569 380072.320783 0.00 0.00
+ LAX NP 3756632.948247 369335.533371 0.00 0.00
+ TCF NP 3772403.408916 352470.414102 0.00 0.00
diff --git a/meshfem3D_examples/socal1D/interfaces.dat b/meshfem3D_examples/socal1D/utm_example/interfaces_utm.dat
similarity index 76%
copy from meshfem3D_examples/socal1D/interfaces.dat
copy to meshfem3D_examples/socal1D/utm_example/interfaces_utm.dat
index 6395dba..3cc5133 100644
--- a/meshfem3D_examples/socal1D/interfaces.dat
+++ b/meshfem3D_examples/socal1D/utm_example/interfaces_utm.dat
@@ -7,16 +7,16 @@
# 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
+ .true. 2 2 122000.d0 3547000.d0 600000.d0 600000.d0
interface1.dat
# interface number 2
- .false. 2 2 -121.d0 32.d0 6.0d0 6.0d0
+ .true. 2 2 122000.d0 3547000.d0 600000.d0 600000.d0
interface2.dat
# interface number 3
- .false. 2 2 -121.d0 32.d0 6.0d0 6.0d0
+ .true. 2 2 122000.d0 3547000.d0 600000.d0 600000.d0
interface3.dat
# interface number 4 (topography, top of the mesh)
- .false. 2 2 -121.d0 32.d0 6.0d0 6.0d0
+ .true. 2 2 122000.d0 3547000.d0 600000.d0 600000.d0
interface4.dat
#
# for each layer, we give the number of spectral elements in the vertical direction
More information about the CIG-COMMITS
mailing list