/** This is a first cut at calculating the great Jupiter-Saturn convergence of 2020. */ use planets.frink lat = 40 deg North long = 105 deg West for d = #2020-12-16# to #2020-12-27# step 30 min { // This can also be angularSeparation or airlessAngularSeparation sep = Planet.angularSeparation[Planet.Jupiter,Planet.Saturn,d,lat,long] [az,alt] = Planet.Jupiter.refractedAzimuthAltitude[d, lat, long] if (alt > 0 degrees) println[d + "\t" + format[sep, degrees, 4] + "\t" + skyDarkness[d,lat,long] + "\t" + format[alt, deg, 3]] }