[cig-commits] [commit] devel,master: added GMT/perl script for plotting wavefield snapshots; added three example figures; tested all examples. (a23beca)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 15:18:40 PDT 2014
Repository : https://github.com/geodynamics/specfem2d
On branches: devel,master
Link : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c
>---------------------------------------------------------------
commit a23beca6580597f111760e053ef0de55d52338ff
Author: Carl Tape <carltape at gi.alaska.edu>
Date: Mon Apr 11 21:39:22 2011 +0000
added GMT/perl script for plotting wavefield snapshots; added three example figures; tested all examples.
>---------------------------------------------------------------
a23beca6580597f111760e053ef0de55d52338ff
M2_UPPA/README | 6 ++++
M2_UPPA/process.sh | 3 ++
M2_UPPA/wavefield_M2_UPPA_PSV.pdf | Bin 0 -> 33883 bytes
Tape2007/README | 5 ++-
Tape2007/process.sh | 3 ++
Tape2007/wavefield_Tape2007_onerec.pdf | Bin 11306 -> 11406 bytes
Tromp2005/Par_file_Tromp2005 | 2 +-
Tromp2005/README | 6 ++++
Tromp2005/process.sh | 3 ++
Tromp2005/wavefield_Tromp2005_PSV.pdf | 59 +++++++++++++++++++++++++++++++++
Tromp2005/wavefield_Tromp2005_SH.pdf | Bin 0 -> 14712 bytes
11 files changed, 85 insertions(+), 2 deletions(-)
diff --git a/M2_UPPA/README b/M2_UPPA/README
index 8b7768f..f6c5383 100644
--- a/M2_UPPA/README
+++ b/M2_UPPA/README
@@ -20,3 +20,9 @@ TO RUN:
./process.sh
5. check out the output files in the local directory OUTPUT_FILES
+
+6. optional: set output_wavefield_snapshot = .true. in Par_file, re-run,
+ then try plotting the wavefield using the script
+ SPECFEM2D/UTILS/visualization/plot_wavefield.pl
+
+----------------------------------------------------------------------
\ No newline at end of file
diff --git a/M2_UPPA/process.sh b/M2_UPPA/process.sh
index 4b4501a..3f6e5a3 100755
--- a/M2_UPPA/process.sh
+++ b/M2_UPPA/process.sh
@@ -41,6 +41,9 @@ ln -s ../../bin/xspecfem2D
# stores setup
cp DATA/Par_file OUTPUT_FILES/
cp DATA/SOURCE OUTPUT_FILES/
+cp DATA/SOURCE_xz.dat OUTPUT_FILES/
+cp DATA/STATIONS OUTPUT_FILES/
+cp DATA/STATIONS_target OUTPUT_FILES/
# runs database generation
echo
diff --git a/M2_UPPA/wavefield_M2_UPPA_PSV.pdf b/M2_UPPA/wavefield_M2_UPPA_PSV.pdf
new file mode 100644
index 0000000..a8c8e88
Binary files /dev/null and b/M2_UPPA/wavefield_M2_UPPA_PSV.pdf differ
diff --git a/Tape2007/README b/Tape2007/README
index dda3fe1..8afd6f7 100644
--- a/Tape2007/README
+++ b/Tape2007/README
@@ -25,7 +25,7 @@ TO RUN:
1. in SPECFEM2D root directory, configure, e.g.,
./configure FC=gfortran
-2. modify constants.h file, setting USER_TO = 48 s
+2. modify setup/constants.h file, setting USER_TO = 48.0 s
3. compile:
make all
@@ -37,4 +37,7 @@ TO RUN:
6. check out the output files in the local directory OUTPUT_FILES
+ optional: try plotting the wavefield using the script
+ SPECFEM2D/UTILS/visualization/plot_wavefield.pl
+
---------------------------
diff --git a/Tape2007/process.sh b/Tape2007/process.sh
index bcb1f4b..fed34ba 100755
--- a/Tape2007/process.sh
+++ b/Tape2007/process.sh
@@ -41,6 +41,9 @@ ln -s ../../bin/xspecfem2D
# stores setup
cp DATA/Par_file OUTPUT_FILES/
cp DATA/SOURCE OUTPUT_FILES/
+cp DATA/SOURCE_xz.dat OUTPUT_FILES/
+cp DATA/STATIONS OUTPUT_FILES/
+cp DATA/STATIONS_target OUTPUT_FILES/
# runs database generation
echo
diff --git a/Tape2007/wavefield_Tape2007_onerec.pdf b/Tape2007/wavefield_Tape2007_onerec.pdf
index 405a60a..8814670 100644
Binary files a/Tape2007/wavefield_Tape2007_onerec.pdf and b/Tape2007/wavefield_Tape2007_onerec.pdf differ
diff --git a/Tromp2005/Par_file_Tromp2005 b/Tromp2005/Par_file_Tromp2005
index 41a7c1d..86b134c 100644
--- a/Tromp2005/Par_file_Tromp2005
+++ b/Tromp2005/Par_file_Tromp2005
@@ -63,7 +63,7 @@ sizemax_arrows = 1.d0 # maximum size of arrows on vec
gnuplot = .false. # generate a GNUPLOT file for the grid
output_grid = .false. # save the grid in a text file or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
+output_wavefield_snapshot = .true. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
diff --git a/Tromp2005/README b/Tromp2005/README
index 3139457..1aa0118 100644
--- a/Tromp2005/README
+++ b/Tromp2005/README
@@ -35,6 +35,9 @@ TO RUN:
The seismograms (S0001.AA.BHX.semd,S0001.AA.BHZ.semd) will be similar to those in Tromp2005, Figure 7b, but with an 8-second time shift due to the difference in origin time convention.
+ optional: try plotting the wavefield using the script
+ SPECFEM2D/UTILS/visualization/plot_wavefield.pl
+
7. mv OUTPUT_FILES OUTPUT_FILES_PSV
8. set p_sv = .false. in Par_file_Tromp2005
@@ -46,4 +49,7 @@ TO RUN:
The seismogram (S0001.AA.BHY.semd) will be similar to those in Tromp2005, Figure 2b, but with an 8-second time shift due to the difference in origin time convention.
+ optional: try plotting the wavefield using the script
+ SPECFEM2D/UTILS/visualization/plot_wavefield.pl
+
---------------------------
diff --git a/Tromp2005/process.sh b/Tromp2005/process.sh
index d48fd42..e4dbd60 100755
--- a/Tromp2005/process.sh
+++ b/Tromp2005/process.sh
@@ -41,6 +41,9 @@ ln -s ../../bin/xspecfem2D
# stores setup
cp DATA/Par_file OUTPUT_FILES/
cp DATA/SOURCE OUTPUT_FILES/
+cp DATA/SOURCE_xz.dat OUTPUT_FILES/
+cp DATA/STATIONS OUTPUT_FILES/
+cp DATA/STATIONS_target OUTPUT_FILES/
# runs database generation
echo
diff --git a/Tromp2005/wavefield_Tromp2005_PSV.pdf b/Tromp2005/wavefield_Tromp2005_PSV.pdf
new file mode 100644
index 0000000..f80c501
--- /dev/null
+++ b/Tromp2005/wavefield_Tromp2005_PSV.pdf
@@ -0,0 +1,59 @@
+%PDF-1.4
+%Çì¢
+4 0 obj
+<</Type/Page/MediaBox [0 0 595 842]
+/Parent 3 0 R
+/Resources<</ProcSet[/PDF]
+>>
+>>
+endobj
+3 0 obj
+<< /Type /Pages /Kids [
+4 0 R
+] /Count 1
+>>
+endobj
+1 0 obj
+<</Type /Catalog /Pages 3 0 R
+/Metadata 5 0 R
+>>
+endobj
+5 0 obj
+<</Type/Metadata
+/Subtype/XML/Length 1319>>stream
+<?xpacket begin='' id='W5M0MpCehiHzreSzNTczkc9d'?>
+<?adobe-xap-filters esc="CRLF"?>
+<x:xmpmeta xmlns:x='adobe:ns:meta/' x:xmptk='XMP toolkit 2.9.1-13, framework 1.6'>
+<rdf:RDF xmlns:rdf='http://www.w3.org/1999/02/22-rdf-syntax-ns#' xmlns:iX='http://ns.adobe.com/iX/1.0/'>
+<rdf:Description rdf:about='6d045576-9c94-11eb-0000-8616d0889b02' xmlns:pdf='http://ns.adobe.com/pdf/1.3/' pdf:Producer='GPL Ghostscript 8.70'/>
+<rdf:Description rdf:about='6d045576-9c94-11eb-0000-8616d0889b02' xmlns:xmp='http://ns.adobe.com/xap/1.0/'><xmp:ModifyDate>2011-04-11T12:11:25-08:00</xmp:ModifyDate>
+<xmp:CreateDate>2011-04-11T12:11:25-08:00</xmp:CreateDate>
+<xmp:CreatorTool>UnknownApplication</xmp:CreatorTool></rdf:Description>
+<rdf:Description rdf:about='6d045576-9c94-11eb-0000-8616d0889b02' xmlns:xapMM='http://ns.adobe.com/xap/1.0/mm/' xapMM:DocumentID='6d045576-9c94-11eb-0000-8616d0889b02'/>
+<rdf:Description rdf:about='6d045576-9c94-11eb-0000-8616d0889b02' xmlns:dc='http://purl.org/dc/elements/1.1/' dc:format='application/pdf'><dc:title><rdf:Alt><rdf:li xml:lang='x-default'>Untitled</rdf:li></rdf:Alt></dc:title></rdf:Description>
+</rdf:RDF>
+</x:xmpmeta>
+
+
+<?xpacket end='w'?>
+endstream
+endobj
+2 0 obj
+<</Producer(GPL Ghostscript 8.70)
+/CreationDate(D:20110411121125-08'00')
+/ModDate(D:20110411121125-08'00')>>endobj
+xref
+0 6
+0000000000 65535 f
+0000000172 00000 n
+0000001631 00000 n
+0000000113 00000 n
+0000000015 00000 n
+0000000236 00000 n
+trailer
+<< /Size 6 /Root 1 0 R /Info 2 0 R
+/ID [<D9361AA852B03FB91A83BBF3BE843AAC><D9361AA852B03FB91A83BBF3BE843AAC>]
+>>
+startxref
+1754
+%%EOF
diff --git a/Tromp2005/wavefield_Tromp2005_SH.pdf b/Tromp2005/wavefield_Tromp2005_SH.pdf
new file mode 100644
index 0000000..16485a1
Binary files /dev/null and b/Tromp2005/wavefield_Tromp2005_SH.pdf differ
More information about the CIG-COMMITS
mailing list