[cig-commits] r12034 - in seismo/3D/automeasure: scripts scripts/matlab_socal user_files/socal_3D
carltape at geodynamics.org
carltape at geodynamics.org
Tue May 27 09:53:23 PDT 2008
Author: carltape
Date: 2008-05-27 09:53:22 -0700 (Tue, 27 May 2008)
New Revision: 12034
Added:
seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T02_m0
seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T02_m1
seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T02_m2
seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06_m0
seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06_m1
seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06_m2
seismo/3D/automeasure/user_files/socal_3D/user_functions_m0.f90
seismo/3D/automeasure/user_files/socal_3D/user_functions_m1.f90
seismo/3D/automeasure/user_files/socal_3D/user_functions_m2.f90
Removed:
seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T02
seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06
seismo/3D/automeasure/user_files/socal_3D/user_functions.f90
Modified:
seismo/3D/automeasure/scripts/extract_event_windowing_stats_carl.sh
seismo/3D/automeasure/scripts/matlab_socal/window_subset.m
seismo/3D/automeasure/scripts/pick_all_windows.pl
seismo/3D/automeasure/scripts/plot_windows_all.pl
seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06_14095540
Log:
Updated socal user files and also run scripts.
Modified: seismo/3D/automeasure/scripts/extract_event_windowing_stats_carl.sh
===================================================================
--- seismo/3D/automeasure/scripts/extract_event_windowing_stats_carl.sh 2008-05-27 15:09:17 UTC (rev 12033)
+++ seismo/3D/automeasure/scripts/extract_event_windowing_stats_carl.sh 2008-05-27 16:53:22 UTC (rev 12034)
@@ -155,6 +155,8 @@
# set output filename
out=${basename}/event_winstats.eps
+out_ps=${basename}/event_winstats.ps
+out_pdf=${basename}/event_winstats.pdf
short_basename=`echo $basename | awk -F"/" '{print $NF}'`
##evlo=`echo $evlo | awk '{printf "%.0f", $1}'`
@@ -244,6 +246,8 @@
# plot window record-section
out2=${basename}/event_recordsection.eps
+out2_ps=${basename}/event_recordsection.ps
+out2_pdf=${basename}/event_recordsection.pdf
gmtset PAPER_MEDIA letter MEASURE_UNIT inch
region=$min_time/$max_time/$min_ddg/$max_ddg
@@ -298,5 +302,11 @@
mv $rectxt $basename # CHT: move pstext file to data directory
+# make pdf files
+cp $out ${out_ps}
+ps2pdf ${out_ps} ${out_pdf}
+cp $out2 ${out2_ps}
+ps2pdf ${out2_ps} ${out2_pdf}
+
rm $t0 $t1 $t2z $t2r $t2t $t3
Modified: seismo/3D/automeasure/scripts/matlab_socal/window_subset.m
===================================================================
--- seismo/3D/automeasure/scripts/matlab_socal/window_subset.m 2008-05-27 15:09:17 UTC (rev 12033)
+++ seismo/3D/automeasure/scripts/matlab_socal/window_subset.m 2008-05-27 16:53:22 UTC (rev 12034)
@@ -32,8 +32,10 @@
win_mat = [Zwin Rwin Twin winall];
% read in list of event IDs
-file_eid = [dir0 'EIDs_fail'];
+%file_eid = [dir0 'EIDs_fail'];
%file_eid = [dir0 'EIDs_pass'];
+file_eid = '/net/sierra/raid1/carltape/results/SOURCES/socal_5/EIDs_all_loc_eid';
+
eids = textread(file_eid,'%s');
neid = length(eids);
@@ -51,7 +53,7 @@
disp([' writing ' ofile]);
% output file
-Wmin = 10;
+Wmin = 500; % KEY COMMAND
ofile = [file_eid '_window_' num2str(Wmin) '_plus'];
fid = fopen(ofile,'w');
for ii = 1:neid
Modified: seismo/3D/automeasure/scripts/pick_all_windows.pl
===================================================================
--- seismo/3D/automeasure/scripts/pick_all_windows.pl 2008-05-27 15:09:17 UTC (rev 12033)
+++ seismo/3D/automeasure/scripts/pick_all_windows.pl 2008-05-27 16:53:22 UTC (rev 12034)
@@ -37,19 +37,22 @@
# prepare_meas_all.pl
#
# EXAMPLES:
-# pick_all_windows.pl 6/40 179/179 1/1/0/0 1 0 # plots only
-# pick_all_windows.pl 6/40 179/179 0/0/1/0 1 0 # measurement code only
+# pick_all_windows.pl m0 6/40 179/179 1/1/0/0 1 0 # plots only
+# pick_all_windows.pl m0 6/40 179/179 0/0/1/0 1 0 # measurement code only
#
-# pick_all_windows.pl 6/40 1/203 1/1/1/0 1 0 # both, T = 6s
-# pick_all_windows.pl 2/40 1/203 1/1/1/1 1 0 # both, T = 2s
+# pick_all_windows.pl m0 6/40 1/203 1/1/1/0 1 0 # both, T = 6s
+# pick_all_windows.pl m0 2/40 1/203 1/1/1/1 1 0 # both, T = 2s
#
-# pick_all_windows.pl 6/40 54/54 1/1/1/0 1 0 # both (54=14095540 : body wave mode for 6s)
+# pick_all_windows.pl m0 6/40 54/54 1/1/1/0 1 0 # both (54=14095540 : body wave mode for 6s)
#
#==========================================================
-if (@ARGV < 5) {die("Usage: pick_all_windows.pl Tmin/Tmax imin/imax idebug/iplot/imeas/ibody idataset iexecute\n");}
-($Ts,$inds,$ibools,$idataset,$iexecute) = @ARGV;
+if (@ARGV < 6) {die("Usage: pick_all_windows.pl model Tmin/Tmax imin/imax idebug/iplot/imeas/ibody idataset iexecute\n");}
+($smodel,$Ts,$inds,$ibools,$idataset,$iexecute) = @ARGV;
+# plot a figure for every record -- this takes a long time (iplot = 1 for this option)
+$iplotall = 0;
+
# split input entries
($Tmin,$Tmax) = split("/",$Ts);
($imin,$imax) = split("/",$inds);
@@ -71,8 +74,8 @@
if($idataset == 1) {
$dir_data = "/net/sierra/raid1/carltape/socal/socal_3D/DATA/FINAL"; # Socal data (Carl)
- $dir_syn = "/net/sierra/raid1/carltape/socal/socal_3D/SYN/model_m0"; # Socal syn
- #$dir_syn = "/net/sierra/raid1/carltape/socal/socal_3D/SYN/model_pre_m0"; # Socal syn
+ $dir_syn = "/net/sierra/raid1/carltape/socal/socal_3D/SYN/model_${smodel}"; # Socal syn
+ #$dir_syn = "/net/sierra/raid1/carltape/socal/socal_3D/SYN/model_pre_${smodel}"; # Socal syn
$sdataset = "socal";
} elsif ($idataset == 2) {
@@ -103,8 +106,8 @@
$dir_user = "${dir_win_code}/user_files";
#$odir = "/home/carltape/UTILS/test";
-#$odir = "/net/sierra/raid1/carltape/results/WINDOWS/model_pre_m0";
-$odir = "/net/sierra/raid1/carltape/results/WINDOWS/model_m0";
+#$odir = "/net/sierra/raid1/carltape/results/WINDOWS/model_pre_${smodel}";
+$odir = "/net/sierra/raid1/carltape/results/WINDOWS/model_${smodel}";
# run directory for windows, measurements, adjoint sources, and kernels
$rundir = "/net/sierra/raid1/carltape/socal/socal_3D/RUNS";
@@ -178,8 +181,9 @@
# copy PAR_FILE into local directory
# directories for data and synthetics
if($idataset == 1) {
- $par_file_in = "${dir_user}/socal_3D/PAR_FILE_T${sTmin}";
- $userfun_file_in = "${dir_user}/socal_3D/user_functions.f90";
+ #$par_file_in = "${dir_user}/socal_3D/PAR_FILE_T06_14095540";
+ $par_file_in = "${dir_user}/socal_3D/PAR_FILE_T${sTmin}_${smodel}";
+ $userfun_file_in = "${dir_user}/socal_3D/user_functions_${smodel}.f90";
} elsif ($idataset == 2) {
$par_file_in = "${dir_user}/japan_3D/PAR_FILE_T${sTmin}";
@@ -189,8 +193,8 @@
if (not -f ${userfun_file_in}) {die("check if PAR_FILE ${userfun_file_in} exist or not\n")}
# copy user files to code directory
-`cp $par_file_in $par_file`;
-`cp $userfun_file_in $userfun_file`;
+`cp ${par_file_in} ${par_file}`;
+`cp ${userfun_file_in} ${userfun_file}`;
# name of executable windowing code (in $dir_win_code)
$win_execute = "test2";
@@ -292,7 +296,7 @@
print CSH "echo Event ID is $eid -- $kk out of $nevent\n";
$eout = "${eid}_T${sTmin}"; # output file tag
- $measdir = "$rundir/$eid/m0/WINDOW_T${sTmin}"; # main run directory
+ $measdir = "$rundir/$eid/$smodel/WINDOW_T${sTmin}"; # 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
@@ -330,10 +334,10 @@
# copy prepared data and synthetic files into a directory for measurement code
if ($imeas == 1) {
- #$measdir = "$rundir/$eid/m0/WINDOW_T${sTmin}";
+ #$measdir = "$rundir/$eid/$smodel/WINDOW_T${sTmin}";
print CSH "echo copying prepared DATA and SYN to ${measdir}...\n";
print CSH "mkdir -p $rundir/$eid\n";
- print CSH "mkdir -p $rundir/$eid/m0/\n";
+ print CSH "mkdir -p $rundir/$eid/$smodel/\n";
print CSH "mkdir -p $measdir\n";
print CSH "\\rm -rf $measdir/SYN ; cp -r ${dir_win_run_syn} $measdir\n";
print CSH "\\rm -rf $measdir/DATA ; cp -r ${dir_win_run_data} $measdir\n";
@@ -367,8 +371,10 @@
print CSH "\\cp ${reclist_dist} ${reclist_out}\n";
# generate composite PDF file using plot_windows_all.pl
- $pdffile = "$odir/${eout}_all.pdf";
- print CSH "${plot_windows_perl} ${dir_win_code} ${dir_win_run_meas} ${reclist_dist} $pdffile\n";
+ if($iplotall == 1) {
+ $pdffile = "$odir/${eout}_all.pdf";
+ print CSH "${plot_windows_perl} ${dir_win_code} ${dir_win_run_meas} ${reclist_dist} $pdffile\n";
+ }
# copy record section to output directory
$rs_in = "${dir_win_run_meas}/event_recordsection.pdf";
Modified: seismo/3D/automeasure/scripts/plot_windows_all.pl
===================================================================
--- seismo/3D/automeasure/scripts/plot_windows_all.pl 2008-05-27 15:09:17 UTC (rev 12033)
+++ seismo/3D/automeasure/scripts/plot_windows_all.pl 2008-05-27 16:53:22 UTC (rev 12034)
@@ -40,7 +40,7 @@
$dir_scripts = "${dir_win}/scripts";
if (not -e ${dir_scripts}) {die("check if ${dir_scripts} exist or not\n")}
- at comps = ("BHZ","BHR","BHT");
+ at comps = ("BHZ","BHR","BHT","HHZ","HHR","HHT");
$ncomp = @comps;
#$ncomp = 2; # PSV only
@@ -92,13 +92,14 @@
#======================
# concatenate the PDF files into one document
-# convert overal stats figures to pdfs
+# convert overall stats figures to pdfs
+# THIS IS NOW DONE IN extract_event_windowing_stats_carl.sh
$allstats = "${dir_win_run_meas}/event_winstats";
$recordsec = "${dir_win_run_meas}/event_recordsection";
-`cp ${allstats}.eps ${allstats}.ps`;
-`cp ${recordsec}.eps ${recordsec}.ps`;
-`ps2pdf ${allstats}.ps ${allstats}.pdf`;
-`ps2pdf ${recordsec}.ps ${recordsec}.pdf`;
+#`cp ${allstats}.eps ${allstats}.ps`;
+#`cp ${recordsec}.eps ${recordsec}.ps`;
+#`ps2pdf ${allstats}.ps ${allstats}.pdf`;
+#`ps2pdf ${recordsec}.ps ${recordsec}.pdf`;
# concatenate pdf files
$k = 1;
Deleted: seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T02
===================================================================
--- seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T02 2008-05-27 15:09:17 UTC (rev 12033)
+++ seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T02 2008-05-27 16:53:22 UTC (rev 12034)
@@ -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 = .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 = 4.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
Added: seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T02_m0
===================================================================
--- seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T02_m0 (rev 0)
+++ seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T02_m0 2008-05-27 16:53:22 UTC (rev 12034)
@@ -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 = .true.
+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 = 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
Added: seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T02_m1
===================================================================
--- seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T02_m1 (rev 0)
+++ seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T02_m1 2008-05-27 16:53:22 UTC (rev 12034)
@@ -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 = .true.
+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 = 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
Added: seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T02_m2
===================================================================
--- seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T02_m2 (rev 0)
+++ seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T02_m2 2008-05-27 16:53:22 UTC (rev 12034)
@@ -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 = .true.
+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 = 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/automeasure/user_files/socal_3D/PAR_FILE_T06
===================================================================
--- seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06 2008-05-27 15:09:17 UTC (rev 12033)
+++ seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06 2008-05-27 16:53:22 UTC (rev 12034)
@@ -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 (0.23)
-STALTA_BASE = 0.22
-
-# -------------------------------------------------------------
-# TSHIFT
-TSHIFT_BASE = 12.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
Modified: seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06_14095540
===================================================================
--- seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06_14095540 2008-05-27 15:09:17 UTC (rev 12033)
+++ seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06_14095540 2008-05-27 16:53:22 UTC (rev 12034)
@@ -29,7 +29,7 @@
# -------------------------------------------------------------
# TSHIFT
-TSHIFT_BASE = 2.0
+TSHIFT_BASE = 5.0
# -------------------------------------------------------------
# limit on CC for window acceptance
Added: seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06_m0
===================================================================
--- seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06_m0 (rev 0)
+++ seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06_m0 2008-05-27 16:53:22 UTC (rev 12034)
@@ -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 = .true.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+WIN_MIN_PERIOD = 6.00
+WIN_MAX_PERIOD = 40.00
+
+# -------------------------------------------------------------
+# E(t) water level (0.23)
+STALTA_BASE = 0.22
+
+# -------------------------------------------------------------
+# TSHIFT
+TSHIFT_BASE = 12.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
Added: seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06_m1
===================================================================
--- seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06_m1 (rev 0)
+++ seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06_m1 2008-05-27 16:53:22 UTC (rev 12034)
@@ -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 = .true.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+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
Added: seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06_m2
===================================================================
--- seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06_m2 (rev 0)
+++ seismo/3D/automeasure/user_files/socal_3D/PAR_FILE_T06_m2 2008-05-27 16:53:22 UTC (rev 12034)
@@ -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 = .true.
+BODY_WAVE_ONLY = .false.
+
+# -------------------------------------------------------------
+# period min/max for filtering
+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/automeasure/user_files/socal_3D/user_functions.f90
===================================================================
--- seismo/3D/automeasure/user_files/socal_3D/user_functions.f90 2008-05-27 15:09:17 UTC (rev 12033)
+++ seismo/3D/automeasure/user_files/socal_3D/user_functions.f90 2008-05-27 16:53:22 UTC (rev 12034)
@@ -1,133 +0,0 @@
-! -------------------------------------------------------------
-! 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.
- 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
- Sw_start = -15.0 + dist_km/3.5
- Sw_end = 35.0 + dist_km/3.1
-
- else
- Pnl_start = -5.0 + dist_km/7.8
- Sw_start = -15.0 + dist_km/3.5
- Sw_end = 35.0 + dist_km/3.1
- S_end = Sw_start
- endif
-
- ! variables for signal to noise ratio criteria.
- signal_end = Sw_end
- noise_end = Pnl_start
- if(DEBUG) then
- write(*,*) 'DEBUG : P_pick = ', sngl(P_pick)
- write(*,*) 'DEBUG : signal_end = ', sngl(sigmal_end)
- write(*,*) 'DEBUG : noise_end = ', sngl(noise_end)
- endif
-
- ! --------------------------------
- ! modulate criteria in time
- do i = 1, npts
- time = b+(i-1)*dt ! time
-
- ! 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.*STALTA_BASE
- endif
-
- ! raises STA/LTA water level after surface wave arrives
- if (BODY_WAVE_ONLY) then
- if(time.gt.S_end) then
- STALTA_W_LEVEL(i) = 10.*STALTA_BASE
- endif
- endif
-
- ! raises STA/LTA water level before P wave arrival.
- if(time.lt.Pnl_start) then
- STALTA_W_LEVEL(i) = 10.*STALTA_BASE
- 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/automeasure/user_files/socal_3D/user_functions_m0.f90
===================================================================
--- seismo/3D/automeasure/user_files/socal_3D/user_functions_m0.f90 (rev 0)
+++ seismo/3D/automeasure/user_files/socal_3D/user_functions_m0.f90 2008-05-27 16:53:22 UTC (rev 12034)
@@ -0,0 +1,135 @@
+! -------------------------------------------------------------
+! 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.
+ 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
+ Sw_start = -15.0 + dist_km/3.5
+ Sw_end = 35.0 + dist_km/3.1
+
+ else
+ Pnl_start = -5.0 + dist_km/7.8
+ Sw_start = -15.0 + dist_km/3.5
+ Sw_end = 35.0 + dist_km/3.1
+ S_end = Sw_start
+ endif
+
+ ! variables for signal to noise ratio criteria.
+ signal_end = Sw_end
+ noise_end = Pnl_start
+ if(DEBUG) then
+ write(*,*) 'DEBUG : P_pick = ', sngl(P_pick)
+ write(*,*) 'DEBUG : signal_end = ', sngl(sigmal_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
+ 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
+
+ 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/automeasure/user_files/socal_3D/user_functions_m1.f90
===================================================================
--- seismo/3D/automeasure/user_files/socal_3D/user_functions_m1.f90 (rev 0)
+++ seismo/3D/automeasure/user_files/socal_3D/user_functions_m1.f90 2008-05-27 16:53:22 UTC (rev 12034)
@@ -0,0 +1,135 @@
+! -------------------------------------------------------------
+! 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.
+ 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
+ Sw_start = -15.0 + dist_km/3.5
+ Sw_end = 35.0 + dist_km/3.1
+
+ else
+ Pnl_start = -5.0 + dist_km/7.8
+ Sw_start = -15.0 + dist_km/3.5
+ Sw_end = 35.0 + dist_km/3.1
+ S_end = Sw_start
+ endif
+
+ ! variables for signal to noise ratio criteria.
+ signal_end = Sw_end
+ noise_end = Pnl_start
+ if(DEBUG) then
+ write(*,*) 'DEBUG : P_pick = ', sngl(P_pick)
+ write(*,*) 'DEBUG : signal_end = ', sngl(sigmal_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
+ 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
+
+ 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/automeasure/user_files/socal_3D/user_functions_m2.f90
===================================================================
--- seismo/3D/automeasure/user_files/socal_3D/user_functions_m2.f90 (rev 0)
+++ seismo/3D/automeasure/user_files/socal_3D/user_functions_m2.f90 2008-05-27 16:53:22 UTC (rev 12034)
@@ -0,0 +1,135 @@
+! -------------------------------------------------------------
+! 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.
+ 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
+ Sw_start = -15.0 + dist_km/3.5
+ Sw_end = 35.0 + dist_km/3.1
+
+ else
+ Pnl_start = -5.0 + dist_km/7.8
+ Sw_start = -15.0 + dist_km/3.5
+ Sw_end = 35.0 + dist_km/3.1
+ S_end = Sw_start
+ endif
+
+ ! variables for signal to noise ratio criteria.
+ signal_end = Sw_end
+ noise_end = Pnl_start
+ if(DEBUG) then
+ write(*,*) 'DEBUG : P_pick = ', sngl(P_pick)
+ write(*,*) 'DEBUG : signal_end = ', sngl(sigmal_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
+ 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
+
+ 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