Python calcStagnationRadius Examples

Example #1

File: select-ben.py Project: hgsteggles/Torch

casa_tot_flux = np.zeros(nstars, dtype=float)
casa_peak_flux = np.zeros(nstars, dtype=float)
g_long = star_data[:,15]
g_lat = star_data[:,16]

for i in range(1, nstars + 1):
	t = age[i - 1] * 1000.0 * YR2S
	m = mcur[i - 1]
	nh = nhs[args.iden]
	logQ = float(logQ_interp(m))
	#dsun = d_sun[i - 1] * 1000.0 * PC2CM
	dsun = 5.0 * 1000.0 * PC2CM

	stro = koo.calcStromgrenRadius(logQ, nh)
	raga = koo.calcSpitzerRadius(logQ, nh, t)
	stag = koo.calcStagnationRadius(logQ, nh)
	#radius = np.minimum(raga, stag)
	radius = raga

	angsize_rad = (2.0 * radius / dsun)
	angsize_as = angsize_rad * 60.0 * 60.0 * 180.0 / math.pi
	ne = nh * (stro / radius)**(1.5)
	T = 8000.0
	f = 5.0e9
	tau = calcTau(T, f, ne, logQ)
	B = calcPlanck2(T, f) * 1.0e26
	S = B * (1.0 - math.exp(-tau))
	F = S * (0.5 * angsize_rad)**2 * 1000

	bmsize = 1.5

Example #2

File: select-harry.py Project: hgsteggles/Torch

for i in range(1, nstars + 1):
	if remaining_set[i]:
		t = t_ms[i - 1] * 1000.0 * YR2S
		m = mcur[i - 1]
		#nh = nhs[args.iden]
		nh = 4.0e4
		logQ = float(logQ_interp(m))
		#dsun = d_sun[i - 1] * 1000.0 * PC2CM
		dsun = 11.0 * 1000.0 * PC2CM

		stro = koo.calcStromgrenRadius(logQ, nh)
		raga = koo.calcSpitzerRadius(logQ, nh, t)
		if args.ben:
			stag = raga
		else:
			stag = koo.calcStagnationRadius(logQ, nh, 25.0, 8000.0)
		radius = np.minimum(raga, stag)

		#print radius * CM2PC

		angsize_rad = (2.0 * radius / dsun)
		angsize_as = angsize_rad * 60.0 * 60.0 * 180.0 / math.pi

		ne = nh * (stro / radius)**(1.5)
		T = 8000.0
		f = 5.0e9
		tau = calcTau(T, f, ne, logQ)
		B = calcPlanck2(T, f) * 1.0e26
		S = B * (1.0 - math.exp(-tau))
		F = S * (0.5 * angsize_rad)*(0.5 * angsize_rad) * 1000