"especially for transformers") kappa = ith or ip net["_options"] = {} _add_ppc_options(net, calculate_voltage_angles=False, trafo_model="pi", check_connectivity=check_connectivity, mode="sc", switch_rx_ratio=2, init_vm_pu="flat", init_va_degree="flat", enforce_q_lims=False, recycle=None) _add_sc_options(net, fault=fault, case=case, lv_tol_percent=lv_tol_percent, tk_s=tk_s, topology=topology, r_fault_ohm=r_fault_ohm, kappa_method=kappa_method, x_fault_ohm=x_fault_ohm, kappa=kappa, ip=ip, ith=ith, branch_results=branch_results, return_all_currents=return_all_currents) init_results(net, "sc") if fault == "3ph": _calc_sc(net)
logger.warning("Branch results are in beta mode and might not always be reliable, " "especially for transformers") # Convert bus to numpy array if bus is None: bus = net.bus.index.values else: bus = np.array([bus]).ravel() kappa = ith or ip net["_options"] = {} _add_ppc_options(net, calculate_voltage_angles=False, trafo_model="pi", check_connectivity=check_connectivity, mode="sc", switch_rx_ratio=2, init_vm_pu="flat", init_va_degree="flat", enforce_q_lims=False, recycle=None) _add_sc_options(net, fault=fault, case=case, lv_tol_percent=lv_tol_percent, tk_s=tk_s, topology=topology, r_fault_ohm=r_fault_ohm, kappa_method=kappa_method, x_fault_ohm=x_fault_ohm, kappa=kappa, ip=ip, ith=ith, branch_results=branch_results, return_all_currents=return_all_currents, inverse_y=inverse_y) init_results(net, "sc") if fault in ("2ph", "3ph"): _calc_sc(net, bus) elif fault == "1ph": _calc_sc_1ph(net, bus) else: raise ValueError("Invalid fault %s" % fault) def _calc_current(net, ppci_orig, bus):