daily_bshg function

Calculate daily sunshine duration