Utilities

using CairoMakie

# Define time range for the entire day (every 5 minutes)
start_time = ZonedDateTime(2023, 6, 21, 0, 0, tz_amsterdam)
end_time = ZonedDateTime(2023, 6, 21, 23, 59, tz_amsterdam)
times = collect(start_time:Minute(5):end_time)

# Compute solar positions for all times
positions = solar_position(obs, times)

# Get key events
events = transit_sunrise_sunset(obs, zdt)
sunrise_elev = solar_position(obs, events.sunrise).elevation
sunset_elev = solar_position(obs, events.sunset).elevation
transit_elev = solar_position(obs, events.transit).elevation

# Convert times to hours for plotting
times_hours = hour.(times) .+ minute.(times) ./ 60

# Create the plot with styling
fig = Figure(backgroundcolor=:transparent, textcolor="#f5ab35", size=(800, 400))
ax = Axis(fig[1, 1],
    xlabel="Time (hours)",
    ylabel="Solar Altitude (°)",
    title="Solar Altitude - Amsterdam, June 21, 2023",
    backgroundcolor=:transparent,
    xticks=0:3:24)

# Plot the solar altitude curve
lines!(ax, times_hours, positions.elevation,
    linewidth=2, color="#f5ab35")

# Add vertical markers and labels for events
sunrise_hour = hour(events.sunrise) + minute(events.sunrise) / 60
transit_hour = hour(events.transit) + minute(events.transit) / 60
sunset_hour = hour(events.sunset) + minute(events.sunset) / 60

vlines!(ax, sunrise_hour, linestyle=:dash, color=:gold, linewidth=1.5)
text!(ax, sunrise_hour, sunrise_elev + 5,
    text=Dates.format(events.sunrise, "HH:MM"),
    align=(:center, :bottom), color=:gold, fontsize=12)

vlines!(ax, transit_hour, linestyle=:dash, color=:red, linewidth=1.5)
text!(ax, transit_hour, transit_elev + 5,
    text=Dates.format(events.transit, "HH:MM"),
    align=(:center, :bottom), color=:red, fontsize=12)

vlines!(ax, sunset_hour, linestyle=:dash, color=:purple, linewidth=1.5)
text!(ax, sunset_hour, sunset_elev + 5,
    text=Dates.format(events.sunset, "HH:MM"),
    align=(:center, :bottom), color=:purple, fontsize=12)

# Add horizontal line at horizon
hlines!(ax, 0, linestyle=:dash, color=:gray, linewidth=1)

Makie.HLines{Tuple{Int64}}

# Generate dates for the entire year 2023 (every day)
year_start = ZonedDateTime(2023, 1, 1, 12, 0, tz_amsterdam)
year_dates = [year_start + Day(i) for i in 0:364]

# Calculate sunrise, sunset, and solar noon for each day
sunrise_times = Float64[]
sunset_times = Float64[]
solar_noon_times = Float64[]

for date in year_dates
    events = transit_sunrise_sunset(obs, date)

    # Convert to hours since midnight
    push!(sunrise_times, hour(events.sunrise) + minute(events.sunrise) / 60)
    push!(sunset_times, hour(events.sunset) + minute(events.sunset) / 60)
    push!(solar_noon_times, hour(events.transit) + minute(events.transit) / 60)
end

# Calculate daylight hours for each day
daylight_hours = sunset_times .- sunrise_times

# Find solstices (longest and shortest days)
summer_solstice_idx = argmax(daylight_hours)
winter_solstice_idx = argmin(daylight_hours)

# Create day of year array for x-axis
day_of_year = 1:365

# Create the sun graph
fig = Figure(backgroundcolor=:transparent, textcolor="#f5ab35", size=(900, 500))
ax = Axis(fig[1, 1],
    xlabel="Month",
    ylabel="Time (hours)",
    title="Sun Graph 2023 - Amsterdam (52.36°N, 4.88°E)",
    backgroundcolor=:transparent,
    yticks=0:2:24,
    xticks=(
        [1, 32, 60, 91, 121, 152, 182, 213, 244, 274, 305, 335],
        ["Jan", "Feb", "Mar", "Apr", "May", "Jun", "Jul", "Aug", "Sep", "Oct", "Nov", "Dec"]
    ))

xlims!(ax, 1, 365)

# Fill night time (top part: midnight to sunrise)
band!(ax, day_of_year, sunrise_times, fill(24.0, length(day_of_year)),
    color=(:darkslategray, 0.8))

# Fill night time (bottom part: sunset to midnight)
band!(ax, day_of_year, fill(0.0, length(day_of_year)), sunset_times,
    color=(:darkslategray, 0.8))

# Fill daylight time
band!(ax, day_of_year, sunrise_times, sunset_times,
    color=(:lightblue, 0.6))

# Plot solar noon line
lines!(ax, day_of_year, solar_noon_times,
    color=:red, linewidth=2, label="Solar Noon")

# Mark solstices
vlines!(ax, summer_solstice_idx, linestyle=:dash, color=:orange, linewidth=2)
text!(ax, summer_solstice_idx, 12,
    text="Summer\nSolstice",
    align=(:center, :bottom), color=:orange, fontsize=10, rotation=π/2)

vlines!(ax, winter_solstice_idx, linestyle=:dash, color=:steelblue, linewidth=2)
text!(ax, winter_solstice_idx, 12,
    text="Winter\nSolstice",
    align=(:center, :bottom), color=:steelblue, fontsize=10, rotation=π/2)

ylims!(ax, 0, 24)

┌ Warning: Assignment to `events` in soft scope is ambiguous because a global variable by the same name exists: `events` will be treated as a new local. Disambiguate by using `local events` to suppress this warning or `global events` to assign to the existing global variable.
└ @ utilities.md:154