function n = valve_flow(Cv, z, rho, DP, MW)
n = Cv*z*sqrt(DP*rho/MW);
