PRO loadabu,file, $
            ION=ion, $
            X=abu, $
            IDX=idx, $
            NION=n, $
            Y=y, $
            Z=z, $
            A=a

@compile_opt

OPENR,unit,file,/GET_LUN

ion=REPLICATE('',1000)
abu=REPLICATE(1.D0,1000)
idx=REPLICATE(0L,1000)

n=0
WHILE NOT EOF(unit) DO BEGIN
    s=''
    readf,unit,s
    IF STRMID(s,0,1) NE ';' THEN BEGIN
        xion=''
        xabu=0.D0
        i=strpos(s,' ')
        READS,STRMID(s,0,i),xion
        READS,STRMID(s,i+1,99),xabu
        ion[n]=xion
        abu[n]=xabu
        idx[n]=ionidx(xion)
        n+=1
    ENDIF

ENDWHILE

CLOSE,unit
FREE_LUN,unit

abu=abu[0:n-1]
ion=ion[0:n-1]
idx=idx[0:n-1]

ionaz,ion,a,z
y=abu/a

END