Changeset 46 for SH_SHM/trunk/util

Timestamp:

10/30/2008 10:54:18 PM (15 years ago)

Author:

marcus

Message:

r31 | svn | 2008-03-03 20:47:24 +0100 (Mo, 03 Mär 2008) | 1 line

first working version for reading ingres sfdb; new check for origin in 1970; fixed ps output of locate_by_ttfit

Location:

SH_SHM/trunk/util

Files:

• check_evt.c (modified) (5 diffs)
• ev2view.c (modified) (1 diff)
• locate_by_ttfit.csh (modified) (3 diffs)

SH_SHM/trunk/util/check_evt.c

@@ -39 +39 @@
 #include BC_SYSBASE
 #include "eventdsc.h"
+#include "tcusrdef.h"
 
 
@@ -68 +69 @@
 static float   efv_magn=0.0;              /* max. magnitude found */
 static int     efv_pkp_count=0;           /* number of PKP phases */
+static BOOLEAN efv_illegal_origin;        /* illegal origin time */
 
 
@@ -99 +101 @@
         /* efv_infile is used for messages, directory will be removed */
         i = strlen( infile ) - 1;
-        while  (infile[i] != '/' && i > 0)  i--;
-                strcpy( efv_infile, infile+i+1 );
+        while  (infile[i] != '/' && i > 0)
+                i--;
+        strcpy( efv_infile, infile+i );
 
         status = EveNOERROR;
@@ -133 +136 @@
         BOOLEAN  is_p_phase;              /* phase is P or (P) */
         BOOLEAN  loc_given;               /* location given in this record */
+        NTIME    ntime;                   /* numeric origin time */
 
         /* executable code */
@@ -181 +185 @@
                 efv_uncorr_loc = TRUE;
                 strcpy( efv_uncorrphase, event->phase );
+        } /*endif*/
+
+        /* is there an origin time from 1970? */
+        if  (loc_given)  {
+                tc_t2n( event->origin_time, &ntime, status );
+                if  (SySevere(status))  {
+                        efv_illegal_origin = TRUE;
+                } else if  (ntime.year == 1970)  {
+                        efv_illegal_origin = TRUE;
+                } /*endif*/
         } /*endif*/
 

SH_SHM/trunk/util/ev2view.c

@@ -240 +240 @@
                 if  (event->mag[EvcMagMs] != EvEMPTY_MAGNITUDE)  {
                         fprintf( view, "  MS: %3.1f", event->mag[EvcMagMs] );
+                } /*endif*/
+                if  (event->mag[EvcMagMbb] != EvEMPTY_MAGNITUDE)  {
+                        fprintf( view, "  Mbb: %3.1f", event->mag[EvcMagMbb] );
                 } /*endif*/
         } /*endif*/

SH_SHM/trunk/util/locate_by_ttfit.csh

@@ -4 +4 @@
 #      ===================
 #
-# version 14, 4-Jun-2007
+# $Revision: 31 $, $Date: 2008-03-03 20:47:24 +0100 (Mo, 03 MÀr 2008) $
 #
 # Locate event via travel time fit.
@@ -52 +52 @@
 
 chdir $SH_SCRATCH
+# set units for plot, want to change this to cm
+gmtdefaults -Du >.gmtdefaults
 
 # delete previously existing result file
@@ -296 +298 @@
 #$PSVIEW -right $PWD/$psfile & # this dies after the shell terminates (on Sun)
 
-set xtmp=`whoami | cut -c1-6`
-
-if  ("$xtmp" == "tmpusr")  then
-
-        $PSVIEW $PWD/$psfile &
-
-else
-
-        set pos=""
-        if  ("`hostname`" == "ersn24")  set pos="-geometry 940x780+0+0"
-
-        rsh -n `hostname` "$PSVIEW -right -display $DISPLAY $pos $PWD/$psfile" &
-
-endif
+$PSVIEW $PWD/$psfile &
 
 #\rm $fitxyz