library("akima")
source("zimm.fn.s")

aaa.999<-newmat(sigma=.999)
aaa.1<-newmat(sigma=.1)
aaa.01<-newmat(sigma=.01)
aaa.001<-newmat(sigma=.001)
aaa.0001<-newmat(sigma=.0001)
aaa.00001<-newmat(sigma=.00001)
aaa.000001<-newmat(sigma=.00001)
aaa.000001<-newmat(sigma=.000001)

cursigma<-0.001