  Using the delta-shell representation we present a successful fit to
neutron-proton and proton-proton scattering data below pion production
threshold. A detailed overview of the theory necessary to calculate observables
with this potential is presented. A new data selection process is used to
obtain the largest mutually consistent data base. The analysis includes data
within the years 1950 to 2013. Using 46 parameters we obtain chi^2/Ndata = 1.04
with Ndata = 6713 including normalization data. Phase shifts with error bars
are provided.
