In this paper, we develop a multi layer method to solve a generalized case of a phytoplankton model introduced in [7]. It is treated by means of a sequence of approximations: the mixed layer is subdivided into a finite number of thin layers within each of which horizontal velocity can be considered constant with respect to depth. Existence, uniqueness and non negativity of solutions are investigated.