docs: def Clift(Re):
      r'''Calculates drag coefficient of a smooth sphere using the method in
      [1]_ as described in [2]_.
  
      .. math::
          C_D = \left\{ \begin{array}{ll}
          \frac{24}{Re} + \frac{3}{16} & \mbox{if $Re < 0.01$}\\
          \frac{24}{Re}(1 + 0.1315Re^{0.82 - 0.05\log Re}) & \mbox{if $0.01 < Re < 20$}\\
          \frac{24}{Re}(1 + 0.1935Re^{0.6305}) & \mbox{if $20 < Re < 260$}\\
          10^{[1.6435 - 1.1242\log Re + 0.1558[\log Re]^2} & \mbox{if $260 < Re < 1500$}\\
          10^{[-2.4571 + 2.5558\log Re - 0.9295[\log Re]^2 + 0.1049[\log Re]^3} & \mbox{if $1500 < Re < 12000$}\\
          10^{[-1.9181 + 0.6370\log Re - 0.0636[\log Re]^2} & \mbox{if $12000 < Re < 44000$}\\
          10^{[-4.3390 + 1.5809\log Re - 0.1546[\log Re]^2} & \mbox{if $44000 < Re < 338000$}\\
          9.78 - 5.3\log Re & \mbox{if $338000 < Re < 400000$}\\
          0.19\log Re - 0.49 & \mbox{if $400000 < Re < 1000000$}\end{array} \right.
  
      Parameters
      ----------
      Re : float
          Reynolds number of the sphere, [-]
  
      Returns
      -------
      Cd : float
          Drag coefficient [-]
  
      Notes
      -----
      Range is Re <= 1E6.
  
      Examples
      --------
      >>> Clift(200)
      0.7756342422322543
  
      References
      ----------
      .. [1] R. Clift, J.R. Grace, M.E. Weber, Bubbles, Drops, and Particles,
         Academic, New York, 1978.
      .. [2] Barati, Reza, Seyed Ali Akbar Salehi Neyshabouri, and Goodarz
         Ahmadi. "Development of Empirical Models with High Accuracy for
         Estimation of Drag Coefficient of Flow around a Smooth Sphere: An
         Evolutionary Approach." Powder Technology 257 (May 2014): 11-19.
         doi:10.1016/j.powtec.2014.02.045.
      '''
      if Re < 0.01:
          Cd = 24./Re + 3/16.
      elif Re < 20:
          Cd = 24./Re*(1 + 0.1315*Re**(0.82 - 0.05*log10(Re)))
      elif Re < 260:
          Cd = 24./Re*(1 + 0.1935*Re**(0.6305))
      elif Re < 1500:
          Cd = 10**(1.6435 - 1.1242*log10(Re) + 0.1558*(log10(Re))**2)
      elif Re < 12000:
          Cd = 10**(-2.4571 + 2.5558*log10(Re) - 0.9295*(log10(Re))**2 + 0.1049*log10(Re)**3)
      elif Re < 44000:
          Cd = 10**(-1.9181 + 0.6370*log10(Re) - 0.0636*(log10(Re))**2)
      elif Re < 338000:
          Cd = 10**(-4.3390 + 1.5809*log10(Re) - 0.1546*(log10(Re))**2)
      elif Re < 400000:
          Cd = 29.78 - 5.3*log10(Re)
      else:
          Cd = 0.19*log10(Re) - 0.49
      return Cd