Download or view manhattanhengesun.frink in plain text format
// Program to calculate sun alignments with Manhattan's East-West streets.
use sun.frink
North = +1
South = -1
West = +1
East = -1
lat = DMS[40,45,261.16] North
long = DMS[73,59,9.02] West
// Approximate direction of the streets.
// see http://en.wikipedia.org/wiki/Manhattanhenge
// http://www.charlespetzold.com/etc/AvenuesOfManhattan/index.html
axis = (270 + 29.008) degrees
// Conversion to coordinate system used by Meeus
axisMeeus = (axis + 180 degrees) mod circle
println[axisMeeus -> degrees]
date = #2016-01-01 12:01 PM#
dateOut = ### yyyy-MM-dd hh:mm:ss a zzz ###
sep = "\t"
while (date <= #2020-12-31#)
{
date = sunSecantAzimuth[date, lat, long, axisMeeus]
[azimuth, altitude] = refractedSunAzimuthAltitude[date, lat, long]
actualSunset = sunset[date, lat, long]
[actualSunsetAz, actualSunsetAlt] = refractedSunAzimuthAltitude[actualSunset, lat, long]
actualSunsetAz = (actualSunsetAz + 180 degrees) mod circle
print[(date -> [dateOut, "Eastern"]) + "$sep"]
print[format[JD[date],day,5] + "$sep"]
print[format[altitude,degrees,2] + "$sep"]
print[format[actualSunsetAz, degrees, 2] + "$sep"]
println[]
date = date + 1 day
}
Download or view manhattanhengesun.frink in plain text format
This is a program written in the programming language Frink.
For more information, view the Frink
Documentation or see More Sample Frink Programs.
Alan Eliasen was born 20143 days, 10 hours, 22 minutes ago.