def calfb(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
b = bdeg * par.degtorad
l = ldeg * par.degtorad
zkpc = dkpc * math.sin(b)
adrc = aplmod(dkpc, b, l) * math.cos(b) #s^-1
errReid = err_Reid14(bdeg, sigb, ldeg, sigl, dkpc, sigd) #s^-1
Excz = MWBHZfo(bdeg, ldeg, dkpc) #s^-1
print("Excess_parallel_Reid2014, Excess_z_galpy(with BH) = ", adrc, ", ",
Excz)
return None
def cald(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
global excpl, exz, errpl, errz
b = bdeg * par.degtorad
l = ldeg * par.degtorad
zkpc = dkpc * math.sin(b)
if zkpc < 0.0:
zkpcm = -zkpc
else:
zkpcm = zkpc
adrc = aplmod(dkpc, b, l) * math.cos(b) #s^-1
errReid = err_Reid14(bdeg, sigb, ldeg, sigl, dkpc, sigd) #s^-1
azbchfh = fhigh(zkpc) * math.sin(b) * 1.08100761142e-19 #s^-1
azbchfl = flow(zkpc) * math.sin(b) * 1.08100761142e-19 #s^-1
errhi = errHFhi(bdeg, sigb, dkpc, sigd) #s^-1
errlo = errHFlo(bdeg, sigb, dkpc, sigd) #s^-1
if Har == 1:
if zkpcm <= 1.5:
print("Excess_parallel_Reid2014, Excess_z_HF04fit = ", adrc, ", ",
azbchfl)
print("Vp/Vs= ", Vprat(Rpkpcfunc(dkpc, b, l, par.Rskpc)))
excpl = adrc
exz = azbchfl
else:
print("Excess_parallel_Reid2014, Excess_z_HF04fit = ", adrc, ", ",
azbchfh)
print("Vp/Vs= ", Vprat(Rpkpcfunc(dkpc, b, l, par.Rskpc)))
excpl = adrc
exz = azbchfh
else:
if zkpcm <= 1.5:
print("Excess_parallel_Reid2014, Excess_z_HF04fit = ", adrc, "+/-",
errReid, ", ", azbchfl, "+/-", errlo)
print("Vp/Vs= ", Vprat(Rpkpcfunc(dkpc, b, l, par.Rskpc)))
excpl = adrc
exz = azbchfl
errpl = errReid
errz = errlo
else:
print("Excess_parallel_Reid2014, Excess_z_HF04fit = ", adrc, "+/-",
errReid, ", ", azbchfh, "+/-", errhi)
print("Vp/Vs= ", Vprat(Rpkpcfunc(dkpc, b, l, par.Rskpc)))
excpl = adrc
exz = azbchfh
errpl = errReid
errz = errhi
return None
def calfa(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
global excpl, exz
b = bdeg * par.degtorad
l = ldeg * par.degtorad
zkpc = dkpc * math.sin(b)
adrc = aplmod(dkpc, b, l) * math.cos(b) #s^-1
errReid = err_Reid14(bdeg, sigb, ldeg, sigl, dkpc, sigd) #s^-1
Excz = MWZfo(bdeg, ldeg, dkpc) #s^-1
print("Excess_parallel_Reid2014, Excess_z_galpy(without BH) = ", adrc,
", ", Excz)
excpl = adrc
exz = Excz
return None
def calc(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
b = bdeg * par.degtorad
l = ldeg * par.degtorad
zkpc = dkpc * math.sin(b)
adrc = aplmod(dkpc, b, l) * math.cos(b) #s^-1
errReid = err_Reid14(bdeg, sigb, ldeg, sigl, dkpc, sigd) #s^-1
azbcnt = g(zkpc) * math.sin(b) #s^-1
errnt = errNT(bdeg, sigb, dkpc, sigd) #s^-1
if Har == 1:
print("Excess_parallel_Reid2014, Excess_z_NT95 = ", adrc, ", ", azbcnt)
else:
print("Excess_parallel_Reid2014, Excess_z_NT95 = ", adrc, "+/-",
errReid, ", ", azbcnt, "+/-", errnt)
return None
def calc(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
b = bdeg*par.degtorad
l = ldeg*par.degtorad
zkpc = dkpc*math.sin(b)
adrc = aplmod(dkpc,b,l)*math.cos(b) #s^-1
errReid = err_Reid14(bdeg, sigb, ldeg, sigl, dkpc, sigd) #s^-1
azbcnt = g(zkpc)*math.sin(b) #s^-1
errnt = errNT(bdeg, sigb, dkpc, sigd) #s^-1
be = (dkpc/par.Rskpc)*math.cos(b) - math.cos(l)#remove afterwards
Vp = par.Vs*Vprat(Rpkpcfunc(dkpc,b,l,par.Rskpc))#km/s
Vpms = 1000.0*Vp
Vsms = 1000.0*par.Vs
Rs = par.kpctom*par.Rskpc #m
Rp = par.kpctom*Rpkpcfunc(dkpc,b,l,par.Rskpc) #m
coslam = be*(par.Rskpc/Rpkpcfunc(dkpc,b,l,par.Rskpc))
Vp2byRp = (Vpms*Vpms)/Rp #m/ss
Vs2byRs = (Vsms*Vsms)/Rs #m/ss
if Har==1:
print ("Excess_parallel_Reid2014, Excess_z_NT95 = ", adrc,", ", azbcnt)
else:
print ("Excess_parallel_Reid2014, Excess_z_NT95 = ", adrc,"+/-",errReid, ", ", azbcnt,"+/-",errnt)
return None;