We describe a procedure to tune the charged hadron multiplicities for deep inelastic events produced by the NEUT neutrino interaction generator. This tuning uses a model based on Koba-Nielsen-Olesen scaling, whose parameters are obtained by fitting multiplicity data from deuterium bubble chamber experiments. After tuning, the multiplicities of the events generated by NEUT are found to be in good agreement with the measurements from the bubble chamber experiments.