We propose Gaussian-beam based Eulerian methods to compute semi-classical solutions of the Schrdinger equation. Traditional Gaussian beam type methods for the Schrdinger equation are based on the Lagrangian ray tracing. Based on the first Eulerian Gaussian beam framework proposed in Leung et al.[S. Leung, J. Qian, R. Burridge, Eulerian Gaussian beams for high frequency wave propagation, Geophysics 72 (2007) SM61SM76], we develop a new Eulerian Gaussian beam method which uses global Cartesian coordinates, level-set based implicit representation and Liouville equations. The resulting method gives uniformly distributed phases and amplitudes in phase space simultaneously. To obtain semi-classical solutions to the Schrdinger equation with different initial wave functions, we only need to slightly modify the summation formula. This yields a very efficient method for computing semi