[cig-commits] [commit] pluggable: Removed phtimes (8374c62)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Apr 9 08:55:57 PDT 2014
Repository : ssh://geoshell/specfem3d_globe
On branch : pluggable
Link : https://github.com/geodynamics/specfem3d_globe/compare/64e1b38f0c5ebb4056cce0b15d41c0b9f94ab6e5...099a4d330d5b173b21e51ad441f9f429e5d37842
>---------------------------------------------------------------
commit 8374c627fa42d5f127bc2f4276da50c98268e12e
Author: Eh Tan <tan2 at earth.sinica.edu.tw>
Date: Thu May 7 20:04:50 2009 +0000
Removed phtimes
>---------------------------------------------------------------
8374c627fa42d5f127bc2f4276da50c98268e12e
UTILS/seis_process/process_data.pl.in | 28 ++++++++++++++--------------
UTILS/seis_process/process_syn.pl.in | 26 +++++++++++++-------------
2 files changed, 27 insertions(+), 27 deletions(-)
diff --git a/UTILS/seis_process/process_data.pl.in b/UTILS/seis_process/process_data.pl.in
index 3d4909a..d9be61d 100755
--- a/UTILS/seis_process/process_data.pl.in
+++ b/UTILS/seis_process/process_data.pl.in
@@ -16,7 +16,6 @@ sub Usage{
-l Start/End -- start and end of record from o
-t Ts/Tl -- shortest and longest period
-i -- transfer record to displacement
- -p -- pick event p and s first arrival into sac headers (t1 and t2)
-x Ext -- extension of new file name
-d Dir -- directory to put the output files (default .)
-s Sps -- sample per second (default 1.0)
@@ -28,7 +27,7 @@ sub Usage{
data_sac_files --- names of the data sac files to be processed
Notice:
- 0. Make sure you have sac, saclst, phtimes in PATH before execution
+ 0. Make sure you have sac, saclst in PATH before execution
1. We require that polezero files are in the same directory as the sac
data files which is generally satisfied. All the needed info are
taken from sac headers
@@ -70,7 +69,7 @@ if ($opt_l) {($lmin,$lmax) = split(/\//,$opt_l);} else {$lmin = 0; $lmax = 3600;
$sac="$sacaux/../bin/sac";
$saclst="$sacaux/../bin/saclst";
-$phtimes="phtimes";
+#$phtimes="phtimes";
if (! -e $sac) {die(" No $sac file\n");}
if (! -e $saclst) {die(" No $saclst file\n");}
#if (! -e $phtimes) {die("No $phtimes file\n");}
@@ -147,17 +146,18 @@ foreach $file (@ARGV) {
printf SAC ("trans from polezero s %20s to none freq %12.6f%12.6f%12.6f%12.6f \n",
$pzfile,$f0,$f1,$f2,$f3);}
- if ($opt_p) { # add p and s arrival info
- print " Add first P and S arrival information\n";
- (undef,$gcarc,$edep)=split(" ",`$saclst gcarc evdp f $outfile`);
- if (abs($gcarc-$undef) < $eps or abs($edep-$undef) < $eps) {die("No gcarc and depth info\n");}
- ($Pph,$Ptime)=split(" ",`$phtimes $edep $gcarc P`);
- ($Sph,$Stime)=split(" ",`$phtimes $edep $gcarc S`);
- if (not $opt_m) {$tshift = 0;}
- print SAC "evaluate to tmp1 $Ptime - $tshift\n";
- print SAC "evaluate to tmp2 $Stime - $tshift\n";
- print SAC "ch t1 %tmp1% t2 %tmp2%\n";
- print SAC "ch kt1 $Pph kt2 $Sph\n";}
+ ### phtimes is not included in this distribution
+ #if ($opt_p) { # add p and s arrival info
+ # print " Add first P and S arrival information\n";
+ # (undef,$gcarc,$edep)=split(" ",`$saclst gcarc evdp f $outfile`);
+ # if (abs($gcarc-$undef) < $eps or abs($edep-$undef) < $eps) {die("No gcarc and depth info\n");}
+ # ($Pph,$Ptime)=split(" ",`$phtimes $edep $gcarc P`);
+ # ($Sph,$Stime)=split(" ",`$phtimes $edep $gcarc S`);
+ # if (not $opt_m) {$tshift = 0;}
+ # print SAC "evaluate to tmp1 $Ptime - $tshift\n";
+ # print SAC "evaluate to tmp2 $Stime - $tshift\n";
+ # print SAC "ch t1 %tmp1% t2 %tmp2%\n";
+ # print SAC "ch kt1 $Pph kt2 $Sph\n";}
if ($opt_y) { # add arrival times to headers
@numbers=split(/\//,$opt_y); $npairs = floor(@numbers/2);
diff --git a/UTILS/seis_process/process_syn.pl.in b/UTILS/seis_process/process_syn.pl.in
index d221dc6..2984633 100755
--- a/UTILS/seis_process/process_syn.pl.in
+++ b/UTILS/seis_process/process_syn.pl.in
@@ -21,7 +21,6 @@ where
-P -- number of poles and passes in butterworth filter(default 4/2)
-i Dir -- convolve synthetics with instrument response in Dir
-A amp_factor -- multiple all files with an amplification factor
- -p -- pick P and S arrivals according to iasp91 model to t1/t2 headers
-a STAFILE -- add station information from STAFILE
-s sps -- resampling rate for synthetics (default is 1.0)
-y t0/V0/t1/V1 -- add to t3/t4 the arrival times of t + dist/V
@@ -31,7 +30,7 @@ where
names of files -- name of syn files to be processed
- Make sure you have sac, saclst, phtimes, asc2sac, ascii2sac.csh
+ Make sure you have sac, saclst, asc2sac
in PATH before execution
Qinya Liu, Caltech, May 2007
@@ -72,7 +71,7 @@ $eps = 0.1;
$sac="$sacaux/../bin/sac";
$saclst="$sacaux/../bin/saclst";
-$phtimes="phtimes";
+#$phtimes="phtimes";
$asc2sac="@prefix@/bin/asc2sac";
if (! -e $sac) {die(" No $sac file\n");}
if (! -e $saclst) {die(" No $saclst file\n");}
@@ -176,16 +175,17 @@ foreach $file (@ARGV) {
print SAC "transfer from none to polezero s ${pz}\n";
print SAC "w over\n";}
- if ($opt_p) { # pick arrival
- print "Picking arrivals...\n";
- (undef,$gcarc)=split(" ",`$saclst gcarc f $outfile`);
- if (abs($gcarc-$undef)< $eps) {die("No gcarc is defined\n");}
- ($Pph,$Ptime)=split(" ",`$phtimes $edep $gcarc P`);
- ($Sph,$Stime)=split(" ",`$phtimes $edep $gcarc S`);
- print SAC "evaluate to tmp1 $Ptime - $tshift\n";
- print SAC "evaluate to tmp2 $Stime - $tshift\n";
- print SAC "ch t1 %tmp1% t2 %tmp2%\n";
- print SAC "ch kt1 $Pph kt2 $Sph\n wh\n";}
+ ### phtimes is not included in this distribution
+ #if ($opt_p) { # pick arrival
+ # print "Picking arrivals...\n";
+ # (undef,$gcarc)=split(" ",`$saclst gcarc f $outfile`);
+ # if (abs($gcarc-$undef)< $eps) {die("No gcarc is defined\n");}
+ # ($Pph,$Ptime)=split(" ",`$phtimes $edep $gcarc P`);
+ # ($Sph,$Stime)=split(" ",`$phtimes $edep $gcarc S`);
+ # print SAC "evaluate to tmp1 $Ptime - $tshift\n";
+ # print SAC "evaluate to tmp2 $Stime - $tshift\n";
+ # print SAC "ch t1 %tmp1% t2 %tmp2%\n";
+ # print SAC "ch kt1 $Pph kt2 $Sph\n wh\n";}
if ($opt_y) { # add arrival time for surface waves
@numbers=split(/\//,$opt_y); $npairs = floor(@numbers/2);
More information about the CIG-COMMITS
mailing list