[cig-commits] r19155 - seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/Visualization/Gnuplot

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2011-11-07 14:49:02 -0800 (Mon, 07 Nov 2011)
New Revision: 19155

Added:
   seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/Visualization/Gnuplot/create_gnuplot_script_to_compare_seismograms.csh
Log:
created a shell script to automatically create a Gnuplot script to compare two sets of seismograms (for instance to validate a small modification or bug fix in the code by checking that the seismograms computed are unchanged)


Added: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/Visualization/Gnuplot/create_gnuplot_script_to_compare_seismograms.csh
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/Visualization/Gnuplot/create_gnuplot_script_to_compare_seismograms.csh	                        (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/Visualization/Gnuplot/create_gnuplot_script_to_compare_seismograms.csh	2011-11-07 22:49:02 UTC (rev 19155)
@@ -0,0 +1,32 @@
+#!/bin/csh
+
+# script by Dimitri Komatitsch, CNRS Marseille, November 2011
+# to compare two sets of seismograms (for instance to make sure that a modification made in the code does not change the seismograms)
+
+# call it with:
+#   ./create_gnuplot_script_to_compare_seismograms.csh OUTPUT_FILES/*.semd > plotall.gnu
+# (assuming that the seismogram files are called *.semd in directory OUTPUT_FILES)
+# and place the reference seismograms to compare to in directory REFERENCE/OUTPUT_FILES
+
+# also make sure that parameter "hdur" is greater than zero in file DATA/CMTSOLUTION
+# because only seismograms convolved with a source time function should be compared
+# (seismograms computed for a Heaviside source time function before convolution cannot be reliably compared
+# because they are too sensitive to roundoff noise)
+
+#echo set term postscript color solid "Helvetica" 22
+#echo set output \"tutu.ps\"
+
+echo set term X11
+
+#echo set xrange \[0:14000\]
+
+foreach file ($*)
+
+# change "w l 1" to "w l lc 1" to use more recent Gnuplot syntax
+# (i.e., draw a line between data points, with line color 1)
+# (same thing with "w l 3")
+echo plot \"REFERENCE/$file\" w l 1, \"$file\" w l 3
+echo pause -1 \"Hit any key to see the next seismogram...\"
+
+end
+


Property changes on: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/Visualization/Gnuplot/create_gnuplot_script_to_compare_seismograms.csh
___________________________________________________________________
Name: svn:executable
   + *