# readselrad.py
# 2026-05-15 masudako
# 熊谷キャンパスの熱収支観測データを読む
# 放射フラックスの値を書き出す
#
# パッケージを予約
import pandas as pd
import datetime
#
# 抽出する日時の範囲
dtstart = datetime.datetime(2017, 10,  1, 12, 0, 0)
dtend   = dtstart + datetime.timedelta(days=8)
#
# データをファイルからデータフレームに読みこむ
indatapath = '/raid/kumagaya_campus/obsdata/'
infilename = 'surface_TBL110_20161103_20251201.dat'
infilepathname = indatapath+infilename
df = pd.read_csv(infilepathname, header=0, skiprows=[0,2,3])
nrec = len(df)
#
# 第0列を日付型にする
for krec in range(nrec):
    df.iloc[krec,0] = datetime.datetime.strptime(df.iloc[krec,0], '%Y-%m-%d %H:%M:%S')
#
# 行を日時で抽出する
dfwork = df[        df.TIMESTAMP <= dtend]
dfsub  = dfwork[dfwork.TIMESTAMP >= dtstart]
#
# 列を抽出する
dfsubsub = dfsub.iloc[:, [0,1,2,3,4,5,6,7]]
#
# 抽出結果を csv ファイルに書き出す
#print(dfsubsub)
dfsubsub.to_csv("readselradout.csv", index=False)
##### END OF PROGRAM #####