use sun.frink
use planets.frink

startdate = #2026-01-03 00:00#
enddate =   #2026-01-05 00:00#
lat = 40 deg North
long = 105 deg West

// Disproves picture of a Moon-Jupiter convergence from Twitter
// https://x.com/M1ONLOVER/status/2007819322333577246
date = #2026-01-03 03:00 PM#
g = new graphics

//for date = startdate to enddate step 10 min
//{
mr = moonRadiusAngle[date]
[maz, malt] = refractedMoonAzimuthAltitude[date, lat, long]
g.add[drawMoonPolygonRelativeToZenith[date, lat, long, maz, -malt, mr, true]]

jr = Planet.Jupiter.radiusAngle[date]
[jaz, jalt] = Planet.Jupiter.refractedAzimuthAltitude[date, lat, long]
g.add[Planet.Jupiter.drawJupiterMoonsRelativeToZenith[date, lat, long, jaz, -jalt, jr, false, true]]
println["Radius ratio: " + (mr/jr)]
println["$date: " + (angularSeparation[maz, malt, jaz, jalt] -> "deg")]
//}

g.invertGrays[].show[]
g.print[]
g.invertGrays[].write["MoonVsJupiter.png", 1000, undef]
g.write["MoonVsJupiter.svg", 1000, undef]