Dear Jeffrey
As far as I know, RGeostats can only run kriging or simulations on the whole sphere.
By using variety.define function, you will activate the mode where all distances are calculated on the sphere (not euclidean).
Please, have a look to the following vignette:
http://rgeostats.free.fr/doc/Simsph.htmlHope this helps.
Fabien