Skip to content

Instantly share code, notes, and snippets.

import math
import matplotlib.pyplot as plt
def calculate_sundial_points(lat, wall_decl=90.0, g=100.0, max_hb=6, declinations=[-23.44, 0.0, 23.44]):
results = {}
dc_labels = ['winter', 'equinox', 'summer']
for hb in range(max_hb + 1):
points = []
for idx, dc in enumerate(declinations):