docs: def get_mask(cls, azim):
          """Linear interpolation between two points of the mask
          """
  
          if cls.mask is None:
              raise ValueError("No mask defined for the station {}".format(cls.name))
  
          azim %= 2 * np.pi
  
          if azim in cls.mask[0, :]:
              return cls.mask[1, np.where(azim == cls.mask[0, :])[0][0]]
  
          for next_i, mask_azim in enumerate(cls.mask[0, :]):
              if mask_azim > azim:
                  break
          else:
              next_i = 0
  
          x0, y0 = cls.mask[:, next_i - 1]
          x1, y1 = cls.mask[:, next_i]
  
          if next_i - 1 == -1:
              x0 = 0
  
          return y0 + (y1 - y0) * (azim - x0) / (x1 - x0)