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