Commit 159928a7 authored by Zoheyr Doctor's avatar Zoheyr Doctor 🔭
Browse files

more orbital frequency fixes

parent da20a0d2
......@@ -48,11 +48,13 @@ def main():
print('evolving ',event)
event_pos_dict = {}
omega = 10./data[event]['mass_1_source']['med']
m1 = data[event]['mass_1_source']['med']
m2 = data[event]['mass_2_source']['med']
omega = 30.*((m1/((m1+m2)**3.))**(1./2))
# create a Binary instance for the event
binary = Binary(
m1 = data[event]['mass_1_source']['med'],
m2 = data[event]['mass_2_source']['med'],
m1 = m1,
m2 = m2,
c1 = 'k',
c2 = 'k',
alpha = 45.,
......
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