By introducing an elliptic vortex ansatz, the 2+1-dimensional two-layer fluid system is reduced to a finite-dimensional nonlinear dynamical system. Time-modulated variables are then introduced and multicomponent Ermakov systems are isolated. The latter is shown to be also Hamiltonian, thereby admitting general solutions in terms of an elliptic integral representation. In particular, a subclass of vortex solutions is obtained and their behaviors are simulated. Such solutions have recently found applications in oceanic and atmospheric dynamics. Moreover, it is proved that the Hamiltonian system is equivalent to the stationary nonlinear cubic Schrödinger equations coupled with a Steen-Ermakov-Pinney equation.