diff options
| author | Aaron LI <aly@aaronly.me> | 2017-08-03 22:10:16 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2017-08-03 22:10:16 +0800 |
| commit | 7bcebf9daf5157ed249d162a75ff61d60fa98226 (patch) | |
| tree | 613a3b563a0c58663dc1632d4b290295a45ae649 /fg21sim/extragalactic/pointsources/starbursting.py | |
| parent | c6f861ee32233733b081209f4578450893816ee7 (diff) | |
| download | fg21sim-7bcebf9daf5157ed249d162a75ff61d60fa98226.tar.bz2 | |
Update against "Fnu_to_Tb_fast()" unit changes
Signed-off-by: Aaron LI <aly@aaronly.me>
Diffstat (limited to 'fg21sim/extragalactic/pointsources/starbursting.py')
| -rw-r--r-- | fg21sim/extragalactic/pointsources/starbursting.py | 11 |
1 files changed, 6 insertions, 5 deletions
diff --git a/fg21sim/extragalactic/pointsources/starbursting.py b/fg21sim/extragalactic/pointsources/starbursting.py index 8b5c597..57fd5a0 100644 --- a/fg21sim/extragalactic/pointsources/starbursting.py +++ b/fg21sim/extragalactic/pointsources/starbursting.py @@ -189,8 +189,9 @@ class StarBursting(BasePointSource): Parameters ------------ - area: `~astropy.units.Quantity` - Area of the PS, e.g., `1.0*au.sr2` + area: float + Area of the PS + Unit: [arcsec^2] freq: `~astropy.units.Quantity` Frequency, e.g., `1.0*au.MHz` @@ -217,8 +218,6 @@ class StarBursting(BasePointSource): Parameters ------------ - area: `~astropy.units.Quantity` - Area of the PS, e.g., `1.0*au.sr` freq: `~astropy.units.Quantity` Frequency, e.g., `1.0*au.MHz` @@ -230,9 +229,11 @@ class StarBursting(BasePointSource): # Tb_list num_ps = self.ps_catalog.shape[0] Tb_list = np.zeros((num_ps,)) + sr_to_arcsec2 = (np.rad2deg(1) * 3600) ** 2 # [sr] -> [arcsec^2] # Iteratively calculate Tb for i in range(num_ps): ps_area = self.ps_catalog['Area (sr)'][i] # [sr] - Tb_list[i] = self.calc_single_Tb(ps_area, freq) + area = ps_area * sr_to_arcsec2 + Tb_list[i] = self.calc_single_Tb(area, freq) return Tb_list |