[cig-commits] r19240 - in seismo/3D/SPECFEM3D_GLOBE/trunk: DATA UTILS/Visualization/Gnuplot

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2011-11-23 05:09:13 -0800 (Wed, 23 Nov 2011)
New Revision: 19240

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/CMTSOLUTION.in
   seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/Visualization/Gnuplot/create_gnuplot_script_to_compare_seismograms.csh
Log:
improved the Gnuplot script to compare seismograms


Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/CMTSOLUTION.in
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/CMTSOLUTION.in	2011-11-22 22:58:23 UTC (rev 19239)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/CMTSOLUTION.in	2011-11-23 13:09:13 UTC (rev 19240)
@@ -10,4 +10,4 @@
 Mpp:      -1.600000e+26
 Mrt:      -2.503000e+28
 Mrp:       4.200000e+26
-Mtp:      -2.480000e+27
\ No newline at end of file
+Mtp:      -2.480000e+27

Modified: 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	2011-11-22 22:58:23 UTC (rev 19239)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/Visualization/Gnuplot/create_gnuplot_script_to_compare_seismograms.csh	2011-11-23 13:09:13 UTC (rev 19240)
@@ -1,32 +1,26 @@
 #!/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)
+#echo set term postscript color solid "Helvetica" 22
+#echo set output \"all_seismograms_comparison.ps\"
 
-# 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
+echo set term pdf color solid
+echo set output \"all_seismograms_comparison.pdf\"
 
-# 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 x11"
+#echo "set term wxt"
 
-#echo set term postscript color solid "Helvetica" 22
-#echo set output \"tutu.ps\"
+echo set xrange \[0:3300\]
 
-echo set term X11
+foreach file ( OUTPUT_FILES/*.sem.ascii )
 
-#echo set xrange \[0:14000\]
+  set newfile = `basename $file`
 
-foreach file ($*)
+  echo plot \"OUTPUT_FILES_reference_OK/$newfile\" w l lc 1, \"$file\" w l lc 3
+# use the line below instead if you have an old install of Gnuplot (the Gnuplot syntax was different back then)
+# echo plot \"OUTPUT_FILES_reference_OK/$newfile\" w l 1, \"$file\" w l 3
 
-# 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...\"
+# uncomment this only when outputting to the screen (X11 or wxt)
+#  echo "pause -1 'hit any key...'"
 
 end