SetPVSat

def setPvsat(P,v):

    data = pd.read_excel('R134a.xlsx', sheet_name='Saturated l-v')
    P_data = np.array(data.iloc[2:,1])
    other_data = np.array(data.iloc[2:,[0,2,3,4,5,6,7,8,9]])

    fun = interp1d(P_data, other_data, axis=0)  # create linear interpolant
    values = fun(P)

    T = values[0]
    x = (v-values[1])/(values[2]-values[1])
    u = (1-x)*values[3] + x*values[4]
    h = (1-x)*values[5] + x*values[6]
    s = (1-x)*values[7] + x*values[8]

    return T, x, u, h, s