p = subprocess.Popen ([cutseismo, start_cmd, str(marginseismo_post)], cwd = d, stdout = subprocess.PIPE) p.wait () out = p.stdout.read () if p.returncode != 0: print "%s: GAKS: Failed to extract seismometer wavefile, skipping." % d continue n = out.find ('file to: ') mseedf = out[n+9:-2] print "%s: GAKS: Moving seismometer wav to 02..: %s" % (d, mseedf) if os.path.exists (os.path.join (ex_02, mseedf)): print "%s: GAKS: mseed exists, deleting old." % d os.remove (os.path.join (ex_02, mseedf)) shutil.move (os.path.join (sdir, mseedf), ex_02) rd = open (os.path.join (ex_02, 'extract_done'), 'w') rd.write ('yes\n') rd.close () else: print "%s: %s already extracted, skipping station." % (d, s) if putautoextracted: print "%s: putting to 04_events_localize.." % d p = Putrevised () p.doevent (d, d, '02')
if not reextract_done: idf = open (os.path.join (d, s, '00', idfiles[s]), 'r') ids = [] for l in idf.readlines (): (ii, sdlag) = l.split (',') ids.append (int(ii)) idf.close () mkm = Makemseed () mkm.station = s mkm.root = os.path.join (d, s) mkm.destdir = os.path.join (d, s, '01') mkm.ids = ids (mseedf, idsf, refsf) = mkm.dorange () rd = open (os.path.join (d, s, '01', 'reextract_done'), 'w') rd.write ('yes\n') rd.close () else: print "%s: %s already reextracted, skipping." % (d, s) print "%s: putting to 04_events_localize.." % d p = Putrevised () p.doevent (d, d, '01') else: print "%s: not fully revised. ignoring." % d
str(marginseismo_post)], cwd=d, stdout=subprocess.PIPE) p.wait() out = p.stdout.read() if p.returncode != 0: print "%s: GAKS: Failed to extract seismometer wavefile, skipping." % d continue n = out.find('file to: ') mseedf = out[n + 9:-2] print "%s: GAKS: Moving seismometer wav to 02..: %s" % (d, mseedf) if os.path.exists(os.path.join(ex_02, mseedf)): print "%s: GAKS: mseed exists, deleting old." % d os.remove(os.path.join(ex_02, mseedf)) shutil.move(os.path.join(sdir, mseedf), ex_02) rd = open(os.path.join(ex_02, 'extract_done'), 'w') rd.write('yes\n') rd.close() else: print "%s: %s already extracted, skipping station." % (d, s) if putautoextracted: print "%s: putting to 04_events_localize.." % d p = Putrevised() p.doevent(d, d, '02')