def METHOD_NAME(start_year, time, season='djf'):
  ''' Get the time index for the given season '''
  # convert time as numpy array
  time = np.asarray(time)
  # getting the datetime values for the time index
  dates_month=[]
  dates_year=[]
  for i_time in time: 
    temp_time = dt.datetime(start_year, 1, 1) + dt.timedelta(days=np.float(i_time)-1)
    dates_month.append(temp_time.month)
    dates_year.append(temp_time.year)
  dates_month = np.asarray(dates_month)
  # getting the time index
  if (season == ''):
    raise Exception('Set which season you want to extract!')
  elif (season == 'djf'):
    time_ind = (dates_month == 12) | (dates_month == 1) | (dates_month == 2)
  elif (season == 'mam'):
    time_ind = (dates_month == 3) | (dates_month == 4) | (dates_month == 5)
  elif (season == 'jja'):
    time_ind = (dates_month == 6) | (dates_month == 7) | (dates_month == 8)
  elif (season == 'son'):
    time_ind = (dates_month == 9) | (dates_month == 10) | (dates_month == 11)
  return time_ind