Create 100_geostat.R