Changeset 4633


Ignore:
Timestamp:
10/25/11 12:22:08 (10 years ago)
Author:
saurel
Message:

P-wave velocity used was expressed in m/s where it should have been in km/s

File:
1 edited

Legend:

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

    r4620 r4633  
    821821                                break; 
    822822                        for(i=0; i<maxEvent ;i++) { 
    823                                 time_diff = Event[i].Origin.ot - Sum.ot - Sum.dmin / 5000; 
     823                                time_diff = Event[i].Origin.ot - Sum.ot - Sum.dmin / 5; 
    824824                                if(ABS(time_diff) < DelayTime) { 
    825825                                        new_event = &Event[i]; 
     
    857857 
    858858/******************** Message Processing Thread ******************* 
    859  *           Take messages from memory queue, create trace        * 
    860  *      snippets, and call the disposal (PutAway) routines        * 
     859 *           Take messages from memory queue,                     * 
     860 *             output message if timeout has been reached         * 
    861861 ******************************************************************/ 
    862862 
     
    986986 
    987987/******************************************************************  
    988  *  Add new line in the Main Courante file hosted in OutputDir    * 
     988 *  Create output HYP2000_ARC message and send it                 * 
    989989 ******************************************************************/ 
    990990int OutputEvent (EVENT *OutEvent) 
     
    10791079    *outptr = *outptr + 8; 
    10801080 
    1081     sprintf(*outptr,"%5.0f   %3d%3d%3d%4.0f",SumP.z*100,SumP.nph,SumP.gap,SumP.dmin,SumP.rms*100); 
    1082     *outptr = *outptr + 21; 
     1081    sprintf(*outptr,"%5.0f",SumP.z*100); 
     1082    *outptr = *outptr + 5; 
     1083    if ((SumP.labelpref == 'L') || (SumP.labelpref == 'X'))     /* If preferred Magnitude is Amplitude, update correct field */ 
     1084        sprintf(*outptr,"%3.0f%3d%3d%3d%4.0f",SumP.Mpref*100,SumP.nph,SumP.gap,SumP.dmin,SumP.rms*100); 
     1085    else 
     1086        sprintf(*outptr,"   %3d%3d%3d%4.0f",SumP.nph,SumP.gap,SumP.dmin,SumP.rms*100); 
     1087    *outptr = *outptr + 16; 
    10831088 
    10841089    sprintf(*outptr,"%3d%2d%4.0f%3d%2d%4.0f",SumP.e0az,SumP.e0dp,SumP.e0*100,SumP.e1az,SumP.e1dp,SumP.e1*100); 
Note: See TracChangeset for help on using the changeset viewer.