[cig-commits] r14513 - seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt
carltape at geodynamics.org
carltape at geodynamics.org
Fri Mar 27 21:58:05 PDT 2009
Author: carltape
Date: 2009-03-27 21:58:04 -0700 (Fri, 27 Mar 2009)
New Revision: 14513
Added:
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/INPUT
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/SAVE
Modified:
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_horz_coverage.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_horz_models.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_vert_models.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_vert_models_basemap.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_xc_seis_all.pl
seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/xc_and_seis.tex
Log:
Modified plotting scripts for tomographic model cross sections.
Added: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/INPUT
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/INPUT (rev 0)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/INPUT 2009-03-28 04:58:04 UTC (rev 14513)
@@ -0,0 +1 @@
+link /home/carltape/ADJOINT_TOMO/ADJOINT_TOMO_OUTPUT/model_plot_matlab_OUTPUT
\ No newline at end of file
Property changes on: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/INPUT
___________________________________________________________________
Name: svn:special
+ *
Added: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/SAVE
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/SAVE (rev 0)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/SAVE 2009-03-28 04:58:04 UTC (rev 14513)
@@ -0,0 +1 @@
+link /home/carltape/ADJOINT_TOMO/ADJOINT_TOMO_OUTPUT/model_plot_gmt_OUTPUT
\ No newline at end of file
Property changes on: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/SAVE
___________________________________________________________________
Name: svn:special
+ *
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_horz_coverage.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_horz_coverage.pl 2009-03-28 03:36:12 UTC (rev 14512)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_horz_coverage.pl 2009-03-28 04:58:04 UTC (rev 14513)
@@ -34,7 +34,7 @@
# KEY COMMAND: HOW DID YOU WEIGHT THE SUMMED KERNELS
$smodeltag = "coverage_sum_abs_nwin";
-$smodeltag = "coverage_sum_abs";
+#$smodeltag = "coverage_sum_abs";
# detailed flags for publication figures
#$ilonlab = 0; # display longitude tick labels at the top
@@ -53,7 +53,7 @@
# USER INPUT
# directory containing the data files
-$dirdat = "/net/sierra/raid1/carltape/results/PLOTTING/SAVE/horz_full";
+$dirdat = "INPUT/horz_01";
# KEY: file showing the cuts and the plotting range for the values
$fcuts = "${dirdat}/horz_xc_elevation_cuts_mod";
@@ -120,7 +120,6 @@
$Wshelf = "-W0.75p,0/0/0,--";
# CMT sources
-$dir_tomo = "/net/denali/scratch1/carltape/svn/cig/seismo/3D/ADJOINT_TOMO";
$dir_sources = "/net/sierra/raid1/carltape/results/SOURCES/socal_16";
$outer_boundary = "${dir_sources}/SPECFEM_outer_points.dat";
$inner_boundary = "${dir_sources}/SPECFEM_inner_points.dat";
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_horz_models.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_horz_models.pl 2009-03-28 03:36:12 UTC (rev 14512)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_horz_models.pl 2009-03-28 04:58:04 UTC (rev 14513)
@@ -11,9 +11,9 @@
# and plots the three cross-sections.
#
# EXAMPLE:
-# plot_horz_models.pl 1/50 m00/m16 1 0.2 # Vs
-# plot_horz_models.pl 1/50 m00/m16 0 0.3 # Vb
-# plot_horz_models.pl 1/6 m00/m01 1 0.2 # Vs
+# plot_horz_models.pl 1/50 m00/m16 1 1/1/1/0.2 # Vs
+# plot_horz_models.pl 1/50 m00/m16 0 1/1/1/0.3 # Vb
+# plot_horz_models.pl 1/6 m00/m01 1 1/1/1/0.2 # Vs
#
# SCIENCE PAPER:
# plot_horz_models.pl 3/3 m00/m16 1 1/0/1/0.2 # Vs
@@ -48,7 +48,7 @@
# USER INPUT
# directory containing the data files
-$dirdat = "/net/sierra/raid1/carltape/results/PLOTTING/SAVE/horz_full";
+$dirdat = "INPUT/horz_01";
# KEY: file showing the cuts and the plotting range for the values
$fcuts = "${dirdat}/horz_xc_elevation_cuts_mod";
@@ -115,7 +115,6 @@
$Wshelf = "-W0.75p,0/0/0,--";
# CMT sources
-$dir_tomo = "/net/denali/scratch1/carltape/svn/cig/seismo/3D/ADJOINT_TOMO";
$dir_sources = "/net/sierra/raid1/carltape/results/SOURCES/socal_16";
$outer_boundary = "${dir_sources}/SPECFEM_outer_points.dat";
$inner_boundary = "${dir_sources}/SPECFEM_inner_points.dat";
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_vert_models.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_vert_models.pl 2009-03-28 03:36:12 UTC (rev 14512)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_vert_models.pl 2009-03-28 04:58:04 UTC (rev 14513)
@@ -43,9 +43,8 @@
# directory containing the data files
$stirun = sprintf("%2.2i",$irun);
-$dirdat = "/net/sierra/raid1/carltape/results/PLOTTING/SAVE/vert_${stirun}";
+$dirdat = "INPUT/vert_${stirun}";
-
#---------------------------------------------------------
# subtitles
@@ -82,7 +81,6 @@
$Wshelf = "-W0.75p,0/0/0,--";
# CMT sources
-$dir_tomo = "/net/denali/scratch1/carltape/svn/cig/seismo/3D/ADJOINT_TOMO";
$dir_sources = "/net/sierra/raid1/carltape/results/SOURCES/socal_16";
$outer_boundary = "${dir_sources}/SPECFEM_outer_points.dat";
$inner_boundary = "${dir_sources}/SPECFEM_inner_points.dat";
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_vert_models_basemap.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_vert_models_basemap.pl 2009-03-28 03:36:12 UTC (rev 14512)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_vert_models_basemap.pl 2009-03-28 04:58:04 UTC (rev 14513)
@@ -76,7 +76,6 @@
$Wshelf = "-W0.75p,0/0/0,--";
# CMT sources
-$dir_tomo = "/net/denali/scratch1/carltape/svn/cig/seismo/3D/ADJOINT_TOMO";
$dir_sources = "/net/sierra/raid1/carltape/results/SOURCES/socal_16";
$outer_boundary = "${dir_sources}/SPECFEM_outer_points.dat";
$inner_boundary = "${dir_sources}/SPECFEM_inner_points.dat";
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_xc_seis_all.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_xc_seis_all.pl 2009-03-28 03:36:12 UTC (rev 14512)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/plot_xc_seis_all.pl 2009-03-28 04:58:04 UTC (rev 14513)
@@ -4,7 +4,7 @@
#
# plot_xc_seis_all.pl
# Carl Tape
-# 07-Dec-2008
+# 27-March-2009
#
# EXAMPLES:
# plot_xc_seis_all.pl T006_T030 1 1
@@ -21,7 +21,7 @@
# make output directory
$otag = "${Ttag}_win${iwin}_1D${i1D}";
-$odir = "./$otag";
+$odir = "SAVE/$otag";
#`rm -rf $odir`;
`mkdir -p $odir`;
@@ -29,9 +29,11 @@
$irun = 1;
$stirun = sprintf("%2.2i",$irun);
-$dirdat = "/net/sierra/raid1/carltape/results/PLOTTING/SAVE/vert_${stirun}";
-$seisdir = "/net/denali/raid1/carltape/svn/cig/seismo/3D/ADJOINT_TOMO/iterate_adj/misfit_plot_work/${otag}";
-$xcdir = "/net/sierra/raid1/carltape/results/PLOTTING/gmt";
+$dir0 = "/home/carltape/ADJOINT_TOMO";
+$dirin = "${dir0}/ADJOINT_TOMO_OUTPUT";
+$dirdat = "${dirin}/model_plot_matlab_OUTPUT/vert_${stirun}";
+$seisdir = "${dirin}/misfit_plot_OUTPUT/${otag}";
+$xcdir = "${dirin}/model_plot_gmt_OUTPUT/VERT_${stirun}";
if (not -e $dirdat) {die("Check if dirdat $dirdat exist or not\n");}
if (not -e $xcdir) {die("Check if xcdir $xcdir exist or not\n");}
if (not -e $seisdir) {die("Check if seisdir $seisdir exist or not\n");}
@@ -52,7 +54,7 @@
$pmin = 1; $pmax = $nump;
#$pmin = 1; $pmax = 10;
-#$pmin = 77; $pmax = $pmin;
+$pmin = 77; $pmax = $pmin;
for ($p = $pmin; $p <= $pmax; $p ++ ) {
$stp = sprintf("%3.3i",$p);
@@ -131,8 +133,8 @@
}
-# move all pdf and ps files to output directory
-`mv xc_and_seis*pdf xc_and_seis*ps $odir`;
-`/home/carltape/bin/pdcat -r $odir/xc_and_seis*pdf $odir/ALL_xc_and_seis_${otag}.pdf`;
+# UNCOMMENT: move all pdf and ps files to output directory
+#`mv xc_and_seis*pdf xc_and_seis*ps $odir`;
+#`/home/carltape/bin/pdcat -r $odir/xc_and_seis*pdf $odir/ALL_xc_and_seis_${otag}.pdf`;
#==================================================
Modified: seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/xc_and_seis.tex
===================================================================
--- seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/xc_and_seis.tex 2009-03-28 03:36:12 UTC (rev 14512)
+++ seismo/3D/ADJOINT_TOMO/iterate_adj/model_plot/gmt/xc_and_seis.tex 2009-03-28 04:58:04 UTC (rev 14513)
@@ -6,9 +6,9 @@
\fcolorbox{white}{white}{
\hspace{-1cm}
\begin{tabular}{c}
-\includegraphics[width=4.5cm,angle=-90]{/net/sierra/raid1/carltape/results/PLOTTING/gmt/vert_01_xc_DSC_CI_214_14155260_vs_m16_m00_086.eps}
+\includegraphics[width=4.5cm,angle=-90]{/home/carltape/ADJOINT_TOMO/ADJOINT_TOMO_OUTPUT/model_plot_gmt_OUTPUT/VERT_01/vert_01_xc_LGU_CI_206_12659440_vs_m16_m00_077.eps}
\\
-\includegraphics[width=3.8cm,angle=-90]{/net/denali/raid1/carltape/svn/cig/seismo/3D/ADJOINT_TOMO/iterate_adj/misfit_plot_work/T002_T030_win0_1D0/14155260_T002_T030_DSC_CI_m16_seis_2.eps}
+\includegraphics[width=3.8cm,angle=-90]{/home/carltape/ADJOINT_TOMO/ADJOINT_TOMO_OUTPUT/misfit_plot_OUTPUT/T006_T030_win1_1D1/12659440_T006_T030_LGU_CI_m16_seis_3.eps}
\end{tabular}
}
\end{slide}
More information about the CIG-COMMITS
mailing list