C-----------------------------------------------------------------
C    NS-BWR Equation
C    This program calculates temperature, volume,entropy and quality
C    of R22-R114 mixture 
C    at P=0.6[MPa], Z=0.5[kg/kg] and H=300[kJ/kg].
C-----------------------------------------------------------------
C
C---------------------------------------------
C    Definition of array PR1,PR2,CP1 and CP2.
C---------------------------------------------
      DIMENSION PR1(5),PR2(5),CP1(9),CP2(9)
C
C----------------------------------------
C    T=[K],P=[Pa],quantity=[kg]
C    standard state=(saturated liquid at 273.15[K])
C-----------------------------------------
      call KPAMES(0,1)
      call STNKAS(1,1)
C
C----------------------------------------------------------------------------
C    Putting fundamental constants into PR1 and PR2 using DATA STATEMENT.
C    DATA(PR1(I),I=1,5)/MW1[KG/KMOL],TC1[K],dummy, VC1[M^3/KMOL],OMEGA1[-]/
C    DATA(PR2(I),I=1,5)/MW2[KG/KMOL],TC2[K],dummy, VC2[M^3/KMOL],OMEGA2[-]/
C----------------------------------------------------------------------------
      DATA(PR1(I),I=1,5)/ 86.469, 369.3,  -1.0, 0.168555, 0.205/
      DATA(PR2(I),I=1,5)/170.922, 418.78, -1.0, 0.296739, 0.234/
C
C----------------------------------------------------------------------------
C    Putting isobaric heat capacity at ideal gas state into CP1 and CP2
C    using DATA STATEMENT.
C    DATA(CP1(I),I=1,8)/EQUATION NUMBER,TMIN[K],TMAX[K], A,B,C,D,E/
C    DATA(CP2(I),I=1,8)/EQUATION NUMBER,TMIN[K],TMAX[K], A,B,C,D,E/
C----------------------------------------------------------------------------
      DATA(CP1(I),I=1,9)
     &  /3.0, 250.0, 350.0,
     &   -917.3,28.755,-2.030E-02,2.5351E-03,-2.185E-06, 7.1E-10/
C        [kJ/kg]
      DATA(CP2(I),I=1,9)
     &  /2.0, 250.0, 400.0,
     &   -0.12916, 4.7009E-03, -8.612E-06, 6.05E-09, 0.0, 0.0 /
C        [kJ/kg]
      do 10 I=4,9
C        [kJ/kg] -> [J/kmol]
         CP1(I) = CP1(I) * PR1(1) * 1000
         CP2(I) = CP2(I) * PR2(1) * 1000
 10   continue
C
C------------------------------------------
C    Setting interaction parametar AKIJ.
C------------------------------------------
      AKIJ=0.976
C
C--------------------------
C    Execution START2.
C--------------------------
      call START2(ECODE, PR1, PR2, CP1, CP2, AKIJ)
C
C--------------------------------------------
C   Calculation of temperature, volume, entropy and quality
C   from pressure, whole composition and enthalpy.
C--------------------------------------------
      P = 0.6E+6
      Z = 0.5
      H = 300.0*1000.0
      CALL SUBMIX(2, J, T, P, Z, V, H, S)
      Q = QMIX(T, P, Z)
      WRITE(*,*) ' P[MPa]     =', P/1e+6
      WRITE(*,*) ' Z[kg/kg]   =', Z
      WRITE(*,*) ' H[kJ/kg]   =', H/1000
      WRITE(*,*) ' T[K]       =', T
      WRITE(*,*) ' Q[kg/kg]   =', Q
      WRITE(*,*) ' V[m3/kg]   =', V
      WRITE(*,*) ' S[kJ/kg]   =', S/1000
      STOP
      END
C     
Cกกกก result
C     
C  P[MPa]     =  0.600000024
C  Z[kg/kg]   =  0.5
C  H[kJ/kg]   =  300.
C  T[K]       =  297.011963
C  Q[kg/kg]   =  0.667366326
C  V[m3/kg]   =  0.00257301494
C  S[kJ/kg]   =  1.34425902