From 27c9bbb96af91df4b19355de7e65443d0f97a29f Mon Sep 17 00:00:00 2001 From: Jeldrik Gaal Date: Thu, 6 May 2021 13:48:44 +0200 Subject: [PATCH] Manualy adding generic Amplitude to .obs files Due to the obspy write function not including the generic Amplitude/A0 in the .obs file (NonLinLoc) even when the Event holds this Data. Obspy Team has been informed by opening an issue on Github --- pylot/core/io/data.py | 10 +++++++++- 1 file changed, 9 insertions(+), 1 deletion(-) diff --git a/pylot/core/io/data.py b/pylot/core/io/data.py index d7faced9..e497783b 100644 --- a/pylot/core/io/data.py +++ b/pylot/core/io/data.py @@ -1,4 +1,4 @@ -#!/usr/bin/env python +#!/usr/bin/env pyth n # -*- coding: utf-8 -*- import copy @@ -353,6 +353,14 @@ class Data(object): header = '# EQEVENT: Label: EQ%s Loc: X 0.00 Y 0.00 Z 10.00 OT 0.00 \n' % evid nllocfile = open(fnout + fnext) l = nllocfile.readlines() + # Adding A0/Generic Amplitude to .obs file + l2 = [] + for li in l: + for amp in evtdata_org.amplitudes: + if amp.waveform_id.station_code == li[0:5].strip(): + li = li[0:64] + '{:0.2e}'.format(amp.generic_amplitude) + li[73:-1] + '\n' + l2.append(li) + l = l2 nllocfile.close() l.insert(0, header) nllocfile = open(fnout + fnext, 'w')