[cig-commits] r11760 - in seismo/3D/automeasure/scripts: matlab_socal prepare_scripts/socal

[carltape at geodynamics.org]

Author: carltape
Date: 2008-04-07 12:08:07 -0700 (Mon, 07 Apr 2008)
New Revision: 11760

Modified:
   seismo/3D/automeasure/scripts/matlab_socal/window_summary.m
   seismo/3D/automeasure/scripts/prepare_scripts/socal/prepare_seis.pl
Log:
Minor edits to run scripts.


Modified: seismo/3D/automeasure/scripts/matlab_socal/window_summary.m
===================================================================
--- seismo/3D/automeasure/scripts/matlab_socal/window_summary.m	2008-04-07 18:47:42 UTC (rev 11759)
+++ seismo/3D/automeasure/scripts/matlab_socal/window_summary.m	2008-04-07 19:08:07 UTC (rev 11760)
@@ -37,7 +37,7 @@
     index = 1;
     dir = [dir0 stper '_' sprintf('%2.2i',index) '/'];
 else
-    dir = '/net/sierra/raid1/carltape/results/WINDOWS/model_pre_m0/';
+    dir = '/net/sierra/raid1/carltape/results/WINDOWS/model_m0/';
     %dir = '/net/sierra/raid1/carltape/results/WINDOWS/final_6s/';
     %dir = '/net/wagholi/scratch1/dchao/automeasure_work/collect_long_biot/socal/';
     %dir = '/net/wagholi/scratch1/dchao/automeasure_work/collect_long/socal/';
@@ -46,13 +46,13 @@
 %--------------------------------------
 % read in list of event IDs
 
-dir_source = '/net/sierra/raid1/carltape/results/SOURCES/socal_4/';
-eid_file = [dir_source 'SOCAL_FINAL_CMT_v4_eid'];
+dir_source = '/net/sierra/raid1/carltape/results/SOURCES/socal_5/';
+eid_file = [dir_source 'SOCAL_FINAL_CMT_v5_eid'];
 eid_all = textread(eid_file,'%s','headerlines',0);
 nevent = length(eid_all);
 
 % read in CMT solutions
-cmt_file_all = [dir_source 'SOCAL_FINAL_CMT_v4'];
+cmt_file_all = [dir_source 'SOCAL_FINAL_CMT_v5'];
 [date,tshift,hdur,slat,slon,dep,M,eid,elabel] = readCMT(cmt_file_all, 13, 1);
 
 %--------------------------------------

Modified: seismo/3D/automeasure/scripts/prepare_scripts/socal/prepare_seis.pl
===================================================================
--- seismo/3D/automeasure/scripts/prepare_scripts/socal/prepare_seis.pl	2008-04-07 18:47:42 UTC (rev 11759)
+++ seismo/3D/automeasure/scripts/prepare_scripts/socal/prepare_seis.pl	2008-04-07 19:08:07 UTC (rev 11760)
@@ -102,10 +102,21 @@
     # e = earliest end time, multiplied by some factor
     if($bd < $bs) {$b0 = $bd} else {$b0 = $bs}
     if($ed < $es) {$e0 = $ed} else {$e0 = $es}
-    $b = $b0;
-    $tlen0 = $e0 - $b;
-    $tlen = $tfac * $tlen0;     # extend record length
-    $e = $b0 + $tlen;
+
+    if (1==1) {
+      $b = $b0;
+      $tlen0 = $e0 - $b;
+      $tlen = $tfac * $tlen0;	# extend record length
+      $e = $b0 + $tlen;
+
+    } else {
+      $e0 = 130;   # cut the records to avoid a later, larger event (14179288)
+      $b = $b0;
+      $tlen0 = $e0 - $b;
+      $tlen = $tlen0;
+      $e = $b0 + $tlen;
+    }
+
     $npt = int( ($e-$b)/$dt );
 
     if (0==1) {
@@ -130,7 +141,7 @@
 
 } 
 `echo quit >> sac.mac`;
-`sac sac.mac`;              # KEY: execute the SAC commands
+`sac sac.mac`;              # KEY: EXECUTE the SAC commands
 #`rm sac.mac`;
 
 print "\n Done with prepare_seis.pl\n\n";