The paper describes a new method of gravity modelling of the density cross section along a DSS profile. We propose a systematic procedure of gravity interpretation using structural information contained both in the velocity distribution of the seismic model and the geometry of its layers. The procedure is based on the optimisation of the density parameters of the model. The values of the parameters are limited by relationships between the seismic velocity and the density (optimisation with inequalities). The formulation of the method allows for discussion of the presence of isostasy and the distribution of compensating masses. The discussion is very important for the considered crust-lithosphere system. The presented calculations are done for one Carpathian profile (CEL01). The analysis of the properties of the fitting procedure provides the characteristic depth of subcrustal compensating density structures, the position of crucial tectonic boundaries along the profiles, locations of the crustal areas having exceptional lithology, the regional density differences of subcrustal matter and the estimation of regional deviations from the isostatic equilibrium.