From 6de8280fec59e4e0c6e3a18c8660ea765c3bedc6 Mon Sep 17 00:00:00 2001
From: Akio Taniguchi <taniguchi@a.phys.nagoya-u.ac.jp>
Date: Thu, 20 Jun 2024 15:43:55 +0000
Subject: [PATCH] #107 Fix temperature units of weather log

---
 demerge/merge/utils.py | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/demerge/merge/utils.py b/demerge/merge/utils.py
index 9ede55c..b5e4ece 100644
--- a/demerge/merge/utils.py
+++ b/demerge/merge/utils.py
@@ -447,6 +447,7 @@ def create_dems(
     readout_hdul = fits.open(readout_path, mode="readonly")
     ddbfits_hdul = fits.open(ddbfits_path, mode="readonly")
     weather_table = ascii.read(weather_path)
+    weather_table["tmperature"] += 273.15
     # 最後の1行は終端を表す意味のないデータが入っているため無視する
     antenna_table = ascii.read(antenna_path)[:-1]
     # 観測スクリプトに含まれているパラメタを抽出する
@@ -489,7 +490,7 @@ def create_dems(
         corresp=corresp,
         to=loadtype,
         T_room=lower_cabin_temp[0],
-        T_amb=np.nanmean(weather_table["tmperature"]) + 273.15,
+        T_amb=np.nanmean(weather_table["tmperature"]),
     )
 
     if loadtype == "Tsignal":