[cig-commits] r14918 - seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/UTILS/seis_process
tan2 at geodynamics.org
tan2 at geodynamics.org
Thu May 7 13:04:50 PDT 2009
Author: tan2
Date: 2009-05-07 13:04:50 -0700 (Thu, 07 May 2009)
New Revision: 14918
Modified:
seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/UTILS/seis_process/process_data.pl.in
seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/UTILS/seis_process/process_syn.pl.in
Log:
Removed phtimes
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/UTILS/seis_process/process_data.pl.in
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/UTILS/seis_process/process_data.pl.in 2009-05-07 18:42:42 UTC (rev 14917)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/UTILS/seis_process/process_data.pl.in 2009-05-07 20:04:50 UTC (rev 14918)
@@ -16,7 +16,6 @@
-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 @@
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 @@
$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 @@
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);
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/UTILS/seis_process/process_syn.pl.in
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/UTILS/seis_process/process_syn.pl.in 2009-05-07 18:42:42 UTC (rev 14917)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/UTILS/seis_process/process_syn.pl.in 2009-05-07 20:04:50 UTC (rev 14918)
@@ -21,7 +21,6 @@
-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 @@
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 @@
$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 @@
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