Added truncation and turbulent initialization.

Added routines to initialize a solenoidal velocity field fitting a desired
spectrum. The mathematics behind the initialization needs to be checked. Also
added routines for mode truncation according to the 3/2 rule.