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;