A non-orthogonal curvilinear systems model of species transport in shallow flow domains