If applyTimes is set, and includes time t=0 in the range, particles may be loaded both during initialization and in the zeroth time step of the model, potentially giving twice the density you want.
Use together with a history of kind=speciesNumPhysical and compare the total population to the number of particles you expected.
Use together with the computePtclNumDensity script to check the overall density of the loaded particles.