def _calendar_month_middles(year):
<EOL>mdays = np.array(calendar.mdays[<NUM_LIT:1>:])<EOL>ydays = <NUM_LIT><EOL>if calendar.isleap(year):<EOL><INDENT>mdays[<NUM_LIT:1>] = mdays[<NUM_LIT:1>] + <NUM_LIT:1><EOL>ydays = <NUM_LIT><EOL><DEDENT>middles = np.concatenate(<EOL>[[-calendar.mdays[-<NUM_LIT:1>] / <NUM_LIT>],  <EOL>np.cumsum(mdays) - np.array(mdays) / <NUM_LIT>,  <EOL>[ydays + calendar.mdays[<NUM_LIT:1>] / <NUM_LIT>]])  <EOL>return middles<EOL>
List of middle day of each month, used by Linke turbidity lookup