[cig-commits] r13084 - in seismo/3D/ADJOINT_TOMO/flexwin: scripts scripts/prepare_scripts/socal user_files/socal_3D
carltape at geodynamics.org
carltape at geodynamics.org
Wed Oct 15 11:06:06 PDT 2008
Author: carltape
Date: 2008-10-15 11:06:06 -0700 (Wed, 15 Oct 2008)
New Revision: 13084
Added:
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m07
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m08
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m09
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m10
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_Daniel
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m00
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m00_orig
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m01
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m02
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m03
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m04
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m05
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m06
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T003_T030_m07
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T003_T030_m08
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T003_T030_m09
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T003_T030_m10
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T004_T040_m07
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m07
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m08
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m09
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m10
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_Daniel
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m00
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m01
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m02
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m03
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m04
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m05
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m06
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m07
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m08.f90
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m09.f90
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m10.f90
Removed:
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_Daniel
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m00
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m00_orig
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m01
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m02
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m03
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m04
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m05
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m06
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m07
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T004_m07
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_14095540
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_Daniel
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m00
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m01
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m02
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m03
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m04
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m05
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m06
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m07
Modified:
seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats_carl.sh
seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows.pl
seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows_local.pl
seismo/3D/ADJOINT_TOMO/flexwin/scripts/plot_seismos_gmt.sh
seismo/3D/ADJOINT_TOMO/flexwin/scripts/plot_windows_all.pl
seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/prepare_input_test
seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/process_data_and_syn.pl
seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m07.f90
Log:
Modified socal user file names. Added some variables to plot_seismos_gmt.sh for flexibility.
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 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats_carl.sh 2008-10-15 18:06:06 UTC (rev 13084)
@@ -193,7 +193,7 @@
proj=X1.75/1.5
## plot window number statistics
-awk '{print $6}' $t0 | pshistogram -J$proj -C -B:"CC":/WS -W0.02 -Lthin -G$grey -O -K -X-4.5 -Y-2. >> $out
+awk '{print $6}' $t0 | pshistogram -J$proj -C -B:"CC":/WS -W0.02 -Lthin -G$grey -O -K -X-4.5 -Y-2.0 >> $out
awk '{print $5}' $t0 | pshistogram -J$proj -C -B:"Tshift":/WS -W1 -Lthin -G$grey -O -K -X2.75 >> $out
awk '{print $7}' $t0 | pshistogram -J$proj -C -B:"dlnA":/WS -W0.1 -Lthin -G$grey -O -X2.75 -K >> $out
@@ -227,6 +227,7 @@
psbasemap -R$region -J$proj -B -O -K >> $out
## plot window number statistics
+## NOTE: IF THERE IS A CONSTANT NUMBER OF WINDOWS FOR A COMPONENT, THEN PSHISTOGRAM CRASHES
grep Z $t1 | awk '{print $2}' | pshistogram -J$proj -C -B:"Windows per Z trace":/WS -W1 -Lthin -G$red -O -K -X-5.5 -Y-2.25 >> $out
grep R $t1 | awk '{print $2}' | pshistogram -J$proj -C -B:"Windows per R trace":/WS -W1 -Lthin -G$green -O -K -X2.75 >> $out
grep T $t1 | awk '{print $2}' | pshistogram -J$proj -C -B:"Windows per T trace":/WS -W1 -Lthin -G$blue -O -X2.75 >> $out # FINISH (no -K)
@@ -310,3 +311,4 @@
rm $t0 $t1 $t2z $t2r $t2t $t3
+#===============================================
Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows.pl 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows.pl 2008-10-15 18:06:06 UTC (rev 13084)
@@ -38,12 +38,12 @@
# prepare_meas_all.pl
#
# EXAMPLES:
-# pick_all_windows.pl m07 0 6/40 1/204 1/1/1/0 1 0 # make plots and WINDOWS file, T = 6-40s
-# pick_all_windows.pl m07 0 4/40 1/204 1/1/1/1 1 0 #
-# pick_all_windows.pl m07 0 2/40 1/204 1/1/1/1 1 0 # make plots and WINDOWS file, T = 2-40s
+# pick_all_windows.pl m10 0 6/30 1/204 1/1/1/0 1 0 # make plots and WINDOWS file, T = 6-40s
+# pick_all_windows.pl m10 0 3/30 1/204 1/1/1/1 1 0 #
+# pick_all_windows.pl m10 0 2/30 1/204 1/1/1/1 1 0 # make plots and WINDOWS file, T = 2-40s
#
-# pick_all_windows.pl m00 0 6/40 179/179 1/1/0/0 1 0 # make plots only, T = 6-40s
-# pick_all_windows.pl m00 0 6/40 179/179 0/0/1/0 1 0 # make WINDOWS file, T = 6-40s
+# pick_all_windows.pl m00 0 6/30 179/179 1/1/0/0 1 0 # make plots only, T = 6-40s
+# pick_all_windows.pl m00 0 6/30 179/179 0/0/1/0 1 0 # make WINDOWS file, T = 6-40s
#
#==========================================================
@@ -143,7 +143,7 @@
$dir_win_run_meas = "${dir_win_run}/MEASURE";
# NEW: EVENT LIST
-$eid_list = "/net/sierra/raid1/carltape/results/EID_LISTS/syn_run_m07";
+$eid_list = "/net/sierra/raid1/carltape/results/EID_LISTS/syn_run_${smodel}";
if (not -f $eid_list) {die("check if eid_list ${eid_list} exist or not\n")}
open(IN,$eid_list); @eids = <IN>; close(IN);
$nevent0 = @eids;
@@ -180,11 +180,11 @@
# copy PAR_FILE into local directory
# directories for data and synthetics
if($idataset == 1) {
- $par_file_in = "${dir_user}/socal_3D/PAR_FILE_${sTmin}_${smodel}";
+ $par_file_in = "${dir_user}/socal_3D/PAR_FILE_${Ttag}_${smodel}";
$userfun_file_in = "${dir_user}/socal_3D/user_functions_${smodel}.f90";
} elsif ($idataset == 2) {
- $par_file_in = "${dir_user}/japan_3D/PAR_FILE_${sTmin}";
+ $par_file_in = "${dir_user}/japan_3D/PAR_FILE_${Ttag}";
$userfun_file_in = "${dir_user}/japan_3D/user_functions.f90";
}
if (not -f ${par_file_in}) {die("check if PAR_FILE ${par_file_in} exist or not\n")}
@@ -265,7 +265,7 @@
# copy PAR_FILE and user_functions.f90 to output directory
print CSH "\\cp ${par_file} ${userfun_file} ${dir_prepare}/${pseis} $odir\n";
-#die("testing");
+#die("TESTING");
#===========================================================================
@@ -305,8 +305,8 @@
print CSH "echo Event ID is $eid -- $kk out of $nevent\n";
- $eout = "${eid}_${sTmin}"; # output file tag
- $measdir = "$rundir/$eid/$smodel/WINDOW_${sTmin}"; # main run directory
+ $eout = "${eid}_${Ttag}_${smodel}"; # output file tag
+ $measdir = "$rundir/$eid/$smodel/WINDOW_${Ttag}"; # main run directory
# If the windowing code has already been run, then move to the next event.
# NOTE: to overwrite, just comment out this option
@@ -321,8 +321,10 @@
# remove folders in the windowing directory
print CSH "\\rm -rf ${dir_win_run_syn} ${dir_win_run_data} ${dir_win_run_meas} \n";
print CSH "mkdir ${dir_win_run_syn} ${dir_win_run_data} ${dir_win_run_meas} \n";
+
# 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
@@ -344,7 +346,7 @@
# copy prepared data and synthetic files into a directory for measurement code
if ($imeas == 1) {
- #$measdir = "$rundir/$eid/$smodel/WINDOW_${sTmin}";
+ #$measdir = "$rundir/$eid/$smodel/WINDOW_${Ttag}";
print CSH "echo copying prepared DATA and SYN to ${measdir}...\n";
print CSH "mkdir -p $rundir/$eid\n";
print CSH "mkdir -p $rundir/$eid/$smodel/\n";
Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows_local.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows_local.pl 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/pick_all_windows_local.pl 2008-10-15 18:06:06 UTC (rev 13084)
@@ -4,7 +4,7 @@
#
# pick_all_windows_local.pl
# Carl Tape
-# 22-August-2007
+# 02-Oct-2008
#
# This is the "local version" of pick_all_windows_local.pl
# This is run from the directory containing MEASURE, DATA, SYN,
@@ -59,8 +59,8 @@
if($ibody==1) {$bbody = ".true."} else {$bbody = ".false."}
# string for PDF
-$sTmin = sprintf("%2.2i",$Tmin);
-$sTmax = sprintf("%2.2i",$Tmax);
+$sTmin = sprintf("%3.3i",$Tmin);
+$sTmax = sprintf("%3.3i",$Tmax);
$sTminD = sprintf("%.2f",$Tmin);
$sTmaxD = sprintf("%.2f",$Tmax);
@@ -71,7 +71,8 @@
$win_execute = "flexwin";
# make command for windowing code
-$make = "make -f make_intel_caltech";
+$make = "make -f make_gfortran"; # change as of 9-30-08
+#$make = "make -f make_intel_caltech";
# location of plot_windows_all.pl
$plot_windows_perl = "${dir_scripts}/plot_windows_all.pl";
@@ -146,7 +147,7 @@
#-----------------------------------------
close(CSH);
-system("csh -f $cshfile");
+#system("csh -f $cshfile");
print "\n";
Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/plot_seismos_gmt.sh
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/plot_seismos_gmt.sh 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/plot_seismos_gmt.sh 2008-10-15 18:06:06 UTC (rev 13084)
@@ -19,15 +19,28 @@
t_step=`echo $t_end $t_start | awk '{print int(($1-$2)/10)}'`
# CHT
-#t_start=-20
+t_start=-20
#t_end=120
-#t_end=200
+t_end=200
# find number of windows
nwin=`grep NUM_WIN ${win} | awk '{print $NF}'`
-gmtset PAPER_MEDIA a4+ MEASURE_UNIT inch HEADER_FONT_SIZE 14p LABEL_FONT_SIZE 16p
+# default settings
+paper="a4+"
+orient="-P"
proj=X8/1.5
+origin="-Y6"
+shift="-Y-2.5"
+
+# carl settings
+#paper="letter"
+#orient=" "
+#proj=X10/1.8
+#origin="-X0.5 -Y6"
+#shift="-Y-2.5"
+
+gmtset PAPER_MEDIA $paper MEASURE_UNIT inch HEADER_FONT_SIZE 14p LABEL_FONT_SIZE 16p
nhead=4
red=255/0/0
black=0/0/0
@@ -40,7 +53,8 @@
#################################
# set output filename
-out=${basename}.seis.eps
+out=${basename}.seis.ps
+outpdf=${basename}.seis.pdf
## for SCSN plots:
#dkm=`grep DKM $info | awk '{print $NF}'`
@@ -53,7 +67,8 @@
region=$t_start/$t_end/-$max/$max
# do plot
-psbasemap -R$region -J$proj -B${t_step}::.${label}:"Time (s)":/S -K -P -Y8 > $out
+
+psbasemap -R$region -J$proj -B${t_step}::.${label}:"Time (s)":/S -K $orient $origin > $out
tail -$nwin $win | awk '{printf "%f %f\n%f %f\n%f %f\n%f %f\n>\n", $2,-1*m,$2,m,$3,m,$3,-1*m}' m=$max | psxy -R$region -J$proj -M -W1 -G$paleblue -O -K >> $out
psxy $obs -R$region -J$proj -H$nhead -W2 -O -K >> $out
psxy $syn -R$region -J$proj -H$nhead -W2/$red -O -K >> $out
@@ -86,11 +101,10 @@
nhead=8
# set region
-proj=X8/1.5
max=`grep STA_LTA_MAX ${stalta} | awk '{print $NF}'`
region=$t_start/$t_end/0/$max
-psbasemap -R$region -J$proj -B${t_step}:"Time (s)":/S -O -K -Y-2.5 >> $out
+psbasemap -R$region -J$proj -B${t_step}:"Time (s)":/S -O -K $shift >> $out
tail -$nwin $win | awk '{printf "%f %f\n%f %f\n%f %f\n%f %f\n>\n", $2,-1*m,$2,m,$3,m,$3,-1*m}' m=$max | psxy -R$region -J$proj -M -W1 -G$paleblue -O -K >> $out
awk '{print $1, $2}' $stalta | psxy -R$region -J$proj -H$nhead -W2/$blue -O -K >> $out
@@ -112,10 +126,9 @@
# set region
max=`grep PLOT_MAX ${envobs} | awk '{print $NF}'`
region=$t_start/$t_end/0/$max
-proj=X8/1.5
# do plot
-psbasemap -R$region -J$proj -B${t_step}:"Time (s)":/S -O -K -Y-2.5 >> $out
+psbasemap -R$region -J$proj -B${t_step}:"Time (s)":/S -O -K $shift >> $out
tail -$nwin $win | awk '{printf "%f %f\n%f %f\n%f %f\n%f %f\n>\n", $2,-1*m,$2,m,$3,m,$3,-1*m}' m=$max | psxy -R$region -J$proj -M -W1 -G$paleblue -O -K >> $out
psxy $envobs -R$region -J$proj -H$nhead -W2 -O -K >> $out
psxy $envsyn -R$region -J$proj -H$nhead -W2/$red -O -K >> $out
@@ -131,6 +144,8 @@
#epstopdf $out
#rm t1.$$
+ps2pdf $out $outpdf
+
#echo $out
#gv $out
Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/plot_windows_all.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/plot_windows_all.pl 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/plot_windows_all.pl 2008-10-15 18:06:06 UTC (rev 13084)
@@ -67,10 +67,11 @@
$filename = "$namerec.$namenet.$comps[$icomp-1]";
$obsfile = "${dir_win_run_meas}/${filename}.obs";
- # generate output PDF files -- make sure that plot_seismos_gmt.sh is doing this
+ # generate output PDF files -- make sure that plot_seismos_gmt.sh is working
if (-f $obsfile) {
print "$obsfile\n";
$ftag = "${dir_win_run_meas}/${filename}";
+ print "${dir_scripts}/plot_seismos_gmt.sh $ftag\n";
`${dir_scripts}/plot_seismos_gmt.sh $ftag`;
}
}
Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/prepare_input_test
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/prepare_input_test 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/prepare_input_test 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,32 +1,35 @@
#!/bin/sh
-if [ -e t1 ] ; then ( rm t1 ) ; fi
-if [ -e t2 ] ; then ( rm t2 ) ; fi
-if [ -e t3 ] ; then ( rm t3 ) ; fi
+# update by Qinya --
+# theoretically t2 should not exist after process_data_and_syn.pl
+rm -f t1 t2 t3
n=0
-for file in DATA_TEST/* ; do
+# only processing ?H[RTZ] components for the windowing code
+for file in DATA_TEST/*.?H[RTZ].* ; do
+ # take network, station, channel and comp info from the sac headers
net=`saclst knetwk f $file | awk '{print $2}'`
sta=`saclst kstnm f $file | awk '{print $2}'`
chan=`saclst kcmpnm f $file | awk '{print $2}'`
- #net=`echo $file | awk -F"." '{print $2}'`
- #sta=`echo $file | awk -F"." '{print $3}'`
- #chan=`echo $file | awk -F"." '{print $4}'`
comp=`echo $chan | awk '{print substr($1,3,1)}'`
#echo $net $sta $chan $comp
# synthetics are always BH_
- synt=SYN_TEST/${sta}.${net}.BH${comp}.semd.sac.d
+ #synt=SYN_TEST/${sta}.${net}.BH${comp}.semd.sac.d
+ #synt=`ls SYN_TEST/${sta}.${net}.BH${comp}*`
+ nsynt=`ls -1 SYN_TEST/${sta}.${net}.BH${comp}* | wc | awk '{print $1}'`
+
out=MEASURE_TEST/${sta}.${net}.${chan}
- if [ -e $synt ] ; then
+ if [ $nsynt == 1 ] ; then # data with unique syn match
- b1=`saclst b f $file | awk '{print $2}'`
- b2=`saclst b f $synt | awk '{print $2}'`
- db=$(echo "sqrt(($b1 - $b2)^2)" | bc)
- dt1=`saclst delta f $file | awk '{print $2}'`
- dt2=`saclst delta f $synt | awk '{print $2}'`
+ synt=`ls SYN_TEST/${sta}.${net}.BH${comp}*`
+ b1=`saclst b f $file | awk '{print $2}'`
+ b2=`saclst b f $synt | awk '{print $2}'`
+ db=$(echo "sqrt(($b1 - $b2)^2)" | bc)
+ dt1=`saclst delta f $file | awk '{print $2}'`
+ dt2=`saclst delta f $synt | awk '{print $2}'`
flag=`echo "$db 2*$dt1" | awk '{if ($1 < $2) print -1; else print 1}'`
#if [ `echo "$dt1 - $dt2"| bc` = 0 ] ; then
@@ -42,7 +45,7 @@
#fi
else
- echo $synt >> t3
+ echo "SYN_TEST/${sta}.${net}.BH${comp}*" >> t3
fi
done
@@ -50,17 +53,17 @@
# input file for window code
echo $n > input_test
cat t1 >> input_test
-rm t1
+rm -f t1
# list of probable bad data files
-rm input_bad_test
-cat t2 >> input_bad_test
-rm t2
+if [ -e "t2" ]; then
+ mv t2 input_test_bad
+fi
# list of missing synthetic files
-rm input_no_syn_test
-cat t3 >> input_no_syn_test
-rm t3
+if [ -e "t3" ]; then
+ mv t3 input_test_no_syn
+fi
@@ -69,4 +72,3 @@
-
Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/process_data_and_syn.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/process_data_and_syn.pl 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/prepare_scripts/socal/process_data_and_syn.pl 2008-10-15 18:06:06 UTC (rev 13084)
@@ -3,7 +3,7 @@
#==========================================================
#
# Carl Tape
-# 24-July-2007
+# 01-Oct-2008
# process_data_and_syn.pl
#
# This script processes data and 3D-Socal-SEM synthetics for Southern California.
@@ -20,37 +20,30 @@
# Trange band-pass filter
#
# ORDER OF OPERATIONS:
-# ~/UTILS/process_data_and_syn.pl 1 m07 1 0 d d 20 6/40 PROCESSED # data, initial pre-processing
-# ~/UTILS/process_data_and_syn.pl 1 m07 0 1 d d 20 6/40 PROCESSED # syn, initial pre-processing
-# ~/UTILS/process_data_and_syn.pl 2 m07 1 1 d d 20 6/40 PROCESSED # both, create cut files
-# ~/UTILS/process_data_and_syn.pl 3 m07 1 0 d d 20 6/40 PROCESSED # data, execute cut file
-# ~/UTILS/process_data_and_syn.pl 3 m07 0 1 d d 20 6/40 PROCESSED # syn, execute cut file
+# ~/UTILS/process_data_and_syn.pl 1 m10 1 0 d d 20 6/30 PROCESSED # data, initial pre-processing
+# ~/UTILS/process_data_and_syn.pl 1 m10 0 1 d d 20 6/30 PROCESSED # syn, initial pre-processing
+# ~/UTILS/process_data_and_syn.pl 2 m10 1 1 d d 20 6/30 PROCESSED # both, create cut files
+# ~/UTILS/process_data_and_syn.pl 3 m10 1 0 d d 20 6/30 PROCESSED # data, execute cut file
+# ~/UTILS/process_data_and_syn.pl 3 m10 0 1 d d 20 6/30 PROCESSED # syn, execute cut file
#
-# ~/UTILS/process_data_and_syn.pl 4 m07 1 0 d d 20 6/40 PROCESSED # data, bandpass T=6-40
-# ~/UTILS/process_data_and_syn.pl 4 m07 0 1 d d 20 6/40 PROCESSED # syn, bandpass T=6-40
+# ~/UTILS/process_data_and_syn.pl 4 m10 1 0 d d 20 6/30 PROCESSED # data, bandpass T=6-30
+# ~/UTILS/process_data_and_syn.pl 4 m10 0 1 d d 20 6/30 PROCESSED # syn, bandpass T=6-30
#
-# ~/UTILS/process_data_and_syn.pl 4 m07 1 0 d d 20 2/40 PROCESSED # data, bandpass T=2-40
-# ~/UTILS/process_data_and_syn.pl 4 m07 0 1 d d 20 2/40 PROCESSED # syn, bandpass T=2-40
+# ~/UTILS/process_data_and_syn.pl 4 m10 1 0 d d 20 2/30 PROCESSED # data, bandpass T=2-30
+# ~/UTILS/process_data_and_syn.pl 4 m10 0 1 d d 20 2/30 PROCESSED # syn, bandpass T=2-30
#
+# Chino Hills synthetics for Hiroo:
+# ~/UTILS/process_data_and_syn.pl 4 m10 0 1 d d 20 50/150 PROCESSED
+#
#==========================================================
if (@ARGV < 9) {die("Usage: process_data_and_syn.pl iprocess smodel idata isyn syn_ext dat_ext sps Tmin/Tmax pdir\n")}
($iprocess,$smodel,$idata,$isyn,$syn_ext,$dat_ext,$sps,$Trange,$pdir) = @ARGV;
-#$smodel = "m05"; # KEY: model iteration index
$iexecute = 0;
$ilist = 1;
-
$iprocess0 = 0;
-# THESE MUST BE DONE IN SEQUENTIAL ORDER
-#$iprocess1 = 0; # only done ONCE, no matter how many band-pass filters you use
-#$iprocess2 = 1; # cut records and pad zeros (and bandpass, if iprocess3=1)
-#$iprocess3 = 0; # bandpass filter (iprocess2=1 also)
-
-#$Lrange = "-l -40/220";
-#$Lrange = " ";
-
# USER PARAMETERS
$tfac = 1.0; # factor to extend lengths of records (should be > 1.0)
$itest = 0; # test directories or not
@@ -62,19 +55,23 @@
# directories
#$CMT_list = "/net/sierra/raid1/carltape/socal/socal_3D/SYN/model_${smodel}";
-$dir_source = "/net/sierra/raid1/carltape/results/SOURCES/socal_9";
-$dirCMT = "${dir_source}/v9_files";
-$CMT_list = "${dir_source}/EIDs_only_eid";
+$dir_source = "/net/sierra/raid1/carltape/results/SOURCES/socal_10";
+$dirCMT = "${dir_source}/v10_files";
+#$CMT_list = "${dir_source}/EIDs_only_eid";
+$CMT_list = "/net/sierra/raid1/carltape/results/EID_LISTS/syn_run_${smodel}";
$dirdat0 = "/net/sierra/raid1/carltape/socal/socal_3D/DATA/FINAL";
$dirsyn0 = "/net/sierra/raid1/carltape/socal/socal_3D/SYN/model_${smodel}";
#$dirsyn0 = "/net/sierra/raid1/carltape/socal/socal_3D/SYN/model_pre_${smodel}";
-# check if the STATIONS file exists
-#$stafile = "/net/denali/home2/carltape/gmt/stations/seismic/Matlab_output/STATIONS";
-$stafile = "/net/denali/home2/carltape/gmt/stations/seismic/Matlab_output/STATIONS_CALIFORNIA_TOMO_INNER_specfem";
-#$stafile = "/net/denali/home2/carltape/gmt/stations/seismic/Matlab_output/STATIONS_CALIFORNIA_TOMO_OUTER_specfem";
-if (not -f $stafile) {print "\n check if $stafile exists -- we will use the default"; $stafile = "none";}
+# STATIONS file
+#$stafile = "/net/denali/home1/carltape/gmt/stations/seismic/Matlab_output/STATIONS";
+$stafile = "/net/denali/home1/carltape/gmt/stations/seismic/Matlab_output/STATIONS_CALIFORNIA_TOMO_INNER_specfem";
+#$stafile = "/net/denali/home1/carltape/gmt/stations/seismic/Matlab_output/STATIONS_CALIFORNIA_TOMO_OUTER_specfem";
+# check if the files exist
+if (not -f $stafile) {print "\n check if $stafile exists";}
+if (not -f ${CMT_list}) {print "\n check if CMT_list ${CMT_list} exists";}
+
# period range
($Tmin,$Tmax) = split("/",$Trange);
$sTmin = sprintf("T%3.3i",$Tmin);
@@ -107,11 +104,11 @@
# write the C-shell script to file
if($idata==1 && $isyn==1) {
- $cshfile = "process_data_and_syn.csh";
+ $cshfile = "process_data_and_syn_${Ttag}.csh";
} elsif($idata==1 && $isyn==0) {
- $cshfile = "process_data.csh";
+ $cshfile = "process_data_${Ttag}.csh";
} elsif($idata==0 && $isyn==1) {
- $cshfile = "process_syn.csh";
+ $cshfile = "process_syn_${Ttag}.csh";
} else {
die("check idata and isyn\n");
}
@@ -124,8 +121,8 @@
}
$imin = 1; $imax = $ncmt; # default
-#$imin = 101; $imax = $ncmt;
-#$imin = 204; $imax = $imin;
+#$imin = 1; $imax = 10;
+#$imin = 95; $imax = $imin;
#----------------------------------------------------------------------
@@ -281,38 +278,24 @@
$dt = $deltad;
}
- # determine the cut for the records
- # b = earliest start time
- # e = earliest end time, multiplied by some factor
- if ($bd < $bs) {
- $b0 = $bd;
- } else {
- $b0 = $bs;
- }
- if ($b0 < $bmin) {
- $b0 = $bmin;
- }
- if ($ed < $es) {
- $e0 = $ed;
- } else {
- $e0 = $es;
- }
+ # determine the cut times for the records
+ $b0 = $bd;
+ #if ($bd < $bs) {$b0 = $bd} else {$b0 = $bs}
+ if ($b0 < $bmin) {$b0 = $bmin}
+ if ($ed < $es) {$e0 = $ed} else {$e0 = $es}
- # avoid events 14263716 and 14179292
- # --> adjust the simulation times, rather than the cut times
- #if($eid == 14179288 || $eid == 14263712) {
- # $e0 = 120;
- #}
-
$b = $b0;
$tlen0 = $e0 - $b;
$tlen = $tfac * $tlen0; # extend record length (if desired)
$e = $b0 + $tlen;
$npt = int( ($e-$b)/$dt );
- #print CUT "$datfile $synfile $b $e $npt $dt\n";
- print CUTDAT "$datfile $b $e $npt $dt\n";
- print CUTSYN "$synfile_base $b $e $npt $dt\n";
+ # print the cut times if the begin time is before the end time
+ if ($b < $e) {
+ #print CUT "$datfile $synfile $b $e $npt $dt\n";
+ print CUTDAT "$datfile $b $e $npt $dt\n";
+ print CUTSYN "$synfile_base $b $e $npt $dt\n";
+ }
if (0==1) {
print "\n Data : $bd $ed $deltad $nptd -- $tlend";
@@ -424,18 +407,22 @@
# KEY: indicate the base directory
$synfile = "${dirsyn_pro_1}/${synfile_base}";
- if (not -f $synfile) {die("synfile $synfile does not exist");}
+ if (not -f $synfile) {
+ print "synfile $synfile does not exist\n";
+ #die("synfile $synfile does not exist");
- print "$j out of $nlines -- $synfile\n";
- #print "-- $datfile -- $synfile -- $b -- $e -- $npt -- $dt -- \n";
+ } else {
+ print "$j out of $nlines -- $synfile\n";
+ #print "-- $datfile -- $synfile -- $b -- $e -- $npt -- $dt -- \n";
- # cut records and fill zeros
- `echo r $synfile >> $sacfile`;
- `echo cuterr fillz >> $sacfile`;
- `echo "cut $b n $npt" >> $sacfile`;
- `echo r $synfile >> $sacfile`;
- `echo cut off >> $sacfile`;
- `echo w over >> $sacfile`;
+ # cut records and fill zeros
+ `echo r $synfile >> $sacfile`;
+ `echo cuterr fillz >> $sacfile`;
+ `echo "cut $b n $npt" >> $sacfile`;
+ `echo r $synfile >> $sacfile`;
+ `echo cut off >> $sacfile`;
+ `echo w over >> $sacfile`;
+ }
}
`echo quit >> $sacfile`;
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_Daniel
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_Daniel 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_Daniel 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,68 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .false.
-BODY_WAVE_ONLY = .true.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-WIN_MIN_PERIOD = 2.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level
-STALTA_BASE = 0.07
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 2.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.85
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.0
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 2.5
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 4.0
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 1.0
-C_1 (small windows) = 4.0
-C_2 (prominence) = 1.0
-C_3a (separation height) = 4.0
-C_3b (separation time) = 2.5
-C_4a (curtail on left) = 2.0
-C_4b (curtail on right) = 6.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 1.0
-WEIGHT_N_WINDOWS = 0.0
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m07 (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m07)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m07 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m07 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 2.00
+WIN_MAX_PERIOD = 30.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.07
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 3.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.85
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.0
+C_1 (small windows) = 5.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 5.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m07
___________________________________________________________________
Name: svn:mergeinfo
+
Added: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m08
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m08 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m08 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 2.00
+WIN_MAX_PERIOD = 30.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.07
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 3.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.85
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.0
+C_1 (small windows) = 5.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 5.0
Added: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m09
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m09 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m09 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 2.00
+WIN_MAX_PERIOD = 30.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.07
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 3.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.85
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.0
+C_1 (small windows) = 5.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 5.0
Added: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m10
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m10 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T030_m10 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 2.00
+WIN_MAX_PERIOD = 30.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.07
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 2.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.85
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.0
+C_1 (small windows) = 5.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 5.0
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_Daniel (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_Daniel)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_Daniel (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_Daniel 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,68 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .false.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+WIN_MIN_PERIOD = 2.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.07
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 2.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.85
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.0
+C_1 (small windows) = 4.0
+C_2 (prominence) = 1.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 1.0
+WEIGHT_N_WINDOWS = 0.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_Daniel
___________________________________________________________________
Name: svn:mergeinfo
+
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m00 (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m00)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m00 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m00 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 2.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.07
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 5.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.85
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.0
+C_1 (small windows) = 5.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 5.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m00
___________________________________________________________________
Name: svn:mergeinfo
+
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m00_orig (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m00_orig)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m00_orig (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m00_orig 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 2.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.07
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 6.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.85
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.0
+C_1 (small windows) = 4.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 5.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m00_orig
___________________________________________________________________
Name: svn:mergeinfo
+
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m01 (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m01)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m01 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m01 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 2.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.07
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 3.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.85
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.0
+C_1 (small windows) = 4.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 5.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m01
___________________________________________________________________
Name: svn:mergeinfo
+
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m02 (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m02)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m02 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m02 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 2.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.07
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 3.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.85
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.0
+C_1 (small windows) = 4.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 5.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m02
___________________________________________________________________
Name: svn:mergeinfo
+
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m03 (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m03)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m03 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m03 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 2.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.07
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 3.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.85
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.0
+C_1 (small windows) = 4.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 5.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m03
___________________________________________________________________
Name: svn:mergeinfo
+
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m04 (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m04)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m04 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m04 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 2.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.07
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 3.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.85
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.0
+C_1 (small windows) = 4.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 5.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m04
___________________________________________________________________
Name: svn:mergeinfo
+
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m05 (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m05)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m05 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m05 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 2.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.07
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 3.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.85
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.0
+C_1 (small windows) = 4.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 5.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m05
___________________________________________________________________
Name: svn:mergeinfo
+
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m06 (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m06)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m06 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m06 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 2.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.07
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 3.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.85
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.0
+C_1 (small windows) = 4.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 5.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_T040_m06
___________________________________________________________________
Name: svn:mergeinfo
+
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m00
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m00 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m00 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .true.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 2.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level
-STALTA_BASE = 0.07
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 5.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.85
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.0
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 2.5
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 4.0
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 1.0
-C_1 (small windows) = 5.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 4.0
-C_3b (separation time) = 2.5
-C_4a (curtail on left) = 2.0
-C_4b (curtail on right) = 6.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 5.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m00_orig
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m00_orig 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m00_orig 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .true.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 2.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level
-STALTA_BASE = 0.07
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 6.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.85
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.0
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 2.5
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 4.0
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 1.0
-C_1 (small windows) = 4.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 4.0
-C_3b (separation time) = 2.5
-C_4a (curtail on left) = 2.0
-C_4b (curtail on right) = 6.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 5.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m01
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m01 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m01 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .true.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 2.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level
-STALTA_BASE = 0.07
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 3.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.85
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.0
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 2.5
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 4.0
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 1.0
-C_1 (small windows) = 4.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 4.0
-C_3b (separation time) = 2.5
-C_4a (curtail on left) = 2.0
-C_4b (curtail on right) = 6.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 5.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m02
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m02 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m02 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .true.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 2.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level
-STALTA_BASE = 0.07
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 3.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.85
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.0
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 2.5
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 4.0
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 1.0
-C_1 (small windows) = 4.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 4.0
-C_3b (separation time) = 2.5
-C_4a (curtail on left) = 2.0
-C_4b (curtail on right) = 6.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 5.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m03
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m03 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m03 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .true.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 2.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level
-STALTA_BASE = 0.07
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 3.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.85
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.0
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 2.5
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 4.0
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 1.0
-C_1 (small windows) = 4.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 4.0
-C_3b (separation time) = 2.5
-C_4a (curtail on left) = 2.0
-C_4b (curtail on right) = 6.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 5.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m04
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m04 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m04 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .true.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 2.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level
-STALTA_BASE = 0.07
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 3.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.85
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.0
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 2.5
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 4.0
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 1.0
-C_1 (small windows) = 4.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 4.0
-C_3b (separation time) = 2.5
-C_4a (curtail on left) = 2.0
-C_4b (curtail on right) = 6.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 5.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m05
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m05 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m05 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .true.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 2.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level
-STALTA_BASE = 0.07
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 3.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.85
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.0
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 2.5
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 4.0
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 1.0
-C_1 (small windows) = 4.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 4.0
-C_3b (separation time) = 2.5
-C_4a (curtail on left) = 2.0
-C_4b (curtail on right) = 6.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 5.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m06
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m06 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m06 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .true.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 2.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level
-STALTA_BASE = 0.07
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 3.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.85
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.0
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 2.5
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 4.0
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 1.0
-C_1 (small windows) = 4.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 4.0
-C_3b (separation time) = 2.5
-C_4a (curtail on left) = 2.0
-C_4b (curtail on right) = 6.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 5.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m07
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m07 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T002_m07 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .true.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 2.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level
-STALTA_BASE = 0.07
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 3.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.85
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.0
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 2.5
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 4.0
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 1.0
-C_1 (small windows) = 5.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 4.0
-C_3b (separation time) = 2.5
-C_4a (curtail on left) = 2.0
-C_4b (curtail on right) = 6.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 5.0
Added: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T003_T030_m07
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T003_T030_m07 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T003_T030_m07 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 3.00
+WIN_MAX_PERIOD = 30.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.11
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 4.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.80
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.3
+C_1 (small windows) = 4.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 0.70
+WEIGHT_AVERAGE_CC = 0.25
+WEIGHT_N_WINDOWS = 0.05
Added: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T003_T030_m08
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T003_T030_m08 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T003_T030_m08 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 3.00
+WIN_MAX_PERIOD = 30.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.11
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 4.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.80
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.3
+C_1 (small windows) = 4.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 0.70
+WEIGHT_AVERAGE_CC = 0.25
+WEIGHT_N_WINDOWS = 0.05
Added: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T003_T030_m09
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T003_T030_m09 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T003_T030_m09 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 3.00
+WIN_MAX_PERIOD = 30.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.11
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 4.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.80
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.3
+C_1 (small windows) = 4.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 0.70
+WEIGHT_AVERAGE_CC = 0.25
+WEIGHT_N_WINDOWS = 0.05
Added: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T003_T030_m10
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T003_T030_m10 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T003_T030_m10 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .true.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 3.00
+WIN_MAX_PERIOD = 30.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.11
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 3.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.80
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 2.5
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 4.0
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 1.3
+C_1 (small windows) = 4.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 4.0
+C_3b (separation time) = 2.5
+C_4a (curtail on left) = 2.0
+C_4b (curtail on right) = 6.0
+
+WEIGHT_SPACE_COVERAGE = 0.70
+WEIGHT_AVERAGE_CC = 0.25
+WEIGHT_N_WINDOWS = 0.05
Added: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T004_T040_m07
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T004_T040_m07 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T004_T040_m07 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 4.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level (0.23)
+STALTA_BASE = 0.22
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 4.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.80
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.5
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 3.0
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 2.5
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 0.7
+C_1 (small windows) = 3.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 3.0
+C_3b (separation time) = 2.0
+C_4a (curtail on left) = 2.5
+C_4b (curtail on right) = 12.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 0.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T004_m07
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T004_m07 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T004_m07 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .false.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 4.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level (0.23)
-STALTA_BASE = 0.22
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 4.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.80
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.5
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 3.0
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 2.5
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 0.7
-C_1 (small windows) = 3.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 3.0
-C_3b (separation time) = 2.0
-C_4a (curtail on left) = 2.5
-C_4b (curtail on right) = 12.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 0.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_14095540
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_14095540 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_14095540 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,68 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .false.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-WIN_MIN_PERIOD = 6.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level
-STALTA_BASE = 0.07
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 5.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.85
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.0
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 2.5
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 4.0
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 1.0
-C_1 (small windows) = 4.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 4.0
-C_3b (separation time) = 2.5
-C_4a (curtail on left) = 2.0
-C_4b (curtail on right) = 6.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 5.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_Daniel
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_Daniel 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_Daniel 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .false.
-BODY_WAVE_ONLY = .false.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 6.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level
-STALTA_BASE = 0.23
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 3.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.7
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.0
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 3.0
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 2.5
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 0.7
-C_1 (small windows) = 1.5
-C_2 (prominence) = 0.3
-C_3a (separation height) = 3.0
-C_3b (separation time) = 2.0
-C_4a (curtail on left) = 1.5
-C_4b (curtail on right) = 4.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 1.0
-WEIGHT_N_WINDOWS = 0.0
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m07 (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m07)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m07 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m07 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 6.00
+WIN_MAX_PERIOD = 30.00
+
+# -------------------------------------------------------------
+# E(t) water level (0.23)
+STALTA_BASE = 0.18
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 4.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.71
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.5
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 3.0
+SNR_MAX_BASE = 2.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 2.5
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 0.7
+C_1 (small windows) = 2.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 3.0
+C_3b (separation time) = 2.0
+C_4a (curtail on left) = 2.5
+C_4b (curtail on right) = 12.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 0.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m07
___________________________________________________________________
Name: svn:mergeinfo
+
Added: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m08
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m08 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m08 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 6.00
+WIN_MAX_PERIOD = 30.00
+
+# -------------------------------------------------------------
+# E(t) water level (0.23)
+STALTA_BASE = 0.18
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 4.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.71
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.5
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 3.0
+SNR_MAX_BASE = 2.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 2.5
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 0.7
+C_1 (small windows) = 2.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 3.0
+C_3b (separation time) = 2.0
+C_4a (curtail on left) = 2.5
+C_4b (curtail on right) = 12.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 0.0
Added: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m09
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m09 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m09 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 6.00
+WIN_MAX_PERIOD = 30.00
+
+# -------------------------------------------------------------
+# E(t) water level (0.23)
+STALTA_BASE = 0.18
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 4.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.71
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.5
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 3.0
+SNR_MAX_BASE = 2.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 2.5
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 0.7
+C_1 (small windows) = 2.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 3.0
+C_3b (separation time) = 2.0
+C_4a (curtail on left) = 2.5
+C_4b (curtail on right) = 12.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 0.0
Added: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m10
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m10 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T030_m10 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 6.00
+WIN_MAX_PERIOD = 30.00
+
+# -------------------------------------------------------------
+# E(t) water level (0.23)
+STALTA_BASE = 0.18
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 4.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.71
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.5
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 3.0
+SNR_MAX_BASE = 2.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 2.5
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 0.7
+C_1 (small windows) = 2.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 3.0
+C_3b (separation time) = 2.0
+C_4a (curtail on left) = 2.5
+C_4b (curtail on right) = 12.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 0.0
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_Daniel (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_Daniel)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_Daniel (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_Daniel 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .false.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 6.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level
+STALTA_BASE = 0.23
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 3.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.7
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.0
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 3.0
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 2.5
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 0.7
+C_1 (small windows) = 1.5
+C_2 (prominence) = 0.3
+C_3a (separation height) = 3.0
+C_3b (separation time) = 2.0
+C_4a (curtail on left) = 1.5
+C_4b (curtail on right) = 4.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 1.0
+WEIGHT_N_WINDOWS = 0.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_Daniel
___________________________________________________________________
Name: svn:mergeinfo
+
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m00 (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m00)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m00 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m00 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 6.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level (0.23)
+STALTA_BASE = 0.18
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 10.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.71
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.5
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 3.0
+SNR_MAX_BASE = 2.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 2.5
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 0.7
+C_1 (small windows) = 2.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 3.0
+C_3b (separation time) = 2.0
+C_4a (curtail on left) = 2.5
+C_4b (curtail on right) = 12.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 0.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m00
___________________________________________________________________
Name: svn:mergeinfo
+
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m01 (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m01)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m01 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m01 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 6.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level (0.23)
+STALTA_BASE = 0.22
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 5.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.74
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.5
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 3.0
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 2.5
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 0.7
+C_1 (small windows) = 2.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 3.0
+C_3b (separation time) = 2.0
+C_4a (curtail on left) = 2.5
+C_4b (curtail on right) = 12.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 0.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m01
___________________________________________________________________
Name: svn:mergeinfo
+
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m02 (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m02)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m02 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m02 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 6.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level (0.23)
+STALTA_BASE = 0.22
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 5.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.74
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.5
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 3.0
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 2.5
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 0.7
+C_1 (small windows) = 2.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 3.0
+C_3b (separation time) = 2.0
+C_4a (curtail on left) = 2.5
+C_4b (curtail on right) = 12.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 0.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m02
___________________________________________________________________
Name: svn:mergeinfo
+
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m03 (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m03)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m03 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m03 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 6.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level (0.23)
+STALTA_BASE = 0.22
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 5.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.74
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.5
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 3.0
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 2.5
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 0.7
+C_1 (small windows) = 2.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 3.0
+C_3b (separation time) = 2.0
+C_4a (curtail on left) = 2.5
+C_4b (curtail on right) = 12.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 0.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m03
___________________________________________________________________
Name: svn:mergeinfo
+
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m04 (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m04)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m04 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m04 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 6.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level (0.23)
+STALTA_BASE = 0.22
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 4.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.74
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.5
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 3.0
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 2.5
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 0.7
+C_1 (small windows) = 2.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 3.0
+C_3b (separation time) = 2.0
+C_4a (curtail on left) = 2.5
+C_4b (curtail on right) = 12.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 0.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m04
___________________________________________________________________
Name: svn:mergeinfo
+
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m05 (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m05)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m05 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m05 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 6.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level (0.23)
+STALTA_BASE = 0.22
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 4.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.74
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.5
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 3.0
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 2.5
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 0.7
+C_1 (small windows) = 2.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 3.0
+C_3b (separation time) = 2.0
+C_4a (curtail on left) = 2.5
+C_4b (curtail on right) = 12.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 0.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m05
___________________________________________________________________
Name: svn:mergeinfo
+
Copied: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m06 (from rev 13083, seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m06)
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m06 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m06 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 6.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level (0.23)
+STALTA_BASE = 0.22
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 4.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.74
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.5
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 3.0
+SNR_MAX_BASE = 3.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 2.5
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 0.7
+C_1 (small windows) = 2.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 3.0
+C_3b (separation time) = 2.0
+C_4a (curtail on left) = 2.5
+C_4b (curtail on right) = 12.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 0.0
Property changes on: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m06
___________________________________________________________________
Name: svn:mergeinfo
+
Added: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m07
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m07 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_T040_m07 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,69 @@
+# -------------------------------------------------------------
+#
+# This is the parameter file for FLEXWIN. It is based on the
+# same syntax as the Par_file for SPECFEM. Variable names are
+# put first, values are placed after the 34th column.
+#
+# Comment lines and blank lines are significant. If you
+# change the layout of this file or add/remove parameters
+# you must also modify the user_variables module and the
+# read_parameter_file subroutine at the start of seismo_subs.f90.
+#
+# -------------------------------------------------------------
+
+# -------------------------------------------------------------
+# boolean parameters
+DEBUG = .true.
+MAKE_SEISMO_PLOTS = .true.
+MAKE_WINDOW_FILES = .true.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+RUN_BANDPASS = .false.
+WIN_MIN_PERIOD = 6.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level (0.23)
+STALTA_BASE = 0.18
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 4.0
+
+# -------------------------------------------------------------
+# limit on CC for window acceptance
+CC_BASE = 0.71
+
+# -------------------------------------------------------------
+# limit on dlnA (dA/A) for window acceptance
+DLNA_BASE = 1.5
+
+# -------------------------------------------------------------
+# boolean switch for check_data_quality
+DATA_QUALITY = .true.
+
+# if DATA_QUALITY = .true. and if two different measurements of
+# signal-to-noise ratios exceeds these two base levels,
+# then the data time series (and syn) is kept
+SNR_INTEGRATE_BASE = 3.0
+SNR_MAX_BASE = 2.5
+
+# -------------------------------------------------------------
+# limit on signal to noise ratio in a particular window.
+WINDOW_AMP_BASE = 2.5
+
+# -------------------------------------------------------------
+# Fine tuning constants
+C_0 (internal minima) = 0.7
+C_1 (small windows) = 2.0
+C_2 (prominence) = 0.0
+C_3a (separation height) = 3.0
+C_3b (separation time) = 2.0
+C_4a (curtail on left) = 2.5
+C_4b (curtail on right) = 12.0
+
+WEIGHT_SPACE_COVERAGE = 1.0
+WEIGHT_AVERAGE_CC = 0.0
+WEIGHT_N_WINDOWS = 0.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m00
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m00 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m00 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .false.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 6.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level (0.23)
-STALTA_BASE = 0.18
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 10.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.71
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.5
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 3.0
-SNR_MAX_BASE = 2.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 2.5
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 0.7
-C_1 (small windows) = 2.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 3.0
-C_3b (separation time) = 2.0
-C_4a (curtail on left) = 2.5
-C_4b (curtail on right) = 12.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 0.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m01
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m01 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m01 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .false.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 6.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level (0.23)
-STALTA_BASE = 0.22
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 5.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.74
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.5
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 3.0
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 2.5
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 0.7
-C_1 (small windows) = 2.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 3.0
-C_3b (separation time) = 2.0
-C_4a (curtail on left) = 2.5
-C_4b (curtail on right) = 12.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 0.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m02
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m02 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m02 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .false.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 6.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level (0.23)
-STALTA_BASE = 0.22
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 5.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.74
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.5
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 3.0
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 2.5
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 0.7
-C_1 (small windows) = 2.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 3.0
-C_3b (separation time) = 2.0
-C_4a (curtail on left) = 2.5
-C_4b (curtail on right) = 12.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 0.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m03
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m03 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m03 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .false.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 6.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level (0.23)
-STALTA_BASE = 0.22
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 5.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.74
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.5
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 3.0
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 2.5
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 0.7
-C_1 (small windows) = 2.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 3.0
-C_3b (separation time) = 2.0
-C_4a (curtail on left) = 2.5
-C_4b (curtail on right) = 12.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 0.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m04
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m04 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m04 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .false.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 6.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level (0.23)
-STALTA_BASE = 0.22
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 4.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.74
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.5
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 3.0
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 2.5
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 0.7
-C_1 (small windows) = 2.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 3.0
-C_3b (separation time) = 2.0
-C_4a (curtail on left) = 2.5
-C_4b (curtail on right) = 12.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 0.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m05
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m05 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m05 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .false.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 6.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level (0.23)
-STALTA_BASE = 0.22
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 4.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.74
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.5
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 3.0
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 2.5
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 0.7
-C_1 (small windows) = 2.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 3.0
-C_3b (separation time) = 2.0
-C_4a (curtail on left) = 2.5
-C_4b (curtail on right) = 12.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 0.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m06
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m06 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m06 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .false.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 6.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level (0.23)
-STALTA_BASE = 0.22
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 4.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.74
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.5
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 3.0
-SNR_MAX_BASE = 3.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 2.5
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 0.7
-C_1 (small windows) = 2.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 3.0
-C_3b (separation time) = 2.0
-C_4a (curtail on left) = 2.5
-C_4b (curtail on right) = 12.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 0.0
Deleted: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m07
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m07 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/PAR_FILE_T006_m07 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,69 +0,0 @@
-# -------------------------------------------------------------
-#
-# This is the parameter file for FLEXWIN. It is based on the
-# same syntax as the Par_file for SPECFEM. Variable names are
-# put first, values are placed after the 34th column.
-#
-# Comment lines and blank lines are significant. If you
-# change the layout of this file or add/remove parameters
-# you must also modify the user_variables module and the
-# read_parameter_file subroutine at the start of seismo_subs.f90.
-#
-# -------------------------------------------------------------
-
-# -------------------------------------------------------------
-# boolean parameters
-DEBUG = .true.
-MAKE_SEISMO_PLOTS = .true.
-MAKE_WINDOW_FILES = .true.
-BODY_WAVE_ONLY = .false.
-
-# -------------------------------------------------------------
-# period min/max for filtering
-RUN_BANDPASS = .false.
-WIN_MIN_PERIOD = 6.00
-WIN_MAX_PERIOD = 40.00
-
-# -------------------------------------------------------------
-# E(t) water level (0.23)
-STALTA_BASE = 0.18
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 4.0
-
-# -------------------------------------------------------------
-# limit on CC for window acceptance
-CC_BASE = 0.71
-
-# -------------------------------------------------------------
-# limit on dlnA (dA/A) for window acceptance
-DLNA_BASE = 1.5
-
-# -------------------------------------------------------------
-# boolean switch for check_data_quality
-DATA_QUALITY = .true.
-
-# if DATA_QUALITY = .true. and if two different measurements of
-# signal-to-noise ratios exceeds these two base levels,
-# then the data time series (and syn) is kept
-SNR_INTEGRATE_BASE = 3.0
-SNR_MAX_BASE = 2.5
-
-# -------------------------------------------------------------
-# limit on signal to noise ratio in a particular window.
-WINDOW_AMP_BASE = 2.5
-
-# -------------------------------------------------------------
-# Fine tuning constants
-C_0 (internal minima) = 0.7
-C_1 (small windows) = 2.0
-C_2 (prominence) = 0.0
-C_3a (separation height) = 3.0
-C_3b (separation time) = 2.0
-C_4a (curtail on left) = 2.5
-C_4b (curtail on right) = 12.0
-
-WEIGHT_SPACE_COVERAGE = 1.0
-WEIGHT_AVERAGE_CC = 0.0
-WEIGHT_N_WINDOWS = 0.0
Modified: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m07.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m07.f90 2008-10-15 17:13:10 UTC (rev 13083)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m07.f90 2008-10-15 18:06:06 UTC (rev 13084)
@@ -1,4 +1,21 @@
! -------------------------------------------------------------
+! edit here to change T0 and T1 on some condition
+! Note, this function is called AFTER the seismogram has been
+! read but before it is filtered.
+! -------------------------------------------------------------
+
+subroutine modify_T0_T1_on_condition
+ use seismo_variables
+
+ ! do nothing
+
+ ! adjust fstart and fend accordingly
+ !FSTART=1./WIN_MAX_PERIOD
+ !FEND=1./WIN_MIN_PERIOD
+
+end subroutine modify_T0_T1_on_condition
+
+! -------------------------------------------------------------
! Edit here to change the time dependent properties of the selection criteria
! Note, this function is called AFTER the seismogram has been read.
! -------------------------------------------------------------
@@ -60,7 +77,10 @@
signal_end = Sw_end
noise_end = Pnl_start
if(DEBUG) then
- write(*,*) 'DEBUG : P_pick = ', sngl(P_pick)
+ if (BODY_WAVE_ONLY) then
+ write(*,*) 'DEBUG : P_pick = ', P_pick
+ write(*,*) 'DEBUG : S_pick = ', S_pick
+ endif
write(*,*) 'DEBUG : signal_end = ', sngl(signal_end)
write(*,*) 'DEBUG : noise_end = ', sngl(noise_end)
endif
@@ -77,12 +97,11 @@
! raises STA/LTA water level after surface wave arrives
if (BODY_WAVE_ONLY) then
-
!if(time.gt.S_end) then
if(time.gt.Sw_end) then
STALTA_W_LEVEL(i) = 10.*STALTA_BASE
endif
-
+
else
!!$ ! set time- and distance-specific Tshift and DlnA to mimic Qinya's criteria
!!$ ! (see Liu et al., 2004, p. 1755; note comment above)
@@ -100,6 +119,11 @@
! STALTA_W_LEVEL(i) = 2.0*STALTA_BASE
!endif
+ ! allow for 100s to possibly capture additional phases
+ if(time.gt. (Sw_end+100.0) ) then
+ STALTA_W_LEVEL(i) = 10.*STALTA_BASE
+ endif
+
endif
enddo
Added: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m08.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m08.f90 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m08.f90 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,161 @@
+! -------------------------------------------------------------
+! edit here to change T0 and T1 on some condition
+! Note, this function is called AFTER the seismogram has been
+! read but before it is filtered.
+! -------------------------------------------------------------
+
+subroutine modify_T0_T1_on_condition
+ use seismo_variables
+
+ ! do nothing
+
+ ! adjust fstart and fend accordingly
+ !FSTART=1./WIN_MAX_PERIOD
+ !FEND=1./WIN_MIN_PERIOD
+
+end subroutine modify_T0_T1_on_condition
+
+! -------------------------------------------------------------
+! Edit here to change the time dependent properties of the selection criteria
+! Note, this function is called AFTER the seismogram has been read.
+! -------------------------------------------------------------
+subroutine set_up_criteria_arrays
+ use seismo_variables
+
+ integer :: i
+ double precision :: time
+
+ ! for qinya's scsn picking
+ double precision :: Pnl_start, S_end, Sw_start, Sw_end
+
+!===========================
+
+! -----------------------------------------------------------------
+! This is the basic version of the subroutine - no variation with time
+! -----------------------------------------------------------------
+ do i = 1, npts
+ time = b+(i-1)*dt
+ DLNA_LIMIT(i) = DLNA_BASE
+ CC_LIMIT(i) = CC_BASE
+ TSHIFT_LIMIT(i) = TSHIFT_BASE ! WIN_MIN_PERIOD/2.0
+ STALTA_W_LEVEL(i) = STALTA_BASE
+ S2N_LIMIT(i) = WINDOW_AMP_BASE
+ enddo
+
+!!$ if (.not. BODY_WAVE_ONLY) then
+!!$ Pnl_start = -5.0 + dist_km/7.8
+!!$ Sw_start = -15.0 + dist_km/3.5
+!!$ Sw_end = 35.0 + dist_km/3.1
+!!$ else
+!!$ Pnl_start = P_pick - 5.0
+!!$ S_end = S_pick + 5.0
+!!$ Sw_start = -15.0 + dist_km/3.5
+!!$ Sw_end = 35.0 + dist_km/3.1
+!!$ endif
+
+ ! regional (Qinya's formulation):
+ ! -------------------------------------------------------------
+ ! see Liu et al. (2004), p. 1755, but note that the PARENTHESES
+ ! that are listed in the publication should not be there
+ ! THESE ARE PROBABLY NOT ACCURATE ENOUGH FOR LONGER PATHS.
+
+ Sw_start = -15.0 + dist_km/3.5
+ Sw_end = 35.0 + dist_km/3.1
+
+ if (BODY_WAVE_ONLY) then
+ !Pnl_start = P_pick - 5.0
+ !S_end = S_pick + 5.0
+ Pnl_start = P_pick - 2.5*WIN_MIN_PERIOD
+ S_end = S_pick + 2.5*WIN_MIN_PERIOD
+
+ else
+ Pnl_start = -5.0 + dist_km/7.8
+ S_end = Sw_start
+ endif
+
+ ! variables for signal to noise ratio criteria.
+ signal_end = Sw_end
+ noise_end = Pnl_start
+ if(DEBUG) then
+ if (BODY_WAVE_ONLY) then
+ write(*,*) 'DEBUG : P_pick = ', P_pick
+ write(*,*) 'DEBUG : S_pick = ', S_pick
+ endif
+ write(*,*) 'DEBUG : signal_end = ', sngl(signal_end)
+ write(*,*) 'DEBUG : noise_end = ', sngl(noise_end)
+ endif
+
+ ! --------------------------------
+ ! modulate criteria in time
+ do i = 1, npts
+ time = b+(i-1)*dt ! time
+
+ ! raises STA/LTA water level before P wave arrival.
+ if(time.lt.Pnl_start) then
+ STALTA_W_LEVEL(i) = 10.*STALTA_BASE
+ endif
+
+ ! raises STA/LTA water level after surface wave arrives
+ if (BODY_WAVE_ONLY) then
+ !if(time.gt.S_end) then
+ if(time.gt.Sw_end) then
+ STALTA_W_LEVEL(i) = 10.*STALTA_BASE
+ endif
+
+ else
+!!$ ! set time- and distance-specific Tshift and DlnA to mimic Qinya's criteria
+!!$ ! (see Liu et al., 2004, p. 1755; note comment above)
+!!$ if(time.ge.Pnl_start .and. time.lt.Sw_start) then
+!!$ !DLNA_LIMIT(i) = 1.5 ! ratio is 2.5, and dlna is ratio-1
+!!$ TSHIFT_LIMIT(i) = 3.0 + dist_km/80.0
+!!$ endif
+!!$ if(time.ge.Sw_start .and. time.le.Sw_end) then
+!!$ !DLNA_LIMIT(i) = 1.5 ! ratio is 2.5, and dlna is ratio-1
+!!$ TSHIFT_LIMIT(i) = 3.0 + dist_km/50.0
+!!$ endif
+
+ ! double the STA/LTA water level after the surface waves
+ !if(time.gt.Sw_end) then
+ ! STALTA_W_LEVEL(i) = 2.0*STALTA_BASE
+ !endif
+
+ ! allow for 100s to possibly capture additional phases
+ if(time.gt. (Sw_end+100.0) ) then
+ STALTA_W_LEVEL(i) = 10.*STALTA_BASE
+ endif
+
+ endif
+
+ enddo
+
+! The following is for check_window quality_s2n
+
+! -----------------------------------------------------------------
+! Start of user-dependent portion
+
+! This is where you modulate the time dependence of the selection
+! criteria. You have access to the following parameters from the
+! seismogram itself:
+!
+! dt, b, kstnm, knetwk, kcmpnm
+! evla, evlo, stla, stlo, evdp, azimuth, backazimuth, dist_deg, dist_km
+! num_phases, ph_names, ph_times
+!
+! Example of modulation:
+!-----------------------
+! To increase s2n limit after arrival of R1 try
+!
+! R_vel=3.2
+! R_time=dist_km/R_vel
+! do i = 1, npts
+! time=b+(i-1)*dt
+! if (time.gt.R_time) then
+! S2N_LIMIT(i)=2*WINDOW_AMP_BASE
+! endif
+! enddo
+!
+! End of user-dependent portion
+! -----------------------------------------------------------------
+
+end subroutine set_up_criteria_arrays
+! -------------------------------------------------------------
Added: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m09.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m09.f90 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m09.f90 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,161 @@
+! -------------------------------------------------------------
+! edit here to change T0 and T1 on some condition
+! Note, this function is called AFTER the seismogram has been
+! read but before it is filtered.
+! -------------------------------------------------------------
+
+subroutine modify_T0_T1_on_condition
+ use seismo_variables
+
+ ! do nothing
+
+ ! adjust fstart and fend accordingly
+ !FSTART=1./WIN_MAX_PERIOD
+ !FEND=1./WIN_MIN_PERIOD
+
+end subroutine modify_T0_T1_on_condition
+
+! -------------------------------------------------------------
+! Edit here to change the time dependent properties of the selection criteria
+! Note, this function is called AFTER the seismogram has been read.
+! -------------------------------------------------------------
+subroutine set_up_criteria_arrays
+ use seismo_variables
+
+ integer :: i
+ double precision :: time
+
+ ! for qinya's scsn picking
+ double precision :: Pnl_start, S_end, Sw_start, Sw_end
+
+!===========================
+
+! -----------------------------------------------------------------
+! This is the basic version of the subroutine - no variation with time
+! -----------------------------------------------------------------
+ do i = 1, npts
+ time = b+(i-1)*dt
+ DLNA_LIMIT(i) = DLNA_BASE
+ CC_LIMIT(i) = CC_BASE
+ TSHIFT_LIMIT(i) = TSHIFT_BASE ! WIN_MIN_PERIOD/2.0
+ STALTA_W_LEVEL(i) = STALTA_BASE
+ S2N_LIMIT(i) = WINDOW_AMP_BASE
+ enddo
+
+!!$ if (.not. BODY_WAVE_ONLY) then
+!!$ Pnl_start = -5.0 + dist_km/7.8
+!!$ Sw_start = -15.0 + dist_km/3.5
+!!$ Sw_end = 35.0 + dist_km/3.1
+!!$ else
+!!$ Pnl_start = P_pick - 5.0
+!!$ S_end = S_pick + 5.0
+!!$ Sw_start = -15.0 + dist_km/3.5
+!!$ Sw_end = 35.0 + dist_km/3.1
+!!$ endif
+
+ ! regional (Qinya's formulation):
+ ! -------------------------------------------------------------
+ ! see Liu et al. (2004), p. 1755, but note that the PARENTHESES
+ ! that are listed in the publication should not be there
+ ! THESE ARE PROBABLY NOT ACCURATE ENOUGH FOR LONGER PATHS.
+
+ Sw_start = -15.0 + dist_km/3.5
+ Sw_end = 35.0 + dist_km/3.1
+
+ if (BODY_WAVE_ONLY) then
+ !Pnl_start = P_pick - 5.0
+ !S_end = S_pick + 5.0
+ Pnl_start = P_pick - 2.5*WIN_MIN_PERIOD
+ S_end = S_pick + 2.5*WIN_MIN_PERIOD
+
+ else
+ Pnl_start = -5.0 + dist_km/7.8
+ S_end = Sw_start
+ endif
+
+ ! variables for signal to noise ratio criteria.
+ signal_end = Sw_end
+ noise_end = Pnl_start
+ if(DEBUG) then
+ if (BODY_WAVE_ONLY) then
+ write(*,*) 'DEBUG : P_pick = ', P_pick
+ write(*,*) 'DEBUG : S_pick = ', S_pick
+ endif
+ write(*,*) 'DEBUG : signal_end = ', sngl(signal_end)
+ write(*,*) 'DEBUG : noise_end = ', sngl(noise_end)
+ endif
+
+ ! --------------------------------
+ ! modulate criteria in time
+ do i = 1, npts
+ time = b+(i-1)*dt ! time
+
+ ! raises STA/LTA water level before P wave arrival.
+ if(time.lt.Pnl_start) then
+ STALTA_W_LEVEL(i) = 10.*STALTA_BASE
+ endif
+
+ ! raises STA/LTA water level after surface wave arrives
+ if (BODY_WAVE_ONLY) then
+ !if(time.gt.S_end) then
+ if(time.gt.Sw_end) then
+ STALTA_W_LEVEL(i) = 10.*STALTA_BASE
+ endif
+
+ else
+!!$ ! set time- and distance-specific Tshift and DlnA to mimic Qinya's criteria
+!!$ ! (see Liu et al., 2004, p. 1755; note comment above)
+!!$ if(time.ge.Pnl_start .and. time.lt.Sw_start) then
+!!$ !DLNA_LIMIT(i) = 1.5 ! ratio is 2.5, and dlna is ratio-1
+!!$ TSHIFT_LIMIT(i) = 3.0 + dist_km/80.0
+!!$ endif
+!!$ if(time.ge.Sw_start .and. time.le.Sw_end) then
+!!$ !DLNA_LIMIT(i) = 1.5 ! ratio is 2.5, and dlna is ratio-1
+!!$ TSHIFT_LIMIT(i) = 3.0 + dist_km/50.0
+!!$ endif
+
+ ! double the STA/LTA water level after the surface waves
+ !if(time.gt.Sw_end) then
+ ! STALTA_W_LEVEL(i) = 2.0*STALTA_BASE
+ !endif
+
+ ! allow for 100s to possibly capture additional phases
+ if(time.gt. (Sw_end+100.0) ) then
+ STALTA_W_LEVEL(i) = 10.*STALTA_BASE
+ endif
+
+ endif
+
+ enddo
+
+! The following is for check_window quality_s2n
+
+! -----------------------------------------------------------------
+! Start of user-dependent portion
+
+! This is where you modulate the time dependence of the selection
+! criteria. You have access to the following parameters from the
+! seismogram itself:
+!
+! dt, b, kstnm, knetwk, kcmpnm
+! evla, evlo, stla, stlo, evdp, azimuth, backazimuth, dist_deg, dist_km
+! num_phases, ph_names, ph_times
+!
+! Example of modulation:
+!-----------------------
+! To increase s2n limit after arrival of R1 try
+!
+! R_vel=3.2
+! R_time=dist_km/R_vel
+! do i = 1, npts
+! time=b+(i-1)*dt
+! if (time.gt.R_time) then
+! S2N_LIMIT(i)=2*WINDOW_AMP_BASE
+! endif
+! enddo
+!
+! End of user-dependent portion
+! -----------------------------------------------------------------
+
+end subroutine set_up_criteria_arrays
+! -------------------------------------------------------------
Added: seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m10.f90
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m10.f90 (rev 0)
+++ seismo/3D/ADJOINT_TOMO/flexwin/user_files/socal_3D/user_functions_m10.f90 2008-10-15 18:06:06 UTC (rev 13084)
@@ -0,0 +1,161 @@
+! -------------------------------------------------------------
+! edit here to change T0 and T1 on some condition
+! Note, this function is called AFTER the seismogram has been
+! read but before it is filtered.
+! -------------------------------------------------------------
+
+subroutine modify_T0_T1_on_condition
+ use seismo_variables
+
+ ! do nothing
+
+ ! adjust fstart and fend accordingly
+ !FSTART=1./WIN_MAX_PERIOD
+ !FEND=1./WIN_MIN_PERIOD
+
+end subroutine modify_T0_T1_on_condition
+
+! -------------------------------------------------------------
+! Edit here to change the time dependent properties of the selection criteria
+! Note, this function is called AFTER the seismogram has been read.
+! -------------------------------------------------------------
+subroutine set_up_criteria_arrays
+ use seismo_variables
+
+ integer :: i
+ double precision :: time
+
+ ! for qinya's scsn picking
+ double precision :: Pnl_start, S_end, Sw_start, Sw_end
+
+!===========================
+
+! -----------------------------------------------------------------
+! This is the basic version of the subroutine - no variation with time
+! -----------------------------------------------------------------
+ do i = 1, npts
+ time = b+(i-1)*dt
+ DLNA_LIMIT(i) = DLNA_BASE
+ CC_LIMIT(i) = CC_BASE
+ TSHIFT_LIMIT(i) = TSHIFT_BASE ! WIN_MIN_PERIOD/2.0
+ STALTA_W_LEVEL(i) = STALTA_BASE
+ S2N_LIMIT(i) = WINDOW_AMP_BASE
+ enddo
+
+!!$ if (.not. BODY_WAVE_ONLY) then
+!!$ Pnl_start = -5.0 + dist_km/7.8
+!!$ Sw_start = -15.0 + dist_km/3.5
+!!$ Sw_end = 35.0 + dist_km/3.1
+!!$ else
+!!$ Pnl_start = P_pick - 5.0
+!!$ S_end = S_pick + 5.0
+!!$ Sw_start = -15.0 + dist_km/3.5
+!!$ Sw_end = 35.0 + dist_km/3.1
+!!$ endif
+
+ ! regional (Qinya's formulation):
+ ! -------------------------------------------------------------
+ ! see Liu et al. (2004), p. 1755, but note that the PARENTHESES
+ ! that are listed in the publication should not be there
+ ! THESE ARE PROBABLY NOT ACCURATE ENOUGH FOR LONGER PATHS.
+
+ Sw_start = -15.0 + dist_km/3.5
+ Sw_end = 35.0 + dist_km/3.1
+
+ if (BODY_WAVE_ONLY) then
+ !Pnl_start = P_pick - 5.0
+ !S_end = S_pick + 5.0
+ Pnl_start = P_pick - 2.5*WIN_MIN_PERIOD
+ S_end = S_pick + 2.5*WIN_MIN_PERIOD
+
+ else
+ Pnl_start = -5.0 + dist_km/7.8
+ S_end = Sw_start
+ endif
+
+ ! variables for signal to noise ratio criteria.
+ signal_end = Sw_end
+ noise_end = Pnl_start
+ if(DEBUG) then
+ if (BODY_WAVE_ONLY) then
+ write(*,*) 'DEBUG : P_pick = ', P_pick
+ write(*,*) 'DEBUG : S_pick = ', S_pick
+ endif
+ write(*,*) 'DEBUG : signal_end = ', sngl(signal_end)
+ write(*,*) 'DEBUG : noise_end = ', sngl(noise_end)
+ endif
+
+ ! --------------------------------
+ ! modulate criteria in time
+ do i = 1, npts
+ time = b+(i-1)*dt ! time
+
+ ! raises STA/LTA water level before P wave arrival.
+ if(time.lt.Pnl_start) then
+ STALTA_W_LEVEL(i) = 10.*STALTA_BASE
+ endif
+
+ ! raises STA/LTA water level after surface wave arrives
+ if (BODY_WAVE_ONLY) then
+ !if(time.gt.S_end) then
+ if(time.gt.Sw_end) then
+ STALTA_W_LEVEL(i) = 10.*STALTA_BASE
+ endif
+
+ else
+!!$ ! set time- and distance-specific Tshift and DlnA to mimic Qinya's criteria
+!!$ ! (see Liu et al., 2004, p. 1755; note comment above)
+!!$ if(time.ge.Pnl_start .and. time.lt.Sw_start) then
+!!$ !DLNA_LIMIT(i) = 1.5 ! ratio is 2.5, and dlna is ratio-1
+!!$ TSHIFT_LIMIT(i) = 3.0 + dist_km/80.0
+!!$ endif
+!!$ if(time.ge.Sw_start .and. time.le.Sw_end) then
+!!$ !DLNA_LIMIT(i) = 1.5 ! ratio is 2.5, and dlna is ratio-1
+!!$ TSHIFT_LIMIT(i) = 3.0 + dist_km/50.0
+!!$ endif
+
+ ! double the STA/LTA water level after the surface waves
+ !if(time.gt.Sw_end) then
+ ! STALTA_W_LEVEL(i) = 2.0*STALTA_BASE
+ !endif
+
+ ! allow for 100s to possibly capture additional phases
+ if(time.gt. (Sw_end+100.0) ) then
+ STALTA_W_LEVEL(i) = 10.*STALTA_BASE
+ endif
+
+ endif
+
+ enddo
+
+! The following is for check_window quality_s2n
+
+! -----------------------------------------------------------------
+! Start of user-dependent portion
+
+! This is where you modulate the time dependence of the selection
+! criteria. You have access to the following parameters from the
+! seismogram itself:
+!
+! dt, b, kstnm, knetwk, kcmpnm
+! evla, evlo, stla, stlo, evdp, azimuth, backazimuth, dist_deg, dist_km
+! num_phases, ph_names, ph_times
+!
+! Example of modulation:
+!-----------------------
+! To increase s2n limit after arrival of R1 try
+!
+! R_vel=3.2
+! R_time=dist_km/R_vel
+! do i = 1, npts
+! time=b+(i-1)*dt
+! if (time.gt.R_time) then
+! S2N_LIMIT(i)=2*WINDOW_AMP_BASE
+! endif
+! enddo
+!
+! End of user-dependent portion
+! -----------------------------------------------------------------
+
+end subroutine set_up_criteria_arrays
+! -------------------------------------------------------------
More information about the cig-commits
mailing list