Initial energy density in heavy ion collisions from a color neutral three-dimensional color glass condensate model of QCD