An algorithm is introduced to generate an initial condition for lattice Boltzmann simulations for which an initial velocity field is prescribed. The basic idea is to run a few initial lattice Boltzmann iterations, but to use the prescribed velocity instead of the moment of particle populations when computing the equilibrium distribution.