docs: def Barati_high(Re):
      r'''Calculates drag coefficient of a smooth sphere using the method in
      [1]_.
  
      .. math::
          C_D = 8\times 10^{-6}\left[(Re/6530)^2 + \tanh(Re) - 8\ln(Re)/\ln(10)\right]
          - 0.4119\exp(-2.08\times10^{43}/[Re + Re^2]^4)
          -2.1344\exp(-\{[\ln(Re^2 + 10.7563)/\ln(10)]^2 + 9.9867\}/Re)
          +0.1357\exp(-[(Re/1620)^2 + 10370]/Re)
          - 8.5\times 10^{-3}\{2\ln[\tanh(\tanh(Re))]/\ln(10) - 2825.7162\}/Re
          + 2.4795
  
      Parameters
      ----------
      Re : float
          Particle Reynolds number of the sphere using the surrounding fluid
          density and viscosity, [-]
  
      Returns
      -------
      Cd : float
          Drag coefficient [-]
  
      Notes
      -----
      Range is Re <= 1E6
      This model is the wider-range model the authors developed.
      At sufficiently low diameters or Re values, drag is no longer a phenomena.
  
      Examples
      --------
      Matching example in [1]_, in a table of calculated values.
  
      >>> Barati_high(200.)
      0.7730544082789523
  
      References
      ----------
      .. [1] 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.
      '''
      Re2 = Re*Re
      t0 = 1.0/Re
      t1 = (Re/6530.)
      t2 = (Re/1620.)
      t3 = log10(Re2 + 10.7563)
      tanhRe = tanh(Re)
      Cd = (8E-6*(t1*t1 + tanhRe - 8.0*log10(Re))
      - 0.4119*exp(-2.08E43/(Re+Re2)**4)
      - 2.1344*exp(-t0*(t3*t3 + 9.9867))
      + 0.1357*exp(-t0*(t2*t2 + 10370.))
      - 8.5E-3*t0*(2.0*log10(tanh(tanhRe)) - 2825.7162) + 2.4795)
      return Cd