Changeset 87 for SH_SHM/trunk/util

Timestamp:

12/10/2008 08:46:37 AM (14 years ago)

Author:

marcus

Message:

r69 | svn | 2008-12-07 15:30:25 +0100 (So, 07 Dez 2008) | 1 line

more tests in check_evt

File:

• SH_SHM/trunk/util/check_evt.c (modified) (10 diffs)

SH_SHM/trunk/util/check_evt.c

@@ -40 +40 @@
 #include "eventdsc.h"
 #include "tcusrdef.h"
+#include "earthloc.h"
 
 
@@ -50 +51 @@
 #define EVTYPE_PKP   3
 
-/* next error number currently 024 */
+/* next error number currently 031 */
 
 
 
 /* global variables */
+static float   efv_reflat=49.0;           /* reference latitude */
+static float   efv_reflon=11.0;           /* reference longitude */
 static char    efv_infile[BC_FILELTH+1];  /* name of input file */
 static char    efv_uncorrphase[EvPHASELTH+1]; /* name of phase on uncorr beam */
 static BOOLEAN efv_first_arrival=FALSE;   /* first arrival phase selected */
 static BOOLEAN efv_only_grf=TRUE;         /* only GRF stations used */
+static BOOLEAN efv_gra1_picked=FALSE;     /* pick on GRA1 found */
 static BOOLEAN efv_uncorr_loc=FALSE;      /* uncorrected beam for location */
 static BOOLEAN efv_telex_phase=FALSE;     /* telex phase found */
@@ -70 +74 @@
 static int     efv_pkp_count=0;           /* number of PKP phases */
 static BOOLEAN efv_illegal_origin;        /* illegal origin time */
-
+static BOOLEAN efv_foundPg=FALSE;         /* Pg phase found */
+static BOOLEAN efv_foundPnSn=FALSE;       /* Pn/Sn found */
+static BOOLEAN efv_foundTele=FALSE;       /* teleseismic phase found */
+static float   efv_latitude=0.0;          /* epicentre latitude */
+static float   efv_longitude=0.0;         /* epicentre longitude */
+static char    efv_analyst[EvANALYSTLTH+1]=""; /* analysts initials */
+static BOOLEAN efv_showinitials=FALSE;    /* show initials */
 
 
@@ -91 +101 @@
         /* executable code */
 
-        if  (argc != 3)  {
-                fprintf( stderr, "Usage: %s <input> <output>\n", argv[0] );
+        if  (argc < 2)  {
+                fprintf( stderr, "Usage: %s <input> [<output>] [<showinitials>]\n",
+                        argv[0] );
                 return 1;
         } /*endif*/
 
         strcpy( infile, argv[1] );
-        strcpy( outfile, argv[2] );
+        if  (argc > 2)
+                strcpy( outfile, argv[2] );
+                if  (*outfile == '\0')
+                        strcpy( outfile, "TT" );
+        else
+                strcpy( outfile, "TT" );
+        efv_showinitials = FALSE;
+        if  (argc > 3)
+                efv_showinitials = (strcmp(argv[3],"showinitials") == 0);
 
         /* efv_infile is used for messages, directory will be removed */
@@ -151 +170 @@
         /* is a location given here ? */
         loc_given = (event->latitude > 0.0 && event->longitude > 0.0);
+        if  (loc_given)  {
+                efv_latitude = event->latitude;
+                efv_longitude = event->longitude;
+        } /*endif*/
 
         /* is this an L phase ? */
@@ -162 +185 @@
         /* is this a PKP phase ? */
         is_pkp_phase = (strstr(event->phase,"PKP") != NULL);
+        if  (is_pkp_phase)
+                efv_foundTele = TRUE;
+
+        /* a depth phase can found only in teleseismic events */
+        if  (event->phase[0] == 'p' || event->phase[0] == 's')
+                efv_foundTele = TRUE;
+
+        /* check for local and regional phases */
+        if  (strcmp(event->phase,"Pg") == 0)
+                efv_foundPg = TRUE;
+        if  (strcmp(event->phase,"Pn") == 0)
+                efv_foundPnSn = TRUE;
+        if  (strcmp(event->phase,"Sn") == 0)
+                efv_foundPnSn = TRUE;
 
         if  (event->event_type != EvcEventTypeUndefined)
@@ -170 +207 @@
 
         /* is it a GRF station ? */
-        if  (strncmp(event->station,"GR",2) != 0)
+        if  (strncmp(event->station,"GR",2) != 0
+                || event->station[2] < 'A' || event->station[2] > 'C')
                 efv_only_grf = FALSE;
+        if  (strcmp(event->station,"GRA1") == 0)
+                efv_gra1_picked = TRUE;
 
         /* is a first arrival phase there ? */
@@ -315 +355 @@
         } /*endif*/
 
+        if  (event->analyst[0] != '\0')
+                strcpy( efv_analyst, event->analyst );
+
 } /* end of EfPutCheck */
 
@@ -365 +408 @@
  */
 {
+        /* local variables */
+        char     expl[cBcLongStrLth+1];   /* explanation string */
+
         /* executable code */
 
         /* fprintf( out, "--> final check\n" ); */
+        if  (efv_showinitials)
+                sprintf( expl, "%s (%s)", efv_infile, efv_analyst );
+        else
+                strcpy( expl, efv_infile );
 
         /* a first arrival must be specified */
         if  (!efv_first_arrival)
-                fprintf( out, "F014 %s: no first arrival phase selected\n", efv_infile );
+                fprintf( out, "F014 %s: no first arrival phase selected\n", expl );
 
         /* beam with uncorrected slowness at GRF is at least suspicious */
         if  (efv_only_grf && efv_uncorr_loc)
                 fprintf( out, "F015 %s: GRF location with uncorrected beam at phase %s\n",
-                        efv_infile, efv_uncorrphase );
+                        expl, efv_uncorrphase );
+
+        /* only GRF picked but no GRA1? */
+        if  (efv_only_grf && !efv_gra1_picked)
+                fprintf( out, "F035 %s: only GRF picks but none at GRA1\n", expl );
 
         /* no telex phase */
         /* disabled 12-Nov-2008 K.S.
         if  (!efv_telex_phase)
-                fprintf( out, "F016 %s: no telex phase found\n", efv_infile );
+                fprintf( out, "F016 %s: no telex phase found\n", expl );
         */
 
         /* mb determined? */
         if  (efv_subtype == EVTYPE_P && !efv_mb_found)
-                fprintf( out, "F017 %s: no Magnitude mb determined for P phase\n",
-                        efv_infile );
+                fprintf( out, "F017 %s: no Magnitude mb determined for P phase\n", expl );
 
         /* Ms determined? */
@@ -393 +446 @@
                 !efv_ms_found)
                 fprintf( out, "F018 %s: no Magnitude Ms determined for magn %3.1f event\n",
-                        efv_infile, efv_magn );
+                        expl, efv_magn );
         if  (efv_subtype == EVTYPE_PKP && efv_depth <= 50.0 && efv_pkp_count > 3 &&
                 !efv_ms_found)
                 fprintf( out, "F019 %s: no Magnitude Ms determined for PKP event\n",
-                        efv_infile );
+                        expl );
         if  (efv_ms_found && efv_depth > 50.0)
                 fprintf( out, "F020 %s: Magnitude Ms determined for depth %5.1f event\n",
-                        efv_infile, efv_depth );
+                        expl, efv_depth );
 
         /* Mbb determined? */
         if  (efv_mb_found && efv_magn > 6.2 && !efv_mbb_found)
                 fprintf( out, "F021 %s: no Broadband Magnitude determined for magn %3.1f event\n",
-                        efv_infile, efv_magn );
+                        expl, efv_magn );
 
         /* depth ok? */
         if  (efv_evtype == EvcEventTypeMining && Abs(efv_depth-1.0) > 0.1)
                 fprintf( out, "F022 %s: depth %5.1f given for mining event\n",
-                        efv_infile, efv_depth );
+                        expl, efv_depth );
         if  (efv_evtype == EvcEventTypeBlast && Abs(efv_depth) > 0.1)
                 fprintf( out, "F023 %s: depth %5.1f given for blast\n",
-                        efv_infile, efv_depth );
+                        expl, efv_depth );
         if  (efv_evtype == EvcEventTypeTeleQuake && (efv_depth < 10.0))
                 fprintf( out, "F026 %s: depth %5.1f for teleseismic event\n",
-                        efv_infile, efv_depth );
+                        expl, efv_depth );
+        if  ((efv_evtype == EvcEventTypeLocalQuake
+                || efv_evtype == EvcEventTypeRegioQuake) && (efv_depth < 2.0))
+                fprintf( out, "F027 %s: depth %5.1f for local/regional event\n",
+                        expl, efv_depth );
+        if  (efv_evtype == EvcEventTypeLocalQuake && efv_depth > 50.0)
+                fprintf( out, "F036 %s: depth %5.1f given for local event\n",
+                        expl, efv_depth );
         if  (efv_evtype == EvcEventTypeUndefined)
-                fprintf( out, "F024 %s: no event type specified\n", efv_infile );
+                fprintf( out, "F024 %s: no event type specified\n", expl );
         if  (efv_depth < -99.0)
-                fprintf( out, "F025 %s: no depth specified\n", efv_infile );
+                fprintf( out, "F025 %s: no depth specified\n", expl );
+
+        /* proper event type? */
+        if  (efv_foundPg && efv_evtype == EvcEventTypeTeleQuake)
+                fprintf( out, "F028 %s: improper teleseismic event type: found Pg\n",
+                        expl );
+        if  (efv_foundPnSn && efv_evtype == EvcEventTypeTeleQuake)
+                fprintf( out, "F029 %s: improper teleseismic event type: found Pn/Sn\n",
+                        expl );
+        if  (efv_foundTele && efv_evtype == EvcEventTypeLocalQuake)
+                fprintf( out, "F030 %s: improper local event type: found tele phase\n",
+                        expl );
+        if  (efv_foundTele && efv_evtype == EvcEventTypeRegioQuake)
+                fprintf( out, "F031 %s: improper regional event type: found tele phase\n",
+                        expl );
+
+        /* check distance and event type */
+        if  (efv_latitude != 0.0 || efv_longitude != 0.0)  {
+                double dist, azim, bazim;   /* result of distance computation */
+                float fdist;    /* single precision distance */
+                mb_locdiff( (double)efv_reflat, (double)efv_reflon, (double)efv_latitude,
+                        (double)efv_longitude, &dist, &azim, &bazim );
+                fdist = dist;
+                if  (dist > 15.0 && efv_evtype == EvcEventTypeLocalQuake)
+                        fprintf( out, "F032 %s: improper local event type (distance %f deg)\n",
+                                expl, dist );
+                if  (dist > 30.0 && efv_evtype == EvcEventTypeRegioQuake)
+                        fprintf( out, "F033 %s: improper regional event type (distance %f deg)\n",
+                                expl, dist );
+                if  (dist < 25.0 && efv_evtype == EvcEventTypeTeleQuake)
+                        fprintf( out, "F034 %s: improper teleseismic event type (distance %f deg)\n",
+                                expl, dist );
+        } /*endif*/
 
 } /* end of EfFInalCheck */