Adjust center of mass position and velocity of simulation for GW150914


In [3]:
import scri
import scri.SpEC
import numpy as np

data_dir = '/Users/boyle/Research/Data/SimulationAnnex/CatalogLinks/SXS:BBH:0317/Lev3/'

In [4]:
w_N2 = scri.SpEC.remove_avg_com_motion(data_dir + 'rhOverM_Asymptotic_GeometricUnits.h5/Extrapolated_N2.dir', file_write_mode='w')
w_N3 = scri.SpEC.remove_avg_com_motion(data_dir + 'rhOverM_Asymptotic_GeometricUnits.h5/Extrapolated_N3.dir', file_write_mode='a')
w_N4 = scri.SpEC.remove_avg_com_motion(data_dir + 'rhOverM_Asymptotic_GeometricUnits.h5/Extrapolated_N4.dir', file_write_mode='a')
w_No = scri.SpEC.remove_avg_com_motion(data_dir + 'rhOverM_Asymptotic_GeometricUnits.h5/OutermostExtraction.dir', file_write_mode='a')


Optimal x_i: [-0.04796393008437783, -0.03749570036059665, 0.0002546296438184966]
Optimal v_i: [-2.3731089790623737e-06, 1.0909676453474665e-05, -4.251262320063484e-07]
t_i, t_f: 165.5, 14917.0
Optimal x_i: [-0.04796393008437783, -0.03749570036059665, 0.0002546296438184966]
Optimal v_i: [-2.3731089790623737e-06, 1.0909676453474665e-05, -4.251262320063484e-07]
t_i, t_f: 165.5, 14917.0
Optimal x_i: [-0.04796393008437783, -0.03749570036059665, 0.0002546296438184966]
Optimal v_i: [-2.3731089790623737e-06, 1.0909676453474665e-05, -4.251262320063484e-07]
t_i, t_f: 165.5, 14917.0
Optimal x_i: [-0.04796393008437783, -0.03749570036059665, 0.0002546296438184966]
Optimal v_i: [-2.3731089790623737e-06, 1.0909676453474665e-05, -4.251262320063484e-07]
t_i, t_f: 165.5, 14917.0

In [ ]: