Changeset 7517


Ignore:
Timestamp:
06/28/18 22:42:14 (3 weeks ago)
Author:
baker
Message:

cleanup warning: variable set but not used [-Wunused-but-set-variable]

Location:
trunk/src/libsrc/earlybird
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/libsrc/earlybird/diskrw.c

    r7219 r7517  
    288288   double  dInt;                       /* Time interval between samps */ 
    289289   int     i, j, k, l, m; 
    290    static  int     iAllRead;           /* 1 -> all needed data files read */ 
    291290   int     iInt;                       /* # samps since start of second*/ 
    292291   static  int     iMaxToRead;         /* Number of files to read */ 
    293    static  int     iRead[MAX_STATIONS];/* Flag which indicates read status 
    294                                           0   -> No data yet for this station 
    295                                           1-X -> Number of files read for stn */ 
    296292   time_t  iTime;                      /* time (1/1/70) at min P time */    
    297293   int32_t lStart;                     /* Starting data index of buffer */ 
     
    318314#endif 
    319315/* Initialize some things */ 
    320    iAllRead = 0; 
    321316   for ( i=0; i<iNumStas; i++ ) 
    322317   { 
    323       iRead[i] = 0; 
    324318      StaArray[i].lRawCircCtr   = 0; 
    325319      StaArray[i].lSampIndexF   = 0; 
     
    519513   static  FILE *hFile;        /* Handle to seismic data file */ 
    520514   static  DISKHEADER dh;      /* DISKHEADER information (see earlybirdlib.h) */ 
    521    static  double dPreTime;    /* Time (s) to look for DC offset */ 
    522515   static  double dStartTime;  /* Nominal File start time in epochal Seconds */ 
    523516   double  dTemp; 
     
    539532   DateToModJulianSec( &stDate, &dSCNLTime ); 
    540533   dSCNLTime -= 3506630400.;  /* Then to epochal seconds */ 
    541 /* For filtering... */ 
    542    if ( iLP == 1 ) dPreTime = PRE_LP_TIME; 
    543    else            dPreTime = PRE_P_TIME; 
    544534 
    545535/* Try to open data file */ 
  • trunk/src/libsrc/earlybird/dummy.c

    r6803 r7517  
    1414#endif 
    1515 
     16#include <fcntl.h> 
     17#include <math.h> 
     18#include <stdio.h> 
    1619#include <stdlib.h> 
    17 #include <stdio.h> 
    18 #include <math.h> 
     20#include <string.h> 
    1921#include <time.h>               
    20 #include <string.h> 
     22#include <sys/stat.h> 
    2123#include <sys/types.h> 
    22 #include <sys/stat.h> 
    23 #include <fcntl.h> 
    24 #include <earthworm.h> 
     24 
     25#include "earthworm.h" 
    2526#include "earlybirdlib.h" 
     27 
    2628#ifdef _WINNT 
    2729 #include <crtdbg.h> 
     
    134136   double  dDum = 0.0; 
    135137   FILE    *hFile;        /* File handle */ 
    136    int     iCnt;          /* Dummy file open counter */ 
    137138   int     iDepth = 0;    /* Quake depth from dummy file */ 
    138139   int     iDum = 0; 
     
    146147         
    147148/* Read Dummy File */ 
    148    iCnt = 0;  
    149149   if ( (hFile = openFile( pszDumFile, "r" )) != NULL ) 
    150150      fscanf( hFile, "%d %lf %lf %d %lf %s %d %d %d %d %d %d " 
  • trunk/src/libsrc/earlybird/logquake.c

    r6803 r7517  
    1212#include <string.h> 
    1313#include <time.h> 
    14 #include <earthworm.h> 
    15 #include <transport.h> 
     14 
     15#include "earthworm.h" 
     16#include "transport.h" 
    1617#include "earlybirdlib.h" 
    1718 
     
    138139   int     iTenths;             /* Tenths of seconds in P-time */ 
    139140   LATLON  ll;                  /* Geographic epicenter coordinates */ 
    140    time_t  lTime,logTime;       /* 1/1/70 time in seconds */ 
     141   time_t  lTime;               /* 1/1/70 time in seconds */ 
    141142   char    *psz;                /* Epicentral region */ 
    142143   char    szHour[3], szMin[3], szSec[3];  /* P-time in string form */ 
     
    178179/* Get current date/time and write on top of output header */ 
    179180   time( &lTime ); 
    180    logTime = lTime;                           
    181181   if ( iType == 1 ) 
    182182   { 
  • trunk/src/libsrc/earlybird/mags.c

    r6803 r7517  
    88     ******************************************************************/ 
    99 
     10#include <math.h> 
     11#include <stdio.h> 
    1012#include <stdlib.h> 
    11 #include <stdio.h> 
    12 #include <math.h> 
     13#include <string.h> 
    1314#include <time.h> 
    14 #include <string.h> 
    15 #include <earthworm.h> 
     15 
     16#include "earthworm.h" 
    1617#include "earlybirdlib.h" 
    1718 
     
    22722273   static  double  x[MAXMWPARRAY], x_next[MAXMWPARRAY], w[MAXMWPARRAY]; 
    22732274   int     ic, index, division_counter; 
    2274    double  y, z, smooth[4], nonsmooth[4]; 
     2275   double  y, z; 
     2276#if 0 
     2277   double  smooth[4]; 
     2278#endif 
    22752279   double  smooth2[6], nonsmooth2[6]; 
    22762280   double  avg_count = 0., avg = 0., two = 2.; 
     
    22782282   double  xmin, tau, a, fc, wi; 
    22792283   double  sum0, sum1, sum2, sum3; 
    2280    static  double  decomp_mag[MAXMWPARRAY]; 
    22812284   double  weight[100]; 
    22822285     
     
    22982301     x_next[n] = Z_in[n]; 
    22992302   } 
    2300    
     2303 
    23012304/* Daubechies order 4 wavelet coeffs */ 
     2305#if 0 
    23022306   smooth[0] = (1. + sqrt(3.0)) / (4. * sqrt(2.)); 
    23032307   smooth[1] = (3. + sqrt(3.0)) / (4. * sqrt(2.)); 
    23042308   smooth[2] = (3. - sqrt(3.0)) / (4. * sqrt(2.)); 
    23052309   smooth[3] = (1. - sqrt(3.0)) / (4. * sqrt(2.)); 
    2306     
     2310#endif 
     2311 
    23072312/* Daubechies order 6 wavelet coeffs */ 
    23082313   smooth2[0] = .332671; 
     
    23172322   { 
    23182323      nonsmooth2[n] = z * smooth2[5 - n]; 
    2319       if (n < 4) nonsmooth[n] = z * smooth[3 - n]; 
    23202324      z *= -1.; 
    23212325   } 
     
    23412345            } 
    23422346            avg += (w[m] * w[m]); 
    2343             decomp_mag[k] = w[m] * w[m]; 
    23442347            avg_count += 1.; 
    23452348                 
     
    23652368            } 
    23662369            avg += (w[m] * w[m]); 
    2367             decomp_mag[k] = w[m] * w[m]; 
    23682370            avg_count += 1.;             
    23692371            m += 1; 
  • trunk/src/libsrc/earlybird/taulib.c

    r7381 r7517  
    1919#endif 
    2020 
    21 #include <stdlib.h> 
     21#include <math.h> 
    2222#include <stddef.h> 
    2323#include <stdio.h> 
     24#include <stdlib.h> 
    2425#include <string.h> 
    25 #include <math.h> 
    26 #include <earthworm.h> 
     26 
     27#include "earthworm.h" 
    2728#include "earlybirdlib.h" 
    2829#include "iasplib.h" 
     
    133134int brnset( int nn, char pcntl[10][PHASE_LENGTH], int prflg[] ) 
    134135{ 
    135    bool    fnd, all; 
     136   bool    all; 
     137#if DEBUGL 
     138   bool    fnd; 
     139#endif 
    136140   static int     nsgpt[JBRN];            /* Based on 0 */ 
    137141   static int     ncmpt[NCMD][2] = { {1, 2}, {1, 7}, {1, 13}, {13, 16} }; 
     
    265269      for ( j=0; j<NCMD; j++ ) if ( !strcmp( phlst[i], cmdcd[j] ) ) goto Match; 
    266270/* If the token isn't a keyword, see if it is a generic branch name. */ 
     271#if DEBUGL 
    267272      fnd = FALSE; 
     273#endif 
    268274      for ( k=0; k<=kseg; k++ ) 
    269275      { 
    270276         if ( strcmp( phlst[i], segcd[k] ) ) continue; 
     277#if DEBUGL 
    271278         fnd        = TRUE; 
     279#endif 
    272280         l          = nsgpt[k]; 
    273281         bSegmsk[l] = TRUE; 
     
    304312   Branch: 
    305313   if ( !prflg[2] ) return( 1 ); 
     314#if DEBUGL 
    306315   fnd = FALSE; 
     316#endif 
    307317   j2  = 0; 
    308318/* Loop over segments. */ 
     
    328338      if ( !fnd )  
    329339         logit( "", "Brnset:  the following phases have been selected:\n" ); 
    330 #endif 
    331340      fnd = TRUE; 
    332 #if DEBUGL 
    333341      logit( "", "     %d  ", i+1 ); 
    334342      for ( j=j1-1; j<j2; j++ ) logit( "", "%10s", segcd[j] ); 
     
    821829   int     le;                                  /* Based on 1 */ 
    822830   double  x, p0, p1, arg, dp, dps = 0.0, delp, tol=3.e-6, ps, deps=1.e-10, 
    823            hsgn, dp0, temp, dsgn, dpn; 
     831           hsgn, temp, dsgn, dpn; 
     832#if DEBUGL 
     833   double  dp0; 
     834#endif 
    824835 
    825836   nph = abs (iDel[jb][2]) - 1; 
     
    848859            if (dps < 0.) dp = dps*dps*(-1.); 
    849860            else        dp = dps*dps; 
     861#if DEBUGL 
    850862            dp0 = dp; 
     863#endif 
    851864            if ( dp < p1-delp || dp > p0+delp ) goto CheckHere; 
    852865            if ( *n >= MAX_PHASES ) goto SayWarning; 
     
    877890               if ( dps < 0. ) dp = dps * dps * (-1.); 
    878891               else            dp = dps * dps; 
     892#if DEBUGL 
    879893               dp0 = dp; 
     894#endif 
    880895            } 
    881896            else 
     
    961976             double coef2[5][JOUT] ) 
    962977{ 
    963    double  a[2][60], ap[3], b[60], alr, g1, gn; 
     978   double  a[2][60], ap[3], b[60], alr, gn; 
    964979   int     i, ie, is, j, n, n1; 
    965980 
     
    10161031      j = i; 
    10171032   } 
    1018    g1 = (x1 - coef2[3][i1-1]*b[0] - coef2[4][i1-1]*b[1]) / coef2[2][i1-1]; 
    10191033 
    10201034   tau2[i1-1][1] = x1; 
     
    19591973double umod( double zs, int *isrc, int nph ) 
    19601974{ 
    1961    double dep; 
    19621975   int    i, j;                /* Index based on 0 */ 
    19631976   double dtol = 1.e-6; 
     
    19651978   for ( i=1; i<iMt[nph]; i++ )  
    19661979      if ( dZm[i][nph] <= zs ) goto ItsOK; 
    1967    dep = (1. - exp (zs)) / dXn; 
    19681980   return 0.; 
    19691981 
Note: See TracChangeset for help on using the changeset viewer.