Calculate maps of azimuth angle in nav

[dahlek]

Some radiative transfer codes (cough cough nemesis) require measurements of azimuth angle (cw angle between the projection of the solar emission angle and north, see fig. 1 in Nemesis manual). I have a python code that does this using various vectors and the zenith and solar zenith angles but it would need to be reformatted to work within pylanetary