wbt_turning_bands_simulation function

Turning bands simulation