Gutzwiller variational wave function for a two-orbitalHubbard model on a square lattice