docs: def round_edge_screen(alpha, Re, angle=0):
      r'''Returns the loss coefficient for a round edged wire screen or bar
      screen, as shown in [1]_. Angle of inclination may be specified as well.
  
      Parameters
      ----------
      alpha : float
          Fraction of screen open to flow [-]
      Re : float
          Reynolds number of flow through screen with D = space between rods, []
      angle : float, optional
          Angle of inclination, with 0 being straight and 90 being parallel to
          flow [degrees]
  
      Returns
      -------
      K : float
          Loss coefficient [-]
  
      Notes
      -----
      Linear interpolation between a table of values. Re table extends
      from 20 to 400, with constant values outside of the table. This behavior
      should be adequate.
      alpha should be between 0.05 and 0.8.
      If angle is over 85 degrees, the value at 85 degrees is used.
  
      Examples
      --------
      >>> round_edge_screen(0.5, 100)
      2.0999999999999996
      >>> round_edge_screen(0.5, 100, 45)
      1.05
  
      References
      ----------
      .. [1] Blevins, Robert D. Applied Fluid Dynamics Handbook. New York, N.Y.:
         Van Nostrand Reinhold Co., 1984.
      '''
      beta = interp(Re, round_Res, round_betas)
      alpha2 = alpha*alpha
      K = beta*(1.0 - alpha2)/alpha2
      if angle:
          if angle <= 45:
              K *= cos(radians(angle))**2
          else:
              K *= interp(angle, round_thetas, round_gammas)
      return K