[cig-commits] [commit] devel, master: done cleaning all the processing and plotting scripts in EXAMPLES/benchmarks/prem_* from Anne Sieminski (0faa0e3)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Thu Nov 6 08:08:12 PST 2014
Repository : https://github.com/geodynamics/specfem3d_globe
On branches: devel,master
Link : https://github.com/geodynamics/specfem3d_globe/compare/bc58e579b3b0838a0968725a076f5904845437ca...be63f20cbb6f462104e949894dbe205d2398cd7f
>---------------------------------------------------------------
commit 0faa0e3f1de49cbd03bd6cb8d0e4476b2365d57a
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date: Sat Sep 28 23:07:20 2013 +0000
done cleaning all the processing and plotting scripts in EXAMPLES/benchmarks/prem_* from Anne Sieminski
>---------------------------------------------------------------
0faa0e3f1de49cbd03bd6cb8d0e4476b2365d57a
...SOLUTION => CMTSOLUTION_to_filter_nonzero_hdur} | 0
.../QMXD/DATA/cmd_line | 2 +-
.../SEMD/DATA/Par_file | 2 +-
...reate_gnuplot_script_to_compare_seismograms.csh | 27 +++++
.../plotall.gnu | 109 +++++++++++++++++++++
.../process_all_seismograms.csh | 29 ++++++
.../processing | 10 --
.../processing_PAS | 10 --
...SOLUTION => CMTSOLUTION_to_filter_nonzero_hdur} | 0
.../QMXD/DATA/cmd_line | 2 +-
.../SEMD/DATA/Par_file | 2 +-
...reate_gnuplot_script_to_compare_seismograms.csh | 27 +++++
.../prem_ani_noocean_Q_Cowling_Vanuatu/plotall.gnu | 109 +++++++++++++++++++++
.../process_all_seismograms.csh | 29 ++++++
.../prem_ani_noocean_Q_Cowling_Vanuatu/processing | 10 --
...SOLUTION => CMTSOLUTION_to_filter_nonzero_hdur} | 0
.../QMXD/DATA/cmd_line | 2 +-
.../SEMD/DATA/Par_file | 2 +-
...reate_gnuplot_script_to_compare_seismograms.csh | 27 +++++
.../prem_iso_noocean_Q_Cowling_Vanuatu/plotall.gnu | 109 +++++++++++++++++++++
.../process_all_seismograms.csh | 29 ++++++
.../prem_iso_noocean_Q_Cowling_Vanuatu/processing | 11 ---
22 files changed, 501 insertions(+), 47 deletions(-)
diff --git a/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/CMTSOLUTION b/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/CMTSOLUTION_to_filter_nonzero_hdur
similarity index 100%
rename from benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/CMTSOLUTION
rename to benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/CMTSOLUTION_to_filter_nonzero_hdur
diff --git a/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/QMXD/DATA/cmd_line b/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/QMXD/DATA/cmd_line
index 0939f36..0b8a3b4 100644
--- a/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/QMXD/DATA/cmd_line
+++ b/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/QMXD/DATA/cmd_line
@@ -1 +1 @@
-./calc_modes_qmxd_no_ocean_prem_ani_noocean_Q_Cowling.pl -m CMTSOLUTION_BOLIVIA -S STATIONS -O SYNT_ANI -C -s 4.768 -n 14200 -c LHZ/LHE/LHN
+./calc_modes_qmxd_no_ocean_prem_ani_noocean_Q_Cowling.pl -m CMTSOLUTION_BOLIVIA -S STATIONS -O SYNT_ANI -C -s 4.768 -n 14200 -c MXZ/MXE/MXN
diff --git a/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA/Par_file b/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA/Par_file
index b2e71fe..8bd6a84 100644
--- a/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA/Par_file
+++ b/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/SEMD/DATA/Par_file
@@ -58,7 +58,7 @@ RECORD_LENGTH_IN_MINUTES = 50.0d0
# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
# but requires a significant amount of disk space for temporary storage.
ATTENUATION_1D_WITH_3D_STORAGE = .true.
-PARTIAL_PHYS_DISPERSION_ONLY = .true.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
UNDO_ATTENUATION = .false.
NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
diff --git a/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/create_gnuplot_script_to_compare_seismograms.csh b/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/create_gnuplot_script_to_compare_seismograms.csh
new file mode 100755
index 0000000..32e4a84
--- /dev/null
+++ b/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/create_gnuplot_script_to_compare_seismograms.csh
@@ -0,0 +1,27 @@
+#!/bin/csh
+
+echo \#set term postscript color solid "Helvetica" 22
+echo \#set output \"all_seismograms_comparison.ps\"
+
+echo set term pdf color solid
+echo set output \"all_seismograms_comparison.pdf\"
+
+#echo set terminal postscript eps color solid
+#echo "set output '| epstopdf --filter > all_seismograms_comparison.pdf'"
+
+#echo "set term x11"
+#echo "set term wxt"
+
+#echo set xrange \[0:2200\]
+
+foreach file ( SEMD/*.sem.ascii.sac.asciinew )
+
+ set newfile = `basename $file .sem.ascii.sac.asciinew`
+
+ echo plot \"QMXD/$newfile.qmxd.sac.asciinew\" w l lc 1, \"$file\" w l lc 3
+
+# uncomment this only when outputting to the screen (X11 or wxt)
+# echo "pause -1 'hit any key...'"
+
+end
+
diff --git a/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/plotall.gnu b/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/plotall.gnu
new file mode 100644
index 0000000..7c1f6f7
--- /dev/null
+++ b/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/plotall.gnu
@@ -0,0 +1,109 @@
+#set term postscript color solid Helvetica 22
+#set output "all_seismograms_comparison.ps"
+set term pdf color solid
+set output "all_seismograms_comparison.pdf"
+plot "QMXD/AIS.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/AIS.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/AIS.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/AIS.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/AIS.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/AIS.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ATD.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/ATD.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ATD.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/ATD.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ATD.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/ATD.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BBSR.IU.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/BBSR.IU.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BBSR.IU.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/BBSR.IU.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BBSR.IU.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/BBSR.IU.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BNG.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/BNG.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BNG.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/BNG.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BNG.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/BNG.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CAN.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/CAN.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CAN.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/CAN.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CAN.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/CAN.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CRZF.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/CRZF.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CRZF.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/CRZF.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CRZF.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/CRZF.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/DRV.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/DRV.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/DRV.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/DRV.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/DRV.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/DRV.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ECH.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/ECH.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ECH.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/ECH.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ECH.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/ECH.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/EVO.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/EVO.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/EVO.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/EVO.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/EVO.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/EVO.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/FDF.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/FDF.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/FDF.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/FDF.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/FDF.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/FDF.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HDC.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/HDC.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HDC.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/HDC.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HDC.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/HDC.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HYB.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/HYB.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HYB.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/HYB.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HYB.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/HYB.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/INU.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/INU.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/INU.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/INU.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/INU.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/INU.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KIP.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/KIP.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KIP.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/KIP.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KIP.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/KIP.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KOG.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/KOG.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KOG.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/KOG.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KOG.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/KOG.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/NOUC.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/NOUC.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/NOUC.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/NOUC.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/NOUC.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/NOUC.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAF.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/PAF.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAF.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/PAF.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAF.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/PAF.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAS.CI.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/PAS.CI.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAS.CI.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/PAS.CI.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAS.CI.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/PAS.CI.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PEL.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/PEL.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PEL.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/PEL.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PEL.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/PEL.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PPT.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/PPT.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PPT.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/PPT.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PPT.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/PPT.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PVC.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/PVC.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PVC.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/PVC.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PVC.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/PVC.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RAO.IU.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/RAO.IU.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RAO.IU.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/RAO.IU.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RAO.IU.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/RAO.IU.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RER.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/RER.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RER.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/RER.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RER.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/RER.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SAML.IU.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SAML.IU.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SAML.IU.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SAML.IU.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SAML.IU.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SAML.IU.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SCZ.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SCZ.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SCZ.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SCZ.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SCZ.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SCZ.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SEY.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SEY.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SEY.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SEY.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SEY.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SEY.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SFJD.IU.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SFJD.IU.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SFJD.IU.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SFJD.IU.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SFJD.IU.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SFJD.IU.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SPB.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SPB.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SPB.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SPB.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SPB.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SPB.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB2.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB2.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB2.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB2.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB2.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB2.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB3.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB3.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB3.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB3.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB3.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB3.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TAM.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/TAM.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TAM.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/TAM.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TAM.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/TAM.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TRIS.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/TRIS.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TRIS.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/TRIS.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TRIS.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/TRIS.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/UNM.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/UNM.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/UNM.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/UNM.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/UNM.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/UNM.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/WUS.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/WUS.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/WUS.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/WUS.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/WUS.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/WUS.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
diff --git a/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/process_all_seismograms.csh b/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/process_all_seismograms.csh
new file mode 100755
index 0000000..3884466
--- /dev/null
+++ b/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/process_all_seismograms.csh
@@ -0,0 +1,29 @@
+#!/bin/csh
+
+# you will find the Perl and "sac2asc" processing scripts below in directory utils/seis_process
+# and you will need the Seismic Analysis Code (SAC) software http://www.iris.edu/dms/nodes/dmc/software/downloads/sac
+# (i.e. the "sac" and "saclst" executable commands) to be installed and executable on your machine
+# for this processing script to run
+
+# Dimitri Komatitsch, CNRS, Marseille, France, September 2013
+
+########### SEM spectral-element results ###########
+
+./process_syn.pl -m CMTSOLUTION_to_filter_nonzero_hdur -a STATIONS -l 0/2900 -t 10/250 SEMD/*.sem.ascii
+
+./rotate.pl -l 0 -L 2900 SEMD/*.MXE.sem.ascii.sac
+
+foreach file ( SEMD/*.MXR.sem.ascii.sac SEMD/*.MXT.sem.ascii.sac SEMD/*.MXZ.sem.ascii.sac )
+ ./sac2asc $file > $file.asciinew
+end
+
+########### QMXD normal-mode results ###########
+
+./process_syn.pl -m CMTSOLUTION_to_filter_nonzero_hdur -a STATIONS -l 0/2900 -t 10/250 QMXD/*.qmxd
+
+./rotate.pl -l 0 -L 2900 QMXD/*.MXE.qmxd.sac
+
+foreach file ( QMXD/*.MXR.qmxd.sac QMXD/*.MXT.qmxd.sac QMXD/*.MXZ.qmxd.sac )
+ ./sac2asc $file > $file.asciinew
+end
+
diff --git a/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/processing b/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/processing
deleted file mode 100755
index b48ae6a..0000000
--- a/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/processing
+++ /dev/null
@@ -1,10 +0,0 @@
-./process_syn.pl -m CMTSOLUTION -a STATIONS -l 0/2900 -t 10/250 QMXD/*.qmxd
-./process_syn.pl -m CMTSOLUTION -a STATIONS -l 0/2900 -t 10/250 SEMD/*.semd.ascii
-./rotate.pl -l 0 -L 2900 QMXD/*.LHE.qmxd.sac
-./rotate.pl -l 0 -L 2900 SEMD/*.LHE.semd.ascii.sac
-mv QMXD/*LHZ*sac QMXD_SAC/
-mv QMXD/*LHT*sac QMXD_SAC/
-mv SEMD/*LHZ*sac SEMD_SAC/
-mv SEMD/*LHT*sac SEMD_SAC/
-tar cvf QMXD_SAC.tar QMXD_SAC/
-tar cvf SEMD_SAC.tar SEMD_SAC/
diff --git a/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/processing_PAS b/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/processing_PAS
deleted file mode 100755
index c4966fd..0000000
--- a/benchmarks/prem_ani_noocean_Q_Cowling_Bolivia_9s/processing_PAS
+++ /dev/null
@@ -1,10 +0,0 @@
-./process_syn.pl -m CMTSOLUTION -a STATIONS -l 0/2900 -t 10/500 QMXD/PAS*.qmxd
-./process_syn.pl -m CMTSOLUTION -a STATIONS -l 0/2900 -t 10/500 SEMD/PAS*.semd.ascii
-./rotate.pl -l 0 -L 2900 QMXD/PAS*.LHE.qmxd.sac
-./rotate.pl -l 0 -L 2900 SEMD/PAS*.LHE.semd.ascii.sac
-mv QMXD/PAS*LHZ*sac QMXD_SAC_PAS/
-mv QMXD/PAS*LHT*sac QMXD_SAC_PAS/
-mv SEMD/PAS*LHZ*sac SEMD_SAC_PAS/
-mv SEMD/PAS*LHT*sac SEMD_SAC_PAS/
-tar cvf QMXD_SAC_PAS.tar QMXD_SAC_PAS/
-tar cvf SEMD_SAC_PAS.tar SEMD_SAC_PAS/
diff --git a/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/CMTSOLUTION b/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/CMTSOLUTION_to_filter_nonzero_hdur
similarity index 100%
rename from benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/CMTSOLUTION
rename to benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/CMTSOLUTION_to_filter_nonzero_hdur
diff --git a/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/QMXD/DATA/cmd_line b/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/QMXD/DATA/cmd_line
index e6cc958..1bd1bbc 100644
--- a/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/QMXD/DATA/cmd_line
+++ b/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/QMXD/DATA/cmd_line
@@ -1 +1 @@
-./calc_modes_qmxd_no_ocean_prem_ani_noocean_Q_Cowling.pl -m CMTSOLUTION_VANUATU -S STATIONS -O SYNT_ANI -C -s 2.0 -n 12100 -c LHZ/LHE/LHN
+./calc_modes_qmxd_no_ocean_prem_ani_noocean_Q_Cowling.pl -m CMTSOLUTION_VANUATU -S STATIONS -O SYNT_ANI -C -s 2.0 -n 12100 -c MXZ/MXE/MXN
diff --git a/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file b/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
index 49a47ec..b0a82b3 100644
--- a/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
+++ b/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
@@ -58,7 +58,7 @@ RECORD_LENGTH_IN_MINUTES = 100.0d0
# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
# but requires a significant amount of disk space for temporary storage.
ATTENUATION_1D_WITH_3D_STORAGE = .true.
-PARTIAL_PHYS_DISPERSION_ONLY = .true.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
UNDO_ATTENUATION = .false.
NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
diff --git a/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/create_gnuplot_script_to_compare_seismograms.csh b/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/create_gnuplot_script_to_compare_seismograms.csh
new file mode 100755
index 0000000..32e4a84
--- /dev/null
+++ b/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/create_gnuplot_script_to_compare_seismograms.csh
@@ -0,0 +1,27 @@
+#!/bin/csh
+
+echo \#set term postscript color solid "Helvetica" 22
+echo \#set output \"all_seismograms_comparison.ps\"
+
+echo set term pdf color solid
+echo set output \"all_seismograms_comparison.pdf\"
+
+#echo set terminal postscript eps color solid
+#echo "set output '| epstopdf --filter > all_seismograms_comparison.pdf'"
+
+#echo "set term x11"
+#echo "set term wxt"
+
+#echo set xrange \[0:2200\]
+
+foreach file ( SEMD/*.sem.ascii.sac.asciinew )
+
+ set newfile = `basename $file .sem.ascii.sac.asciinew`
+
+ echo plot \"QMXD/$newfile.qmxd.sac.asciinew\" w l lc 1, \"$file\" w l lc 3
+
+# uncomment this only when outputting to the screen (X11 or wxt)
+# echo "pause -1 'hit any key...'"
+
+end
+
diff --git a/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/plotall.gnu b/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/plotall.gnu
new file mode 100644
index 0000000..7c1f6f7
--- /dev/null
+++ b/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/plotall.gnu
@@ -0,0 +1,109 @@
+#set term postscript color solid Helvetica 22
+#set output "all_seismograms_comparison.ps"
+set term pdf color solid
+set output "all_seismograms_comparison.pdf"
+plot "QMXD/AIS.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/AIS.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/AIS.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/AIS.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/AIS.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/AIS.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ATD.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/ATD.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ATD.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/ATD.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ATD.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/ATD.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BBSR.IU.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/BBSR.IU.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BBSR.IU.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/BBSR.IU.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BBSR.IU.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/BBSR.IU.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BNG.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/BNG.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BNG.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/BNG.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BNG.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/BNG.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CAN.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/CAN.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CAN.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/CAN.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CAN.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/CAN.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CRZF.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/CRZF.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CRZF.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/CRZF.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CRZF.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/CRZF.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/DRV.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/DRV.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/DRV.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/DRV.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/DRV.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/DRV.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ECH.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/ECH.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ECH.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/ECH.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ECH.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/ECH.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/EVO.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/EVO.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/EVO.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/EVO.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/EVO.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/EVO.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/FDF.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/FDF.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/FDF.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/FDF.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/FDF.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/FDF.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HDC.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/HDC.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HDC.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/HDC.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HDC.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/HDC.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HYB.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/HYB.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HYB.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/HYB.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HYB.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/HYB.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/INU.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/INU.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/INU.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/INU.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/INU.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/INU.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KIP.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/KIP.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KIP.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/KIP.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KIP.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/KIP.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KOG.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/KOG.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KOG.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/KOG.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KOG.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/KOG.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/NOUC.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/NOUC.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/NOUC.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/NOUC.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/NOUC.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/NOUC.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAF.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/PAF.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAF.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/PAF.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAF.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/PAF.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAS.CI.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/PAS.CI.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAS.CI.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/PAS.CI.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAS.CI.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/PAS.CI.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PEL.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/PEL.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PEL.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/PEL.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PEL.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/PEL.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PPT.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/PPT.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PPT.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/PPT.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PPT.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/PPT.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PVC.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/PVC.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PVC.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/PVC.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PVC.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/PVC.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RAO.IU.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/RAO.IU.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RAO.IU.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/RAO.IU.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RAO.IU.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/RAO.IU.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RER.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/RER.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RER.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/RER.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RER.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/RER.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SAML.IU.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SAML.IU.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SAML.IU.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SAML.IU.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SAML.IU.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SAML.IU.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SCZ.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SCZ.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SCZ.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SCZ.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SCZ.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SCZ.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SEY.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SEY.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SEY.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SEY.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SEY.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SEY.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SFJD.IU.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SFJD.IU.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SFJD.IU.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SFJD.IU.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SFJD.IU.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SFJD.IU.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SPB.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SPB.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SPB.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SPB.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SPB.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SPB.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB2.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB2.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB2.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB2.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB2.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB2.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB3.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB3.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB3.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB3.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB3.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB3.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TAM.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/TAM.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TAM.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/TAM.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TAM.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/TAM.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TRIS.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/TRIS.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TRIS.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/TRIS.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TRIS.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/TRIS.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/UNM.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/UNM.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/UNM.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/UNM.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/UNM.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/UNM.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/WUS.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/WUS.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/WUS.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/WUS.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/WUS.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/WUS.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
diff --git a/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/process_all_seismograms.csh b/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/process_all_seismograms.csh
new file mode 100755
index 0000000..c513497
--- /dev/null
+++ b/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/process_all_seismograms.csh
@@ -0,0 +1,29 @@
+#!/bin/csh
+
+# you will find the Perl and "sac2asc" processing scripts below in directory utils/seis_process
+# and you will need the Seismic Analysis Code (SAC) software http://www.iris.edu/dms/nodes/dmc/software/downloads/sac
+# (i.e. the "sac" and "saclst" executable commands) to be installed and executable on your machine
+# for this processing script to run
+
+# Dimitri Komatitsch, CNRS, Marseille, France, September 2013
+
+########### SEM spectral-element results ###########
+
+./process_syn.pl -m CMTSOLUTION_to_filter_nonzero_hdur -a STATIONS -l 0/6000 -t 150/500 SEMD/*.sem.ascii
+
+./rotate.pl -l 0 -L 6000 SEMD/*.MXE.sem.ascii.sac
+
+foreach file ( SEMD/*.MXR.sem.ascii.sac SEMD/*.MXT.sem.ascii.sac SEMD/*.MXZ.sem.ascii.sac )
+ ./sac2asc $file > $file.asciinew
+end
+
+########### QMXD normal-mode results ###########
+
+./process_syn.pl -m CMTSOLUTION_to_filter_nonzero_hdur -a STATIONS -l 0/6000 -t 150/500 QMXD/*.qmxd
+
+./rotate.pl -l 0 -L 6000 QMXD/*.MXE.qmxd.sac
+
+foreach file ( QMXD/*.MXR.qmxd.sac QMXD/*.MXT.qmxd.sac QMXD/*.MXZ.qmxd.sac )
+ ./sac2asc $file > $file.asciinew
+end
+
diff --git a/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/processing b/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/processing
deleted file mode 100755
index 8f7244a..0000000
--- a/benchmarks/prem_ani_noocean_Q_Cowling_Vanuatu/processing
+++ /dev/null
@@ -1,10 +0,0 @@
-./process_syn.pl -m CMTSOLUTION -a STATIONS -l 0/6000 -t 150/500 QMXD/*.qmxd
-./process_syn.pl -m CMTSOLUTION -a STATIONS -l 0/6000 -t 150/500 SEMD/*.semd.ascii
-./rotate.pl -l 0 -L 6000 QMXD/*.LHE.qmxd.sac
-./rotate.pl -l 0 -L 6000 SEMD/*.LHE.semd.ascii.sac
-mv QMXD/*LHZ*sac QMXD_SAC/
-mv QMXD/*LHT*sac QMXD_SAC/
-mv SEMD/*LHZ*sac SEMD_SAC/
-mv SEMD/*LHT*sac SEMD_SAC/
-tar cvf QMXD_SAC.tar QMXD_SAC/
-tar cvf SEMD_SAC.tar SEMD_SAC/
diff --git a/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/CMTSOLUTION b/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/CMTSOLUTION_to_filter_nonzero_hdur
similarity index 100%
rename from benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/CMTSOLUTION
rename to benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/CMTSOLUTION_to_filter_nonzero_hdur
diff --git a/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/QMXD/DATA/cmd_line b/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/QMXD/DATA/cmd_line
index 4acf096..cd38c8f 100644
--- a/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/QMXD/DATA/cmd_line
+++ b/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/QMXD/DATA/cmd_line
@@ -1 +1 @@
-./calc_modes_qmxd_no_ocean_prem_iso_noocean_Q_Cowling.pl -m CMTSOLUTION_VANUATU -S STATIONS -O SYNT_ISO -C -s 2.0 -n 12100 -c LHZ/LHE/LHN
+./calc_modes_qmxd_no_ocean_prem_iso_noocean_Q_Cowling.pl -m CMTSOLUTION_VANUATU -S STATIONS -O SYNT_ISO -C -s 2.0 -n 12100 -c MXZ/MXE/MXN
diff --git a/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file b/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
index 8422f4d..184f0d9 100644
--- a/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
+++ b/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/SEMD/DATA/Par_file
@@ -58,7 +58,7 @@ RECORD_LENGTH_IN_MINUTES = 100.0d0
# use one (and only one) of the two flags below. UNDO_ATTENUATION is much better (it is exact)
# but requires a significant amount of disk space for temporary storage.
ATTENUATION_1D_WITH_3D_STORAGE = .true.
-PARTIAL_PHYS_DISPERSION_ONLY = .true.
+PARTIAL_PHYS_DISPERSION_ONLY = .false.
UNDO_ATTENUATION = .false.
NT_DUMP_ATTENUATION = 100 # how often we dump restart files to undo attenuation, only needed when using UNDO_ATTENUATION
diff --git a/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/create_gnuplot_script_to_compare_seismograms.csh b/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/create_gnuplot_script_to_compare_seismograms.csh
new file mode 100755
index 0000000..32e4a84
--- /dev/null
+++ b/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/create_gnuplot_script_to_compare_seismograms.csh
@@ -0,0 +1,27 @@
+#!/bin/csh
+
+echo \#set term postscript color solid "Helvetica" 22
+echo \#set output \"all_seismograms_comparison.ps\"
+
+echo set term pdf color solid
+echo set output \"all_seismograms_comparison.pdf\"
+
+#echo set terminal postscript eps color solid
+#echo "set output '| epstopdf --filter > all_seismograms_comparison.pdf'"
+
+#echo "set term x11"
+#echo "set term wxt"
+
+#echo set xrange \[0:2200\]
+
+foreach file ( SEMD/*.sem.ascii.sac.asciinew )
+
+ set newfile = `basename $file .sem.ascii.sac.asciinew`
+
+ echo plot \"QMXD/$newfile.qmxd.sac.asciinew\" w l lc 1, \"$file\" w l lc 3
+
+# uncomment this only when outputting to the screen (X11 or wxt)
+# echo "pause -1 'hit any key...'"
+
+end
+
diff --git a/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/plotall.gnu b/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/plotall.gnu
new file mode 100644
index 0000000..7c1f6f7
--- /dev/null
+++ b/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/plotall.gnu
@@ -0,0 +1,109 @@
+#set term postscript color solid Helvetica 22
+#set output "all_seismograms_comparison.ps"
+set term pdf color solid
+set output "all_seismograms_comparison.pdf"
+plot "QMXD/AIS.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/AIS.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/AIS.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/AIS.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/AIS.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/AIS.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ATD.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/ATD.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ATD.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/ATD.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ATD.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/ATD.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BBSR.IU.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/BBSR.IU.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BBSR.IU.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/BBSR.IU.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BBSR.IU.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/BBSR.IU.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BNG.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/BNG.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BNG.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/BNG.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/BNG.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/BNG.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CAN.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/CAN.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CAN.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/CAN.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CAN.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/CAN.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CRZF.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/CRZF.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CRZF.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/CRZF.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/CRZF.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/CRZF.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/DRV.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/DRV.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/DRV.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/DRV.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/DRV.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/DRV.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ECH.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/ECH.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ECH.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/ECH.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/ECH.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/ECH.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/EVO.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/EVO.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/EVO.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/EVO.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/EVO.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/EVO.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/FDF.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/FDF.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/FDF.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/FDF.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/FDF.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/FDF.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HDC.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/HDC.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HDC.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/HDC.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HDC.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/HDC.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HYB.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/HYB.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HYB.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/HYB.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/HYB.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/HYB.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/INU.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/INU.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/INU.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/INU.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/INU.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/INU.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KIP.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/KIP.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KIP.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/KIP.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KIP.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/KIP.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KOG.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/KOG.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KOG.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/KOG.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/KOG.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/KOG.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/NOUC.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/NOUC.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/NOUC.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/NOUC.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/NOUC.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/NOUC.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAF.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/PAF.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAF.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/PAF.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAF.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/PAF.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAS.CI.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/PAS.CI.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAS.CI.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/PAS.CI.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PAS.CI.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/PAS.CI.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PEL.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/PEL.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PEL.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/PEL.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PEL.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/PEL.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PPT.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/PPT.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PPT.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/PPT.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PPT.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/PPT.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PVC.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/PVC.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PVC.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/PVC.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/PVC.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/PVC.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RAO.IU.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/RAO.IU.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RAO.IU.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/RAO.IU.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RAO.IU.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/RAO.IU.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RER.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/RER.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RER.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/RER.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/RER.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/RER.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SAML.IU.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SAML.IU.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SAML.IU.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SAML.IU.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SAML.IU.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SAML.IU.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SCZ.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SCZ.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SCZ.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SCZ.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SCZ.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SCZ.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SEY.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SEY.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SEY.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SEY.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SEY.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SEY.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SFJD.IU.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SFJD.IU.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SFJD.IU.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SFJD.IU.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SFJD.IU.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SFJD.IU.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SPB.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SPB.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SPB.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SPB.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SPB.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SPB.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB2.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB2.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB2.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB2.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB2.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB2.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB3.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB3.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB3.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB3.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB3.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB3.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/SSB.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/SSB.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TAM.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/TAM.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TAM.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/TAM.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TAM.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/TAM.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TRIS.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/TRIS.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TRIS.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/TRIS.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/TRIS.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/TRIS.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/UNM.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/UNM.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/UNM.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/UNM.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/UNM.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/UNM.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/WUS.G.MXR.qmxd.sac.asciinew" w l lc 1, "SEMD/WUS.G.MXR.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/WUS.G.MXT.qmxd.sac.asciinew" w l lc 1, "SEMD/WUS.G.MXT.sem.ascii.sac.asciinew" w l lc 3
+plot "QMXD/WUS.G.MXZ.qmxd.sac.asciinew" w l lc 1, "SEMD/WUS.G.MXZ.sem.ascii.sac.asciinew" w l lc 3
diff --git a/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/process_all_seismograms.csh b/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/process_all_seismograms.csh
new file mode 100755
index 0000000..c513497
--- /dev/null
+++ b/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/process_all_seismograms.csh
@@ -0,0 +1,29 @@
+#!/bin/csh
+
+# you will find the Perl and "sac2asc" processing scripts below in directory utils/seis_process
+# and you will need the Seismic Analysis Code (SAC) software http://www.iris.edu/dms/nodes/dmc/software/downloads/sac
+# (i.e. the "sac" and "saclst" executable commands) to be installed and executable on your machine
+# for this processing script to run
+
+# Dimitri Komatitsch, CNRS, Marseille, France, September 2013
+
+########### SEM spectral-element results ###########
+
+./process_syn.pl -m CMTSOLUTION_to_filter_nonzero_hdur -a STATIONS -l 0/6000 -t 150/500 SEMD/*.sem.ascii
+
+./rotate.pl -l 0 -L 6000 SEMD/*.MXE.sem.ascii.sac
+
+foreach file ( SEMD/*.MXR.sem.ascii.sac SEMD/*.MXT.sem.ascii.sac SEMD/*.MXZ.sem.ascii.sac )
+ ./sac2asc $file > $file.asciinew
+end
+
+########### QMXD normal-mode results ###########
+
+./process_syn.pl -m CMTSOLUTION_to_filter_nonzero_hdur -a STATIONS -l 0/6000 -t 150/500 QMXD/*.qmxd
+
+./rotate.pl -l 0 -L 6000 QMXD/*.MXE.qmxd.sac
+
+foreach file ( QMXD/*.MXR.qmxd.sac QMXD/*.MXT.qmxd.sac QMXD/*.MXZ.qmxd.sac )
+ ./sac2asc $file > $file.asciinew
+end
+
diff --git a/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/processing b/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/processing
deleted file mode 100755
index ac426b1..0000000
--- a/benchmarks/prem_iso_noocean_Q_Cowling_Vanuatu/processing
+++ /dev/null
@@ -1,11 +0,0 @@
-./process_syn.pl -m CMTSOLUTION -a STATIONS -l 0/6000 -t 150/500 QMXD/*.qmxd
-./process_syn.pl -m CMTSOLUTION -a STATIONS -l 0/6000 -t 150/500 SEMD/*.semd.ascii
-./rotate.pl -l 0 -L 6000 QMXD/*.LHE.qmxd.sac
-./rotate.pl -l 0 -L 6000 SEMD/*.LHE.semd.ascii.sac
-mv QMXD/*LHZ*sac QMXD_SAC/
-mv QMXD/*LHT*sac QMXD_SAC/
-mv SEMD/*LHZ*sac SEMD_SAC/
-mv SEMD/*LHT*sac SEMD_SAC/
-tar cvf QMXD_SAC.tar QMXD_SAC/
-tar cvf SEMD_SAC.tar SEMD_SAC/
-bzip2 *.tar
More information about the CIG-COMMITS
mailing list