Changeset 64

Timestamp:

10/30/2008 11:07:23 PM (14 years ago)

Author:

marcus

Message:

r49 | svn | 2008-10-14 22:36:40 +0200 (Di, 14 Okt 2008) | 1 line

autoevt: filter depending on distance; appropriate variables in globalparams.h

Location:

SH_SHM/trunk

Files:

• inputs/shm-config.txt (modified) (1 diff)
• source/globalparams.c (modified) (7 diffs)
• source/globalparams.h (modified) (4 diffs)
• source/motif/mfexec.c (modified) (3 diffs)
• source/motif/mfexec.h (modified) (1 diff)
• source/motif/shm.c (modified) (2 diffs)

SH_SHM/trunk/inputs/shm-config.txt

@@ -384 +384 @@
 keysym_arrow_right           268762963
 
+!                            automatic evt files (used in detection lists)
+!                            ---------------------------------------------
+!                            when an autoevt contains a location, should the
+!                            trace automatically be sorted by distance?
+autosort_by_distance         true
+!                            translation of :LOC expression at station buttons:
+!                            if a autoevt file contains a location the station
+!                            set read in is determined from it: the location
+!                            is compared with the following list of lat,lon
+!                            areas specified in v$loc* variables (order:
+!                            <lat1>,<lon1>,<lat2>,<lon2>).  The first matching
+!                            area sets the station list.  The station list
+!                            itself may be a variable.  If no area matches
+!                            the list of v$loc is taken.
+!                            Vogtland area use station list v$e from above
+v$loc1                       50.0,11.0,52.0,13.0 $e
+!                            Saarbruecken area use station list v$sw
+v$loc2                       48.9,6.4,49.9,7.4 $sw
+!                            an arbitrary list of v$loc* lines may follow
+!                            the escape is list v$loc (if no area matches)
+v$loc                        $grsn
+!                            if an autoevt file contains a location, should the
+!                            default filter (on Read New) be determined from
+!                            the distance between the epicentre and the
+!                            reference station?
+autofilter_by_distance       true
+!                            if yes, what is the separation distance in degrees
+autofilter_sepdist           5.0
+!                            and what are the local and teleseismic filters,
+!                            respectively
+autofilter_local             SHM_BP_1HZ_8HZ_4
+autofilter_tele              S+G_WWSSN_SP
+
 !!                            autofilter definitions
 !v$autofilter                 cll-bh-z  SHM_HP_1HZ_3

SH_SHM/trunk/source/globalparams.c

@@ -3 +3 @@
  *      ==============
  *
- * $Revision: 42 $, $Date: 2008-10-13 22:51:21 +0200 (Mo, 13 Okt 2008) $
+ * $Revision: 49 $, $Date: 2008-10-14 22:36:40 +0200 (Di, 14 Okt 2008) $
  *
  * Management of global parameters
@@ -211 +211 @@
         "idphases_tol_azim",
         "axis_label_rounding",
+        "autofilter_sepdist",
         "f_last_parameter"
 };
@@ -275 +276 @@
         20.0,      /* idphases_tol_azim */
         0.01,      /* axis_label_rounding */
+        5.0,       /* autofilter_sepdist */
         0.0        /* last parameter */
 };
@@ -334 +336 @@
         "sfdb_exec_qual",
         "sfdb_trailer",
+        "autofilter_local",
+        "autofilter_tele",
         "s_last_parameter"
 };
@@ -373 +377 @@
         "-B -e",                         /* sfdb_exec_qual */
         "",                              /* sfdb_trailer */
+        "SHM_BP_1HZ_8HZ_4",              /* autofilter_local */
+        "S+G_WWSSN_SP",                  /* autofilter_tele */
         "s_last_parameter"               /* last parameter */
 };
@@ -396 +402 @@
         "reads_invhdr",
         "autosort_by_distance",
+        "autofilter_by_distance",
         "b_last_parameter"
 };
@@ -405 +412 @@
         (1 << cGpB_own_accelerators)       |
         (1 << cGpB_autosort_by_distance)   |
+        (1 << cGpB_autofilter_by_distance) |
         (1 << cGpB_reads_swap)
 ;

SH_SHM/trunk/source/globalparams.h

@@ -3 +3 @@
  *      ==============
  *
- * $Revision: 42 $, $Date: 2008-10-13 22:51:21 +0200 (Mo, 13 Okt 2008) $
+ * $Revision: 49 $, $Date: 2008-10-14 22:36:40 +0200 (Di, 14 Okt 2008) $
  *
  * Management of global parameters
@@ -159 +159 @@
         cGpF_idphases_tol_azim,
         cGpF_axis_label_rounding,
+        cGpF_autofilter_sepdist,
         cGpF_last
 } TGpFloatIndex;
@@ -213 +214 @@
         cGpS_sfdb_exec_qual,
         cGpS_sfdb_trailer,
+        cGpS_autofilter_local,
+        cGpS_autofilter_tele,
         cGpS_last
 } TGpStringIndex;
@@ -236 +239 @@
         cGpB_reads_invhdr,
         cGpB_autosort_by_distance,
+        cGpB_autofilter_by_distance,
         cGpB_last
 } TGpBooleanIndex;

SH_SHM/trunk/source/motif/mfexec.c

@@ -3 +3 @@
  *      ========
  *
- * version 63, 21-Nov-2007
+ * $Revision: 49 $, $Date: 2008-10-14 22:36:40 +0200 (Di, 14 Okt 2008) $
  *
  * execution of SH commands
@@ -57 +57 @@
 #include "shvars.h"
 #include "glusrdef.h"
+#include "earthloc.h"
 #include "phaseinf.h"
 #include "phasemgr.h"
@@ -436 +437 @@
 
 } /* end of mxh_translate_button_name */
+
+
+
+/*--------------------------------------------------------------------------*/
+
+
+
+void mx_filter_from_autoevt( char autoevt[], char filtername[], float *autocut )
+
+/* Determines filter using distance between autoevt location and reference
+ * station
+ *
+ * parameters of routine
+ * char autoevt[]; input;      name of autoevt file
+ * char filtername[]; output;  name of filter to use
+ * float *autocut; output;     cut off that many seconds after filtering
+ */
+{
+        /* local variables */
+        float    epilat, epilon;         /* location from autoevt */
+        TSyBoolean epifound;             /* epicentre found? */
+        GLT_STATINF *statinf;            /* station info of reference station */
+        TSyStatus locstat;               /* local status variable */
+        double dist, az, baz;            /* for distance computation */
+
+        /* executable code */
+
+        /* preset output vars in case of error returns */
+        *filtername = '\0';
+        *autocut = 5.0;
+
+        /* read the automatic location from the evt text file */
+        mxh_read_location_from_evt( autoevt, &epilat, &epilon, &epifound );
+        if  (!epifound)  return;
+
+        /* get the station location of the reference station */
+        locstat = cBcNoError;
+        statinf = gl_store_station( GpGetString(cGpS_refstation), TRUE, &locstat );
+        if  (SySevere(&locstat))  return;
+
+        /* distance of automatic ocation from reference station */
+        mb_locdiff( statinf->lat, statinf->lon, epilat, epilon, &dist, &az, &baz );
+
+        /* take one of two possible filters, depending on the distance */
+        if  (dist < GpGetFloat(cGpF_autofilter_sepdist))  {
+                strcpy( filtername, GpGetString(cGpS_autofilter_local) );
+        } else {
+                strcpy( filtername, GpGetString(cGpS_autofilter_tele) );
+        } /*endif*/
+
+} /* end of mx_filter_from_autoevt */
 
 

SH_SHM/trunk/source/motif/mfexec.h

@@ -332 +332 @@
 
 /*--------------------------------------------------------------------------*/
+
+
+void mx_filter_from_autoevt( char autoevt[], char filtername[],
+        float *autocut );
+
+/* Determines filter using distance between autoevt location and reference
+ * station
+ *
+ * parameters of routine
+ * char autoevt[]; input;      name of autoevt file
+ * char filtername[]; output;  name of filter to use
+ * float *autocut; output;     cut off that many seconds after filtering
+ */
+
+
+/*--------------------------------------------------------------------------*/

SH_SHM/trunk/source/motif/shm.c

@@ -3 +3 @@
  *      =====
  *
- * $Revision: 48 $, $Date: 2008-10-14 18:19:33 +0200 (Di, 14 Okt 2008) $
+ * $Revision: 49 $, $Date: 2008-10-14 22:36:40 +0200 (Di, 14 Okt 2008) $
  *
  * main module of shm
@@ -1444 +1444 @@
                 strcpy( xmv_phase.name, GpGetString(cGpS_auto_phase) );
                 cu_phase_box_defaults( xmv_w, &xmv_phase );
-                xmv_cmd_filter.autocut = 5.0;  /* shv_globals not yet implemented */
-                strcpy( xmv_cmd_filter.name, GpGetString(cGpS_default_filter) );
+                if  (TRUE)  {
+                        mx_filter_from_autoevt( autoevt, xmv_cmd_filter.name,
+                                &xmv_cmd_filter.autocut );
+                        if  (*xmv_cmd_filter.name == '\0')  {
+                                xmv_cmd_filter.autocut = 5.0;  /* shv_globals not yet implemented */
+                                strcpy( xmv_cmd_filter.name, GpGetString(cGpS_default_filter) );
+                        } /*endif*/
+                } else {
+                        xmv_cmd_filter.autocut = 5.0;  /* shv_globals not yet implemented */
+                        strcpy( xmv_cmd_filter.name, GpGetString(cGpS_default_filter) );
+                } /*endif*/
                 mx_filter_box_defaults( xmv_w, &xmv_cmd_filter );
                 xmv_dspctrl.zoom = 1.0;