code: def vel_horizontal(HeightWaterCritical):
      
      #Checking input validity
      ut.check_range([HeightWaterCritical, ">0", "Critical height of water"])
      return np.sqrt(gravity.magnitude * HeightWaterCritical)
docs: Return the horizontal velocity.