Ignore:
Timestamp:
06/13/19 15:31:08 (3 months ago)
Author:
kevin
Message:

Added Arias Intensity and ColorPGA option

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/libsrc/util/rw_strongmotionII.c

    r7513 r8019  
    4040#define SM_RSA  1<<6 
    4141#define SM_QID  1<<7 
    42 #define SM_ALL (SM_SCNL | SM_TIME | SM_ALT | SM_PGA | SM_PGV | SM_PGD | SM_RSA | SM_QID) 
     42#define SM_AI   1<<8 
     43#define SM_ALL (SM_SCNL | SM_TIME | SM_ALT | SM_PGA | SM_PGV | SM_PGD | SM_RSA | SM_QID | SM_AI) 
    4344#define SM_REQ (SM_SCNL | SM_TIME) 
    4445 
     
    625626           } 
    626627           /* End of QID line processing */ 
     628       } 
     629       else if (strcmp(token, "AI:") == 0) { 
     630           if (HAVE(SM_SCNL) ) { 
     631               token = strtok(NULL, " \n"); 
     632               if (token == NULL) { 
     633                   if (logErr)  
     634                       logit("t", "error parsing SM AI line <%s>\n", line); 
     635                   status = -1; 
     636                   goto NextLine; 
     637               } 
     638               if( sscanf( token, "%lf", &sm->ai) != 1) { 
     639                   if (logErr)  
     640                       logit("t", "error parsing AI from <%s>\n", line); 
     641                   status = -1; 
     642                   goto NextLine; 
     643               } 
     644               smState |= SM_AI; 
     645           } else { 
     646               /* still looking for SCNL line to start message */ 
     647               goto NextLine; 
     648           } 
     649           /* End of AI line processing */ 
    627650       } else { 
    628651           /* unknown line; skip it for now */ 
     
    704727   if( strappend( buf, buflen, tmp ) ) return( -1 ); 
    705728 
     729/* Print the AI (Arias Intensity) */ 
     730   if( sm->ai != 0. ) 
     731   { 
     732     sprintf( tmp, "AI: %.6lf\n", sm->ai ); 
     733     if( strappend( buf, buflen, tmp ) ) return( -1 ); 
     734   } 
     735 
    706736   return( 0 ); 
    707737} 
     
    755785/* Print eventid & author */ 
    756786   logit("","\nQID: %s %s\n", sm->qid, sm->qauthor); 
     787 
     788/* Print the AI (Arias Intensity) */ 
     789   if( sm->ai != 0. ) 
     790   { 
     791     logit("","AI: %.6lf\n", sm->ai ); 
     792   } 
    757793 
    758794   return; 
Note: See TracChangeset for help on using the changeset viewer.