can be used to perform sequential Gaussian simulation using the same arguments as “krig”. Only extra options is that the number of realizations computed can be set using

options.nsim = 10;