A spatial multivariate turning bands method is developed for non-LMC.
The method allows different anisotropies on the variables or model components.
The method is fast and easy to implement.
One billion points are simulated in less than one hour on a desktop computer.
Anisotropic covariances to simulate along the lines are derived.