get_monolayer_distance_map() should explicitly require construct_zone_axes()
construct_zone_axes() only gives a TypeError: 'NoneType' object is not iterable
if the sublattice.construct_zone_axes() has not been called. It would be nice if it informed the user of what to do.