Changeset 58


Ignore:
Timestamp:
30.10.2008 23:04:47 (11 years ago)
Author:
marcus
Message:

r43 | svn | 2008-10-14 11:58:13 +0200 (Di, 14 Okt 2008) | 1 line

improved :LOC support

Location:
SH_SHM/trunk/source/motif
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • SH_SHM/trunk/source/motif/mfexec.c

    r57 r58  
    111111static void mxh_read_location_from_evt( char autoevt[], float *lat, float *lon, 
    112112        TSyBoolean *found ); 
     113static void mxh_translate_button_name( char code[], char autoevt[], 
     114        char addlist[], TSyStatus *status ); 
    113115 
    114116 
     
    137139        char     line[cBcVeryLongStrLth+1]; /* command line */ 
    138140        char     statlist[cBcVeryLongStrLth+1]; /* station list */ 
    139         char     varname[MXC_STATSTR_LTH+2]; /* variable name */ 
    140141        char     addlist[cBcLongStrLth+1];  /* parsed station string */ 
    141142        int      i;                         /* bit count */ 
    142143        char     *eptr;                     /* pointer to char for getenv */ 
    143144        TSyBoolean redraw;                  /* dummy */ 
    144         TSyBoolean ok;                      /* return code */ 
    145145        TSyStatus  locstat;                 /* local status */ 
    146146        char     gsename[cBcFileLth+1];     /* name of GSE file */ 
    147         float    epilat, epilon;            /* epicentre location rea from evt */ 
    148         float    vlat1, vlat2, vlon1, vlon2;/* lat and lon area */ 
    149         TSyBoolean epifound;                /* epicentre found in evt */ 
    150         int      locvarcnt;                 /* location variable counter */ 
    151147 
    152148        /* executable code */ 
     
    196192                if  ((1<<i) & par->stations)  { 
    197193                        if  (*statlist != '\0')  strcat( statlist, "," ); 
    198                         if  (par->sl.code[i][0] == '$')  { 
    199                                 varname[0] = 'v'; 
    200                                 strcpy( varname+1, par->sl.code[i] ); 
    201                                 GpReadParameter( varname, cBcLongStrLth, addlist, &ok ); 
    202                                 if  (!ok)  { 
    203                                         *status = MXE_ILL_SUBSET; 
    204                                         err_setcontext( " ## variable: " ); 
    205                                         err_setcontext( varname ); 
    206                                         return; 
    207                                 } /*endif*/ 
    208                                 if  (strlen(statlist)+strlen(addlist) < cBcVeryLongStrLth)  { 
    209                                         strcat( statlist, addlist ); 
    210                                 } else { 
    211                                         *status = MXE_STROVFL; 
    212                                         return; 
    213                                 } /*endif*/ 
    214                         } else if  (strcmp(par->sl.code[i],":AE") == 0)  { 
    215                                 if  (autoevt == NULL)  { 
    216                                         printf( "*SHM: cannot translate :AE, autoevt is NULL\n" ); 
    217                                 } else if  (*autoevt == '\0')  { 
    218                                         printf( "*SHM: cannot translate :AE, autoevt is empty\n" ); 
    219                                 } else { 
    220                                         printf( "get statlist\n" ); 
    221                                         mxh_read_stations_from_evt( autoevt, addlist ); 
    222                                         if  (strlen(statlist)+strlen(addlist) < cBcVeryLongStrLth)  { 
    223                                                 strcat( statlist, addlist ); 
    224                                         } else { 
    225                                                 *status = MXE_STROVFL; 
    226                                                 return; 
    227                                         } /*endif*/ 
    228                                         printf( "statlist: %s\n", statlist ); 
    229                                 } /*endif*/ 
    230                         } else if  (strcmp(par->sl.code[i],":LOC") == 0)  { 
    231                                 if  (autoevt == NULL)  { 
    232                                         printf( "*SHM: cannot translate :LOC, autoevt is NULL\n" ); 
    233                                 } else if  (*autoevt == '\0')  { 
    234                                         printf( "*SHM: cannot translate :LOC, autoevt is empty\n" ); 
    235                                 } else { 
    236                                         mxh_read_location_from_evt( autoevt, &epilat, &epilon, 
    237                                                 &epifound ); 
    238                                         if  (!epifound)  { 
    239                                                 printf( "*SHM: cannot translate :LOC, no epicentre found\n" ); 
    240                                                 return; 
    241                                         } /*endif*/ 
    242                                         if  (GpGetInt(cGpI_debug_level) > 2) 
    243                                                 printf( "SHM-dbg3: :LOC found epi %f,%f\n", epilat, epilon ); 
    244                                 } /*endif*/ 
    245                                 epifound = FALSE; 
    246                                 locvarcnt = 1; 
    247                                 FOREVER { 
    248                                         /* read through all v$loc* variables */ 
    249                                         sprintf( varname, "v$loc%d", locvarcnt ); 
    250                                         GpReadParameter( varname, cBcLongStrLth, addlist, &ok ); 
    251                                         if  (!ok)  break; 
    252                                         if  (GpGetInt(cGpI_debug_level) > 4) 
    253                                                 printf( "SHM-dbg4: %s translated to %s\n", varname, addlist ); 
    254                                         if  (sscanf(addlist,"%f,%f,%f,%f %s", 
    255                                                 &vlat1,&vlon1,&vlat2,&vlon2,statlist) != 5) 
    256                                                 continue; 
    257                                         if  (vlat1 <= epilat && epilat <= vlat2 && vlon1 <= epilon 
    258                                                 && epilon <= vlon2)  { 
    259                                                 epifound = TRUE; 
    260                                                 break; 
    261                                         } /*endif*/ 
    262                                         locvarcnt++; 
    263                                 } /*endfor*/ 
    264                                 if  (!epifound)  { 
    265                                         GpReadParameter("v$loc", cBcLongStrLth, statlist, &ok ); 
    266                                         if  (!ok)  { 
    267                                                 printf( "*SHM: no v$loc definition\n" ); 
    268                                                 return; 
    269                                         } /*endif*/ 
    270                                 } /*endif*/ 
     194                        mxh_translate_button_name( par->sl.code[i], autoevt, addlist, status ); 
     195                        if  (strlen(statlist)+strlen(addlist) < cBcVeryLongStrLth)  { 
     196                                strcat( statlist, addlist ); 
    271197                        } else { 
    272                                 if  (strlen(statlist)+strlen(par->sl.code[i]) < cBcVeryLongStrLth) { 
    273                                         strcat( statlist, par->sl.code[i] ); 
    274                                 } else { 
    275                                         *status = MXE_STROVFL; 
    276                                         return; 
    277                                 } /*endif*/ 
     198                                *status = MXE_STROVFL; 
     199                                return; 
    278200                        } /*endif*/ 
    279201                } /*endif*/ 
     
    386308 
    387309} /* end of mx_readg */ 
     310 
     311 
     312 
     313/*--------------------------------------------------------------------------*/ 
     314 
     315 
     316static void mxh_translate_button_name( char code[], char autoevt[], 
     317        char addlist[], TSyStatus *status ) 
     318 
     319/* Translates button name into a list of stations, returned by addlist 
     320 * 
     321 * parameters of routine 
     322 * char       code[]; input; button name (station code or list abbrev) 
     323 * char       autoevt[]; input; name of evt file with automatic picks and loc 
     324 * char       addlist[]; output; resulting station list; max lth cBcLongStrLth 
     325 * TSyStatus  *status; output; return status 
     326 */ 
     327{ 
     328        /* local variables */ 
     329        char     varname[MXC_STATSTR_LTH+2]; /* variable name */ 
     330        TSyBoolean ok;                       /* worked ok? */ 
     331        float    epilat, epilon;            /* epicentre location rea from evt */ 
     332        float    vlat1, vlat2, vlon1, vlon2;/* lat and lon area */ 
     333        TSyBoolean epifound;                /* epicentre found in evt */ 
     334        int      locvarcnt;                 /* location variable counter */ 
     335        char     vstr[cBcLongStrLth+1];     /* scratch for variable evaluation */ 
     336 
     337        /* executable code */ 
     338 
     339        *addlist = '\0'; 
     340        if  (code[0] == '$')  { 
     341                /* find list of stations in config file */ 
     342                varname[0] = 'v'; 
     343                strcpy( varname+1, code ); 
     344                GpReadParameter( varname, cBcLongStrLth, addlist, &ok ); 
     345                if  (!ok)  { 
     346                        *status = MXE_ILL_SUBSET; 
     347                        err_setcontext( " ## variable: " ); 
     348                        err_setcontext( varname ); 
     349                        *addlist = '\0'; 
     350                        return; 
     351                } /*endif*/ 
     352        } else if  (strcmp(code,":AE") == 0)  { 
     353                if  (autoevt == NULL)  { 
     354                        printf( "*SHM: cannot translate :AE, autoevt is NULL\n" ); 
     355                } else if  (*autoevt == '\0')  { 
     356                        printf( "*SHM: cannot translate :AE, autoevt is empty\n" ); 
     357                } else { 
     358                        mxh_read_stations_from_evt( autoevt, addlist ); 
     359                } /*endif*/ 
     360        } else if  (strcmp(code,":LOC") == 0)  { 
     361                if  (autoevt == NULL)  { 
     362                        /* if no autoevt file there, then use v$loc setting */ 
     363                        GpReadParameter( "v$loc", cBcLongStrLth, addlist, &ok ); 
     364                        if  (!ok)  { 
     365                                printf( "*SHM: cannot translate :LOC, autoevt is NULL and no v$loc found\n" ); 
     366                                *status = MXE_CONFIG_SETUP; 
     367                                err_setcontext( " ## variable: V$LOC not found" ); 
     368                                return; 
     369                        } /*endif*/ 
     370                } else if  (*autoevt == '\0')  { 
     371                        /* same as above */ 
     372                        GpReadParameter( "v$loc", cBcLongStrLth, addlist, &ok ); 
     373                        if  (!ok)  { 
     374                                printf( "*SHM: cannot translate :LOC, autoevt is NULL and no v$loc found\n" ); 
     375                                *status = MXE_CONFIG_SETUP; 
     376                                err_setcontext( " ## variable: V$LOC not found" ); 
     377                                return; 
     378                        } /*endif*/ 
     379                } /*endif*/ 
     380                /* get location from evt file */ 
     381                mxh_read_location_from_evt( autoevt, &epilat, &epilon, &epifound ); 
     382                if  (!epifound)  { 
     383                        /* no epicentre in evt file, behave like :AE */ 
     384                        printf( "*SHM: cannot translate :LOC, no epicentre found -> :AE\n" ); 
     385                        mxh_translate_button_name( ":AE", autoevt, addlist, status ); 
     386                        return; 
     387                } /*endif*/ 
     388                if  (GpGetInt(cGpI_debug_level) > 2) 
     389                        printf( "SHM-dbg3: :LOC found epi %f,%f\n", epilat, epilon ); 
     390                epifound = FALSE; 
     391                locvarcnt = 1; 
     392                FOREVER { 
     393                        /* read through all v$loc* variables and compare with curr. location */ 
     394                        sprintf( varname, "v$loc%d", locvarcnt ); 
     395                        GpReadParameter( varname, cBcLongStrLth, vstr, &ok ); 
     396                        if  (!ok)  break; 
     397                        if  (GpGetInt(cGpI_debug_level) > 3) 
     398                                printf( "SHM-dbg4: %s translated to %s\n", varname, vstr ); 
     399                        if  (sscanf(vstr,"%f,%f,%f,%f %s", 
     400                                &vlat1,&vlon1,&vlat2,&vlon2,addlist) != 5) 
     401                                continue; 
     402                        if  (vlat1 <= epilat && epilat <= vlat2 && vlon1 <= epilon 
     403                                && epilon <= vlon2)  { 
     404                                epifound = TRUE; 
     405                                break; 
     406                        } /*endif*/ 
     407                        locvarcnt++; 
     408                } /*endfor*/ 
     409                if  (!epifound)  { 
     410                        GpReadParameter("v$loc", cBcLongStrLth, addlist, &ok ); 
     411                        if  (!ok)  { 
     412                                printf( "*SHM: no v$loc definition\n" ); 
     413                                *status = MXE_CONFIG_SETUP; 
     414                                err_setcontext( " ## variable: V$LOC not found" ); 
     415                                return; 
     416                        } /*endif*/ 
     417                } /*endif*/ 
     418                /* if addlist is a variable, translate it */ 
     419                if  (addlist[0] == '$')  { 
     420                        strcpy( vstr, addlist ); 
     421                        GpReadParameter( vstr, cBcLongStrLth, addlist, &ok ); 
     422                        if  (!ok)  { 
     423                                printf( "*SHM: no %s definition\n", vstr ); 
     424                                *status = MXE_CONFIG_SETUP; 
     425                                err_setcontext( " ## variable: not found: " ); 
     426                                err_setcontext( vstr ); 
     427                                return; 
     428                        } /*endif*/ 
     429                } /*endif*/ 
     430                if  (GpGetInt(cGpI_debug_level) > 3) 
     431                        printf( "SHM-dbg4: :LOC addlist: %s\n", addlist ); 
     432        } else { 
     433                strcpy( addlist, code ); 
     434        } /*endif*/ 
     435 
     436} /* end of mxh_translate_button_name */ 
    388437 
    389438 
  • SH_SHM/trunk/source/motif/mfexec.h

    r16 r58  
    5050#define MXE_NOT_FOUND    (MXE_OFFSET+15)     /* not found */ 
    5151#define MXE_ILL_SUBSET   (MXE_OFFSET+16)     /* undefined station subset */ 
     52#define MXE_CONFIG_SETUP (MXE_OFFSET+17)     /* conig file setup error */ 
    5253 
    5354 
Note: See TracChangeset for help on using the changeset viewer.