We study numerically the semiclassical limit for the nonlinearSchrödinger equation thanks to a modification of the Madelungtransform due to Grenier. This approach allows for the presence ofvacuum. Even if the mesh size and the time step do not depend on the Planck constant, we recover the position and current densities in thesemiclassical limit, with a numerical rate of convergence inaccordance with the theoreticalresults, before shocks appear in the limiting Eulerequation. By using simple projections, the mass and the momentum ofthe solution are well preserved by the numerical scheme,while the variation of the energy is not negligiblenumerically. Experiments suggest that beyond the critical time for theEuler equation, Grenier's approach yields smooth but highlyoscillatory terms.
