[cig-commits] r18973 - in seismo/3D/ADJOINT_TOMO/flexwin/scripts: . prepare_scripts/socal
carltape at geodynamics.org
carltape at geodynamics.org
Mon Sep 26 14:38:15 PDT 2011
Author: carltape
Date: 2011-09-26 14:38:14 -0700 (Mon, 26 Sep 2011)
New Revision: 18973
Modified:
seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats_carl.sh
seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats_min.sh
seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows.pl
seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/prepare_input
Log:
updated socal scripts to reflect new channel labels for SPECFEM3D synthetics
Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats_carl.sh
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats_carl.sh 2011-09-26 21:36:50 UTC (rev 18972)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats_carl.sh 2011-09-26 21:38:14 UTC (rev 18973)
@@ -165,7 +165,7 @@
#echo $proj
#psbasemap -Rg -J$proj -Bg30/g30 -K -P -Y8 -V > $out
#pscoast -Rg -J$proj -Dl -A5000/0 -W1 -G$grey -O -K -V >> $out
-#psxy $t3 -Rg -J$proj -: -M -W2 -O -K -V >> $out
+#psxy $t3 -Rg -J$proj -: -m -W2 -O -K -V >> $out
# CHT
# set projection for map plotting
@@ -175,7 +175,7 @@
echo $proj
psbasemap $bounds $proj $tick -K -Y7.5 -X0.8 -P -V > $out
pscoast $bounds $proj -Df -A100 -W1p -Na/1.0p -S150/200/255 -G200/255/150 -O -K -V >> $out
-psxy $t3 $bounds $proj -: -M -W2 -O -K -V >> $out
+psxy $t3 $bounds $proj -: -m -W2 -O -K -V >> $out
psxy $proj $bounds -W1p -Sa0.2 -G250/0/0 -O -K >> $out <<EFX
$evlo $evla
EFX
@@ -255,13 +255,13 @@
proj=X3.0/-6.8
psbasemap -R$region -J$proj -B${timestep}:"Time(s) -- Z":/${degstep}:"Distance (deg)":WN -K -X1.0 -Y0.75 -U/0/-0.25/${basename} > $out2
-pswiggle $t2z -R$region -J$proj -Z20 -M -W1 -G$red -O -K >> $out2
+pswiggle $t2z -R$region -J$proj -Z20 -m -W1 -G$red -O -K >> $out2
psbasemap -R$region -J$proj -B${timestep}:"Time(s) -- R":/${degstep}:"Distance (deg)":wN -O -K -X3.0 >> $out2
-pswiggle $t2r -R$region -J$proj -Z20 -M -W1 -G$green -O -K >> $out2
+pswiggle $t2r -R$region -J$proj -Z20 -m -W1 -G$green -O -K >> $out2
psbasemap -R$region -J$proj -B${timestep}:"Time(s) -- T":/${degstep}:"Distance (deg)":wN -O -K -X3.0 >> $out2
-pswiggle $t2t -R$region -J$proj -Z20 -M -W1 -G$blue -O -K >> $out2
+pswiggle $t2t -R$region -J$proj -Z20 -m -W1 -G$blue -O -K >> $out2
pstext $rectxt -R$region -J$proj -N -O -G0 >> $out2 # FINISH (no -K)
@@ -273,13 +273,13 @@
#proj=X5/-2
#psbasemap -R$region -J$proj -B${timestep}:"Time(s)":/${degstep}:"Distance / degree":WS -K -P -Y8 > $out2
-#pswiggle $t2z -R$region -J$proj -Z20 -M -P -W1/$black -G$red -O -K >> $out2
+#pswiggle $t2z -R$region -J$proj -Z20 -m -P -W1/$black -G$red -O -K >> $out2
#psbasemap -R$region -J$proj -B:"Vertical component":/N -O -K >> $out2
#psbasemap -R$region -J$proj -B${timestep}:"Time(s)":/${degstep}:"Distance / degree":WS -O -K -Y-3.5 >> $out2
-#pswiggle $t2r -R$region -J$proj -Z20 -M -P -W1/$black -G$green -O -K >> $out2
+#pswiggle $t2r -R$region -J$proj -Z20 -m -P -W1/$black -G$green -O -K >> $out2
#psbasemap -R$region -J$proj -B:"Radial component":/N -O -K >> $out2
#psbasemap -R$region -J$proj -B${timestep}:"Time(s)":/${degstep}:"Distance / degree":WS -O -K -Y-3.5 >> $out2
-#pswiggle $t2t -R$region -J$proj -Z20 -M -P -W1/$black -G$blue -O -K >> $out2
+#pswiggle $t2t -R$region -J$proj -Z20 -m -P -W1/$black -G$blue -O -K >> $out2
#psbasemap -R$region -J$proj -B:"Transverse component":/N -O >> $out2
#echo "%!PS-Adobe-3.0" >t1
Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats_min.sh
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats_min.sh 2011-09-26 21:36:50 UTC (rev 18972)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats_min.sh 2011-09-26 21:38:14 UTC (rev 18973)
@@ -185,7 +185,7 @@
echo $proj
psbasemap $bounds $proj $tick -K -Y12 -X0.8 -P -V > $out
pscoast $bounds $proj -Dh -A100 -W0.5 -S150/200/255 -W2 -G200/255/150 -O -K -V >> $out
-psxy $t3 $bounds $proj -: -M -W2 -O -K -V >> $out
+psxy $t3 $bounds $proj -: -m -W2 -O -K -V >> $out
psxy $proj $bounds -W1 -Sa0.2 -G250/0/0 -O -K -: >> $out <<EFX
$evla $evlo
EFX
@@ -214,11 +214,11 @@
proj=X1.4/-2.5
psbasemap -R$region -J$proj -B${timestep}:"Time(s)":/${degstep}:"Distance / degree":WN -O -K -X-3.8 -Y-3.7 -U/0/-0.25/${basename} >> $out
-pswiggle $t2z -R$region -J$proj -Z20 -M -W1 -G$red -O -K >> $out
+pswiggle $t2z -R$region -J$proj -Z20 -m -W1 -G$red -O -K >> $out
psbasemap -R$region -J$proj -B${timestep}:"Time(s)":/${degstep}WN -O -K -X1.9>> $out
-pswiggle $t2r -R$region -J$proj -Z20 -M -W1 -G$green -O -K >> $out
+pswiggle $t2r -R$region -J$proj -Z20 -m -W1 -G$green -O -K >> $out
psbasemap -R$region -J$proj -B${timestep}:"Time(s)":/${degstep}WN -O -K -X1.9 >> $out
-pswiggle $t2t -R$region -J$proj -Z20 -M -W1 -G$blue -O >> $out
+pswiggle $t2t -R$region -J$proj -Z20 -m -W1 -G$blue -O >> $out
bbox=`grep "%%BoundingBox: " $out | tail -1 `
sed s/"$bbox"/"%"/ $out | sed s/"%%BoundingBox: (atend)"/"$bbox"/ > t1
@@ -237,13 +237,13 @@
proj=X3.0/-6.8 # negative sign to flip the y-axis
psbasemap -R$region -J$proj -B${timestep}:"Time(s) -- Z":/${degstep}:"Distance (deg)":WN -K -X1.0 -Y0.75 -U/0/-0.25/${basename} > $out
-pswiggle $t2z -R$region -J$proj -Z20 -M -W1 -G$red -O -K >> $out
+pswiggle $t2z -R$region -J$proj -Z20 -m -W1 -G$red -O -K >> $out
psbasemap -R$region -J$proj -B${timestep}:"Time(s) -- R":/${degstep}:"Distance (deg)":wN -O -K -X3.0 >> $out
-pswiggle $t2r -R$region -J$proj -Z20 -M -W1 -G$green -O -K >> $out
+pswiggle $t2r -R$region -J$proj -Z20 -m -W1 -G$green -O -K >> $out
psbasemap -R$region -J$proj -B${timestep}:"Time(s) -- T":/${degstep}:"Distance (deg)":wN -O -K -X3.0 >> $out
-pswiggle $t2t -R$region -J$proj -Z20 -M -W1 -G$blue -O -K >> $out
+pswiggle $t2t -R$region -J$proj -Z20 -m -W1 -G$blue -O -K >> $out
pstext $rectxt -R$region -J$proj -N -O -G0 >> $out # FINISH (no -K)
Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows.pl 2011-09-26 21:36:50 UTC (rev 18972)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows.pl 2011-09-26 21:38:14 UTC (rev 18973)
@@ -74,11 +74,17 @@
if($ibody==1) {$bbody = ".true."} else {$bbody = ".false."}
if($ibp==1) {$bbp = ".true."} else {$bbp = ".false."}
+#-------------------------------------
+# USER INPUT -- CHANGE FOR EACH USER
+
# directories for data and synthetics
# 1=socal; 2=Japan
if($idataset == 1) {
- $dir_data = "/home/carltape/SOCAL_ADJOINT/DATA/FINAL";
- $dir_syn = "/home/carltape/SOCAL_ADJOINT/SYN/model_${smodel}";
+ #$dir_data = "/home/carltape/SOCAL_ADJOINT/DATA/FINAL";
+ #$dir_syn = "/home/carltape/SOCAL_ADJOINT/SYN/model_${smodel}";
+ #$dir_data = "/data2/data/TEST";
+ $dir_data = "/data2/data/calif/FINAL";
+ $dir_syn = "/data2/syn/socal/model_${smodel}";
$sdataset = "socal";
} elsif ($idataset == 2) {
$dir_data = "/home/mchen/DATA/TEST";
@@ -88,30 +94,38 @@
die("\n idataset must be 1 (socal) or 2 (japan)\n");
}
-#-------------------------------------
-
# specify various directories (MUST BE MODIFIED FOR EACH USER)
-$dir0 = "/data1/cig/seismo/3D"; # SVN path
+# note: may want a second copy to run simultaneously
+$dir0 = "/data2/SVN/seismo/3D"; # path to CIG SVN repository
$dir_win_code = "$dir0/ADJOINT_TOMO/flexwin_work";
$dir_win_run = "$dir0/flexwin_run";
#$dir_win_code = "$dir0/ADJOINT_TOMO/flexwin_work_copy2";
#$dir_win_run = "$dir0/flexwin_run_copy2";
$dir_scripts = "${dir_win_code}/scripts";
-# scripts for preparing data and synthetics
-$dir_prepare = "${dir_scripts}/prepare_scripts/${sdataset}";
-
-# scripts for user files
-$dir_user = "${dir_win_code}/user_files";
-
# directory to collect COPIES of various output files
$odir = "/home/carltape/results/WINDOWS/model_${smodel}";
# run directory for windows, measurements, adjoint sources, and kernels
-$rundir = "/home/carltape/SOCAL_ADJOINT/RUNS";
+#$rundir = "/home/carltape/SOCAL_ADJOINT/RUNS";
+$rundir = "/data2/RUNS";
+# channels to use for data and synthetics
+# NOTE: update prepare_input also ($pinput)
+#@chans = ("BHZ","BHT","BHR","HHZ","HHR","HHT");
+ at dchans = ("BHZ","BHT","BHR","HHZ","HHR","HHT");
+#@schans = ("BHZ","BHT","BHR","BHZ","BHR","BHT"); # original
+ at schans = ("HXZ","HXT","HXR","HXZ","HXR","HXT");
+$nchan = @dchans;
+
#-------------------------------------
+# scripts for preparing data and synthetics
+$dir_prepare = "${dir_scripts}/prepare_scripts/${sdataset}";
+
+# scripts for user files
+$dir_user = "${dir_win_code}/user_files";
+
# check that directories exist
if (not -e ${dir_data}) {die("check if ${dir_data} exist or not\n")}
if (not -e ${dir_syn}) {die("check if ${dir_syn} exist or not\n")}
@@ -133,10 +147,6 @@
$suffix_dat = "bp12_150s";
}
-# channels to use for data and synthetics
- at chans = ("BHZ","BHT","BHR","HHZ","HHR","HHT");
-$nchan = @chans;
-
# directories for windowing code
$dir_win_run_syn = "${dir_win_run}/SYN";
$dir_win_run_data = "${dir_win_run}/DATA";
@@ -328,9 +338,10 @@
# copy data and synthetics into windowing directory (BHZ,BHR,BHT)
print CSH "echo copying data and syn files to ${dir_win_run}\n";
for ($j = 1; $j <= $nchan; $j++) {
- $chan = $chans[$j-1];
- print CSH "cp $syndir1/*${chan}.${suffix_syn} ${dir_win_run_syn}\n"; # synthetics
- print CSH "cp $datadir1/*${chan}.${suffix_dat} ${dir_win_run_data}\n"; # data
+ $dchan = $dchans[$j-1];
+ $schan = $schans[$j-1];
+ print CSH "cp $syndir1/*${schan}.${suffix_syn} ${dir_win_run_syn}\n"; # synthetics
+ print CSH "cp $datadir1/*${dchan}.${suffix_dat} ${dir_win_run_data}\n"; # data
}
# print CSH "cp $syndir1/$sfile1 ${dir_win_run_syn}\n";
@@ -343,7 +354,7 @@
# prepare sac files for windowing code (BHZ,BHR,BHT)
# NOTE: the file extensions must be checked within prepare_seis.pl and prepare_input
print CSH "cd ${dir_win_run}\n";
- #print CSH "${pseis}\n"; # this is done in PRE-PROCESSING
+ #print CSH "${pseis}\n"; # now this is done in PRE-PROCESSING
print CSH "${pinput}\n";
# copy prepared data and synthetic files into a directory for measurement code
Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/prepare_input
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/prepare_input 2011-09-26 21:36:50 UTC (rev 18972)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/prepare_input 2011-09-26 21:38:14 UTC (rev 18973)
@@ -15,16 +15,16 @@
comp=`echo $chan | awk '{print substr($1,3,1)}'`
#echo $net $sta $chan $comp
- # synthetics are always BH_ (for socal problem)
- #synt=SYN/${sta}.${net}.BH${comp}.semd.sac.d
- #synt=`ls SYN/${sta}.${net}.BH${comp}*`
- nsynt=`ls -1 SYN/${sta}.${net}.BH${comp}* | wc | awk '{print $1}'`
+ # synthetics are always HX_ (for socal problem)
+ #synt=SYN/${sta}.${net}.HX${comp}.semd.sac.d
+ #synt=`ls SYN/${sta}.${net}.HX${comp}*`
+ nsynt=`ls -1 SYN/${sta}.${net}.HX${comp}* | wc | awk '{print $1}'`
out=MEASURE/${sta}.${net}.${chan}
if [ $nsynt -eq 1 ] ; then # data with unique syn match
- synt=`ls SYN/${sta}.${net}.BH${comp}*`
+ synt=`ls SYN/${sta}.${net}.HX${comp}*`
b1=`saclst b f $file | awk '{print $2}'`
b2=`saclst b f $synt | awk '{print $2}'`
db=$(echo "sqrt(($b1 - $b2)^2)" | bc)
@@ -45,7 +45,7 @@
#fi
else
- echo "SYN/${sta}.${net}.BH${comp}*" >> t3
+ echo "SYN/${sta}.${net}.HX${comp}*" >> t3
fi
done
More information about the CIG-COMMITS
mailing list