Plotting a solution may give you a hint:
With[{p = 1},
{xSol, oSol} =
NDSolveValue[{x'[t] == p1[x, o, t], o'[t] == n1[x, o, t],
x[0] == initial[[p]], o[0] == Pi},
{x, o}, {t, 0, time}]]
{Plot[xSol[t], {t, 0, 2.5*^6}, PlotRange -> {-1, 1}],
Plot[oSol[t], {t, 0, 2.5*^6}, PlotRange -> {-1, 6}]}
