Dear all,
I am trying to run gaussian unconditional simulation with spatial structure
using the grf function of geoR, with the following code:

library(geoR)
a=c(1:611)
b=c(1:634)
gr=expand.grid(a,b)
im=grf(661*634, grid=gr, nugget=nug, cov.par=c(sill,range), mean=mean)

When using geoR version 1.6.32 (or older) with version 1.3.45 (or older) of
RandomFields the results present an evident spatial structure. While using
the same code with version 1.6.33 or 1.6.34 of GeoR and version 2.0.4x of
RandomFields the resulting pattern is like pure random noise.
I tried with different values for sill, range, nugget and mean, but the
results are always without spatial structure. If needed, I can provide more
information and figures for more visual examples.

Any hints would be greatly appreciated.

Thank you very much in advance

Laura

	[[alternative HTML version deleted]]

