SET_SONSET.SHC

! set the p-onset time

 ! pick S-onset automatically


default 1 1         first trace
default 2 $dsptrcs  last seismogram

sdef cnt 1               ! counter
sdef trav                ! travel time of S wave
 
statevent_start:
   if  "cnt gti #2  goto/forward statevent_exit:

   call travel s ^distance("cnt) ^depth("cnt) &trav   
   calc t &g3 = ^origin("cnt) tadd "trav
   set/file "cnt s-onset "g3
!   set/file "cnt s-onset "trav

   calc i &cnt = "cnt + 1
   goto statevent_start:
statevent_exit:

return