def calja(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
ExcRforce = MWRfo(bdeg,ldeg,dkpc) #s^-1
if zkpcm<=1.5:
print ("Excess_parallel(galpy-Rforce,without BH), Excess_z_HF04fit = ", ExcRforce,", ", azbchfl)
excpl = ExcRforce
exz = azbchfl
else:
print ("Excess_parallel(galpy-Rforce,without BH), Excess_z_HF04fit = ", ExcRforce,", ", azbchfh)
excpl = ExcRforce
exz = azbchfh
return None;
def calla(bdeg, sigb, ldeg, sigl, dkpc, sigd, Har):
b = bdeg * par.degtorad
l = ldeg * par.degtorad
zkpc = dkpc * math.sin(b)
ExcRforce = MWRfo(bdeg, ldeg, dkpc) #s^-1
Excz = MWZfo(bdeg, ldeg, dkpc) #s^-1
print(
"Excess_parallel(galpy-Rforce,without BH), Excess_z_galpy(without BH) = ",
ExcRforce, ", ", Excz)
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 Expl(ldeg, bdeg, dkpc):
excpl = MWRfo(ldeg, bdeg, dkpc) #s^-1
return -excpl