monthly_bshg function

Calculate monthly sunshine duration