mcgf_rs_sim function

Simulate regime-switching Markov chain Gaussian field