From 71c183da0193040bd89ea6be1f52666ddbda5425 Mon Sep 17 00:00:00 2001
From: Aaron LI <aly@aaronly.me>
Date: Sat, 25 Nov 2017 17:12:48 +0800
Subject: Add astro/eor_window.py: Calc total power within the defined EoR
 window

---
 astro/ps2d.py | 5 ++---
 1 file changed, 2 insertions(+), 3 deletions(-)

(limited to 'astro/ps2d.py')

diff --git a/astro/ps2d.py b/astro/ps2d.py
index 718b59a..0bd0ac4 100755
--- a/astro/ps2d.py
+++ b/astro/ps2d.py
@@ -339,9 +339,8 @@ class PS2D:
         Reference: Ref.[liu2014].Eq.(A9)
         """
         dfreq = self.dfreq  # [MHz]
-        c = ac.c.si.value  # [m/s]
-        Ez = cosmo.efunc(self.zc)
-        Hz = Ez * H0 * 1000.0  # [m/s/Mpc]
+        c = ac.c.to("km/s").value  # [km/s]
+        Hz = cosmo.H(self.zc).value  # [km/s/Mpc]
         d_z = c * (1+self.zc)**2 * dfreq / Hz / freq21cm
         return d_z
 
-- 
cgit v1.2.2