added z_ids and fixed z_n in gapfill_georegions() with @bbest