FUNCTION logpf2y, data0, SOLABUFILE=solabufile

IF N_ELEMENTS(solabufile) EQ 0 THEN solabufile=''

@solabu_var
solabu, FILE=solabufile

;make local copy
IF N_ELEMENTS(data0) NE 0 THEN data=data0

data.abu=10.D0^(data.abu)*sol_el_mol[data.Z-1]
data.sigma=(10.D0^(ABS(data.sigma))-1.D0)*data.abu
data.sigma*=SGN(data0.sigma)

RETURN, data
END