simulate_secondary function

Simulate secondary observations from primary observations