code: def flow_orifice_vert(Diam, Height, RatioVCOrifice):
      
      #Checking input validity
      ut.check_range([RatioVCOrifice, "0-1", "VC orifice ratio"])
      if Height > -Diam / 2:
          flow_vert = integrate.quad(lambda z: (Diam * np.sin(np.arccos(z/(Diam/2)))
                                                     * np.sqrt(Height - z)
                                                     ),
                                                     - Diam / 2,
                                                     min(Diam/2, Height))
          return flow_vert[0] * RatioVCOrifice * np.sqrt(2 * gravity.magnitude)
      else:
          return 0
docs: Return the vertical flow rate of the orifice.