Changeset 17 for SH_SHM/trunk

Timestamp:

10/30/2008 08:28:00 PM (14 years ago)

Author:

marcus

Message:

internal: r2 | svn | 2008-01-04 07:52:19 +0100 (Fr, 04 Jan 2008) | 1 line

implemented measurement of amplitudes (vel/disp) on traces filtered with bandpass filters

File:

• SH_SHM/trunk/source/motif/cbutil.c (modified) (5 diffs)

SH_SHM/trunk/source/motif/cbutil.c

@@ -1937 +1937 @@
         float      frq;                    /* frequency of signal in Hz */
         float      tffac;                  /* transfer function factor */
+        float      tffac2;                 /* second factor of this kind */
         char       lfilname[BC_LINELTH+1]; /* local filter name */
         TIME       ttime;                  /* trace time */
@@ -1957 +1958 @@
         frq = 1/period;
 
+        /* if a filter is applied the comment info entry contains its name */
         locstat = BC_NOERROR;
         db_gets( trcptr, ES_COMMENT, BC_LINELTH, str, &locstat );
         if  (locstat == BC_NOERROR)  {
+                /* we have to look for a substring starting with 'FLT:' */
                 filptr = strstr( str, "FLT:" );
                 if  (filptr == NULL)  filptr = strstr( str, "flt:" );
@@ -1968 +1971 @@
                 } /*endif*/
         } else {
+                /* if no filter found in the comment string it is an unfiltered trace */
                 *filter = '\0';
         } /*endif*/
         strcpy( lfilname, filter );
 
+        /* if it is an unfiltered trace, take the transfer function of the recording
+         * instrument.  This is returned by cu_lookup_filter. */
         if  (*lfilname == '\0')  {
                 locstat = cBcNoError;
@@ -1982 +1988 @@
         } /*endif*/
 
-        /* build filter name & read amplitude */
+        /* if the trace is unfiltered or filtered with an simulation filter, the
+         * transfer function is found by applying a prefix 'TF_VEL_' or 'TF_DSP_'.
+         * Then this this contains the complete information needed here.
+         * Build filter name & read amplitude */
         strcpy( str, "TF_" );
         strcat( str, kind );
@@ -1988 +1997 @@
         strcat( str, lfilname );
         tffac = ff_filter_amplitude( str, frq, status );
+
+        /* if the filter cannot be found with a TF_VEL/DSP prefix it is not a
+         * simulation filter, so we need the transfer function of the filter
+         * and the transfer function of the recording instrument.
+         */
+        if  (SySevere(status))  {
+                /* at an unfiltered trace this means that the transfer function of the
+                 * recording instrument has not been found before.  So nothing can
+                 * be done here. */
+                if  (*lfilname == '\0')  return;
+                *status = cBcNoError;
+                /* get amplification of applied filter */
+                tffac = ff_filter_amplitude( lfilname, frq, status );
+                printf( "SHM-dbg2: fac1: filter applied: %s -> %f\n", lfilname, tffac );
+                if  (SySevere(status))  return;
+                /* add transfer function of recording instrument (like code above) */
+                locstat = cBcNoError;
+                db_gett( trcptr, ET_START, &ttime, &locstat );
+                tc_a2t( &ttime, strtime, &locstat );
+                cu_get_stream_string( trcptr, str, status );
+                if  (SySevere(status))  return;
+                cu_lookup_filter( str, strtime, lfilname, &trcflags, status );
+                if  (SySevere(status))  return;
+                strcpy( str, "TF_" );
+                strcat( str, kind );
+                strcat( str, "_" );
+                strcat( str, lfilname );
+                tffac2 = ff_filter_amplitude( str, frq, status );
+                printf( "SHM-dbg2: fac2: transfer function of rec. instr. %s -> %f\n",
+                        str, tffac2 );
+                tffac *= tffac2;
+        } /*endif*/
 
         if  (GpGetInt(cGpI_debug_level) > 1)