We model molecular outflows produced by the time dependent interaction between a stellar wind and a rotating cloud envelope in gravitational collapse, studied by Ulrich. We consider spherical and anisotropic stellar winds. We assume that the bipolar outflow is a thin shocked shell, with axial symmetry around the cloud rotation axis and obtain the mass and momentum fluxes into the shell. We solve numerically a set of partial differential equations in space and time, and obtain the shape of the shell, the mass surface density, the velocity field, and the angular momentum of the material in the shell. We find that there is a critical value of the ratio between the wind and the accretion flow momentum rates $\beta$ that allows the shell to expand. As expected, the elongation of the shells increase with the stellar wind anisotropy. In our models, the rotation velocity of the shell is the order to 0.1 - 0.2 km s$^{-1}$, a factor of 5-10 lower than the values measured in several sources. We compare our models with those of Wilkin and Stahler for early evolutionary times and find that our shells have the same sizes at the pole, although we use different boundary conditions at the equator.