In [79]:

    

# Setup ipython environment
%load_ext autoreload
%autoreload 2
# %matplotlib auto
%matplotlib inline
# Import useful things
from nrutils import scsearch,gwylm
# Setup plotting backend
import matplotlib as mpl
from mpl_toolkits.mplot3d import axes3d
mpl.rcParams['lines.linewidth'] = 0.8
mpl.rcParams['font.family'] = 'serif'
mpl.rcParams['font.size'] = 12
mpl.rcParams['axes.labelsize'] = 20
mpl.rcParams['axes.titlesize'] = 20
from matplotlib.pyplot import *
from mpl_toolkits.mplot3d import Axes3D
#
from numpy import *


    

    




The autoreload extension is already loaded. To reload it, use:
  %reload_ext autoreload

In [89]:

    

#
# A = scsearch(keyword='silures',verbose=True,nonspinning=True)
# A = scsearch(keyword='q1.2_base',verbose=True)
A = scsearch(q=[9.9,11],verbose=True,nonspinning=True)


    

    




[scsearch]>> Found nonspinning (=True) keyword.
[scsearch]>> Found q (=[9.9, 11]) keyword.
[scsearch]>> Found validate_remnant (=False) keyword.
[scsearch]>> Found verbose (=True) keyword.
## Found 3 possibly degenerate simulations:
[0001][silures] q10c25eC_T_80_420: ns-q10.00	(q10c25eC_T_80_420)
[0002][athena] HR-series: ns-q10.00	(D8.4_q10.00_a0.0_m400)
[0003][sxs] SXS0185: qc-ns-q9.99	(SXS0185)

In [90]:

    

b = A[0]
y = gwylm(b,lm=[2,2],verbose=True)


    

    




(gwylm)>> Found lm (=[2, 2]) keyword.
(gwylm)>> Found load (=True) keyword.
(gwylm)>> Found scentry_obj (=<nrutils.core.nrsc.scentry instance at 0x1078b5368>) keyword.
(gwylm)>> Found verbose (=True) keyword.
(gwylm)>> The (extraction_parameter,level) is (6,6), which differs from the config values of (5,5). You have either manually input the non-config values, or the handler has set them by looking at the contents of the simulation directory. 
(gwylm.__make_lmlist__)>> The following spherical multipoles will be loaded:[(2, 2)]
(load)>> Loading: psi3col.r6.l6.l2.m2.gz
(gwylm)>> Using w22 from a PN estimate to calculate strain multipoles [see pnw0 in basics.py, and/or arxiv:1310.1528v4].
* w0(w22) = 0.063014 (this is the lower frequency used for FFI method [arxiv:1006.1632v3])
(gwylm.calchlm)>> The user should note that there is no minus sign used in front of the double time integral for strain (i.e. Eq 4 of arxiv:1006.1632). This differs from Eq 3.4 of arxiv:0707.4654v3. The net effect is a rotation of the overall polarization of pi degrees. The user should also note that there is no minus sign applied to h_cross meaning that the user must be mindful to write h_pluss-1j*h_cross when appropriate.

In [85]:

    

y.plot()


    

    
Out[85]:





[<matplotlib.axes._subplots.AxesSubplot at 0x10f2a0990>,
 <matplotlib.axes._subplots.AxesSubplot at 0x10776bbd0>,
 <matplotlib.axes._subplots.AxesSubplot at 0x10f2a0810>]