  We simulate the creation and evolution of non-boost-invariant Glasma in the
early stages of heavy ion collisions within the color glass condensate
framework. This is accomplished by extending the McLerran-Venugopalan model to
include a parameter for the Lorentz-contracted but finite width of the nucleus
in the beam direction. We determine the rapidity profile of the Glasma energy
density, which shows deviations from the boost-invariant result. Varying the
parameters both broad and narrow profiles can be produced. We compare our
results to experimental data from RHIC and find surprising agreement.
