import numpy as np
import pandas as pd


def buck_equation(temperature): # temp in Celsius
    saturation_vapor_pressure = 0.61121 * np.exp((18.678-temperature/234.5)*(temperature/(257.14+temperature)))
    return saturation_vapor_pressure * 1000 # KPa -> Pa

def absolute_humidity(relative_humidity, temperature):
    relative_humidity = np.array(relative_humidity)
    temperature = np.array(temperature)
    saturation_vapor_pressure = buck_equation(temperature)
    # 461.5/Kg Kelvin is specific gas constatnt
    return saturation_vapor_pressure * relative_humidity * 0.01 /(461.5 * (temperature + 273.15) ) # g/m^3

if __name__ == "__main__":
    df = pd.read_csv('data/202007010000_202308010000.csv')

# 절대 습도 구하기
# interpolation을 통해 시간 사이의 온습도 구하기