code: def _season_from_months(months):
      
      # TODO: Move "season" accessor upstream into pandas
      seasons = np.array(['DJF', 'MAM', 'JJA', 'SON'])
      months = np.asarray(months)
      return seasons[(months // 3) % 4]
docs: Compute season (DJF, MAM, JJA, SON) from month ordinal