Model potentials in numerical local density functional methods