Changeset 57

Timestamp:

10/30/2008 11:04:21 PM (15 years ago)

Author:

marcus

Message:

r42 | svn | 2008-10-13 22:51:21 +0200 (Mo, 13 Okt 2008) | 1 line

implementation of :LOC expression as station button name

Location:

SH_SHM/trunk/source

Files:

• globalparams.c (modified) (3 diffs)
• globalparams.h (modified) (2 diffs)
• motif/mfexec.c (modified) (4 diffs)
• motif/shm.c (modified) (2 diffs)

SH_SHM/trunk/source/globalparams.c

@@ -3 +3 @@
  *      ==============
  *
- * $Revision: 29 $, $Date: 2008-02-23 23:02:22 +0100 (Sa, 23 Feb 2008) $
+ * $Revision: 42 $, $Date: 2008-10-13 22:51:21 +0200 (Mo, 13 Okt 2008) $
  *
  * Management of global parameters
@@ -395 +395 @@
         "reads_swap",
         "reads_invhdr",
+        "autosort_by_distance",
         "b_last_parameter"
 };
@@ -403 +404 @@
         (1 << cGpB_overwrite_string_lists) |
         (1 << cGpB_own_accelerators)       |
+        (1 << cGpB_autosort_by_distance)   |
         (1 << cGpB_reads_swap)
 ;

SH_SHM/trunk/source/globalparams.h

@@ -3 +3 @@
  *      ==============
  *
- * $Revision: 29 $, $Date: 2008-02-23 23:02:22 +0100 (Sa, 23 Feb 2008) $
+ * $Revision: 42 $, $Date: 2008-10-13 22:51:21 +0200 (Mo, 13 Okt 2008) $
  *
  * Management of global parameters
@@ -235 +235 @@
         cGpB_reads_swap,
         cGpB_reads_invhdr,
+        cGpB_autosort_by_distance,
         cGpB_last
 } TGpBooleanIndex;

SH_SHM/trunk/source/motif/mfexec.c

@@ -109 +109 @@
         float *azim, TSyStatus *status );
 static void mxh_read_stations_from_evt( char autoevt[], char addlist[] );
+static void mxh_read_location_from_evt( char autoevt[], float *lat, float *lon,
+        TSyBoolean *found );
 
 
@@ -143 +145 @@
         TSyStatus  locstat;                 /* local status */
         char     gsename[cBcFileLth+1];     /* name of GSE file */
+        float    epilat, epilon;            /* epicentre location rea from evt */
+        float    vlat1, vlat2, vlon1, vlon2;/* lat and lon area */
+        TSyBoolean epifound;                /* epicentre found in evt */
+        int      locvarcnt;                 /* location variable counter */
 
         /* executable code */
@@ -221 +227 @@
                                         } /*endif*/
                                         printf( "statlist: %s\n", statlist );
+                                } /*endif*/
+                        } else if  (strcmp(par->sl.code[i],":LOC") == 0)  {
+                                if  (autoevt == NULL)  {
+                                        printf( "*SHM: cannot translate :LOC, autoevt is NULL\n" );
+                                } else if  (*autoevt == '\0')  {
+                                        printf( "*SHM: cannot translate :LOC, autoevt is empty\n" );
+                                } else {
+                                        mxh_read_location_from_evt( autoevt, &epilat, &epilon,
+                                                &epifound );
+                                        if  (!epifound)  {
+                                                printf( "*SHM: cannot translate :LOC, no epicentre found\n" );
+                                                return;
+                                        } /*endif*/
+                                        if  (GpGetInt(cGpI_debug_level) > 2)
+                                                printf( "SHM-dbg3: :LOC found epi %f,%f\n", epilat, epilon );
+                                } /*endif*/
+                                epifound = FALSE;
+                                locvarcnt = 1;
+                                FOREVER {
+                                        /* read through all v$loc* variables */
+                                        sprintf( varname, "v$loc%d", locvarcnt );
+                                        GpReadParameter( varname, cBcLongStrLth, addlist, &ok );
+                                        if  (!ok)  break;
+                                        if  (GpGetInt(cGpI_debug_level) > 4)
+                                                printf( "SHM-dbg4: %s translated to %s\n", varname, addlist );
+                                        if  (sscanf(addlist,"%f,%f,%f,%f %s",
+                                                &vlat1,&vlon1,&vlat2,&vlon2,statlist) != 5)
+                                                continue;
+                                        if  (vlat1 <= epilat && epilat <= vlat2 && vlon1 <= epilon
+                                                && epilon <= vlon2)  {
+                                                epifound = TRUE;
+                                                break;
+                                        } /*endif*/
+                                        locvarcnt++;
+                                } /*endfor*/
+                                if  (!epifound)  {
+                                        GpReadParameter("v$loc", cBcLongStrLth, statlist, &ok );
+                                        if  (!ok)  {
+                                                printf( "*SHM: no v$loc definition\n" );
+                                                return;
+                                        } /*endif*/
                                 } /*endif*/
                         } else {
@@ -2141 +2188 @@
 
 /*----------------------------------------------------------------------------*/
+
+
+
+static void mxh_read_location_from_evt( char autoevt[], float *lat, float *lon,
+        TSyBoolean *found )
+
+/* Reads latitude and longitude from evt file
+ *
+ * parameters of routine
+ * char       autoevt[]; input; name of evt file
+ * float      *lat, *lon; output; location read from above file
+ * TSyBoolean *found; output; TRUE if location is found
+ */
+{
+        /* local variables */
+        FILE     *fp;                  /* pointer to file */
+        char     line[cBcLineLth+1];   /* current line of file */
+        TSyBoolean latfound, lonfound; /* latitude and longitude found */
+        int      i;                    /* counter */
+
+        /* executable code */
+
+        /* open evt file */
+        *found = TRUE;
+        fp = fopen( autoevt, "r" );
+        if  (fp == NULL)  {
+                *found = FALSE;
+                return;
+        } /*endif*/
+
+        /* read through file */
+        latfound = lonfound = FALSE;
+        while (fgets(line,cBcLineLth,fp) != NULL)  {
+                if  (strncmp(line,"Latitude ",9) == 0)  {
+                        i = 9;  /* length of Latitude string */
+                        while (line[i] == ' ' || line[i] == ':')
+                                i++;
+                        if  (sscanf(line+i,"%f",lat) == 1)
+                                latfound = TRUE;
+                } /*endif*/
+                if  (strncmp(line,"Longitude ",10) == 0)  {
+                        i = 10;  /* length of Latitude string */
+                        while (line[i] == ' ' || line[i] == ':')
+                                i++;
+                        if  (sscanf(line+i,"%f",lon) == 1)
+                                lonfound = TRUE;
+                } /*endif*/
+                if  (latfound && lonfound)
+                        break;
+        } /*endwhile*/
+
+        if  (!latfound || !lonfound)  {
+                *found = FALSE;
+                return;
+        } /*endif*/
+
+        fclose( fp );
+
+} /* end of mxh_read_location_from_evt */
+
+
+
+/*----------------------------------------------------------------------------*/

SH_SHM/trunk/source/motif/shm.c

@@ -3 +3 @@
  *      =====
  *
- * $Revision: 33 $, $Date: 2008-05-19 20:03:26 +0200 (Mo, 19 Mai 2008) $
+ * $Revision: 42 $, $Date: 2008-10-13 22:51:21 +0200 (Mo, 13 Okt 2008) $
  *
  * main module of shm
@@ -1551 +1551 @@
                                 /*mg_tracedisplay( xmv_w[k_widget_draw], &xmv_dspctrl, &status );*/
                         } /*endif*/
+                        /* sort by distance automatically if requested */
+                        if  (GpGetBoolean(cGpB_autosort_by_distance))  {
+                                if  (GpGetInt(cGpI_debug_level) > 3)
+                                        printf( "SHM-dbg4: auto sort by distance\n" );
+                                strcpy( xmv_sortstr, "DISTANCE" );
+                                strcpy( xmv_sortinfo, "sort:d" );
+                                /*cu_get_param_values( xmv_w, xmv_cpar );*/
+                                cl4_sort_traces( xmv_w[k_widget_draw], &xmv_dspctrl, xmv_sortstr,
+                                        xmv_cpar, &status );
+                                if  (Severe(&status))  {
+                                        cu_alert( status );
+                                        status = cBcNoError;
+                                } /*endif*/
+                                mg_print_sortinfo( xmv_w[k_widget_draw], xmv_sortinfo );
+                        } /*endif*/
                 } /*endif*/
                 cl2_display_busy( xmv_w, FALSE );