pro chtab,nz1
    common cdat,x,y,z,nx,ny,nz,mp,ntmax,date0,time0
    common ctable,ttmean,rhm,drhm,eemin,eemax,dee,mtab,itab,tab

    nz=nz1
    rhm=rhm(0:n)
    eemin=eemin(0:n)
    eemax=eemax(0:n)
    dee=dee(0:n)
    mtab=mtab(0:n)
    itab=itab(0:n)
end

PRO eos_table,file,sstab=sstab,quiet=quiet, _extra=extra
;
;  Load Table
;
    common cdat,x,y,z,nx,ny,nz,mp,ntmax,date0,time0
    common ctable,ttmean,rhm,drhm,eemin,eemax,dee,mtab,itab,tab
    common mycom,cray,one,two

    if n_params() eq 0 then begin
    print,'eos_table,file'
    return
    endif

    close,2
    openr,2,file,_extra=extra,/f77,/swap_if_big_endian
    mrho=0L
    iupdte=0L
    nvar=0L
    mbox=0L
    xmin=0.
    dbox=0.
    ul=0.
    ur=0.
    ut=0.
    eps=0.
    tff=0.
    grv=0.
    abnd=0.
    readu,2,mrho,iupdte,nvar,mbox,xmin,dbox,ul,ut,ur,eps,tff,grv,abnd
    print,'mrho,mbox=',mrho,mbox,nvar
    a=fltarr(mrho)
    tmean=a
    tamp=a
    rhm=a
    xcorr=a
    thmin=a
    thmax=a
    dth=a
    eemin=a
    eemax=a
    dee=a
    itab=lonarr(mrho)
    mtab=itab
    readu,2,tmean,tamp,rhm,xcorr,thmin,thmax,dth,eemin,eemax,dee,itab,mtab
    ntab=0L
    readu,2,ntab
    tab=fltarr(ntab)
    readu,2,tab
    close,2

    rhm1=alog(rhm(0))
    rhm2=alog(rhm(mrho-1))
    drhm=(rhm2-rhm1)/(mrho-1)

    if n_elements(sstab) ne 0 then sstab
end