From cb839cb55a39bbd1342c5b03e85e79757ab6163b Mon Sep 17 00:00:00 2001 From: Aaron LI 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(-) (limited to 'fg21sim') 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