From cb839cb55a39bbd1342c5b03e85e79757ab6163b Mon Sep 17 00:00:00 2001
From: Aaron LI <aaronly.me@outlook.com>
Date: Sat, 22 Oct 2016 15:50:58 +0800
Subject: galactic/snr.py: Fix wrong astropy quantity operation
When two quantities of same unit class divides, to get the correct
numerical value, whether use `.decompose().value`, or just convert to
plain Python object, e.g., `float(...)`, which will fail if the input
quantity is not dimensionless.
---
fg21sim/galactic/snr.py | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/fg21sim/galactic/snr.py b/fg21sim/galactic/snr.py
index 1230f58..380f61a 100644
--- a/fg21sim/galactic/snr.py
+++ b/fg21sim/galactic/snr.py
@@ -220,7 +220,7 @@ class SuperNovaRemnants:
"""
freq = frequency * self.freq_unit
flux = flux * self.units["flux"]
- Fnu = flux * (freq / self.catalog_flux_freq).value ** (-specindex)
+ Fnu = flux * float(freq / self.catalog_flux_freq) ** (-specindex)
omega = size[0]*self.units["size"] * size[1]*self.units["size"]
Tb = Fnu_to_Tb(Fnu, omega, freq)
return Tb.value
--
cgit v1.2.2