Stable numerical schemes for nonlinear dispersive equations with counter-propagation and gain dynamics