[cig-commits] r17173 - seismo/3D/ADJOINT_TOMO/flexwin/scripts

[liuqy at geodynamics.org]

Author: liuqy
Date: 2010-09-06 09:21:27 -0700 (Mon, 06 Sep 2010)
New Revision: 17173

Modified:
   seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats.sh
   seismo/3D/ADJOINT_TOMO/flexwin/scripts/plot_flexwin.py
   seismo/3D/ADJOINT_TOMO/flexwin/scripts/plot_socal_ds.pl
   seismo/3D/ADJOINT_TOMO/flexwin/scripts/write_flexwin_out.py
Log:
Update scripts for processing and plotting flexwin outputs



Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats.sh
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats.sh	2010-09-06 05:20:44 UTC (rev 17172)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/extract_event_windowing_stats.sh	2010-09-06 16:21:27 UTC (rev 17173)
@@ -1,5 +1,9 @@
 #!/bin/bash
 
+if [ $# != 1 ]; then
+  echo "extract_event_windowing_stats.sh basename(MEASURE)"; exit
+fi
+
 basename=$1
 figname=window_stats.fig
 
@@ -165,5 +169,5 @@
 
 echo $out2
 
-rm $t0 $t1 $t2z $t2r $t2t $t3
+rm -f $t0 $t1 $t2z $t2r $t2t $t3
 

Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/plot_flexwin.py
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/plot_flexwin.py	2010-09-06 05:20:44 UTC (rev 17172)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/plot_flexwin.py	2010-09-06 16:21:27 UTC (rev 17173)
@@ -6,28 +6,28 @@
 import sys,os
 
 if (len(sys.argv) != 2):
-  print "plot_flexwin.py measure_dir(MEASURE)"; exit()
+  sys.exit('plot_flexwin.py measure_dir(MEASURE)')
 
 dir=sys.argv[1]
 if not os.path.isdir(dir):
-  print "no such dir "+dir; exit()
+  sys.exit("no such dir "+dir)
 script1='scripts/plot_seismos_gmt.sh'
 script2='scripts/extract_event_windowing_stats.sh'
 if not os.path.isfile(script1) or not os.path.isfile(script2):
-  print 'no '+script1+' or '+script2; exit()
+  sys.exit('no '+script1+' or '+script2)
 
 if (os.system(script2+' '+dir+' > /dev/null') != 0):
-  print 'Error executing '+ script2; exit()
+  sys.exit('Error executing '+ script2)
 
 ps=dir+'/event_winstats.pdf'+' '+dir+'/event_recordsection.pdf'
 
 for basename in os.popen('grep DDG '+dir+"/*.info | awk '{print $1,$4}'| sort -k 2 -g | awk '{print $1}' | sed 's/\.info:#//g'").readlines():
   input = basename.rstrip()  
   if (os.system(script1+' '+input) != 0):
-    print "Error plotting individual seismograms"; exit()
+    sys.exit("Error plotting individual seismograms")
   ps = ps + ' '+input+'.seis.pdf'
 
 if (os.system('pdcat -r '+ps+' flexwin_seis.pdf') != 0):
-  print "Error concatenate all seismogram plots"; exit()
+  sys.exit("Error concatenate all seismogram plots")
 
 

Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/plot_socal_ds.pl
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/plot_socal_ds.pl	2010-09-06 05:20:44 UTC (rev 17172)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/plot_socal_ds.pl	2010-09-06 16:21:27 UTC (rev 17173)
@@ -1,11 +1,9 @@
 #!/usr/bin/perl -w
 
-# this program plots data and syn waveforms according to either
+# this program plots sophisticated data and syn waveforms according to either
 # given data and syn dir
 # or flexwin output file 
 # sort by distance and then by azimuth, a regional SC map is also supplied
-# scaled by magnitude 5 earthquake
-# Qinya Liu, Oct 2009
 
 use lib '/home/liuqy/local/perl-lib';
 use GMT_PLOT;
@@ -16,7 +14,7 @@
 use List::Util qw[min max];
 
 
-if (@ARGV == 0) {die("Usage: plot_socal_ds.pl -d data_dir,data_ext -s syn_dir,syn_ext -Snew_syn_dir,new_syn_ext -M flexwin_file -m CMTSOLUTION -A Z/R/T(0.03/0.03/0.04)\n");}
+if (@ARGV == 0) {die("Usage: plot_data_and_syn.pl -d data_dir,data_ext -s syn_dir,syn_ext -Snew_syn_dir,new_syn_ext -M flexwin_file -m CMTSOLUTION -A Z/R/T(0.03/0.03/0.01)\n");}
 
 if (!getopts('d:s:S:m:M:A:')) {die("Check input options\n");}
 

Modified: seismo/3D/ADJOINT_TOMO/flexwin/scripts/write_flexwin_out.py
===================================================================
--- seismo/3D/ADJOINT_TOMO/flexwin/scripts/write_flexwin_out.py	2010-09-06 05:20:44 UTC (rev 17172)
+++ seismo/3D/ADJOINT_TOMO/flexwin/scripts/write_flexwin_out.py	2010-09-06 16:21:27 UTC (rev 17173)
@@ -12,24 +12,26 @@
 import os,sys,glob
 
 if (len(sys.argv) != 3 and len(sys.argv) != 4):
-  print "write_flexwin_out.py measure_dir out_filename [manual-pick-file]"; exit()
+  sys.exit("write_flexwin_out.py measure_dir out_filename [manual-pick-file]")
 
 dir=sys.argv[1]; outfile=sys.argv[2]
 if (not os.path.isdir(dir)):
-  print 'check if '+dir+' exists or not'; exit()
+  sys.exit('check if '+dir+' exists or not')
 
 manual_pick=False
 if (len(sys.argv) == 4):
   manual_pick=True
   manual_file=sys.argv[3]
   if (not os.path.isfile(manual_file)):
-    print 'Check if '+manual_file+' exists or not'; exit()
+    sys.exit('Check if '+manual_file+' exists or not')
   manual_lines=open(manual_file,'r').readlines()
   flexwin_input = manual_lines[0].rstrip('\n');
   print '*** Adding manual pick info from '+ manual_file+ ' ***'
   
 output=''
-files=glob.glob(dir+'/*mt*');
+files=glob.glob(dir+'/*mt*')
+if len(files) <= 0:
+  sys.exit('No such file as '+dir+'/*mt*')
 nfiles=0
 
 for file in files: