  In the present work we propose a new initial state model for hydrodynamic
simulation of relativistic heavy ion collisions based on Bjorken-like solution
applied streak by streak in the transverse plane. Previous fluid dynamical
calculations in Cartesian coordinates with an initial state based on a streak
by streak Yang-Mills field led for peripheral higher energy collisions to large
angular momentum, initial shear flow and significant local vorticity. Recent
experiments verified the existence of this vorticity via the resulting
polarization of emitted $\Lambda$ and $\bar{\Lambda}$ particles. At the same
time parton cascade models indicated the existence of more compact initial
state configurations, which we are going to simulate in our approach.
  The proposed model satisfies all the conservation laws including conservation
of a strong initial angular momentum which is present in non-central
collisions. As a consequence of this large initial angular momentum we observe
the rotation of the whole system as well as the fluid shear in the initial
state, which leads to large flow vorticity. Another advantage of the proposed
model is that the initial state can be given in both [t,x,y,z] and $[\tau, x,
y, \eta]$ coordinates, and thus can be tested by all 3+1D hydrodynamical codes
which exist in the field.
