renamed ha to hour angle for clarity
added geometric mean longitude of sun in deg calculation
added geometric mean anomoly of sun in deg