Changeset 4902


Ignore:
Timestamp:
07/12/12 09:51:26 (6 years ago)
Author:
paulf
Message:

upgraded to EB public release version May 2012, compile tested only

Location:
trunk/src/seismic_processing/loc_wcatwc
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/seismic_processing/loc_wcatwc/loc_wcatwc.c

    r4895 r4902  
    124124 
    125125   FullTablePath[0] = 0; 
     126 
    126127/* Check command line arguments 
    127128   ****************************/ 
     
    360361            logit( "", "reset StaArray Nsta = %d\n", Nsta );      
    361362            for ( i=0; i<Gparm.NumPBuffs; i++ ) free( StaArray[i] ); 
    362             StaArray[0] = (STATION *) calloc( MAX_STATIONS, sizeof( STATION ) ); 
    363 /* First time through to get size */ 
    364             Nsta = ReadLineupFile( Gparm.ATPLineupFileBB, StaArray[0] ); 
    365             free( StaArray[0] ); 
    366             if ( Nsta < 1 ) 
    367             { 
    368                logit( "", "Bad Lineup File read-%s\n", Gparm.ATPLineupFileBB ); 
    369                continue; 
    370             }                           
    371             StaArray[0] = (STATION *) calloc( Nsta, sizeof( STATION ) ); 
    372             Nsta = ReadLineupFile( Gparm.ATPLineupFileBB, StaArray[0] ); 
     363            Nsta = ReadLineupFile( Gparm.ATPLineupFileBB, &StaArray[0],  
     364                                   MAX_STATIONS ); 
    373365            logit( "", "New StaArray Nsta = %d\n", Nsta );        
    374366            for ( j=0; j<Nsta; j++ ) InitP( &StaArray[0][j] ); 
     
    387379               } 
    388380               memcpy( StaArray[i], StaArray[0], Nsta*sizeof( STATION ) ); 
    389                for ( j=0; j<Nsta; j++ ) InitP( &StaArray[i][j] ); 
    390381               iPBufCnt[i] = 0; 
    391382               iNumPBufRem[i] = 0; 
     
    545536      for ( j=0; j<iNSta; j++ ) 
    546537      {                            
    547          if ( !strcmp(  szChn[j], "SHZ" ) ||   /* Only use vertical in near stn table */ 
    548               !strcmp(  szChn[j], "HHZ" ) ||  
    549               !strcmp(  szChn[j], "EHZ" ) ||  
    550               !strcmp(  szChn[j], "SZ" ) ||  
    551               !strncmp( szChn[j], "BHZ", 3 ) ) 
     538//         if ( !strcmp(  szChn[j], "SHZ" ) ||   /* Only use vertical in near stn table */ 
     539//              !strcmp(  szChn[j], "HHZ" ) ||  
     540//              !strcmp(  szChn[j], "EHZ" ) ||  
     541//              !strcmp(  szChn[j], "SZ" ) ||  
     542//              !strncmp( szChn[j], "BHZ", 3 ) ) 
     543//            GetDistanceAz( &ll[i], &ll[j], &azidelt ); 
     544         if ( szChn[j][2] == 'Z' )  /* Only use vertical in near stn table */ 
    552545            GetDistanceAz( &ll[i], &ll[j], &azidelt ); 
    553546         else 
  • trunk/src/seismic_processing/loc_wcatwc/locate.c

    r4895 r4902  
    5555 
    5656     /************************************************************** 
    57       *                        AddInMw()                           * 
     57      *                        AddInMwL()                          * 
    5858      *                                                            * 
    5959      * This function reads the Mw data file, determines if the    * 
     
    6767      **************************************************************/ 
    6868       
    69 void AddInMw( char *pszFile, HYPO *pHypo ) 
     69void AddInMwL( char *pszFile, HYPO *pHypo ) 
    7070{ 
    7171   double  dMw[MAX_STATIONS];           /* Mw from data file */ 
    7272   double  dMwSumMod;                   /* Mw summation */ 
     73   double  dDum;                        /* Temp for spectral info */ 
    7374   FILE    *hFile;                      /* Mw File handle */ 
    74    int     i; 
     75   int     i, j, iNumPers; 
    7576   int     iMwCountMod;                 /* # Mws from file(within .6 of avg.) */ 
    7677   HYPO    HypoT;                       /* Quake params. in Mw file */ 
     
    101102      while ( !feof( hFile ) ) 
    102103      {                                 /* Yes, so add in Mw */ 
    103          fscanf( hFile, "%s %s %s %s %lf\n", szStn, szChn, szNet, szLoc, 
    104                  &dMw[pHypo->iNumMw] ); 
     104         fscanf( hFile, "%s %s %s %s %lf %d \n", szStn, szChn, szNet, szLoc, 
     105                 &dMw[pHypo->iNumMw], &iNumPers ); 
    105106         pHypo->dMwAvg += dMw[pHypo->iNumMw]; 
    106          pHypo->iNumMw += 1; 
     107         pHypo->iNumMw += 1; 
     108         for ( j=0; j<iNumPers; j++ )  
     109            fscanf( hFile, "%lf %lf %lf %lf %lf\n",  
     110                   &dDum, &dDum, &dDum, &dDum, &dDum ); 
    107111      } 
    108112      if ( pHypo->iNumMw ) pHypo->dMwAvg /= (double) pHypo->iNumMw; 
     
    703707   ZeroMagnitudes( pSta, iNum ); 
    704708   ComputeMagnitudes( iNum, pSta, pHypo ); 
    705    AddInMw( Gparm->szMwFile, pHypo ); 
     709   AddInMwL( Gparm->szMwFile, pHypo ); 
    706710   GetPreferredMag( pHypo ); 
    707711    
Note: See TracChangeset for help on using the changeset viewer.