[1]:

import CactusTool

name = 'BBH0001_D15_Rin200'
sim = CactusTool.load(name, '/Volumes/simulations/HUST_Catalog')
qlm = sim.ThornOutput('quasilocalmeasures')
bbh = sim.ThornOutput('TwoPunctures')

[4]:

qlm.spin.preview()

[3]:

qlm.Spin1.norm.preview()

[7]:

print(bbh)

m1 = 0.5538461538
m2 = 0.4461538462
ADMMass = 0.9899366929086094
Mass Ratio = 1.2413793101129609
eta = 0.24710059172094676
separation = 10.0
spin1 = [0.         0.         0.09509112]
spin2 = [ 0.         0.        -0.0915645]
chi1 = [0.   0.   0.31]
chi2 = [ 0.    0.   -0.46]
chi_eff = -0.033538461502796824
chi_p = 0.0
Omega_orb = 0.027784006478250793
CutoffFrequency = 0.041676009717376186


[16]:

import matplotlib.pyplot as plt

plt.xlabel('time [M]')
plt.ylabel('qlm_spin')
plt.plot(qlm.spin.t, qlm.spin.y[:, 0], alpha = 0.7, label='qlm_spin[0]')
plt.plot(qlm.spin.t, qlm.spin.y[:, 1], alpha = 0.7, label='qlm_spin[1]')
plt.savefig('/Users/liuyu/Downloads/qlm_spin.png')
../_images/Tutorial_BBH_quasilocalmeasures_4_0.png 
[12]:

import numpy as np
a = [[1,1], [1,1], [1,1]]
b = np.array(a)
b.shape

[12]:

(3, 2)

[ ]: