python调用refprop9.1
时间: 2023-11-03 16:06:43 浏览: 128
Refprop9.1
5星 · 资源好评率100%
要在Python中调用Refprop 9.1,您需要使用Python Refprop接口。这个接口允许您在Python中调用Refprop子例程。
以下是一些步骤,以便在Python中调用Refprop 9.1:
1.下载Python Refprop接口。您可以从官方网站https://github.com/usnistgov/REFPROP-cmake下载接口。
2.下载并安装Refprop 9.1。您可以从官方网站https://www.nist.gov/srd/refprop下载和安装Refprop 9.1软件。
3.将Python Refprop接口解压缩到您的计算机上。
4.使用pip安装numpy和ctypes库。
5.在Python中导入ctypes库。
6.使用ctypes库中的cdll函数加载Refprop库。例如,如果您的Refprop库在C盘的文件夹中,您可以使用以下代码加载它:
```python
import ctypes
RP = ctypes.cdll.LoadLibrary('C:/REFPROP/REFPRP64.DLL')
```
7.在Python中调用Refprop子例程。例如,如果您想用Refprop计算水的密度和压力,您可以使用以下代码:
```python
import numpy as np
# 设置输入变量
x = np.zeros(20)
x[0] = 0.0 # 组分1的摩尔分数
x[1] = 1.0 # 组分2的摩尔分数
x[2] = 0.0 # 组分3的摩尔分数
x[3] = 0.0 # 组分4的摩尔分数
x[4] = 0.0 # 组分5的摩尔分数
x[5] = 0.0 # 组分6的摩尔分数
x[6] = 0.0 # 组分7的摩尔分数
x[7] = 0.0 # 组分8的摩尔分数
x[8] = 0.0 # 组分9的摩尔分数
x[9] = 0.0 # 组分10的摩尔分数
x[10] = 0.0 # 组分11的摩尔分数
x[11] = 0.0 # 组分12的摩尔分数
x[12] = 0.0 # 组分13的摩尔分数
x[13] = 0.0 # 组分14的摩尔分数
x[14] = 0.0 # 组分15的摩尔分数
x[15] = 0.0 # 组分16的摩尔分数
x[16] = 0.0 # 输入压力单位(0-所需单位,1-kPa,2-bar,3-psi,4-kg/cm2)
x[17] = 0.0 # 输入温度单位(0-所需单位,1-K,2-℃,3-℉,4-R)
x[18] = 0.0 # 输入参数1(根据所需功能而定)
x[19] = 0.0 # 输入参数2(根据所需功能而定)
# 设置输出变量
y = np.zeros(96)
# 调用Refprop子例程
RP.SETUP(x)
RP.TPFLSH(y, x[16], x[17])
# 输出结果
print('密度:', y[0])
print('压力:', y[1])
```
这是一个基本的例子,您可以根据需要调整输入变量和输出变量,并使用不同的Refprop子例程。
阅读全文