1. Critical curve

Calculate critical curve of H2O-NaCl system based on IAPS84 and IAPWS95 EOS, and compare result with [Driesner, 2007].

Key code snippet

import numpy as np
from xThermo import H2O
from xThermo import H2ONaCl
sw = H2ONaCl.cH2ONaCl("IAPS84") #IAPS84, IAPWS95
T=np.linspace(sw.Tmin(), sw.Tmax(), 200) # Temperature array
P,X=np.array(sw.P_X_Critical(T)) # calculate critical pressure and composition for given temperature

See also

The C++ member function P_X_Critical and python wrapper function xThermo.H2ONaCl.cH2ONaCl.P_X_Critical().

Python code

import numpy as np
import time
import linecache
import matplotlib.pyplot as plt
import matplotlib as mpl
from matplotlib import patches
import matplotlib.ticker as ticker
from matplotlib.ticker import MultipleLocator
from tabulate import tabulate
import warnings
warnings.filterwarnings("ignore", category=RuntimeWarning)
# 3d plot
import helpfunc
mpl.rcParams['font.family'] = 'Arial'  # default font family
mpl.rcParams['mathtext.fontset'] = 'cm'  # font for math
fmt_figs = ['pdf']  # ['svg','pdf']
figpath = '.'
result_path='../../../gallery_H2ONaCl/pT'
def savefig(figname):
    for fmt_fig in fmt_figs:
        figname_full = '%s/%s.%s' % (figpath, figname, fmt_fig)
        plt.savefig(figname_full, bbox_inches='tight')
        print('figure saved: ', figname_full)
compare = lambda a,b : float(str('%.6e'%(a)))-float(str('%.6e'%(b)))
# Import package of xThermo
from xThermo import H2O
from xThermo import H2ONaCl
sw_84 = H2ONaCl.cH2ONaCl("IAPS84")
sw_95 = H2ONaCl.cH2ONaCl("IAPWS95")

def plot_3d():
    fig=plt.figure(figsize=(14,14))
    ax = fig.add_subplot(111,projection='3d',facecolor='None')
    helpfunc.set_axis_diagram_3D(ax)
    # calculate critical curve
    sw=H2ONaCl.cH2ONaCl('IAPS84')
    T=np.linspace(sw.Tmin(), sw.Tmax(), 200)
    P,X=np.array(sw.P_X_Critical(T))
    # plot
    ax.plot(X*100, T-273.15, P/1E5, label='Critical curve', color='r')

    savefig('CriticalCurve_3D')

def benchmark_CriticalCurve(sw,mmc2='../Driesner2007b/1-s2.0-S0016703707002955-mmc2.txt'):
    T=np.linspace(sw.Tmin(), sw.Tmax(), 200)
    P,X=np.array(sw.P_X_Critical(T))
    T0,P0,X0=[],[],[]
    # try to open Driesner's result file
    data=np.loadtxt(mmc2,skiprows=6)
    T0,P0,X0,Rho0,H0=data[:,0],data[:,1],data[:,2],data[:,3],data[:,4]
    P_,X_wt_=np.array(sw.P_X_Critical(T0+273.15))
    X_=np.array(sw.Wt2Mol(X_wt_))
    Rho_=np.array(sw.Rho_phase(T0+273.15, P_, X_wt_, H2ONaCl.Liquid))
    H_=np.array(sw.H_phase(T0+273.15, P_, X_wt_, H2ONaCl.Liquid))
    Data0 = {'p':P0,'X':X0,'rho':Rho0,'h':H0}
    Data_={'p':P_/1E5,'X':X_,'rho':Rho_,'h':H_}
    Err,RErr={},{}
    for key in Data0.keys(): Err[key],RErr[key] = Data0[key]-Data_[key], np.abs(Data0[key]-Data_[key])/(Data0[key])*100.0
    # print to file and compare
    fpout = open('%s/mmc2_%s.csv'%(result_path,sw.name_backend()),'w')
    fpout.write('T[deg.C],P(Driesner)[bar],P(xThermo)[bar],P(diff)[bar],X(Driesner)[mol],X(xThermo)[mol],X(diff)[mol],Rho(Driesner)[kg/m3],Rho(xThermo),Rho(err),H(Driesner)[J/kg],H(xThermo),H(err)\n')
    for i in range(0,len(T0)):
        fpout.write('%.6e'%(T0[i]))
        for key in Data0.keys():
            fpout.write(',%.6e,%.6e,%.6e'%(Data0[key][i], Data_[key][i],compare(Data0[key][i],Data_[key][i])))
        fpout.write('\n')
    fpout.close()

    # plot
    fig,axes=plt.subplots(1,3, figsize=(12,3),gridspec_kw={'width_ratios':[1, 1, 1],'wspace':0.3})
    axes_err=['']*len(axes)
    for i,ax in enumerate(axes):
        axes_err[i]=ax.inset_axes([0,1,1,0.3])
        axes_err[i].set_yscale('log')
        # axes_err[i].set_ylim(1E-6,5E-3)
        axes_err[i].set_ylabel('RE(%)')
        axes_err[i].text(-0.08,1.0,'(%s)'%(chr(97+i)),ha='right',va='bottom',fontsize=10,fontweight='bold',transform=axes_err[i].transAxes)
    # 1. critical curve
    ax,ax_err=axes[0],axes_err[0]
    l,=ax.plot(T-273.15,P/1E5, color='r',label='xThermo',lw=4)
    if(len(T0)>0): ax.plot(T0,P0,ls='dashed',color='lightgray',label='Driesner(2007b)')
    ax.set_xlabel('Temperature ($^{\circ}$C)')
    ax.set_ylabel('Critical pressure (bar)',color=l.get_color())
    ax_err.bar(T0+2.5, RErr['p'],color=l.get_color(),width=5)
    # ax.yaxis.set_minor_locator(MultipleLocator(100))
    ax.legend(loc='upper left')
    # ax.text(0.005,0.98,'(a)',transform=ax.transAxes,ha='left',va='top',fontsize=12,fontweight='bold')
    ax_X = ax.twinx()
    l,=ax_X.plot(T-273.15,X*100,label='blue',lw=4)
    ax_X.set_ylim(0, ax_X.get_ylim()[1])
    if(len(T0)>0): ax_X.plot(T0,np.array(sw.Mol2Wt(X0))*100,ls='dashed',color='lightgray')
    ax_X.set_ylabel('Critical composition (wt.% NaCl)',color=l.get_color())
    ax_err.bar(T0-2.5, RErr['X'],color=l.get_color(),width=5)
    # 2. density
    ax,ax_err=axes[1],axes_err[1]
    ax.plot(T0, Rho_,lw=4,label='xThermo')
    ax.plot(T0, Rho0, ls='dashed',color='lightgray',label='Driesner(2007b)')
    ax_err.bar(T0, RErr['rho'],width=5)
    ax.yaxis.set_ticks_position('right')
    ax.yaxis.set_label_position('right')
    ax.set_ylabel('Density (km/m$^{\mathregular{3}}$)')
    ax.set_xlabel('Temperature ($^{\circ}$C)')
    ax.legend()
    # 3. specific enthalpy
    ax,ax_err=axes[2],axes_err[2]
    ax.plot(T0, H_/1E6,lw=4,label='xThermo')
    ax.plot(T0, H0/1E6, ls='dashed',color='lightgray',label='Driesner(2007b)')
    ax.yaxis.set_ticks_position('right')
    ax.yaxis.set_label_position('right')
    ax.set_ylabel('Specific enthalpy (MJ/kg)')
    ax.set_xlabel('Temperature ($^{\circ}$C)')
    ax_err.bar(T0, RErr['h'],width=5)
    ax.legend()
    # set axes
    for i,ax in enumerate(axes):
        axes_err[i].set_xlim(ax.get_xlim())
        axes_err[i].xaxis.set_ticks([])
        ax.xaxis.set_minor_locator(MultipleLocator(40))
        ax.grid(which='major',lw=0.02,color='k')
        ax.grid(which='minor',lw=0.01,color='gray')

    # statistics of the difference
    table=[]
    for key,name in zip(list(Err.keys()),['Pressure (bar)','Composition (mole fraction)','Density (kg/m3)','Specific enthalpy (J/kg)']):
        RErr[key] = RErr[key][~(np.isnan(RErr[key]) | np.isinf(RErr[key]))]
        table.append([name,Err[key].min(),Err[key].max(),RErr[key].min(),RErr[key].max()])
    print(tabulate(table, headers=['Critical property', 'Err. Min', 'Err. Max','RE. Min(%)','RE. Max(%)'],numalign="right",floatfmt=".6f"))
    savefig('H2ONaCl_CriticalCurve')
    return T0,P0

def plot_critpT_H2O(T_crit,P_crit,cmap='GnBu'):
    T,p = np.linspace(274,T_crit.max(),100), np.linspace(1E5,P_crit.max(),100)
    TT,PP = np.meshgrid(T,p)
    prop = np.zeros_like(TT)
    iaps84 = H2O.cIAPS84()
    iapws95 = H2O.cIAPWS95_CoolProp()
    for i in range(0,TT.shape[0]):
        for j in range(0,TT.shape[1]):
            props = iapws95.UpdateState_TPX(TT[i][j], PP[i][j])
            prop[i][j] = props.Rho
    fig,axes=plt.subplots(1,2,figsize=(12,4),gridspec_kw={'wspace':0.25})
    ax=axes[0]
    CS = ax.contourf(TT-273.15,PP/1E5, prop, levels=50, cmap=cmap)
    ax_cb = ax.inset_axes([1.01, 0, 0.05, 1])
    plt.colorbar(CS, cax=ax_cb, orientation='vertical',label='%s (%s)'%('Density','kg/$m^3$'))
    l,=ax.plot(T_crit-273.15, P_crit/1E5,color='r',lw=2,label='Critical p,T of $H_{\mathregular{2}}$O-NaCl system')
    x0,y0=T_crit[int(len(T_crit)/3*2)]-273.15, P_crit[int(len(P_crit)/3*2)]/1E5
    ax.annotate("Critical p,T of\nH$_{\mathregular{2}}$O-NaCl system",
                xy=(x0,y0),xytext=(x0-100,y0), ha='right',va='center',bbox={'fc':'None','ec':l.get_color()},fontsize=14,fontweight='bold',
                arrowprops=dict(arrowstyle="->",connectionstyle="arc3"),)
    # labels
    # ax.text(0.98,0.98,water.name(),ha='right',va='top',bbox={'fc':'w','ec':'gray'}, transform=ax.transAxes)
    ax.set_xlabel('Temperature ($^{\circ}$C)')
    ax.set_ylabel('Pressure (bar)')
    # IAPS84 - IAPWS95
    rho_84,rho_95,h_84,h_95=np.zeros_like(T_crit),np.zeros_like(T_crit),np.zeros_like(T_crit),np.zeros_like(T_crit)
    for i in range(0,len(T_crit)):
        props_84 = iaps84.UpdateState_TPX(T_crit[i],P_crit[i])
        props_95 = iapws95.UpdateState_TPX(T_crit[i],P_crit[i])
        rho_84[i],rho_95[i],h_84[i],h_95[i] = props_84.Rho, props_95.Rho, props_84.H, props_95.H
    ax=axes[1]
    barwidth=5
    ax.set_yscale('log')
    ax.bar(T_crit-273.15-barwidth/2,np.abs(rho_84-rho_95)/rho_84*100,width=barwidth,label='Density')
    ax.bar(T_crit-273.15+barwidth/2,np.abs(h_84-h_95)/h_84*100,width=barwidth,label='Specific enthalpy')
    ax.legend(ncol=2,title='Relative difference: $\\frac{|IAPS84-IAPWS95|}{IAPS84}\\times100$')
    ax.set_ylabel('Relative difference (%)')
    ax.set_xlabel('Temperature ($^{\circ}$C)')
    ax.yaxis.set_label_position('right')
    ax.yaxis.set_ticks_position('right')
    savefig('CriticalCurve_on_H2O')

General 3D view

plot_3d()

Out:

figure saved:  ./CriticalCurve_3D.pdf

Based on IAPS84 EOS

T_crit,P_crit=benchmark_CriticalCurve(sw_84)

Out:

Critical property              Err. Min    Err. Max    RE. Min(%)    RE. Max(%)
---------------------------  ----------  ----------  ------------  ------------
Pressure (bar)                -0.000322    0.000713      0.000000      0.000043
Composition (mole fraction)   -0.000000    0.000000      0.000000      0.000044
Density (kg/m3)               -0.000203    0.003298      0.000000      0.000505
Specific enthalpy (J/kg)      -1.270571  114.473760      0.000000      0.005488
figure saved:  ./H2ONaCl_CriticalCurve.pdf

Based on IAPWS95 EOS

T_crit,P_crit=benchmark_CriticalCurve(sw_95)

Out:

Critical property                 Err. Min     Err. Max    RE. Min(%)    RE. Max(%)
---------------------------  -------------  -----------  ------------  ------------
Pressure (bar)                   -0.193065     0.000713      0.000000      0.073312
Composition (mole fraction)      -0.000013     0.000000      0.000000    125.840340
Density (kg/m3)                  -5.513025     2.734699      0.002108      1.589922
Specific enthalpy (J/kg)     -33523.582505  7377.508972      0.000176      1.607075
figure saved:  ./H2ONaCl_CriticalCurve.pdf

H2O properties at critical p,T condition of H2O-NaCl

plot_critpT_H2O(T_crit+273.15,P_crit*1E5)

Out:

figure saved:  ./CriticalCurve_on_H2O.pdf

Result table

See also

Electronic Annex EA-2 in [Driesner, 2007] and Fig. 1,5,6 in [Driesner & Heinrich, 2007].

Result data calculated by xThermo based on both water EOS of IAPS84 and IAPWS95.

IAPS84

Table 14 Comparison between result of [Driesner, 2007] and result calculated by xThermo.

T[deg.C]	P(Driesner)[bar]	P(xThermo)[bar]	P(diff)[bar]	X(Driesner)[mol]	X(xThermo)[mol]	X(diff)[mol]	Rho(Driesner)[kg/m3]	Rho(xThermo)	Rho(err)	H(Driesner)[J/kg]	H(xThermo)	H(err)
3.739760e+02	2.205491e+02	2.205491e+02	0.000000e+00	0.000000e+00	0.000000e+00	0.000000e+00	3.218900e+02	3.218900e+02	0.000000e+00	2.086000e+06	2.085886e+06	1.140000e+02
3.740000e+02	2.206058e+02	2.206058e+02	0.000000e+00	1.925758e-06	1.925758e-06	0.000000e+00	3.220024e+02	3.220025e+02	-1.000000e-04	2.090401e+06	2.090400e+06	1.000000e+00
3.741000e+02	2.208416e+02	2.208416e+02	0.000000e+00	1.007350e-05	1.007350e-05	0.000000e+00	3.226283e+02	3.226277e+02	6.000000e-04	2.111710e+06	2.111711e+06	-1.000000e+00
3.742000e+02	2.210772e+02	2.210772e+02	0.000000e+00	1.842022e-05	1.842022e-05	0.000000e+00	3.234210e+02	3.234210e+02	0.000000e+00	2.127722e+06	2.127722e+06	0.000000e+00
3.743000e+02	2.213126e+02	2.213126e+02	0.000000e+00	2.696511e-05	2.696511e-05	0.000000e+00	3.243026e+02	3.243028e+02	-2.000000e-04	2.133951e+06	2.133949e+06	2.000000e+00
3.744000e+02	2.215478e+02	2.215478e+02	0.000000e+00	3.570734e-05	3.570734e-05	0.000000e+00	3.252234e+02	3.252234e+02	0.000000e+00	2.136659e+06	2.136659e+06	0.000000e+00
3.745000e+02	2.217828e+02	2.217828e+02	0.000000e+00	4.464610e-05	4.464610e-05	0.000000e+00	3.261537e+02	3.261538e+02	-1.000000e-04	2.137989e+06	2.137987e+06	2.000000e+00
3.746000e+02	2.220176e+02	2.220176e+02	0.000000e+00	5.378059e-05	5.378059e-05	0.000000e+00	3.270754e+02	3.270755e+02	-1.000000e-04	2.138636e+06	2.138636e+06	0.000000e+00
3.747000e+02	2.222522e+02	2.222522e+02	0.000000e+00	6.310998e-05	6.310998e-05	0.000000e+00	3.279778e+02	3.279778e+02	0.000000e+00	2.138899e+06	2.138899e+06	0.000000e+00
3.748000e+02	2.224866e+02	2.224866e+02	0.000000e+00	7.263348e-05	7.263348e-05	0.000000e+00	3.288547e+02	3.288547e+02	0.000000e+00	2.138922e+06	2.138922e+06	0.000000e+00
3.749000e+02	2.227209e+02	2.227209e+02	0.000000e+00	8.235027e-05	8.235027e-05	0.000000e+00	3.297035e+02	3.297030e+02	5.000000e-04	2.138788e+06	2.138789e+06	-1.000000e+00
3.750000e+02	2.229550e+02	2.229550e+02	0.000000e+00	9.225956e-05	9.225956e-05	0.000000e+00	3.305218e+02	3.305219e+02	-1.000000e-04	2.138547e+06	2.138547e+06	0.000000e+00
3.751000e+02	2.231890e+02	2.231890e+02	0.000000e+00	1.023605e-04	1.023605e-04	0.000000e+00	3.313108e+02	3.313108e+02	0.000000e+00	2.138230e+06	2.138230e+06	0.000000e+00
3.756000e+02	2.243566e+02	2.243566e+02	0.000000e+00	1.557130e-04	1.557130e-04	0.000000e+00	3.348582e+02	3.348582e+02	0.000000e+00	2.136103e+06	2.136103e+06	0.000000e+00
3.761000e+02	2.255209e+02	2.255209e+02	0.000000e+00	2.137389e-04	2.137389e-04	0.000000e+00	3.379003e+02	3.379003e+02	0.000000e+00	2.133750e+06	2.133750e+06	0.000000e+00
3.766000e+02	2.266822e+02	2.266822e+02	0.000000e+00	2.763407e-04	2.763407e-04	0.000000e+00	3.406247e+02	3.406246e+02	1.000000e-04	2.131521e+06	2.131521e+06	0.000000e+00
3.771000e+02	2.278411e+02	2.278411e+02	0.000000e+00	3.434222e-04	3.434222e-04	0.000000e+00	3.431574e+02	3.431574e+02	0.000000e+00	2.129505e+06	2.129505e+06	0.000000e+00
3.776000e+02	2.289976e+02	2.289976e+02	0.000000e+00	4.148884e-04	4.148884e-04	0.000000e+00	3.455768e+02	3.455768e+02	0.000000e+00	2.127703e+06	2.127703e+06	0.000000e+00
3.781000e+02	2.301523e+02	2.301523e+02	0.000000e+00	4.906457e-04	4.906457e-04	0.000000e+00	3.479295e+02	3.479295e+02	0.000000e+00	2.126086e+06	2.126086e+06	0.000000e+00
3.786000e+02	2.313052e+02	2.313052e+02	0.000000e+00	5.706021e-04	5.706021e-04	0.000000e+00	3.502430e+02	3.502430e+02	0.000000e+00	2.124621e+06	2.124621e+06	0.000000e+00
3.791000e+02	2.324566e+02	2.324566e+02	0.000000e+00	6.546664e-04	6.546664e-04	0.000000e+00	3.525330e+02	3.525330e+02	0.000000e+00	2.123277e+06	2.123277e+06	0.000000e+00
3.796000e+02	2.336069e+02	2.336069e+02	0.000000e+00	7.427492e-04	7.427492e-04	0.000000e+00	3.548078e+02	3.548078e+02	0.000000e+00	2.122029e+06	2.122029e+06	0.000000e+00
3.810000e+02	2.368226e+02	2.368226e+02	0.000000e+00	1.010012e-03	1.010012e-03	0.000000e+00	3.611200e+02	3.611199e+02	1.000000e-04	2.118897e+06	2.118897e+06	0.000000e+00
3.820000e+02	2.391168e+02	2.391168e+02	0.000000e+00	1.218734e-03	1.218734e-03	0.000000e+00	3.655758e+02	3.655758e+02	0.000000e+00	2.116880e+06	2.116879e+06	1.000000e+00
3.830000e+02	2.414102e+02	2.414102e+02	0.000000e+00	1.441552e-03	1.441552e-03	0.000000e+00	3.699761e+02	3.699761e+02	0.000000e+00	2.114985e+06	2.114985e+06	0.000000e+00
3.840000e+02	2.437038e+02	2.437038e+02	0.000000e+00	1.677822e-03	1.677822e-03	0.000000e+00	3.743100e+02	3.743100e+02	0.000000e+00	2.113184e+06	2.113184e+06	0.000000e+00
3.850000e+02	2.459988e+02	2.459988e+02	0.000000e+00	1.926920e-03	1.926920e-03	0.000000e+00	3.785689e+02	3.785689e+02	0.000000e+00	2.111460e+06	2.111460e+06	0.000000e+00
3.860000e+02	2.482960e+02	2.482960e+02	0.000000e+00	2.188239e-03	2.188239e-03	0.000000e+00	3.827470e+02	3.827469e+02	1.000000e-04	2.109803e+06	2.109803e+06	0.000000e+00
3.870000e+02	2.505964e+02	2.505964e+02	0.000000e+00	2.461191e-03	2.461191e-03	0.000000e+00	3.868403e+02	3.868403e+02	0.000000e+00	2.108209e+06	2.108209e+06	0.000000e+00
3.880000e+02	2.529007e+02	2.529007e+02	0.000000e+00	2.745207e-03	2.745207e-03	0.000000e+00	3.908469e+02	3.908469e+02	0.000000e+00	2.106674e+06	2.106674e+06	0.000000e+00
3.890000e+02	2.552097e+02	2.552097e+02	0.000000e+00	3.039736e-03	3.039736e-03	0.000000e+00	3.947660e+02	3.947659e+02	1.000000e-04	2.105198e+06	2.105198e+06	0.000000e+00
3.900000e+02	2.575241e+02	2.575241e+02	0.000000e+00	3.344241e-03	3.344241e-03	0.000000e+00	3.985976e+02	3.985975e+02	1.000000e-04	2.103781e+06	2.103781e+06	0.000000e+00
3.910000e+02	2.598445e+02	2.598445e+02	0.000000e+00	3.658206e-03	3.658206e-03	0.000000e+00	4.023425e+02	4.023424e+02	1.000000e-04	2.102421e+06	2.102421e+06	0.000000e+00
3.920000e+02	2.621714e+02	2.621714e+02	0.000000e+00	3.981129e-03	3.981129e-03	0.000000e+00	4.060019e+02	4.060018e+02	1.000000e-04	2.101120e+06	2.101120e+06	0.000000e+00
3.930000e+02	2.645056e+02	2.645056e+02	0.000000e+00	4.312524e-03	4.312524e-03	0.000000e+00	4.095772e+02	4.095771e+02	1.000000e-04	2.099877e+06	2.099877e+06	0.000000e+00
3.940000e+02	2.668474e+02	2.668474e+02	0.000000e+00	4.651924e-03	4.651924e-03	0.000000e+00	4.130703e+02	4.130702e+02	1.000000e-04	2.098693e+06	2.098693e+06	0.000000e+00
3.950000e+02	2.691974e+02	2.691974e+02	0.000000e+00	4.998874e-03	4.998874e-03	0.000000e+00	4.164830e+02	4.164829e+02	1.000000e-04	2.097568e+06	2.097568e+06	0.000000e+00
3.960000e+02	2.715561e+02	2.715561e+02	0.000000e+00	5.352936e-03	5.352936e-03	0.000000e+00	4.198173e+02	4.198172e+02	1.000000e-04	2.096501e+06	2.096501e+06	0.000000e+00
3.970000e+02	2.739237e+02	2.739237e+02	0.000000e+00	5.713687e-03	5.713687e-03	0.000000e+00	4.230752e+02	4.230750e+02	2.000000e-04	2.095492e+06	2.095492e+06	0.000000e+00
3.980000e+02	2.763009e+02	2.763009e+02	0.000000e+00	6.080717e-03	6.080717e-03	0.000000e+00	4.262587e+02	4.262585e+02	2.000000e-04	2.094542e+06	2.094542e+06	0.000000e+00
3.990000e+02	2.786878e+02	2.786878e+02	0.000000e+00	6.453632e-03	6.453632e-03	0.000000e+00	4.293698e+02	4.293697e+02	1.000000e-04	2.093650e+06	2.093650e+06	0.000000e+00
4.000000e+02	2.810849e+02	2.810849e+02	0.000000e+00	6.832051e-03	6.832051e-03	0.000000e+00	4.324106e+02	4.324105e+02	1.000000e-04	2.092815e+06	2.092815e+06	0.000000e+00
4.100000e+02	3.056780e+02	3.056780e+02	0.000000e+00	1.084455e-02	1.084455e-02	0.000000e+00	4.593637e+02	4.593635e+02	2.000000e-04	2.087488e+06	2.087488e+06	0.000000e+00
4.200000e+02	3.315379e+02	3.315379e+02	0.000000e+00	1.508056e-02	1.508056e-02	0.000000e+00	4.811573e+02	4.811570e+02	3.000000e-04	2.087125e+06	2.087125e+06	0.000000e+00
4.300000e+02	3.587626e+02	3.587626e+02	0.000000e+00	1.934403e-02	1.934403e-02	0.000000e+00	4.990526e+02	4.990521e+02	5.000000e-04	2.090916e+06	2.090916e+06	0.000000e+00
4.400000e+02	3.873578e+02	3.873578e+02	0.000000e+00	2.351927e-02	2.351927e-02	0.000000e+00	5.139825e+02	5.139819e+02	6.000000e-04	2.098093e+06	2.098093e+06	0.000000e+00
4.500000e+02	4.172612e+02	4.172612e+02	0.000000e+00	2.754975e-02	2.754975e-02	0.000000e+00	5.266585e+02	5.266579e+02	6.000000e-04	2.107963e+06	2.107963e+06	0.000000e+00
4.600000e+02	4.483573e+02	4.483573e+02	0.000000e+00	3.142030e-02	3.142030e-02	0.000000e+00	5.376330e+02	5.376322e+02	8.000000e-04	2.119912e+06	2.119912e+06	0.000000e+00
4.700000e+02	4.804883e+02	4.804883e+02	0.000000e+00	3.514237e-02	3.514237e-02	0.000000e+00	5.473332e+02	5.473324e+02	8.000000e-04	2.133417e+06	2.133417e+06	0.000000e+00
4.800000e+02	5.134606e+02	5.134606e+02	0.000000e+00	3.874218e-02	3.874218e-02	0.000000e+00	5.560821e+02	5.560812e+02	9.000000e-04	2.148044e+06	2.148044e+06	0.000000e+00
4.900000e+02	5.470509e+02	5.470509e+02	0.000000e+00	4.225142e-02	4.225142e-02	0.000000e+00	5.641113e+02	5.641103e+02	1.000000e-03	2.163459e+06	2.163459e+06	0.000000e+00
5.000000e+02	5.810101e+02	5.810101e+02	0.000000e+00	4.570045e-02	4.570045e-02	0.000000e+00	5.715720e+02	5.715709e+02	1.100000e-03	2.179422e+06	2.179422e+06	0.000000e+00
5.100000e+02	6.150162e+02	6.150161e+02	1.000000e-04	4.911344e-02	4.911344e-02	0.000000e+00	5.785204e+02	5.785191e+02	1.300000e-03	2.195824e+06	2.195824e+06	0.000000e+00
5.200000e+02	6.489246e+02	6.489245e+02	1.000000e-04	5.250551e-02	5.250551e-02	0.000000e+00	5.850483e+02	5.850469e+02	1.400000e-03	2.212499e+06	2.212499e+06	0.000000e+00
5.300000e+02	6.827353e+02	6.827352e+02	1.000000e-04	5.588146e-02	5.588146e-02	0.000000e+00	5.912500e+02	5.912485e+02	1.500000e-03	2.229295e+06	2.229295e+06	0.000000e+00
5.400000e+02	7.164483e+02	7.164483e+02	0.000000e+00	5.923584e-02	5.923584e-02	0.000000e+00	5.971601e+02	5.971585e+02	1.600000e-03	2.246163e+06	2.246163e+06	0.000000e+00
5.500000e+02	7.500637e+02	7.500637e+02	0.000000e+00	6.255422e-02	6.255422e-02	0.000000e+00	6.027708e+02	6.027691e+02	1.700000e-03	2.263132e+06	2.263133e+06	-1.000000e+00
5.600000e+02	7.835814e+02	7.835814e+02	0.000000e+00	6.581530e-02	6.581530e-02	0.000000e+00	6.080472e+02	6.080455e+02	1.700000e-03	2.280284e+06	2.280284e+06	0.000000e+00
5.700000e+02	8.170015e+02	8.170014e+02	1.000000e-04	6.899362e-02	6.899362e-02	0.000000e+00	6.129411e+02	6.129393e+02	1.800000e-03	2.297724e+06	2.297724e+06	0.000000e+00
5.800000e+02	8.503239e+02	8.503238e+02	1.000000e-04	7.206277e-02	7.206277e-02	0.000000e+00	6.174027e+02	6.174008e+02	1.900000e-03	2.315566e+06	2.315567e+06	-1.000000e+00
5.900000e+02	8.835486e+02	8.835485e+02	1.000000e-04	7.499864e-02	7.499864e-02	0.000000e+00	6.213913e+02	6.213894e+02	1.900000e-03	2.333912e+06	2.333912e+06	0.000000e+00
6.000000e+02	9.166756e+02	9.166755e+02	1.000000e-04	7.777610e-02	7.777610e-02	0.000000e+00	6.248636e+02	6.248616e+02	2.000000e-03	2.352868e+06	2.352868e+06	0.000000e+00
6.100000e+02	9.497050e+02	9.497049e+02	1.000000e-04	8.043811e-02	8.043811e-02	0.000000e+00	6.279859e+02	6.279838e+02	2.100000e-03	2.372155e+06	2.372155e+06	0.000000e+00
6.200000e+02	9.826367e+02	9.826366e+02	1.000000e-04	8.301677e-02	8.301677e-02	0.000000e+00	6.308758e+02	6.308737e+02	2.100000e-03	2.391575e+06	2.391575e+06	0.000000e+00
6.300000e+02	1.015471e+03	1.015471e+03	0.000000e+00	8.551363e-02	8.551363e-02	0.000000e+00	6.335503e+02	6.335481e+02	2.200000e-03	2.411109e+06	2.411110e+06	-1.000000e+00
6.400000e+02	1.048207e+03	1.048207e+03	0.000000e+00	8.793025e-02	8.793025e-02	0.000000e+00	6.360240e+02	6.360217e+02	2.300000e-03	2.430747e+06	2.430747e+06	0.000000e+00
6.500000e+02	1.080846e+03	1.080846e+03	0.000000e+00	9.026816e-02	9.026816e-02	0.000000e+00	6.383099e+02	6.383075e+02	2.400000e-03	2.450477e+06	2.450477e+06	0.000000e+00
6.600000e+02	1.113387e+03	1.113387e+03	0.000000e+00	9.252891e-02	9.252891e-02	0.000000e+00	6.404195e+02	6.404171e+02	2.400000e-03	2.470289e+06	2.470289e+06	0.000000e+00
6.700000e+02	1.145830e+03	1.145830e+03	0.000000e+00	9.471406e-02	9.471406e-02	0.000000e+00	6.423631e+02	6.423606e+02	2.500000e-03	2.490177e+06	2.490177e+06	0.000000e+00
6.800000e+02	1.178176e+03	1.178176e+03	0.000000e+00	9.682516e-02	9.682516e-02	0.000000e+00	6.441500e+02	6.441475e+02	2.500000e-03	2.510134e+06	2.510134e+06	0.000000e+00
6.900000e+02	1.210424e+03	1.210424e+03	0.000000e+00	9.886374e-02	9.886374e-02	0.000000e+00	6.457886e+02	6.457861e+02	2.500000e-03	2.530154e+06	2.530154e+06	0.000000e+00
7.000000e+02	1.242574e+03	1.242574e+03	0.000000e+00	1.008314e-01	1.008314e-01	0.000000e+00	6.472867e+02	6.472841e+02	2.600000e-03	2.550234e+06	2.550234e+06	0.000000e+00
7.100000e+02	1.274627e+03	1.274627e+03	0.000000e+00	1.027296e-01	1.027296e-01	0.000000e+00	6.486513e+02	6.486487e+02	2.600000e-03	2.570368e+06	2.570368e+06	0.000000e+00
7.200000e+02	1.306582e+03	1.306582e+03	0.000000e+00	1.045599e-01	1.045599e-01	0.000000e+00	6.498890e+02	6.498863e+02	2.700000e-03	2.590554e+06	2.590554e+06	0.000000e+00
7.300000e+02	1.338439e+03	1.338439e+03	0.000000e+00	1.063240e-01	1.063240e-01	0.000000e+00	6.510059e+02	6.510031e+02	2.800000e-03	2.610789e+06	2.610789e+06	0.000000e+00
7.400000e+02	1.370199e+03	1.370199e+03	0.000000e+00	1.080232e-01	1.080232e-01	0.000000e+00	6.520076e+02	6.520048e+02	2.800000e-03	2.631070e+06	2.631070e+06	0.000000e+00
7.500000e+02	1.401861e+03	1.401861e+03	0.000000e+00	1.096593e-01	1.096593e-01	0.000000e+00	6.528996e+02	6.528968e+02	2.800000e-03	2.651395e+06	2.651395e+06	0.000000e+00
7.600000e+02	1.433426e+03	1.433425e+03	1.000000e-03	1.112337e-01	1.112337e-01	0.000000e+00	6.536868e+02	6.536840e+02	2.800000e-03	2.671762e+06	2.671762e+06	0.000000e+00
7.700000e+02	1.464892e+03	1.464892e+03	0.000000e+00	1.127479e-01	1.127479e-01	0.000000e+00	6.543741e+02	6.543713e+02	2.800000e-03	2.692168e+06	2.692168e+06	0.000000e+00
7.800000e+02	1.496261e+03	1.496261e+03	0.000000e+00	1.142036e-01	1.142036e-01	0.000000e+00	6.549661e+02	6.549632e+02	2.900000e-03	2.712612e+06	2.712612e+06	0.000000e+00
7.900000e+02	1.527533e+03	1.527533e+03	0.000000e+00	1.156023e-01	1.156023e-01	0.000000e+00	6.554670e+02	6.554641e+02	2.900000e-03	2.733092e+06	2.733092e+06	0.000000e+00
8.000000e+02	1.558706e+03	1.558706e+03	0.000000e+00	1.169455e-01	1.169455e-01	0.000000e+00	6.558812e+02	6.558782e+02	3.000000e-03	2.753607e+06	2.753607e+06	0.000000e+00
8.100000e+02	1.589782e+03	1.589782e+03	0.000000e+00	1.182348e-01	1.182348e-01	0.000000e+00	6.562124e+02	6.562094e+02	3.000000e-03	2.774154e+06	2.774154e+06	0.000000e+00
8.200000e+02	1.620761e+03	1.620761e+03	0.000000e+00	1.194717e-01	1.194717e-01	0.000000e+00	6.564648e+02	6.564618e+02	3.000000e-03	2.794732e+06	2.794732e+06	0.000000e+00
8.300000e+02	1.651642e+03	1.651641e+03	1.000000e-03	1.206577e-01	1.206577e-01	0.000000e+00	6.566419e+02	6.566389e+02	3.000000e-03	2.815339e+06	2.815339e+06	0.000000e+00
8.400000e+02	1.682425e+03	1.682424e+03	1.000000e-03	1.217945e-01	1.217945e-01	0.000000e+00	6.567475e+02	6.567444e+02	3.100000e-03	2.835973e+06	2.835973e+06	0.000000e+00
8.500000e+02	1.713110e+03	1.713110e+03	0.000000e+00	1.228836e-01	1.228836e-01	0.000000e+00	6.567850e+02	6.567819e+02	3.100000e-03	2.856633e+06	2.856633e+06	0.000000e+00
8.600000e+02	1.743698e+03	1.743697e+03	1.000000e-03	1.239265e-01	1.239265e-01	0.000000e+00	6.567579e+02	6.567548e+02	3.100000e-03	2.877317e+06	2.877317e+06	0.000000e+00
8.700000e+02	1.774188e+03	1.774187e+03	1.000000e-03	1.249247e-01	1.249247e-01	0.000000e+00	6.566696e+02	6.566664e+02	3.200000e-03	2.898022e+06	2.898022e+06	0.000000e+00
8.800000e+02	1.804580e+03	1.804580e+03	0.000000e+00	1.258798e-01	1.258798e-01	0.000000e+00	6.565232e+02	6.565201e+02	3.100000e-03	2.918747e+06	2.918747e+06	0.000000e+00
8.900000e+02	1.834875e+03	1.834874e+03	1.000000e-03	1.267934e-01	1.267934e-01	0.000000e+00	6.563221e+02	6.563189e+02	3.200000e-03	2.939489e+06	2.939489e+06	0.000000e+00
9.000000e+02	1.865072e+03	1.865071e+03	1.000000e-03	1.276671e-01	1.276671e-01	0.000000e+00	6.560693e+02	6.560661e+02	3.200000e-03	2.960246e+06	2.960246e+06	0.000000e+00
9.100000e+02	1.895171e+03	1.895171e+03	0.000000e+00	1.285022e-01	1.285022e-01	0.000000e+00	6.557679e+02	6.557647e+02	3.200000e-03	2.981015e+06	2.981015e+06	0.000000e+00
9.200000e+02	1.925173e+03	1.925173e+03	0.000000e+00	1.293005e-01	1.293005e-01	0.000000e+00	6.554209e+02	6.554177e+02	3.200000e-03	3.001794e+06	3.001794e+06	0.000000e+00
9.300000e+02	1.955077e+03	1.955077e+03	0.000000e+00	1.300635e-01	1.300635e-01	0.000000e+00	6.550314e+02	6.550282e+02	3.200000e-03	3.022581e+06	3.022581e+06	0.000000e+00
9.400000e+02	1.984883e+03	1.984883e+03	0.000000e+00	1.307926e-01	1.307926e-01	0.000000e+00	6.546022e+02	6.545989e+02	3.300000e-03	3.043371e+06	3.043371e+06	0.000000e+00
9.500000e+02	2.014592e+03	2.014592e+03	0.000000e+00	1.314895e-01	1.314895e-01	0.000000e+00	6.541362e+02	6.541329e+02	3.300000e-03	3.064163e+06	3.064163e+06	0.000000e+00
9.600000e+02	2.044203e+03	2.044203e+03	0.000000e+00	1.321557e-01	1.321557e-01	0.000000e+00	6.536362e+02	6.536330e+02	3.200000e-03	3.084952e+06	3.084952e+06	0.000000e+00
9.700000e+02	2.073716e+03	2.073716e+03	0.000000e+00	1.327928e-01	1.327928e-01	0.000000e+00	6.531051e+02	6.531018e+02	3.300000e-03	3.105735e+06	3.105735e+06	0.000000e+00
9.800000e+02	2.103132e+03	2.103132e+03	0.000000e+00	1.334022e-01	1.334022e-01	0.000000e+00	6.525456e+02	6.525423e+02	3.300000e-03	3.126509e+06	3.126509e+06	0.000000e+00
9.900000e+02	2.132450e+03	2.132450e+03	0.000000e+00	1.339856e-01	1.339856e-01	0.000000e+00	6.519604e+02	6.519571e+02	3.300000e-03	3.147269e+06	3.147269e+06	0.000000e+00
1.000000e+03	2.161670e+03	2.161670e+03	0.000000e+00	1.345445e-01	1.345445e-01	0.000000e+00	6.513522e+02	6.513489e+02	3.300000e-03	3.168012e+06	3.168012e+06	0.000000e+00

IAPWS95

Table 15 Comparison between result of [Driesner, 2007] and result calculated by xThermo.

T[deg.C]	P(Driesner)[bar]	P(xThermo)[bar]	P(diff)[bar]	X(Driesner)[mol]	X(xThermo)[mol]	X(diff)[mol]	Rho(Driesner)[kg/m3]	Rho(xThermo)	Rho(err)	H(Driesner)[J/kg]	H(xThermo)	H(err)
3.739760e+02	2.205491e+02	2.207108e+02	-1.617000e-01	0.000000e+00	2.408996e-06	-2.408996e-06	3.218900e+02	3.262793e+02	-4.389300e+00	2.086000e+06	2.119524e+06	-3.352400e+04
3.740000e+02	2.206058e+02	2.207674e+02	-1.616000e-01	1.925758e-06	4.349138e-06	-2.423380e-06	3.220024e+02	3.265468e+02	-4.544400e+00	2.090401e+06	2.123328e+06	-3.292700e+04
3.741000e+02	2.208416e+02	2.210032e+02	-1.616000e-01	1.007350e-05	1.255666e-05	-2.483160e-06	3.226283e+02	3.275549e+02	-4.926600e+00	2.111710e+06	2.129138e+06	-1.742800e+04
3.742000e+02	2.210772e+02	2.212387e+02	-1.615000e-01	1.842022e-05	2.096292e-05	-2.542700e-06	3.234210e+02	3.285209e+02	-5.099900e+00	2.127722e+06	2.131079e+06	-3.357000e+03
3.743000e+02	2.213126e+02	2.214740e+02	-1.614000e-01	2.696511e-05	2.956709e-05	-2.601980e-06	3.243026e+02	3.294588e+02	-5.156200e+00	2.133951e+06	2.131924e+06	2.027000e+03
3.744000e+02	2.215478e+02	2.217091e+02	-1.613000e-01	3.570734e-05	3.836837e-05	-2.661030e-06	3.252234e+02	3.303702e+02	-5.146800e+00	2.136659e+06	2.132277e+06	4.382000e+03
3.745000e+02	2.217828e+02	2.219441e+02	-1.613000e-01	4.464610e-05	4.736593e-05	-2.719830e-06	3.261537e+02	3.312564e+02	-5.102700e+00	2.137989e+06	2.132362e+06	5.627000e+03
3.746000e+02	2.220176e+02	2.221788e+02	-1.612000e-01	5.378059e-05	5.655898e-05	-2.778390e-06	3.270754e+02	3.321182e+02	-5.042800e+00	2.138636e+06	2.132286e+06	6.350000e+03
3.747000e+02	2.222522e+02	2.224134e+02	-1.612000e-01	6.310998e-05	6.594669e-05	-2.836710e-06	3.279778e+02	3.329567e+02	-4.978900e+00	2.138899e+06	2.132105e+06	6.794000e+03
3.748000e+02	2.224866e+02	2.226478e+02	-1.612000e-01	7.263348e-05	7.552826e-05	-2.894780e-06	3.288547e+02	3.337727e+02	-4.918000e+00	2.138922e+06	2.131853e+06	7.069000e+03
3.749000e+02	2.227209e+02	2.228820e+02	-1.611000e-01	8.235027e-05	8.530288e-05	-2.952610e-06	3.297035e+02	3.345674e+02	-4.863900e+00	2.138788e+06	2.131551e+06	7.237000e+03
3.750000e+02	2.229550e+02	2.231161e+02	-1.611000e-01	9.225956e-05	9.526976e-05	-3.010200e-06	3.305218e+02	3.353417e+02	-4.819900e+00	2.138547e+06	2.131214e+06	7.333000e+03
3.751000e+02	2.231890e+02	2.233500e+02	-1.610000e-01	1.023605e-04	1.054281e-04	-3.067600e-06	3.313108e+02	3.360968e+02	-4.786000e+00	2.138230e+06	2.130852e+06	7.378000e+03
3.756000e+02	2.243566e+02	2.245174e+02	-1.608000e-01	1.557130e-04	1.590637e-04	-3.350700e-06	3.348582e+02	3.396164e+02	-4.758200e+00	2.136103e+06	2.128859e+06	7.244000e+03
3.761000e+02	2.255209e+02	2.256815e+02	-1.606000e-01	2.137389e-04	2.173670e-04	-3.628100e-06	3.379003e+02	3.427881e+02	-4.887800e+00	2.133750e+06	2.126806e+06	6.944000e+03
3.766000e+02	2.266822e+02	2.268427e+02	-1.605000e-01	2.763407e-04	2.802403e-04	-3.899600e-06	3.406247e+02	3.456979e+02	-5.073200e+00	2.131521e+06	2.124837e+06	6.684000e+03
3.771000e+02	2.278411e+02	2.280014e+02	-1.603000e-01	3.434222e-04	3.475875e-04	-4.165300e-06	3.431574e+02	3.484092e+02	-5.251800e+00	2.129505e+06	2.122999e+06	6.506000e+03
3.776000e+02	2.289976e+02	2.291578e+02	-1.602000e-01	4.148884e-04	4.193137e-04	-4.425300e-06	3.455768e+02	3.509683e+02	-5.391500e+00	2.127703e+06	2.121303e+06	6.400000e+03
3.781000e+02	2.301523e+02	2.303123e+02	-1.600000e-01	4.906457e-04	4.953256e-04	-4.679900e-06	3.479295e+02	3.534089e+02	-5.479400e+00	2.126086e+06	2.119747e+06	6.339000e+03
3.786000e+02	2.313052e+02	2.314652e+02	-1.600000e-01	5.706021e-04	5.755310e-04	-4.928900e-06	3.502430e+02	3.557560e+02	-5.513000e+00	2.124621e+06	2.118318e+06	6.303000e+03
3.791000e+02	2.324566e+02	2.326165e+02	-1.599000e-01	6.546664e-04	6.598389e-04	-5.172500e-06	3.525330e+02	3.580283e+02	-5.495300e+00	2.123277e+06	2.117003e+06	6.274000e+03
3.796000e+02	2.336069e+02	2.337667e+02	-1.598000e-01	7.427492e-04	7.481600e-04	-5.410800e-06	3.548078e+02	3.602399e+02	-5.432100e+00	2.122029e+06	2.115787e+06	6.242000e+03
3.810000e+02	2.368226e+02	2.369823e+02	-1.597000e-01	1.010012e-03	1.016062e-03	-6.050000e-06	3.611200e+02	3.661929e+02	-5.072900e+00	2.118897e+06	2.112803e+06	6.094000e+03
3.820000e+02	2.391168e+02	2.392765e+02	-1.597000e-01	1.218734e-03	1.225217e-03	-6.483000e-06	3.655758e+02	3.702891e+02	-4.713300e+00	2.116880e+06	2.110955e+06	5.925000e+03
3.830000e+02	2.414102e+02	2.415698e+02	-1.596000e-01	1.441552e-03	1.448448e-03	-6.896000e-06	3.699761e+02	3.742926e+02	-4.316500e+00	2.114985e+06	2.109279e+06	5.706000e+03
3.840000e+02	2.437038e+02	2.438635e+02	-1.597000e-01	1.677822e-03	1.685112e-03	-7.290000e-06	3.743100e+02	3.782229e+02	-3.912900e+00	2.113184e+06	2.107736e+06	5.448000e+03
3.850000e+02	2.459988e+02	2.461585e+02	-1.597000e-01	1.926920e-03	1.934585e-03	-7.665000e-06	3.785689e+02	3.820925e+02	-3.523600e+00	2.111460e+06	2.106297e+06	5.163000e+03
3.860000e+02	2.482960e+02	2.484559e+02	-1.599000e-01	2.188239e-03	2.196261e-03	-8.022000e-06	3.827470e+02	3.859089e+02	-3.161900e+00	2.109803e+06	2.104944e+06	4.859000e+03
3.870000e+02	2.505964e+02	2.507563e+02	-1.599000e-01	2.461191e-03	2.469553e-03	-8.362000e-06	3.868403e+02	3.896763e+02	-2.836000e+00	2.108209e+06	2.103665e+06	4.544000e+03
3.880000e+02	2.529007e+02	2.530608e+02	-1.601000e-01	2.745207e-03	2.753893e-03	-8.686000e-06	3.908469e+02	3.933965e+02	-2.549600e+00	2.106674e+06	2.102450e+06	4.224000e+03
3.890000e+02	2.552097e+02	2.553699e+02	-1.602000e-01	3.039736e-03	3.048728e-03	-8.992000e-06	3.947660e+02	3.970696e+02	-2.303600e+00	2.105198e+06	2.101294e+06	3.904000e+03
3.900000e+02	2.575241e+02	2.576845e+02	-1.604000e-01	3.344241e-03	3.353525e-03	-9.284000e-06	3.985976e+02	4.006947e+02	-2.097100e+00	2.103781e+06	2.100192e+06	3.589000e+03
3.910000e+02	2.598445e+02	2.600050e+02	-1.605000e-01	3.658206e-03	3.667766e-03	-9.560000e-06	4.023425e+02	4.042696e+02	-1.927100e+00	2.102421e+06	2.099141e+06	3.280000e+03
3.920000e+02	2.621714e+02	2.623322e+02	-1.608000e-01	3.981129e-03	3.990950e-03	-9.821000e-06	4.060019e+02	4.077921e+02	-1.790200e+00	2.101120e+06	2.098138e+06	2.982000e+03
3.930000e+02	2.645056e+02	2.646666e+02	-1.610000e-01	4.312524e-03	4.322592e-03	-1.006800e-05	4.095772e+02	4.112593e+02	-1.682100e+00	2.099877e+06	2.097182e+06	2.695000e+03
3.940000e+02	2.668474e+02	2.670087e+02	-1.613000e-01	4.651924e-03	4.662225e-03	-1.030100e-05	4.130703e+02	4.146686e+02	-1.598300e+00	2.098693e+06	2.096270e+06	2.423000e+03
3.950000e+02	2.691974e+02	2.693589e+02	-1.615000e-01	4.998874e-03	5.009395e-03	-1.052100e-05	4.164830e+02	4.180173e+02	-1.534300e+00	2.097568e+06	2.095401e+06	2.167000e+03
3.960000e+02	2.715561e+02	2.717178e+02	-1.617000e-01	5.352936e-03	5.363664e-03	-1.072800e-05	4.198173e+02	4.213032e+02	-1.485900e+00	2.096501e+06	2.094574e+06	1.927000e+03
3.970000e+02	2.739237e+02	2.740858e+02	-1.621000e-01	5.713687e-03	5.724608e-03	-1.092100e-05	4.230752e+02	4.245244e+02	-1.449200e+00	2.095492e+06	2.093788e+06	1.704000e+03
3.980000e+02	2.763009e+02	2.764632e+02	-1.623000e-01	6.080717e-03	6.091821e-03	-1.110400e-05	4.262587e+02	4.276795e+02	-1.420800e+00	2.094542e+06	2.093042e+06	1.500000e+03
3.990000e+02	2.786878e+02	2.788504e+02	-1.626000e-01	6.453632e-03	6.464906e-03	-1.127400e-05	4.293698e+02	4.307678e+02	-1.398000e+00	2.093650e+06	2.092335e+06	1.315000e+03
4.000000e+02	2.810849e+02	2.812479e+02	-1.630000e-01	6.832051e-03	6.843485e-03	-1.143400e-05	4.324106e+02	4.337888e+02	-1.378200e+00	2.092815e+06	2.091667e+06	1.148000e+03
4.100000e+02	3.056780e+02	3.058445e+02	-1.665000e-01	1.084455e-02	1.085705e-02	-1.250000e-05	4.593637e+02	4.604769e+02	-1.113200e+00	2.087488e+06	2.087161e+06	3.270000e+02
4.200000e+02	3.315379e+02	3.317083e+02	-1.704000e-01	1.508056e-02	1.509339e-02	-1.283000e-05	4.811573e+02	4.818504e+02	-6.931000e-01	2.087125e+06	2.086847e+06	2.780000e+02
4.300000e+02	3.587626e+02	3.589372e+02	-1.746000e-01	1.934403e-02	1.935673e-02	-1.270000e-05	4.990526e+02	4.993654e+02	-3.128000e-01	2.090916e+06	2.090644e+06	2.720000e+02
4.400000e+02	3.873578e+02	3.875365e+02	-1.787000e-01	2.351927e-02	2.353160e-02	-1.233000e-05	5.139825e+02	5.140217e+02	-3.920000e-02	2.098093e+06	2.097911e+06	1.820000e+02
4.500000e+02	4.172612e+02	4.174436e+02	-1.824000e-01	2.754975e-02	2.756160e-02	-1.185000e-05	5.266585e+02	5.265205e+02	1.380000e-01	2.107963e+06	2.107889e+06	7.400000e+01
4.600000e+02	4.483573e+02	4.485431e+02	-1.858000e-01	3.142030e-02	3.143168e-02	-1.138000e-05	5.376330e+02	5.373904e+02	2.426000e-01	2.119912e+06	2.119916e+06	-4.000000e+00
4.700000e+02	4.804883e+02	4.806769e+02	-1.886000e-01	3.514237e-02	3.515334e-02	-1.097000e-05	5.473332e+02	5.470389e+02	2.943000e-01	2.133417e+06	2.133454e+06	-3.700000e+01
4.800000e+02	5.134606e+02	5.136514e+02	-1.908000e-01	3.874218e-02	3.875282e-02	-1.064000e-05	5.560821e+02	5.557736e+02	3.085000e-01	2.148044e+06	2.148078e+06	-3.400000e+01
4.900000e+02	5.470509e+02	5.472432e+02	-1.923000e-01	4.225142e-02	4.226185e-02	-1.043000e-05	5.641113e+02	5.638149e+02	2.964000e-01	2.163459e+06	2.163464e+06	-5.000000e+00
5.000000e+02	5.810101e+02	5.812032e+02	-1.931000e-01	4.570045e-02	4.571074e-02	-1.029000e-05	5.715720e+02	5.713054e+02	2.666000e-01	2.179422e+06	2.179387e+06	3.500000e+01
5.100000e+02	6.150162e+02	6.150161e+02	1.000000e-04	4.911344e-02	4.912364e-02	-1.020000e-05	5.785204e+02	5.781991e+02	3.213000e-01	2.195824e+06	2.195884e+06	-6.000000e+01
5.200000e+02	6.489246e+02	6.489245e+02	1.000000e-04	5.250551e-02	5.251566e-02	-1.015000e-05	5.850483e+02	5.847792e+02	2.691000e-01	2.212499e+06	2.212506e+06	-7.000000e+00
5.300000e+02	6.827353e+02	6.827352e+02	1.000000e-04	5.588146e-02	5.589156e-02	-1.010000e-05	5.912500e+02	5.910346e+02	2.154000e-01	2.229295e+06	2.229258e+06	3.700000e+01
5.400000e+02	7.164483e+02	7.164483e+02	0.000000e+00	5.923584e-02	5.924586e-02	-1.002000e-05	5.971601e+02	5.969976e+02	1.625000e-01	2.246163e+06	2.246094e+06	6.900000e+01
5.500000e+02	7.500637e+02	7.500637e+02	0.000000e+00	6.255422e-02	6.256411e-02	-9.890000e-06	6.027708e+02	6.026587e+02	1.121000e-01	2.263132e+06	2.263042e+06	9.000000e+01
5.600000e+02	7.835814e+02	7.835814e+02	0.000000e+00	6.581530e-02	6.582497e-02	-9.670000e-06	6.080472e+02	6.079818e+02	6.540000e-02	2.280284e+06	2.280184e+06	1.000000e+02
5.700000e+02	8.170015e+02	8.170014e+02	1.000000e-04	6.899362e-02	6.900300e-02	-9.380000e-06	6.129411e+02	6.129176e+02	2.350000e-02	2.297724e+06	2.297625e+06	9.900000e+01
5.800000e+02	8.503239e+02	8.503238e+02	1.000000e-04	7.206277e-02	7.207178e-02	-9.010000e-06	6.174027e+02	6.174157e+02	-1.300000e-02	2.315566e+06	2.315480e+06	8.600000e+01
5.900000e+02	8.835486e+02	8.835485e+02	1.000000e-04	7.499864e-02	7.500723e-02	-8.590000e-06	6.213913e+02	6.214352e+02	-4.390000e-02	2.333912e+06	2.333847e+06	6.500000e+01
6.000000e+02	9.166756e+02	9.166755e+02	1.000000e-04	7.777610e-02	7.777610e-02	0.000000e+00	6.248636e+02	6.249066e+02	-4.300000e-02	2.352868e+06	2.352882e+06	-1.400000e+01
6.100000e+02	9.497050e+02	9.497049e+02	1.000000e-04	8.043811e-02	8.043811e-02	0.000000e+00	6.279859e+02	6.280498e+02	-6.390000e-02	2.372155e+06	2.372207e+06	-5.200000e+01
6.200000e+02	9.826367e+02	9.826366e+02	1.000000e-04	8.301677e-02	8.301677e-02	0.000000e+00	6.308758e+02	6.309546e+02	-7.880000e-02	2.391575e+06	2.391672e+06	-9.700000e+01
6.300000e+02	1.015471e+03	1.015471e+03	0.000000e+00	8.551363e-02	8.551363e-02	0.000000e+00	6.335503e+02	6.336379e+02	-8.760000e-02	2.411109e+06	2.411261e+06	-1.520000e+02
6.400000e+02	1.048207e+03	1.048207e+03	0.000000e+00	8.793025e-02	8.793025e-02	0.000000e+00	6.360240e+02	6.361144e+02	-9.040000e-02	2.430747e+06	2.430961e+06	-2.140000e+02
6.500000e+02	1.080846e+03	1.080846e+03	0.000000e+00	9.026816e-02	9.026816e-02	0.000000e+00	6.383099e+02	6.383970e+02	-8.710000e-02	2.450477e+06	2.450759e+06	-2.820000e+02
6.600000e+02	1.113387e+03	1.113387e+03	0.000000e+00	9.252891e-02	9.252891e-02	0.000000e+00	6.404195e+02	6.404973e+02	-7.780000e-02	2.470289e+06	2.470646e+06	-3.570000e+02
6.700000e+02	1.145830e+03	1.145830e+03	0.000000e+00	9.471406e-02	9.471406e-02	0.000000e+00	6.423631e+02	6.424256e+02	-6.250000e-02	2.490177e+06	2.490614e+06	-4.370000e+02
6.800000e+02	1.178176e+03	1.178176e+03	0.000000e+00	9.682516e-02	9.682516e-02	0.000000e+00	6.441500e+02	6.441914e+02	-4.140000e-02	2.510134e+06	2.510656e+06	-5.220000e+02
6.900000e+02	1.210424e+03	1.210424e+03	0.000000e+00	9.886374e-02	9.886374e-02	0.000000e+00	6.457886e+02	6.458031e+02	-1.450000e-02	2.530154e+06	2.530767e+06	-6.130000e+02
7.000000e+02	1.242574e+03	1.242574e+03	0.000000e+00	1.008314e-01	1.008314e-01	0.000000e+00	6.472867e+02	6.472687e+02	1.800000e-02	2.550234e+06	2.550941e+06	-7.070000e+02
7.100000e+02	1.274627e+03	1.274627e+03	0.000000e+00	1.027296e-01	1.027296e-01	0.000000e+00	6.486513e+02	6.485953e+02	5.600000e-02	2.570368e+06	2.571174e+06	-8.060000e+02
7.200000e+02	1.306582e+03	1.306582e+03	0.000000e+00	1.045599e-01	1.045599e-01	0.000000e+00	6.498890e+02	6.497897e+02	9.930000e-02	2.590554e+06	2.591463e+06	-9.090000e+02
7.300000e+02	1.338439e+03	1.338439e+03	0.000000e+00	1.063240e-01	1.063240e-01	0.000000e+00	6.510059e+02	6.508581e+02	1.478000e-01	2.610789e+06	2.611803e+06	-1.014000e+03
7.400000e+02	1.370199e+03	1.370199e+03	0.000000e+00	1.080232e-01	1.080232e-01	0.000000e+00	6.520076e+02	6.518063e+02	2.013000e-01	2.631070e+06	2.632193e+06	-1.123000e+03
7.500000e+02	1.401861e+03	1.401861e+03	0.000000e+00	1.096593e-01	1.096593e-01	0.000000e+00	6.528996e+02	6.526400e+02	2.596000e-01	2.651395e+06	2.652630e+06	-1.235000e+03
7.600000e+02	1.433426e+03	1.433425e+03	1.000000e-03	1.112337e-01	1.112337e-01	0.000000e+00	6.536868e+02	6.533643e+02	3.225000e-01	2.671762e+06	2.673110e+06	-1.348000e+03
7.700000e+02	1.464892e+03	1.464892e+03	0.000000e+00	1.127479e-01	1.127479e-01	0.000000e+00	6.543741e+02	6.539842e+02	3.899000e-01	2.692168e+06	2.693633e+06	-1.465000e+03
7.800000e+02	1.496261e+03	1.496261e+03	0.000000e+00	1.142036e-01	1.142036e-01	0.000000e+00	6.549661e+02	6.545045e+02	4.616000e-01	2.712612e+06	2.714195e+06	-1.583000e+03
7.900000e+02	1.527533e+03	1.527533e+03	0.000000e+00	1.156023e-01	1.156023e-01	0.000000e+00	6.554670e+02	6.549298e+02	5.372000e-01	2.733092e+06	2.734796e+06	-1.704000e+03
8.000000e+02	1.558706e+03	1.558706e+03	0.000000e+00	1.169455e-01	1.169455e-01	0.000000e+00	6.558812e+02	6.552643e+02	6.169000e-01	2.753607e+06	2.755432e+06	-1.825000e+03
8.100000e+02	1.589782e+03	1.589782e+03	0.000000e+00	1.182348e-01	1.182348e-01	0.000000e+00	6.562124e+02	6.555123e+02	7.001000e-01	2.774154e+06	2.776103e+06	-1.949000e+03
8.200000e+02	1.620761e+03	1.620761e+03	0.000000e+00	1.194717e-01	1.194717e-01	0.000000e+00	6.564648e+02	6.556779e+02	7.869000e-01	2.794732e+06	2.796806e+06	-2.074000e+03
8.300000e+02	1.651642e+03	1.651641e+03	1.000000e-03	1.206577e-01	1.206577e-01	0.000000e+00	6.566419e+02	6.557651e+02	8.768000e-01	2.815339e+06	2.817540e+06	-2.201000e+03
8.400000e+02	1.682425e+03	1.682424e+03	1.000000e-03	1.217945e-01	1.217945e-01	0.000000e+00	6.567475e+02	6.557775e+02	9.700000e-01	2.835973e+06	2.838302e+06	-2.329000e+03
8.500000e+02	1.713110e+03	1.713110e+03	0.000000e+00	1.228836e-01	1.228836e-01	0.000000e+00	6.567850e+02	6.557191e+02	1.065900e+00	2.856633e+06	2.859092e+06	-2.459000e+03
8.600000e+02	1.743698e+03	1.743697e+03	1.000000e-03	1.239265e-01	1.239265e-01	0.000000e+00	6.567579e+02	6.555933e+02	1.164600e+00	2.877317e+06	2.879906e+06	-2.589000e+03
8.700000e+02	1.774188e+03	1.774187e+03	1.000000e-03	1.249247e-01	1.249247e-01	0.000000e+00	6.566696e+02	6.554037e+02	1.265900e+00	2.898022e+06	2.900742e+06	-2.720000e+03
8.800000e+02	1.804580e+03	1.804580e+03	0.000000e+00	1.258798e-01	1.258798e-01	0.000000e+00	6.565232e+02	6.551537e+02	1.369500e+00	2.918747e+06	2.921599e+06	-2.852000e+03
8.900000e+02	1.834875e+03	1.834874e+03	1.000000e-03	1.267934e-01	1.267934e-01	0.000000e+00	6.563221e+02	6.548468e+02	1.475300e+00	2.939489e+06	2.942475e+06	-2.986000e+03
9.000000e+02	1.865072e+03	1.865071e+03	1.000000e-03	1.276671e-01	1.276671e-01	0.000000e+00	6.560693e+02	6.544863e+02	1.583000e+00	2.960246e+06	2.963366e+06	-3.120000e+03
9.100000e+02	1.895171e+03	1.895171e+03	0.000000e+00	1.285022e-01	1.285022e-01	0.000000e+00	6.557679e+02	6.540753e+02	1.692600e+00	2.981015e+06	2.984270e+06	-3.255000e+03
9.200000e+02	1.925173e+03	1.925173e+03	0.000000e+00	1.293005e-01	1.293005e-01	0.000000e+00	6.554209e+02	6.536172e+02	1.803700e+00	3.001794e+06	3.005185e+06	-3.391000e+03
9.300000e+02	1.955077e+03	1.955077e+03	0.000000e+00	1.300635e-01	1.300635e-01	0.000000e+00	6.550314e+02	6.531149e+02	1.916500e+00	3.022581e+06	3.026108e+06	-3.527000e+03
9.400000e+02	1.984883e+03	1.984883e+03	0.000000e+00	1.307926e-01	1.307926e-01	0.000000e+00	6.546022e+02	6.525716e+02	2.030600e+00	3.043371e+06	3.047035e+06	-3.664000e+03
9.500000e+02	2.014592e+03	2.014592e+03	0.000000e+00	1.314895e-01	1.314895e-01	0.000000e+00	6.541362e+02	6.519904e+02	2.145800e+00	3.064163e+06	3.067964e+06	-3.801000e+03
9.600000e+02	2.044203e+03	2.044203e+03	0.000000e+00	1.321557e-01	1.321557e-01	0.000000e+00	6.536362e+02	6.513742e+02	2.262000e+00	3.084952e+06	3.088891e+06	-3.939000e+03
9.700000e+02	2.073716e+03	2.073716e+03	0.000000e+00	1.327928e-01	1.327928e-01	0.000000e+00	6.531051e+02	6.507259e+02	2.379200e+00	3.105735e+06	3.109813e+06	-4.078000e+03
9.800000e+02	2.103132e+03	2.103132e+03	0.000000e+00	1.334022e-01	1.334022e-01	0.000000e+00	6.525456e+02	6.500485e+02	2.497100e+00	3.126509e+06	3.130726e+06	-4.217000e+03
9.900000e+02	2.132450e+03	2.132450e+03	0.000000e+00	1.339856e-01	1.339856e-01	0.000000e+00	6.519604e+02	6.493447e+02	2.615700e+00	3.147269e+06	3.151626e+06	-4.357000e+03
1.000000e+03	2.161670e+03	2.161670e+03	0.000000e+00	1.345445e-01	1.345445e-01	0.000000e+00	6.513522e+02	6.486175e+02	2.734700e+00	3.168012e+06	3.172509e+06	-4.497000e+03

Tip

The help function for 3D plot can be downloaded at here: helpfunc.py