added restart_ic = true for NEUTRINOS in cubepm.par