def _define_forest(ns):
global my_sp
my_name = 'my_source'
# create some node stubs for the sixpack
# first some parameters
ns.f << Meq.Parm(meq.polclog([1, 0.1, 0.01]))
ns.t << Meq.Parm(meq.polclog([0.01, 0.1, 1]))
# next the node stubs
stubI = ns['Istub'] << 1.1 * Meq.Sin(ns.f + ns.t)
stubQ = ns['Qstub'] << 2.0 * Meq.Cos(ns.f)
stubU = ns['Ustub'] << 2.1 * Meq.Sin(ns.f - 2)
stubV = ns['Vstub'] << 2.1 * Meq.Cos(ns.f - 2)
stubRA = ns['RAstub'] << 2.1 * Meq.Cos(ns.f - 2) * Meq.Sin(ns.t)
stubDec = ns['Decstub'] << 2.1 * Meq.Cos(ns.f - 2) * Meq.Sin(ns.t)
# now create the sixpack
my_sp=TDL_Sixpack.Sixpack(label=my_name,\
ns=ns, ra=stubRA,dec=stubDec,stokesI=stubI,\
stokesQ=stubQ,stokesU=stubU,stokesV=stubV)
my_sp.display()
# resolve the forest
ns.Resolve()
def _define_forest (ns): global my_sp my_name='my_source' # create some node stubs for the sixpack # first some parameters ns.f<<Meq.Parm(meq.polclog([1,0.1,0.01])) ns.t<<Meq.Parm(meq.polclog([0.01,0.1,1])) # next the node stubs stubI=ns['Istub']<<1.1*Meq.Sin(ns.f+ns.t) stubQ=ns['Qstub']<<2.0*Meq.Cos(ns.f) stubU=ns['Ustub']<<2.1*Meq.Sin(ns.f-2) stubV=ns['Vstub']<<2.1*Meq.Cos(ns.f-2) stubRA=ns['RAstub']<<2.1*Meq.Cos(ns.f-2)*Meq.Sin(ns.t) stubDec=ns['Decstub']<<2.1*Meq.Cos(ns.f-2)*Meq.Sin(ns.t) # now create the sixpack my_sp=TDL_Sixpack.Sixpack(label=my_name,\ ns=ns, ra=stubRA,dec=stubDec,stokesI=stubI,\ stokesQ=stubQ,stokesU=stubU,stokesV=stubV) my_sp.display() # resolve the forest ns.Resolve()
sixpack = Sixpack(label=punit, **pp)
if pp.has_key('parmtable'):
sixpack.parmtable(pp['parmtable'])
return sixpack
#############################################################################################
#############################################################################################
if __name__ == '__main__':
ns = NodeScope()
my_name = 'my_sixpack'
# create some node stubs for the sixpack
# first some parameters
ns.f << Meq.Parm(meq.polclog([1, 0.1, 0.01]))
ns.t << Meq.Parm(meq.polclog([0.01, 0.1, 1]))
# next the node stubs
stubI = ns['Istub'] << 1.1 * Meq.Sin(ns.f + ns.t)
stubQ = ns['Qstub'] << 2.0 * Meq.Cos(ns.f)
stubU = ns['Ustub'] << 2.1 * Meq.Sin(ns.f - 2)
stubV = ns['Vstub'] << 2.1 * Meq.Cos(ns.f - 2)
stubRA = ns['RAstub'] << 2.1 * Meq.Cos(ns.f - 2) * Meq.Sin(ns.t)
stubDec = ns['Decstub'] << 2.1 * Meq.Cos(ns.f - 2) * Meq.Sin(ns.t)
# now create the sixpack
my_sp = Sixpack(label=my_name,
ns=ns,
ra=stubRA,
dec=stubDec,
stokesI=stubI,
sixpack.parmtable(pp['parmtable'])
return sixpack
#############################################################################################
#############################################################################################
if __name__=='__main__':
ns=NodeScope()
my_name='my_sixpack'
# create some node stubs for the sixpack
# first some parameters
ns.f<<Meq.Parm(meq.polclog([1,0.1,0.01]))
ns.t<<Meq.Parm(meq.polclog([0.01,0.1,1]))
# next the node stubs
stubI=ns['Istub']<<1.1*Meq.Sin(ns.f+ns.t)
stubQ=ns['Qstub']<<2.0*Meq.Cos(ns.f)
stubU=ns['Ustub']<<2.1*Meq.Sin(ns.f-2)
stubV=ns['Vstub']<<2.1*Meq.Cos(ns.f-2)
stubRA=ns['RAstub']<<2.1*Meq.Cos(ns.f-2)*Meq.Sin(ns.t)
stubDec=ns['Decstub']<<2.1*Meq.Cos(ns.f-2)*Meq.Sin(ns.t)
# now create the sixpack
my_sp=Sixpack(label=my_name, ns=ns,
ra=stubRA, dec=stubDec,
stokesI=stubI,
stokesQ=stubQ,
stokesU=stubU,
