Аннотация:The article is devoted to the description of a new method of numerical solutionfor hydrostatic approximation of incompressible hydrodynamic problems with free surfaces andvariable density. This technique develops a model of multilayer shallow water, taking intoaccount the mass and momentum exchange between the layers.