# manhattanhengesun.frink

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 19355 days, 10 hours, 25 minutes ago.