sim_gns function

Simulate General Nested Model