From 85fc8baf982226f1e9119b119893709b521f99ec Mon Sep 17 00:00:00 2001
From: hapipapijuris <“yamanaka.juri.h0@s.mail.nagoya-u.ac.jp”>
Date: Mon, 29 Jul 2024 01:55:10 +0000
Subject: [PATCH] #179 Add atmospheric transmission correction to daisy command

---
 decode/qlook.py | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/decode/qlook.py b/decode/qlook.py
index 0bfa305..83c2cb9 100644
--- 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)