;
pro log10f, x,spval
;
for i=0,n_elements(x)-1 do begin
  if x[i] ne spval then $
    if (x[i] le 1.e-20) then x[i] = spval else x[i] = alog10(x[i])
endfor
end