Update 100_geostat.R