Improvement of numerical solution smoothness for the hydrodynamics problems modeling on rectangular grids