(gwylm)>> Found clean (=True) keyword.
(gwylm)>> Found dt (=0.5) keyword.
(gwylm)>> Found lm (=[2, 1]) keyword.
(gwylm)>> Found scentry_obj (=<nrutils.core.nrsc.scentry instance at 0x13c23a5f0>) keyword.
(gwylm)>> Found verbose (=True) keyword.
(gwylm)>> The (extraction_parameter,level) is (2,5), 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!)>> The l=m=2 multipole will be loaded in order to determine important characteristice of all modes such as noise floor and junk radiation location.
(__make_lmlist__)>> The following spherical multipoles will be loaded:[(2, 1), (2, 2)]
(load)>> Loading: psi3col.r2.l5.l2.m1.gz
(3,)
[ 0. 0. 0.]
---------------------------------------------------------------------------
ValueError Traceback (most recent call last)
<ipython-input-181-22900deeb0ec> in <module>()
1 a = A[0]
2 ll,mm = 2,1
----> 3 Y = gwylm( a, lm=[ll,mm], clean=True, dt=0.5, verbose=True )
4 print a.mf,a.xf
5 a.mf,a.xf = Mf14067295(e.m1,e.m2,e.S1[-1]/(e.m1**2),e.S2[-1]/(e.m2**2)),jf14067295(e.m1,e.m2,e.S1[-1]/(e.m1**2),e.S2[-1]/(e.m2**2))
/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, calcstrain, verbose)
1725
1726 # Load the waveform data
-> 1727 if load==True: this.__load__(lmax=lmax,lm=lm,dt=dt)
1728
1729 # Characterize the waveform's start and store related information to this.preinspiral
/Users/book/JOKI/Libs/KOALA/nrutils_dev/nrutils/core/nrsc.pyc in __load__(this, lmax, lm, extraction_parameter, level, dt, verbose)
1880 # Load all values in __lmlist__
1881 for lm in this.__lmlist__:
-> 1882 this.load(lm=lm,dt=dt,extraction_parameter=extraction_parameter,level=level,verbose=verbose)
1883
1884 # Ensuer that all modes are the same length
/Users/book/JOKI/Libs/KOALA/nrutils_dev/nrutils/core/nrsc.pyc in load(this, lm, file_location, dt, extraction_parameter, level, output, verbose)
2029
2030 #
-> 2031 y_ = mkgwf(wfarr)
2032
2033 # ---------------------------------------------------- #
/Users/book/JOKI/Libs/KOALA/nrutils_dev/nrutils/core/nrsc.pyc in mkgwf(wfarr_)
2026 label = this.label,
2027 ref_scentry = this.__scentry__,
-> 2028 kind='$rM\psi_{%i%i}$'%(l,m) )
2029
2030 #
/Users/book/JOKI/Libs/KOALA/nrutils_dev/nrutils/core/nrsc.pyc in __init__(this, wfarr, dt, ref_scentry, l, m, extraction_parameter, kind, friend, mf, xf, m1, m2, label, preinspiral, postringdown, verbose)
891
892 # Fix nans, nonmonotinicities and jumps in time series waveform array
--> 893 wfarr = straighten_wfarr( wfarr, verbose=this.verbose )
894
895 # use the raw waveform data to define all fields
/Users/book/JOKI/Libs/KOALA/nrutils_dev/nrutils/core/basics.py in straighten_wfarr(wfarr, verbose)
1835 print wfarr.shape
1836 print wfarr
-> 1837 finite_mask = isfinite( sum( wfarr, 1 ) )
1838 if sum(finite_mask)!=len(finite_mask):
1839 if verbose: alert('Non-finite values found in waveform array. Corresponding rows will be removed.',thisfun)
/Library/Python/2.7/site-packages/numpy/core/fromnumeric.pyc in sum(a, axis, dtype, out, keepdims)
1838 else:
1839 return _methods._sum(a, axis=axis, dtype=dtype,
-> 1840 out=out, keepdims=keepdims)
1841
1842
/Library/Python/2.7/site-packages/numpy/core/_methods.pyc in _sum(a, axis, dtype, out, keepdims)
30
31 def _sum(a, axis=None, dtype=None, out=None, keepdims=False):
---> 32 return umr_sum(a, axis, dtype, out, keepdims)
33
34 def _prod(a, axis=None, dtype=None, out=None, keepdims=False):
ValueError: 'axis' entry is out of bounds