Kriging wiht anisotropy modeling using great circle distance in gstat [SEC=UNCLASSIFIED]

Hi All,

I tried to do kriging with anisotropy using great circle distance as bellow.

library(gstat)
loadMeuse()

require(rgdal)
  proj4string(meuse) = CRS("+init=epsg:28992")
  meuse.ll = spTransform(meuse, CRS("+proj=longlat"))

data(meuse.grid)
  coordinates(meuse.grid) = ~x+y
  proj4string(meuse.grid) = CRS("+init=epsg:28992")
  meuse.grid.ll = spTransform(meuse.grid, CRS("+proj=longlat"))

v =  variogram(zinc~1,meuse, alpha=c(0, 45, 90, 135), tol.hor=45)  
m1<-fit.variogram(v, vgm(1, "Exp", 3,  anis=c(90,0.1)))
  ok<-krige(log(zinc)~1, meuse.ll, meuse.grid.ll, m1)

[using ordinary kriging]
Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim,  : 
  gstat: value not allowed for: for long/lat data, anisotropy cannot be defined

It produced an error message as above. If the anisotropy specification is removed, it works.

This was done in:

Sys.info()

sysname                      release                      version 
  "Windows"                    "XP"         "build 2600, Service Pack 3" 
                                      
Is there anything wrong in my code or it is a bug? Any suggestions are appreciated!

Regards,
Jin