From 7bcebf9daf5157ed249d162a75ff61d60fa98226 Mon Sep 17 00:00:00 2001
From: Aaron LI <aly@aaronly.me>
Date: Thu, 3 Aug 2017 22:10:16 +0800
Subject: Update against "Fnu_to_Tb_fast()" unit changes
Signed-off-by: Aaron LI <aly@aaronly.me>
---
fg21sim/extragalactic/clusters/clusters.py | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
(limited to 'fg21sim/extragalactic/clusters')
diff --git a/fg21sim/extragalactic/clusters/clusters.py b/fg21sim/extragalactic/clusters/clusters.py
index 6886d84..94f00be 100644
--- a/fg21sim/extragalactic/clusters/clusters.py
+++ b/fg21sim/extragalactic/clusters/clusters.py
@@ -25,6 +25,7 @@ from ...utils.wcs import make_wcs
from ...utils.random import spherical_uniform
from ...utils.convert import Fnu_to_Tb_fast
from ...utils.grid import make_ellipse
+from ...utils.units import UnitConversions as AUC
logger = logging.getLogger(__name__)
@@ -480,7 +481,7 @@ class GalaxyClusters:
# Conversion coefficient: [ W/Hz/Mpc^2 ] => [ Jy ]
coef = 1.0502650403056097e-19
Fnu = coef * Lnu / (4*np.pi * distance**2) # [ Jy ]
- omega = size[0] * size[1] # [ deg^2 ]
+ omega = size[0] * size[1] * AUC.arcsec2deg**2 # [ arcsec^2 ]
Tb = Fnu_to_Tb_fast(Fnu, omega, freq)
return Tb
--
cgit v1.2.2