annual_bshg function

Calculate annual sunshine duration