def calgb(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
b = bdeg * par.degtorad
l = ldeg * par.degtorad
zkpc = dkpc * math.sin(b)
azbcnt = g(zkpc) * math.sin(b) #s^-1
errnt = errNT(bdeg, sigb, dkpc, sigd) #s^-1
Excplbh = MWBHpl(bdeg, ldeg, dkpc) #s^-1
print("Excess_parallel_galpy(from Vp/Vs, with BH), Excess_z_NT95 = ",
Excplbh, ", ", azbcnt)
return None
def calhb(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
b = bdeg*par.degtorad
l = ldeg*par.degtorad
zkpc = dkpc*math.sin(b)
azbcnt = g(zkpc)*math.sin(b) #s^-1
errnt = errNT(bdeg, sigb, dkpc, sigd) #s^-1
ExcRforceBH = MWBHRfo(bdeg,ldeg,dkpc) #s^-1
print ("Excess_parallel(galpy-Rforce,with BH), Excess_z_NT95 = ", ExcRforceBH,", ", azbcnt)
return None;
def calha(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
global excpl, exz
b = bdeg * par.degtorad
l = ldeg * par.degtorad
zkpc = dkpc * math.sin(b)
azbcnt = g(zkpc) * math.sin(b) #s^-1
errnt = errNT(bdeg, sigb, dkpc, sigd) #s^-1
ExcRforce = MWRfo(bdeg, ldeg, dkpc)
print("Excess_parallel(galpy-Rforce,without BH), Excess_z_NT95 = ",
ExcRforce, ", ", azbcnt)
excpl = ExcRforce
exz = azbcnt
return None
def calga(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
global excpl,exz
b = bdeg*par.degtorad
l = ldeg*par.degtorad
zkpc = dkpc*math.sin(b)
azbcnt = g(zkpc)*math.sin(b) #s^-1
errnt = errNT(bdeg, sigb, dkpc, sigd) #s^-1
Excpl = MWpl(bdeg,ldeg,dkpc) #s^-1
print ("Excess_parallel_galpy(from Vp/Vs, without BH), Excess_z_NT95 = ", Excpl,", ", azbcnt)
excpl = Excpl
exz = azbcnt
return None;
def cala(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
b = bdeg * par.degtorad
l = ldeg * par.degtorad
zkpc = dkpc * math.sin(b)
adrcold = aplold(dkpc, b, l) * math.cos(b) #s^-1
errDT91 = err_DT91(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_DT91, Excess_z_NT95 = ", adrcold, ", ", azbcnt)
else:
print("Excess_parallel_DT91, Excess_z_NT95 = ", adrcold, "+/-",
errDT91, ", ", 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
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;
