""" Define LISA here eta0 --- true anomaly of LISA's guiding centre xi0 --- rotational phase of the LISA array sw --- sw < 0 swaps spacecraft so LISA becomes 1->3->2->1 viewed from above t0 --- the time at which LISA has the above """ eta0 , xi0 , sw , t0 = 0 , 0 , 1 , 0 lisa = synthlisa.EccentricInclined( eta0 , xi0 , sw , t0 ) "Define the sky" nlon , nlat = 40 , 21 sky = mlisar.mySpharmt( nlon , nlat ) "Which first & last days?" days = [ 180 ] "Specify the function over the sky" whichfunc = 'tdiORF' tdiI , tdiJ = ( 'Michelson' , 'G2' , 'A' , '1' ) , ( 'Michelson' , 'G2' , 'E' , '1' ) "Which frequencies?" f0 , df , Nf = 1e-5 , 1e-5 , 99999 f = f0 + df * np.arange( Nf ) "lmax" lmax = nlat - 1
import synthlisa
import myLISAmodule as mlisar
parser = optparse.OptionParser('usage: ./x_tdiORF_SpHs.py SETUP.pkl')
(options, args) = parser.parse_args()
if len(args) < 1:
parser.error('You must specify SETUP.pkl containing input parameters!')
file = open(args[0], 'rb')
setup = cpkl.load(file)
file.close()
if setup['lisa']['type'] == 'EccentricInclined':
eta0, xi0, sw, lisa_t0 = setup['lisa']['initial conditions']
lisa = synthlisa.EccentricInclined(eta0, xi0, sw, lisa_t0)
sky = mlisar.mySpharmt(setup['sky']['nlon'], setup['sky']['nlat'])
t0 = setup['t0']
tdiI, tdiJ = setup['tdiI'], setup['tdiJ']
f0, df, Nf = setup['f0'], setup['df'], setup['Nf']
lmax = setup['lmax']
orfdir = setup['orfdir']
orfpath = orfdir + 'orf_t0_%f.pkl' % t0
if orfpath not in glob.glob(orfpath):
f = f0 + df * np.arange(Nf)
lisky = mlisar.LISA_in_the_Sky(lisa, sky)
# SpHreal , SpHimag = lisky.get_SpHs( t = t0 , *( lmax , 'tdiORF' , tdiI , tdiJ , f ) )
print lmax, tdiI, tdiJ, f, t0
SpHreal, SpHimag = lisky.get_SpHs(lmax, 'tdiORF', tdiI, tdiJ, f, t=t0)
print SpHreal
orfdict = {
import synthlisa
import myLISAmodule as mlisar
parser = optparse.OptionParser( 'usage: ./x_tdiORF_SpHs.py SETUP.pkl' )
( options , args ) = parser.parse_args()
if len( args ) < 1 :
parser.error( 'You must specify SETUP.pkl containing input parameters!' )
file = open( args[0] , 'rb' ) ; setup = cpkl.load( file ) ; file.close()
if setup['lisa']['type'] == 'EccentricInclined' :
eta0 , xi0 , sw , lisa_t0 = setup['lisa']['initial conditions']
lisa = synthlisa.EccentricInclined( eta0 , xi0 , sw , lisa_t0 )
sky = mlisar.mySpharmt( setup['sky']['nlon'] , setup['sky']['nlat'] )
t0 = setup['t0']
tdiI , tdiJ = setup['tdiI'] , setup['tdiJ']
f0 , df , Nf = setup['f0'] , setup['df'] , setup['Nf']
lmax = setup['lmax']
orfdir = setup['orfdir']
orfpath = orfdir + 'orf_t0_%f.pkl' % t0
if orfpath not in glob.glob( orfpath ) :
f = f0 + df * np.arange( Nf )
lisky = mlisar.LISA_in_the_Sky( lisa , sky )
# SpHreal , SpHimag = lisky.get_SpHs( t = t0 , *( lmax , 'tdiORF' , tdiI , tdiJ , f ) )
print lmax , tdiI , tdiJ , f , t0
SpHreal , SpHimag = lisky.get_SpHs( lmax , 'tdiORF' , tdiI , tdiJ , f , t = t0 )
print SpHreal
orfdict = {'OrfMultipleMoments':
#from LISAresponse import *
import synthlisa
import myLISAmodule as mlisar
"""
Define LISA here
eta0 --- true anomaly of LISA's guiding centre
xi0 --- rotational phase of the LISA array
sw --- sw < 0 swaps spacecraft so LISA becomes 1->3->2->1 viewed from above
t0 --- the time at which LISA has the above
"""
eta0, xi0, sw, t0 = 0, 0, 1, 0
lisa = synthlisa.EccentricInclined(eta0, xi0, sw, t0)
"Define the sky"
nlon, nlat = 40, 21
sky = mlisar.mySpharmt(nlon, nlat)
"Which first & last days?"
days = [180]
"Specify the function over the sky"
whichfunc = 'tdiORF'
tdiI, tdiJ = ('Michelson', 'G2', 'A', '1'), ('Michelson', 'G2', 'E', '1')
"Which frequencies?"
f0, df, Nf = 1e-5, 1e-5, 99999
f = f0 + df * np.arange(Nf)
"lmax"
lmax = nlat - 1