  We perform two high resolution direct numerical simulations of stratified
turbulence for Reynolds number equal to Re~25000 and Froude number respectively
of Fr~0.1 and Fr~0.03. The flows are forced at large scale and discretized on
an isotropic grid of 2048^3 points. Stratification makes the flow anisotropic
and introduces two extra characteristic scales with respect to homogeneous
isotropic turbulence: the buoyancy scale, L_B, and the Ozmidov scale, l_{oz}.
The former is related to the number of layers that the flow develops in the
direction of gravity, the latter is regarded as the scale at which isotropy is
recovered. The values of L_B and l_{oz} depend on the Froude number and their
absolute and relative size affect the repartition of energy among Fourier
modes. By contrasting the behavior of the two simulated flows we identify some
surprising similarities: after an initial transient the two flows evolve
towards comparable values of the kinetic and potential enstrophy, and energy
dissipation rate. Further similarities emerge at large scales: the same ratio
between potential and total energy (~0.1) is spontaneously selected by the
flows, and slow modes grow monotonically in both regimes causing a slow
increase of the total energy in time. The axisymmetric total energy spectrum
shows a wide variety of spectral slopes as a function of the angle between the
imposed stratification and the wave vector. One-dimensional energy spectra
computed in the direction parallel to gravity are flat from the forcing up to
buoyancy scale. At intermediate scales a ~ k^{-3} parallel spectrum develops
for the Fr ~ 0.03 run, whereas for weaker stratification, the saturation
spectrum does not have enough scales to develop and instead one observes a
power law compatible with Kolmogorov scaling. Finally, the spectrum of helicity
is flat until L_B, as observed in the nocturnal planetary boundary layer.
