Does anyones knows about Kriging interpolation method? Is there a example fortran code ?

Any help is appreciated.