Abstract. The nonlinear Schrödinger (NLS) equation describing the propagation of weakly rotational wave packets in an infinitely deep fluid in Lagrangian coordinates has been derived. The vorticity is assumed to be an arbitrary function of Lagrangian coordinates and quadratic in the small parameter proportional to the wave steepness. The vorticity effects manifest themselves in a shift of the wave number in the carrier wave and in variation in the coefficient multiplying the nonlinear term. In the case of vorticity dependence on the vertical Lagrangian coordinate only (Gouyon waves), the shift of the wave number and the respective coefficient are constant. When the vorticity is dependent on both Lagrangian coordinates, the shift of the wave number is horizontally inhomogeneous. There are special cases (e.g., Gerstner waves) in which the vorticity is proportional to the squared wave amplitude and nonlinearity disappears, thus making the equations for wave packet dynamics linear. It is shown that the NLS solution for weakly rotational waves in the Eulerian variables may be obtained from the Lagrangian solution by simply changing the horizontal coordinates.