In [66]:

    

# Setup enviroment to autoreload packages
%load_ext autoreload
%autoreload 2
%matplotlib inline
# Use a linear interpolator, and a reverse greedy process
from matplotlib.pyplot import *
#
from nrutils import scsearch,gwylm


    

    




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

In [67]:

    

A = scsearch(keyword='athena',nonspinning=True,q=2,unique=True,verbose=True)
y = gwylm( A[0], lm = [2,2], verbose=True, dt = 0.5, lowpass=True )


    

    




[scsearch]>> Found keyword (='athena') keyword.
[scsearch]>> Found nonspinning (=True) keyword.
[scsearch]>> Found q (=2) keyword.
[scsearch]>> Found unique (=True) keyword.
[scsearch]>> Found verbose (=True) keyword.
(scsearch)>> List of keywords or string keyword found: ALL scentry objects matching will be passed. To pass ANY entries matching the keywords, input the keywords using an iterable of not of type list.
## Found 1 unique simulations:
[0001][athena] HR-series: ns-q2.00

(gwylm)>> Found clean (=False) keyword.
(gwylm)>> Found dt (=0.5) keyword.
(gwylm)>> Found lm (=[2, 2]) keyword.
(gwylm)>> Found load (=True) keyword.
(gwylm)>> Found lowpass (=True) keyword.
(gwylm)>> Found scentry_obj (=<nrutils.core.nrsc.scentry instance at 0x12fe7c128>) keyword.
(gwylm)>> Found verbose (=True) keyword.
(load)>> Loading: Ylm_WEYLSCAL4::Psi4r_l2_m2_r75.00.asc
(gwylm.setfields)>> Interpolating data to dt=0.500000
(gwylm.setfields)>> Interpolating data to dt=0.500000
(load)>> Re-orienting waveform phase to be consistent with internal sign convention for Psi4, where sign(dPhi/dt)=1*sign(m). Note that the internal sign convention is defined in ... nrutils/core/__init__.py as "M_RELATIVE_SIGN_CONVENTION". This message has appeared becuase the waveform is determioned to obey and sign convention: sign(dPhi/dt)=-1*sign(m).
(gwylm)>> Using w22 from a PN estimate to calculate strain multipoles [see pnw0 in basics.py, and/or arxiv:1310.1528v4].






    




---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
<ipython-input-67-ad5c72149cef> in <module>()
      1 A = scsearch(keyword='athena',nonspinning=True,q=2,unique=True,verbose=True)
----> 2 y = gwylm( A[0], lm = [2,2], verbose=True, dt = 0.5, lowpass=True )

/Users/book/JOKI/Libs/KOALA/nrutils_dev/nrutils/core/nrsc.pyc in __init__(this, scentry_obj, lm, lmax, dt, load, clean, extraction_parameter, level, w22, lowpass, verbose)
   1570         this.__lowpassfiltered__ = False
   1571         if lowpass:
-> 1572             this.lowpass()
   1573 
   1574         # Calculate strain

/Users/book/JOKI/Libs/KOALA/nrutils_dev/nrutils/core/nrsc.pyc in lowpass(this)
   2208         #
   2209         msg = 'Howdy, partner! This function is experimental and should NOT be used.'
-> 2210         error(msg,'lowpass')
   2211 
   2212         #

/Users/book/JOKI/Libs/KOALA/nrutils_dev/nrutils/core/basics.py in error(msg, fname)
    911         fname = 'error'
    912 
--> 913     raise ValueError( '('+red(fname+'!!')+')>> '+msg )
    914 
    915 

ValueError: (lowpass!!)>> Howdy, partner! This function is experimental and should NOT be used.

In [65]:

    

y.hlm[0].plot(show=True,domain='freq')


    

    
Out[65]:





([<matplotlib.axes._subplots.AxesSubplot at 0x125535f90>,
  <matplotlib.axes._subplots.AxesSubplot at 0x140bad090>,
  <matplotlib.axes._subplots.AxesSubplot at 0x1350c7c10>],
 <matplotlib.figure.Figure at 0x142110510>)

In [ ]: