Changeset 3169

Timestamp:

08/05/08 16:05:04 (5 years ago)

Author:

given

Message:

First checkin of hypophasedump

File:

• PP/trunk/perl_utils/hypophasedump.pl (copied) (copied from PP/trunk/perl_utils/hypophasedump) (4 diffs)

PP/trunk/perl_utils/hypophasedump.pl

@@ -7 +7 @@
 
 use DBI;            # Load the DBI module
-use lib "/home/tpp/bin/perlmodules";use DefNames;   #dbase access info
+use strict;
 
-  if ($#ARGV < 0) {
-   print "Syntax: $0 [<evid>] [<dbasename>]\n";
-   exit;
-  }
+# get masterdb and username/password
+use lib "$ENV{TPP_BIN_HOME}/perlmodules";
+use DbConn;
+use MasterDbs;
 
- if ($#ARGV >= 1) {
-   $masterdb = $ARGV[1];
- } else {
-   # lookup current masterdb
-   $masterdb = `/home/tpp/bin/masterdb`; chomp $masterdb;
- }
+# parse optional switches (N.B. this 'shift's the @ARGV array) 
+use Getopt::Std;
+use vars qw/ $opt_d $opt_h /;  #prevent "only used once" error
 
-$dbase = "dbi:Oracle:$masterdb";
-$evid  = $ARGV[0];
+getopts('d:h');  # "d:" means a value must follow
+if ($opt_h) {&usage(); exit;};
 
-  my $dbconn = DBI->connect( $dbase, $dbusername, $dbpassword, {PrintError => 0}  )
-         or die "Can't connect to Oracle database $masterdb: $DBI::errstr\n";
+my $dbase = $masterdb;
+if ($opt_d) { $dbase = $opt_d };  
+
+if ($#ARGV == -1) {&usage(); exit;};
+
+my $evid = $ARGV[0];
+if (!($evid =~ /^\d+$/) ) { print"<id> must be a number: $evid\n"; exit; } 
+
+my $dbconn = DbConn->new($dbase)->getConn();
 
 # get a list of candidate phases
-     $sql = $dbconn->prepare(" Select 
+    my $sql = $dbconn->prepare(" Select 
         net, sta, seedchan, datetime, AssocArO.iphase, Arrival.subsource, 
         fm, quality, Arrival.rflag, qual
@@ -42 +46 @@
 print "\n";
 
-while (@data = $sql->fetchrow_array()) {
+my ($sec,$min,$hour,$mday,$mon,$year,$frac);
+
+while (my @data = $sql->fetchrow_array()) {
     my $net       = $data[0] || ' ';
     my $sta       = $data[1] || ' ';
@@ -73 +79 @@
     }
 
-# Parse & compose time string
-my $frac = ($phasetime - int($phasetime) + 0.005) *100;
-($sec,$min,$hour,$mday,$mon,$year) = gmtime($phasetime);
-$timstr1 = sprintf "%4d%2.2d%2.2d%2.2d%2.2d" , $year+1900, $mon+1, $mday, $hour, $min;
-$timstr2 = sprintf "%2.2d%2.2d" ,$sec,$frac;
+    # Parse & compose time string
+    $frac = ($phasetime - int($phasetime) + 0.005) *100;
+    ($sec,$min,$hour,$mday,$mon,$year) = gmtime($phasetime);
+    my $timstr1 = sprintf "%4d%2.2d%2.2d%2.2d%2.2d" , $year+1900, $mon+1, $mday, $hour, $min;
+    my $timstr2 = sprintf "%2.2d%2.2d" ,$sec,$frac;
+    
+    my $str = sprintf ("%-5s%2s %1s%3s ",$sta, $net , substr($chan, 2,1), $chan );
+    
+    #         1         2         3         4         5         6         7         8         9        10        11
+    #12345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890
+    #STATNnt cCMP pr  yearmodyhrmnseconPresPwtS-secSr sSrespek2pekmmSwtPdlySdlyDistAngADPerRCodaAzmDurMagPimpSimpSFX
+    #MCV  NC VVHZ  PU0199806262007 4310  -9110    0   0   0      0 0  0   0   0  1616400  0  159297404  0 502   0WD 
+    #HMT  CI VVHZ     199811102059    0   0  0 5632IS 0 -27    151 0  0   0   0 708 9000  0    0228  0  0   0
+    
+    if ($phase eq "P" ) {
+     $str = $str.$eiw.$phase.$fm.$wt . $timstr1 ." " . $timstr2;
+    } else {   # S-wave
+     $str = $str ."   4" . $timstr1 ."             " . $timstr2 . $eiw. $phase . " " . $wt;
+    }
+    
+    print "$str\n";
 
-$str = sprintf ("%-5s%2s %1s%3s ",$sta, $net , substr($chan, 2,1), $chan );
+} # end of while loop
 
-#         1         2         3         4         5         6         7         8         9        10        11
-#12345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890
-#STATNnt cCMP pr  yearmodyhrmnseconPresPwtS-secSr sSrespek2pekmmSwtPdlySdlyDistAngADPerRCodaAzmDurMagPimpSimpSFX
-#MCV  NC VVHZ  PU0199806262007 4310  -9110    0   0   0      0 0  0   0   0  1616400  0  159297404  0 502   0WD 
-#HMT  CI VVHZ     199811102059    0   0  0 5632IS 0 -27    151 0  0   0   0 708 9000  0    0228  0  0   0
+# #########################################################################
+# Add the terminator line with ID#
+#                            
+#get the current orid
 
-if ($phase eq "P" ) {
- $str = $str.$eiw.$phase.$fm.$wt . $timstr1 ." " . $timstr2;
-} else {   # S-wave
- $str = $str ."   4" . $timstr1 ."             " . $timstr2 . $eiw. $phase . " " . $wt;
-}
+    $sql = $dbconn->prepare( "Select
+        Event.evid, Origin.datetime, Origin.lat, Origin.lon, Origin.depth,
+        Origin.fdepth, Origin.fepi
+         from Event, Origin, NetMag 
+         WHERE selectflag = 1 and 
+        (Event.prefor = Origin.orid(+))
+         and Event.evid = ? 
+       ");
+
+ $sql->execute($evid);
+
+ my @data = $sql->fetchrow_array();
+
+    my $evidx = $data[0] || 'null';
+    my $ot   = $data[1] || 0;
+    my $lat  = $data[2] || 0;
+    my $lon  = $data[3] || 0;
+    my $z    = $data[4] || 0;
+    my $fixz = $data[5] || 'n';        # y/n or null
+    my $fixloc = $data[6] || 'n';      # y/n
+
+# convert, parse & format the time string (doesn't use yr,mo,day)
+ $frac = ($ot - int($ot) + 0.005) *100;
+#($sec,$min,$hour,$mday,$mon,$year) = gmtime($ot);
+($sec,$min,$hour) = gmtime($ot);
+my $otstr = sprintf "%2.2d%2.2d%2.2d%2.2d" , $hour, $min ,$sec, $frac;
+
+# convert decimal lat/lon to Hypo phase format
+$lat = abs($lat);
+$lon = abs($lon);
+
+my $latdeg = int($lat);
+my $latmin = ($lat - $latdeg+0.00005) * 6000;    # * 600000 dec -> min & shift decimal
+my $londeg = int($lon);
+my $lonmin = ($lon - $londeg+0.00005) * 6000; 
+my $zx = ($z + 0.005) * 100;
+
+# handle fix flags (hypo can't do fixed origin time)
+my $fixFlag = " ";
+if ($fixloc eq "y") { $fixFlag = "X";}
+if ($fixz eq "y") { $zx = -$zx;}
+
+#         1         2         3         4         5         6         7
+#1234567890123456789012345678901234567890123456789012345678901234567890
+#    xxi2i2f4.2i2 f4.2i3  f4.2 f5.2                            i10
+#^^^^xxhrmnsecsLT minsLON minszzzzz                            IDnumberxx
+#      0218152533 3050117 4512  600                            9876543
+#      2102210036 0108117 4638 1.36X                           9082441
+#      2149378533 4929117 3098 00830                          10065513
+my $str = sprintf ("      %7s%2.2d %4.2d%3.3d %4.4d%5d%s                         %10d", 
+          $otstr, $latdeg, $latmin,  $londeg, $lonmin, $zx, $fixFlag, $evid);
 
 print "$str\n";
-
-}
 
 ### Disconnect from the database
@@ -101 +166 @@
 $dbconn->disconnect
   or warn "Disconnection failed: $DBI::errstr\n";
+  
+exit;
 
-#
-# Add the terminator line with ID#
-#                            
-print `/home/tpp/bin/hypoterminator ${evid} ${masterdb}`;
+# -----------------------------------------------------------------
+sub usage() {
+print <<"EOF";
+
+    Dump event phase info in hypoinverse format.
+    
+    usage: $0 -[h] [-d dbase] <evid>
+     
+     -h        : print usage info
+     -d dbase  : use this particular dbase (defaults to "MasterDB")
+    
+    example: $0 -d k2db  9277430
+EOF
+exit;
+}
+