Changeset 303

Timestamp:

01/11/2011 11:21:30 PM (13 years ago)

Author:

marcus

Message:

r145 | walther | 2011-01-11 23:18:02 +0100 (Di, 11 Jan 2011) | 12 lines

• Extended sfdline for output of additional information about network and location code (on by default). Setting the new option "-nonetloc" will reveal the old behaviour. This change should be transparent for any Seismic Handler installations, since this extra information is ignored by default. To access a code using network and location code, simply change for station MOX (network GR and empty location code):

reads . 10-APR-2010_00:00:19 5 mox z bh

to

reads . 10-APR-2010_00:00:19 5 mox.gr.. z bh

Please note: SHM doesn't yet support these schema.

• Added python version of sfdline (only subset of options supported).

Location:

SH_SHM/trunk/source/seed_io

Files:

• sfdline.c (modified) (7 diffs)
• sfdline.py (added)

SH_SHM/trunk/source/seed_io/sfdline.c

@@ -27 +27 @@
 
 #define TIMECORRUNIT 1000.0
-
 
 
@@ -46 +45 @@
         int      sfd_recno;                 /* number of records in file */
         char     sfd_stream[BC_LINELTH+1];  /* stream string */
+        char     sfd_netloc[5];             /* network and location code string */
         BOOLEAN  sfd_swap_hdr;              /* swap header necessary */
         int      i;                         /* counter */
@@ -58 +58 @@
         TSyBoolean include_bad;             /* include bad data files */
         TSyBoolean use_timecorr;            /* use time correction */
+        TSyBoolean netloc;                  /* display network/location code */
 
         /* executable code */
@@ -109 +110 @@
         if  (pa_qspecified("-timecorr"))
                 use_timecorr = TRUE;
+
+        /* output network and location code */
+        netloc = TRUE;
+        if (pa_qspecified("-nonetloc"))
+                netloc = FALSE;
 
         /* open seed file */
@@ -204 +210 @@
         if  (chanstr[1] <= ' ' || chanstr[1] > 'z')  chanstr[1] = '\0';
         ut_uncap( chanstr );
-#ifdef XXX
-        if  (strcmp(chanstr,"bh") == 0)  {         strcat( sfd_stream, "vbb-" );
-        } else if  (strcmp(chanstr,"hh") == 0)  {  strcat( sfd_stream, "vsp-" );
-        } else if  (strcmp(chanstr,"lh") == 0)  {  strcat( sfd_stream, "lp-" );
-        } else {
-                strcat( sfd_stream, chanstr );
-                strcat( sfd_stream, "-" );
-        } /*endif*/
-#endif
-        strcat( sfd_stream, chanstr );
+
+        strcat( sfd_stream, chanstr );
         strcat( sfd_stream, "-" );
         i = (int)strlen( sfd_stream );
@@ -221 +219 @@
         sfd_stream[i+1] = '\0';
         ut_uncap( sfd_stream );
+
+        /* network and location code */
+        sfd_netloc[0] = Cap(seedhdr->network[0]);
+        sfd_netloc[1] = Cap(seedhdr->network[1]);
+        sfd_netloc[2] = Cap(seedhdr->locid[0]);
+        sfd_netloc[3] = Cap(seedhdr->locid[1]);
+        sfd_netloc[4] = '\0';
+
+        /* blanks are masked by dots */
+        for (i=0; i<4; i++)
+            if (sfd_netloc[i] == ' ')
+                sfd_netloc[i] = '.';
 
         /* get sample distance */
@@ -305 +315 @@
         fclose( fp );
 
-        fprintf( out, "%c>%s %c>%s %c>%s %c>%s %c>%d %c>%d %c>%d %c>%d\n",
+        fprintf( out, "%c>%s %c>%s %c>%s %c>%s %c>%d %c>%d %c>%d %c>%d",
                 Seed_C_SfdStream, sfd_stream, Seed_C_SfdName, seedfile,
                 Seed_C_SfdTStart, sfd_t_start, Seed_C_SfdTEnd, sfd_t_end,
                 Seed_C_SfdRecno, sfd_recno, Seed_C_SfdSwapH, sfd_swap_hdr,
                 Seed_C_SfdReclth, recsize, Seed_C_SfdOffset, byteoff );
+        
+        if (netloc)
+            fprintf( out, " %c>%s",
+                Seed_C_SfdAddinf, sfd_netloc);
+
+        fprintf(out, "\n");
 
         if  (out != stdout)  fclose( out );