# -*- coding: utf-8 -*-
"""
Created on Tue Dec 16 15:46:51 2025

@author: LJZ
"""

import ctREFPROP.ctREFPROP as ct

path = "C:/Program Files (x86)/REFPROP 10.0+/REFPROP"
R = ct.REFPROPFunctionLibrary(path)
R.SETPATHdll(path)
fluid = "CO2.FLD"
hfm = "HMX.BNC"
hrf = "DEF"
z = [1.0]
ierr, herr = R.SETUPdll(1, fluid, hfm, hrf)
R.SETUPdll(2, 'SI', 'SI', 'DEF')
R.SETREFdll("DEF",1,[1.0],0,0,0,0)
mw = R.WMOLdll(z)

p = 20000     # kpa
T = 391.05100419252705     # K


result = R.TPFLSHdll(T,p,z)
s_SI = result.s / mw
s = result.s
result2 = R.PSFLSHdll(p, s, z).T
rho = result.D * mw
h = result.h / mw

x = min(1, 50)