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