Added function domain.azimuth_up() to calculate bearing (azimuth) of the upwards orientation of the domain. Replaces upwards_azimuth_direction which only calculates bearing in the center of the domain.