From 8a843428838dd56460056ee5dfd88c88ec206a85 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 10:27:42 +0100
Subject: [PATCH] Implementation of module /loc/nll finished, running for
 single and multiple event processing.

---
 autoPyLoT.py | 52 ++++++++++++++--------------------------------------
 1 file changed, 14 insertions(+), 38 deletions(-)

diff --git a/autoPyLoT.py b/autoPyLoT.py
index d0c94689..85d65a48 100755
--- a/autoPyLoT.py
+++ b/autoPyLoT.py
@@ -78,12 +78,11 @@ def autoPyLoT(inputfile):
             phasef = parameter.getParam('phasefile')
             phasefile = '%s/obs/%s' % (nllocroot, phasef)
             # get name of NLLoc-control file
-            locf = parameter.getParam('locfile')
-            locfile = '%s/run/%s' % (nllocroot, locf)
-            # patter of NLLoc ttimes from location grid
+            ctrf = parameter.getParam('ctrfile')
+            ctrfile = '%s/run/%s' % (nllocroot, ctrf)
+            # pattern of NLLoc ttimes from location grid
             ttpat = parameter.getParam('ttpatter')
-            ttpatter = '%s/time/%s' % (nllocroot, ttpat)
-            # patter of NLLoc-output file
+            # pattern of NLLoc-output file
             nllocoutpatter = parameter.getParam('outpatter')
         else:
             locflag = 0
@@ -112,27 +111,15 @@ def autoPyLoT(inputfile):
                     picksExport(picks, 'NLLoc', phasefile)
 
                     # For locating the event the NLLoc-control file has to be modified!
-                    # create comment line for NLLoc-control file
-                    # NLLoc-output file
                     evID = event[string.rfind(event, "/") + 1 : len(events) - 1]
-                    nllocout = '%s/loc/%s_%s' % (nllocroot, evID, nllocoutpatter)
-                    locfiles = 'LOCFILES %s NLLOC_OBS %s %s 0' % (phasefile, ttpatter, nllocout)
-                    print ("Modifying  NLLoc-control file %s ..." % locfile)
-                    # modification of NLLoc-control file
-                    filedata = None
-                    nllfile = open(locfile, 'r')
-                    filedata = nllfile.read()
-                    if filedata.find(locfiles) < 0:
-                        # replace old command
-                        filedata = filedata.replace('LOCFILES', locfiles)
-                        nllfile = open(locfile, 'w')
-                        nllfile.write(filedata)
-                        nllfile.close()
+                    nllocout = '%s_%s' % (evID, nllocoutpatter)
+                    # create comment line for NLLoc-control file
+                    modifyInputFile(ctrf, nllocroot, nllocout, phasef, ttpat)
 
                     # locate the event
-                    subprocess.call([nlloccall, locfile])
+                    locate(nlloccall, ctrfile)
 
-                    # !iterative picking if traces remained unpicked or with bad picks!
+                    # !iterative picking if traces remained unpicked or occupied with bad picks!
                     # get theoretical onset times for picks with weights >= 4
                     # in order to reprocess them using smaller time windows
                 ##########################################################
@@ -167,25 +154,14 @@ def autoPyLoT(inputfile):
                 picksExport(picks, 'NLLoc', phasefile)
 
                 # For locating the event the NLLoc-control file has to be modified!
-                # create comment line for NLLoc-control file NLLoc-output file
-                nllocout = '%s/loc/%s_%s' % (nllocroot, parameter.getParam('eventID'), nllocoutpatter)
-                locfiles = 'LOCFILES %s NLLOC_OBS %s %s 0' % (phasefile, ttpatter, nllocout)
-                print ("Modifying  NLLoc-control file %s ..." % locfile)
-                # modification of NLLoc-control file
-                filedata = None
-                nllfile = open(locfile, 'r')
-                filedata = nllfile.read()
-                if filedata.find(locfiles) < 0:
-                    # replace old command
-                    filedata = filedata.replace('LOCFILES', locfiles)
-                    nllfile = open(locfile, 'w')
-                    nllfile.write(filedata)
-                    nllfile.close()
+                nllocout = '%s_%s' % (parameter.getParam('eventID'), nllocoutpatter)
+                # create comment line for NLLoc-control file
+                modifyInputFile(ctrf, nllocroot, nllocout, phasef, ttpat)
 
                 # locate the event
-                subprocess.call([nlloccall, locfile])
+                locate(nlloccall, ctrfile)
 
-                # !iterative picking if traces remained unpicked or with bad picks!
+                # !iterative picking if traces remained unpicked or occupied with bad picks!
                 # get theoretical onset times for picks with weights >= 4
                 # in order to reprocess them using smaller time windows
             ##########################################################