From e1c2629f97513fa0b139a4822118c6800a44a313 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Ludger=20K=C3=BCperkoch?= <kueperkoch@bestec-for-nature.com>
Date: Wed, 18 Nov 2015 11:11:10 +0100
Subject: [PATCH] Modified writephases for NLLoc: writes dummy onset times
 (start of seismic trace) into phase file in order to derive theoretical onset
 times later needed for iterative picking.

---
 pylot/core/pick/utils.py | 92 +++++++++++++++++++++++-----------------
 1 file changed, 52 insertions(+), 40 deletions(-)

diff --git a/pylot/core/pick/utils.py b/pylot/core/pick/utils.py
index 9887c78d..f4eb2282 100644
--- a/pylot/core/pick/utils.py
+++ b/pylot/core/pick/utils.py
@@ -961,46 +961,58 @@ def writephases(arrivals, fformat, filename):
         # write header
         fid.write('# EQEVENT:  Label: EQ001  Loc:  X 0.00  Y 0.00  Z 10.00  OT 0.00 \n')
         for key in arrivals:
-            if arrivals[key]['P']['weight'] < 4:
-                fm = arrivals[key]['P']['fm']
-                if fm == None:
-                   fm = '?'
-                onset = arrivals[key]['P']['mpp']
-                year = onset.year
-                month = onset.month
-                day = onset.day
-                hh = onset.hour
-                mm = onset.minute
-                ss = onset.second
-                ms = onset.microsecond
-                ss_ms = ss + ms / 1000000.0
-                fid.write('%s ? ? ? P   %s %d%02d%02d %02d%02d %7.4f GAU 0 0 0 0 1 \n' % (key,
-                 fm,
-                 year,
-                 month,
-                 day,
-                 hh,
-                 mm,
-                 ss_ms))
-            if arrivals[key]['S']['weight'] < 4:
-                fm = '?'
-                onset = arrivals[key]['S']['mpp']
-                year = onset.year
-                month = onset.month
-                day = onset.day
-                hh = onset.hour
-                mm = onset.minute
-                ss = onset.second
-                ms = onset.microsecond
-                ss_ms = ss + ms / 1000000.0
-                fid.write('%s ? ? ? S   %s %d%02d%02d %02d%02d %7.4f GAU 0 0 0 0 1 \n' % (key,
-                 fm,
-                 year,
-                 month,
-                 day,
-                 hh,
-                 mm,
-                 ss_ms))
+             # P onsets
+             if arrivals[key]['P']:
+                 fm = arrivals[key]['P']['fm']
+                 if fm == None:
+                     fm = '?'
+                 onset = arrivals[key]['P']['mpp']
+                 year = onset.year
+                 month = onset.month
+                 day = onset.day
+                 hh = onset.hour
+                 mm = onset.minute
+                 ss = onset.second
+                 ms = onset.microsecond
+                 ss_ms = ss + ms / 1000000.0
+                 if arrivals[key]['P']['weight'] < 4:
+                     pweight = 1 # use pick
+                 else:
+                     pweight = 0 # do not use pick
+                 fid.write('%s ? ? ? P   %s %d%02d%02d %02d%02d %7.4f GAU 0 0 0 0 %d \n' % (key,
+                  fm,
+                  year,
+                  month,
+                  day,
+                  hh,
+                  mm,
+                  ss_ms,
+                  pweight))
+             # S onsets
+             if arrivals[key]['S']:
+                 fm = '?'
+                 onset = arrivals[key]['S']['mpp']
+                 year = onset.year
+                 month = onset.month
+                 day = onset.day
+                 hh = onset.hour
+                 mm = onset.minute
+                 ss = onset.second
+                 ms = onset.microsecond
+                 ss_ms = ss + ms / 1000000.0
+                 if arrivals[key]['S']['weight'] < 4:
+                     sweight = 1 # use pick
+                 else:
+                     sweight = 0 # do not use pick
+                 fid.write('%s ? ? ? S   %s %d%02d%02d %02d%02d %7.4f GAU 0 0 0 0 %d \n' % (key,
+                  fm,
+                  year,
+                  month,
+                  day,
+                  hh,
+                  mm,
+                  ss_ms,
+                  sweight))
 
         fid.close()