Calculate distance/azimuth in SAC output

iris-edu · Jul 3, 2023 · ea0c05f · ea0c05f · adam-iris · Sep 7, 2023
1 parent 3304089
commit ea0c05f
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/pyweed/waveforms_handler.py b/pyweed/waveforms_handler.py
@@ -455,6 +455,7 @@ def to_sac_stream(self, st, waveform):
         tr = SACTrace.from_obspy_trace(st[0])
         tr.kevnm = waveform.event_description[:16]
         event = waveform.event_ref()
+        origin = None
         if not event:
             LOGGER.warn("Lost reference to event %s", waveform.event_description)
         else:
@@ -482,6 +483,9 @@ def to_sac_stream(self, st, waveform):
             tr.cmpaz = channel.azimuth
             tr.cmpinc = channel.dip + 90
             tr.kinst = channel.sensor.description[:8]
+        if origin and channel:
+            # Calculate distances/azimuths
+            tr.lcalda = True
         return Stream([tr.to_obspy_trace()])