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')
