set airmass < 0 = 0
clean up. add solar_constant parameter. suggest solpos/pyephem