from math import sqrt

# Fungsi untuk menghitung kompensasi flow
def kompensasi_flow(Qd, Pd, Td, P, T):
    """ 
    Menghitung flow aktual berdasarkan kondisi operasi
    
    Qd : Flow desain (Nm³/h)
    Pd : Tekanan desain (barA)
    Td : Suhu desain (K)
    P  : Tekanan operasi (barA)
    T  : Suhu operasi (K)
    
    Output: Flow aktual yang telah dikompensasi (Nm³/h)
    """
    Qa = Qd * sqrt((P / Pd) * (Td / T))
    return Qa

# Contoh penggunaan
Qd = 5000  # Flow desain dalam Nm³/h
Pd = 1 + 1.01325   # Tekanan desain dalam bar absolut
Td = 25 + 273.15   # Suhu desain dalam Kelvin (25°C)

P_actual = 1.2 + 1.01325   # Tekanan operasi dalam bar absolut
T_actual = 30 + 273.15   # Suhu operasi dalam Kelvin (30°C)

Qa = kompensasi_flow(Qd, Pd, Td, P_actual, T_actual)

print(f"Flow yang telah dikompensasi: {Qa:.2f} Nm³/h")