code: def max_linear_flow(Diam, HeadlossCDC, Ratio_Error, KMinor):
      
      flow = (pc.area_circle(Diam)).magnitude * np.sqrt((2 * Ratio_Error * HeadlossCDC * pc.gravity)/ KMinor)
      return flow.magnitude
docs: Return the maximum flow that will meet the linear requirement.
      Maximum flow that can be put through a tube of a given diameter without
      exceeding the allowable deviation from linear head loss behavior