xref: /phasta/phSolver/common/gettab.f (revision 595995161822a203c8467e0e4a253d7bd7d6df32)
1*59599516SKenneth E. Jansen        subroutine gettab  (mut,  rhot,  xst)
2*59599516SKenneth E. Jansenc
3*59599516SKenneth E. Jansenc-----------------------------------------------------------------------
4*59599516SKenneth E. Jansenc
5*59599516SKenneth E. Jansenc  This subroutine reads the three tables for equilibrium chemistry.
6*59599516SKenneth E. Jansenc
7*59599516SKenneth E. Jansenc
8*59599516SKenneth E. Jansenc output:
9*59599516SKenneth E. Jansenc
10*59599516SKenneth E. Jansenc    mut  (71,451)      : specific chemical potential function of (p,T)
11*59599516SKenneth E. Jansenc    rhot (71,451)      : density function of (p,T)
12*59599516SKenneth E. Jansenc    xst  (5,71,451)    : mole fractions functions of (p,T)
13*59599516SKenneth E. Jansenc
14*59599516SKenneth E. Jansenc Note: These three arrays are always in double precision.
15*59599516SKenneth E. Jansenc
16*59599516SKenneth E. Jansenc Frederic Chalot and Zdenek Johan, Fall 1990.
17*59599516SKenneth E. Jansenc-----------------------------------------------------------------------
18*59599516SKenneth E. Jansenc
19*59599516SKenneth E. Jansen        include "common.h"
20*59599516SKenneth E. Jansenc
21*59599516SKenneth E. Jansen        real*8  mut(71,451),  rhot(71,451),  xst(5,71,451)
22*59599516SKenneth E. Jansenc
23*59599516SKenneth E. Jansenc.... open table file
24*59599516SKenneth E. Jansenc
25*59599516SKenneth E. Jansen        open (unit=itable, file=ftable, form='unformatted',
26*59599516SKenneth E. Jansen     &                                  status='unknown')
27*59599516SKenneth E. Jansenc
28*59599516SKenneth E. Jansenc.... read tables
29*59599516SKenneth E. Jansenc
30*59599516SKenneth E. Jansen        read (itable) mut
31*59599516SKenneth E. Jansenc
32*59599516SKenneth E. Jansen        read (itable) rhot
33*59599516SKenneth E. Jansenc
34*59599516SKenneth E. Jansen        read (itable) xst
35*59599516SKenneth E. Jansenc
36*59599516SKenneth E. Jansenc.... close table file
37*59599516SKenneth E. Jansenc
38*59599516SKenneth E. Jansen        close(unit=itable)
39*59599516SKenneth E. Jansenc
40*59599516SKenneth E. Jansenc.... end
41*59599516SKenneth E. Jansenc
42*59599516SKenneth E. Jansen        return
43*59599516SKenneth E. Jansen        end
44