Changeset 64 for SH_SHM/trunk/source/motif/mfexec.c

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

File:

• SH_SHM/trunk/source/motif/mfexec.c (modified) (3 diffs)

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 */