Commit 609ea060 authored by Roland Haas's avatar Roland Haas
Browse files

power: fix code when not having TwoPunctures.bbh

parent 3c7bcc06
......@@ -491,7 +491,8 @@ def POWER(sim_path, psi4_glob, radii, modes, f0 = FROM_TWOPUNCTURES,
"""
simdirs = os.path.join(sim_path, OUTPUT_DIR_GLOB)
meta_name = glob.glob(os.path.join(simdirs, "TwoPunctures.bbh"))[0]
if f0 == FROM_TWOPUNCTURES or ADMMass == FROM_TWOPUNCTURES:
meta_name = glob.glob(os.path.join(simdirs, "TwoPunctures.bbh"))[0]
if f0 == FROM_TWOPUNCTURES:
f0 = getCutoffFrequencyFromTwoPuncturesBBH(meta_name)
#Get simulation total mass
......@@ -674,7 +675,8 @@ def NakanoKerr(sim_path, psi4_glob, radii, modes, f0 = FROM_TWOPUNCTURES,
a_final = a_M[0]
if M_final == FROM_QUASILOCALMEASURES:
M_final = a_M[1]
meta_name = glob.glob(os.path.join(simdirs, "TwoPunctures.bbh"))[0]
if f0 == FROM_TWOPUNCTURES or ADMMass == FROM_TWOPUNCTURES:
meta_name = glob.glob(os.path.join(simdirs, "TwoPunctures.bbh"))[0]
if f0 == FROM_TWOPUNCTURES:
f0 = getCutoffFrequencyFromTwoPuncturesBBH(meta_name)
if ADMMass == FROM_TWOPUNCTURES:
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment