  We study numerically the nonlinear stage of modulational instability (MI) of
cnoidal waves, in the framework of the focusing one-dimensional Nonlinear
Schrodinger (NLS) equation. Cnoidal waves are the exact periodic solutions of
the NLS equation and can be represented as a lattice of overlapping solitons.
MI of these lattices lead to development of "integrable turbulence" [Zakharov
V.E., Stud. Appl. Math. 122, 219-234 (2009)]. We study the major
characteristics of the turbulence for dn-branch of cnoidal waves and
demonstrate how these characteristics depend on the degree of "overlapping"
between the solitons within the cnoidal wave.
  Integrable turbulence, that develops from MI of dn-branch of cnoidal waves,
asymptotically approaches to it's stationary state in oscillatory way. During
this process kinetic and potential energies oscillate around their asymptotic
values. The amplitudes of these oscillations decay with time as t^{-a},
1<a<1.5, the phases contain nonlinear phase shift decaying as t^{-1/2}, and the
frequency of the oscillations is equal to the double maximal growth rate of the
MI, s=2g_{max}. In the asymptotic stationary state the ratio of potential to
kinetic energy is equal to -2. The asymptotic PDF of wave amplitudes is close
to Rayleigh distribution for cnoidal waves with strong overlapping, and is
significantly non-Rayleigh one for cnoidal waves with weak overlapping of
solitons. In the latter case the dynamics of the system reduces to two-soliton
collisions, which occur with exponentially small rate and provide up to
two-fold increase in amplitude compared with the original cnoidal wave.
