Add argument 'read_xy_from_geo_dat' and set to TRUE