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