mcgf_sim function

Simulate Markov chain Gaussian field