def calib(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
global excpl, exz
b = bdeg * par.degtorad
l = ldeg * par.degtorad
zkpc = dkpc * math.sin(b)
if zkpc < 0.0:
zkpcm = -zkpc
else:
zkpcm = zkpc
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
Excplbh = MWBHpl(bdeg, ldeg, dkpc) #s^-1
if zkpcm <= 1.5:
print(
"Excess_parallel_galpy(from Vp/Vs, with BH), Excess_z_HF04fit = ",
Excplbh, ", ", azbchfl)
excpl = Excplbh
exz = azbchfl
else:
print(
"Excess_parallel_galpy(from Vp/Vs, with BH), Excess_z_HF04fit = ",
Excplbh, ", ", azbchfh)
excpl = Excplbh
exz = azbchfh
return None
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 calkb(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
b = bdeg * par.degtorad
l = ldeg * par.degtorad
zkpc = dkpc * math.sin(b)
Excplbh = MWBHpl(bdeg, ldeg, dkpc) #s^-1
Exczbh = MWBHZfo(bdeg, ldeg, dkpc) #s^-1
print(
"Excess_parallel_galpy(from Vp/Vs, with BH), Excess_z_galpy(with BH) = ",
Excplbh, ", ", Exczbh)
return None