[cig-commits] r18111 - seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS

carltape at geodynamics.org carltape at geodynamics.org
Mon Mar 14 20:16:09 PDT 2011 

Author: carltape
Date: 2011-03-14 20:16:09 -0700 (Mon, 14 Mar 2011)
New Revision: 18111

Modified:
   seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/9818433_T006_T030_MPM_CI_m16_iker07_win_adj.pdf
   seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/make_pdf_by_event.pl
   seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/make_pdf_by_station.pl
   seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/plot_win_adj.pl
   seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/plot_win_adj_all.pl
   seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/plot_win_stats.pl
   seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/plot_win_stats_all.pl
Log:
ran default example listed in README


Modified: seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/9818433_T006_T030_MPM_CI_m16_iker07_win_adj.pdf
===================================================================
(Binary files differ)

Modified: seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/make_pdf_by_event.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/make_pdf_by_event.pl	2011-03-15 03:11:04 UTC (rev 18110)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/make_pdf_by_event.pl	2011-03-15 03:16:09 UTC (rev 18111)
@@ -1,9 +1,7 @@
 #!/usr/bin/perl -w
-
+#
 #==========================================================
 #  make_pdf_by_event.pl
-#  Carl Tape
-#  20-Oct-2009
 #
 #  This sorts the output PDF files in PLOTS into order by
 #    1  distance

Modified: seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/make_pdf_by_station.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/make_pdf_by_station.pl	2011-03-15 03:11:04 UTC (rev 18110)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/make_pdf_by_station.pl	2011-03-15 03:16:09 UTC (rev 18111)
@@ -1,10 +1,8 @@
 #!/usr/bin/perl -w
-
+#
 #==========================================================
 #
 #  make_pdf_by_station.pl
-#  Carl Tape
-#  01-Oct-2008
 #
 #  This sorts the output PDF files in PLOTS into order by
 #    1  distance

Modified: seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/plot_win_adj.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/plot_win_adj.pl	2011-03-15 03:11:04 UTC (rev 18110)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/plot_win_adj.pl	2011-03-15 03:16:09 UTC (rev 18111)
@@ -1,8 +1,6 @@
 #!/usr/bin/perl -w
-
+#
 #-----------------------------------
-# Min Chen, 09-July-2007
-# Carl Tape, 20-Oct-2009
 # plot_win_adj.pl
 #
 # This script plots a single figure for a single source-receiver pair
@@ -12,7 +10,7 @@
 #       (Search for "USER".)
 # 
 # EXAMPLE:
-#   plot_win_adj.pl -m ../CMTSOLUTION_9818433 -n MPM/CI/BH -b 0 -l -10/200 -k 7/1 -a STATIONS_ADJOINT -d DATA -s SYN -c RECON -w MEASUREMENT.WINDOWS -i m16 -j 6/30
+#   plot_win_adj.pl -M ../CMTSOLUTION_9818433 -n MPM/CI/BH -b 0 -l -10/200 -k 7/1 -a STATIONS_ADJOINT -d DATA -s SYN -c RECON -w MEASUREMENT.WINDOWS -i m16 -j 6/30
 #
 #-----------------------------------
 
@@ -22,10 +20,10 @@
 sub Usage{
   print STDERR <<END;
 
-  plot_win_adj.pl -m CMTFILE -n station/network -b iboth -r -a STATION_FILE -d data_dir -s syn_dir -c recon_dir -w winfile
+  plot_win_adj.pl -M CMTFILE -n station/network -b iboth -r -a STATION_FILE -d data_dir -s syn_dir -c recon_dir -w winfile
   with
 
-       -m -- CMTSOLUTION file to plot event focal mechanism (beach ball) on the map
+       -M -- CMTSOLUTION file to plot event focal mechanism (beach ball) on the map
        -l -- Start/End of trace to be cut
        -b -- plot MT and CC adjoint sources (=1)
        -k -- type of measurement (1-8) / plot adjoint source(0 or 1)
@@ -138,7 +136,7 @@
 $jpg_file = "$name.jpg";
 
 #$saclst="/opt/seismo-util/bin/saclst";
-$saclst="/opt/sac/bin/saclst";
+$saclst="/usr/local/sac/bin/saclst";
 
 #------------
 $xmin = -122; $xmax = -114;
@@ -161,14 +159,14 @@
    $fault_file = "$dir0/faults/jennings.xy";
    $kcf_file   = "$dir0/faults/kcf.xy";
    $breck_file = "$dir0/faults/breck.xy";
-   $fault_infoK = "-M -W0.5p,0/0/0";
+   $fault_infoK = "-m -W0.5p,0/0/0";
    `psxy $fault_file $bounds $proj $fault_infoK -O -K >> $ps_file`;
    `psxy $kcf_file $bounds $proj $fault_infoK -O -K >> $ps_file`;
    `psxy $breck_file $bounds $proj $fault_infoK -O -K >> $ps_file`;
 }
 
 # plot the stations
-`awk \'{print \$4,\$3,0,4,B,C}\' $opt_a | psxy $proj $bounds -M -St0.15 -W2 -G200/200/200 -N -O -K >>$ps_file`;
+`awk \'{print \$4,\$3,0,4,B,C}\' $opt_a | psxy $proj $bounds -m -St0.15 -W2 -G200/200/200 -N -O -K >>$ps_file`;
 
 #------------
 # data, synthetics, and (optional) reconstructed synthetics

Modified: seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/plot_win_adj_all.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/plot_win_adj_all.pl	2011-03-15 03:11:04 UTC (rev 18110)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/plot_win_adj_all.pl	2011-03-15 03:16:09 UTC (rev 18111)
@@ -1,8 +1,9 @@
 #!/usr/bin/perl -w
-
+#
 #----------------------------------------------------
-# Min Chen, 06-July-2007
-# plot_all_win_adj.pl
+# plot_win_adj_all.pl
+#
+# This calls plot_win_adj.pl to make plots of data, syn, windows, measurements, and adjoint sources.
 # 
 # EXAMPLE:
 #   plot_win_adj_all.pl -l -10/200 -m ../CMTSOLUTION_9818433 -n BH -b 0 -k 7/1 -a STATIONS_ADJOINT -d DATA -s SYN -c RECON -w MEASUREMENT.WINDOWS -i m16 -j 6/30

Modified: seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/plot_win_stats.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/plot_win_stats.pl	2011-03-15 03:11:04 UTC (rev 18110)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/plot_win_stats.pl	2011-03-15 03:16:09 UTC (rev 18111)
@@ -1,7 +1,6 @@
 #!/usr/bin/perl -w
-
+#
 #-----------------------------------
-# Carl Tape, 13-Nov-2007
 # plot_win_stats.pl
 #
 # This script plots a single figure for a single source-receiver pair
@@ -64,14 +63,14 @@
    $fault_file = "${dir_fault}/jennings_more.xy";
    #$kcf_file   = "${dir_fault}/kcf.xy";
    #$breck_file = "${dir_fault}/breck.xy";
-   $fault_infoK = "-M -W0.5p,0/0/0";
+   $fault_infoK = "-m -W0.5p,0/0/0";
    print CSH "psxy $fault_file $R $J $fault_infoK -O -K -V >> $psfile\n";
    #print CSH "psxy $kcf_file $R $J $fault_infoK -O -K -V >> $psfile\n";
    #print CSH "psxy $breck_file $R $J $fault_infoK -O -K -V >> $psfile\n";
 }
 
 # plot the stations
-print CSH "awk \'{print \$4,\$3,0,4,B,C}\' ${sta_file} | psxy $J $R -M -Si0.15 -W2 -G200/200/200 -N -O -K -V >>$psfile\n";
+print CSH "awk \'{print \$4,\$3,0,4,B,C}\' ${sta_file} | psxy $J $R -m -Si0.15 -W2 -G200/200/200 -N -O -K -V >>$psfile\n";
 
 # plot labels next to the map
 open(IN,"${eid_text}"); @mapstrings = <IN>; close(IN);

Modified: seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/plot_win_stats_all.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/plot_win_stats_all.pl	2011-03-15 03:11:04 UTC (rev 18110)
+++ seismo/3D/ADJOINT_TOMO/measure_adj/PLOTS/plot_win_stats_all.pl	2011-03-15 03:16:09 UTC (rev 18111)
@@ -1,7 +1,6 @@
 #!/usr/bin/perl -w
-
+#
 #-----------------------------------
-# Carl Tape, 06-April-2008
 # plot_win_stats_all.pl
 #
 # This script calls plot_win_stats.pl to plot histograms of measurements for each event.

More information about the CIG-COMMITS mailing list