docs: def K_globe_stop_check_valve_Crane(D1, D2, fd=None, style=0):
      r'''Returns the loss coefficient for a globe stop check valve as shown in 
      [1]_.
          
      If β = 1:
          
      .. math::
          K = K_1 = K_2 = N\cdot f_d
          
      Otherwise:
      
      .. math::
          K_2 = \frac{K + \left[0.5(1-\beta^2) + (1-\beta^2)^2\right]}{\beta^4}
          
      Style 0 is the standard form; style 1 is angled, with a restrition to force
      the flow up through the valve; style 2 is also angled but with a smaller
      restriction forcing the flow up. N is 400, 300, and 55 for those cases
      respectively.
      
      Parameters
      ----------
      D1 : float
          Diameter of the valve seat bore (must be smaller or equal to `D2`), [m]
      D2 : float
          Diameter of the pipe attached to the valve, [m]
      fd : float, optional
          Darcy friction factor calculated for the actual pipe flow in clean 
          steel (roughness = 0.0018 inch) in the fully developed turbulent 
          region; do not specify this to use the original Crane friction factor!,
          [-]
      style : int, optional
          One of 0, 1, or 2; refers to three different types of angle valves
          as shown in [1]_ [-]
          
      Returns
      -------
      K : float
          Loss coefficient with respect to the pipe inside diameter [-]
  
      Notes
      -----    
      This method is not valid in the laminar regime and the pressure drop will
      be underestimated in those conditions.
  
      Examples
      --------
      >>> K_globe_stop_check_valve_Crane(.1, .02, style=1)
      4.5235076518969795
      
      References
      ----------
      .. [1] Crane Co. Flow of Fluids Through Valves, Fittings, and Pipe. Crane,
         2009.
      '''
      if fd is None:
          fd = ft_Crane(D2)
      try:
          K = globe_stop_check_valve_Crane_coeffs[style]*fd
      except KeyError:
          raise KeyError('Accepted valve styles are 0, 1, and 2 only')
      beta = D1/D2
      if beta == 1:
          return K
      else:
          return (K + beta*(0.5*(1 - beta**2) + (1 - beta**2)**2))/beta**4