#179 Add atmospheric transmission correction to daisy command

deshima-dev · Jul 29, 2024 · 85fc8ba · 85fc8ba
1 parent da3a29d
commit 85fc8ba
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/decode/qlook.py b/decode/qlook.py
@@ -208,7 +208,8 @@ def daisy(
                 kwargs={"fill_value": "extrapolate"},
             )
         )
-        da_sub = da_on - da_base.values
+        t_atm = da_on.temperature
+        da_sub = t_atm * (da_on - da_base) / (t_atm - da_base)
 
         # make continuum series
         weight = get_chan_weight(da_off, method=chan_weight, pwv=pwv)