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,