I'm trying to get the position of the sun in the sky at a given time and location. Similar to this: https://discourse.julialang.org/t/calculate-positions-of-astronomical-objects-from-a-terrestrial-observatory/107772

sunpos returns right ascension and declanation and there exists hadec2altaz, but I don't want to manually convert from right ascension to hour angle.

I'd like to be able to call radec2altaz(sunpos(time)[1:2], lat, lon, time) and get an altitude and azimuth.