Real large-scale electrical power systems offer several difficulties when performing network analysis. The AC power flow model is a non-convex and non-linear problem which is difficult to solve quickly and reliably. Additionally, real power system often suffers from voltage instability and are composed with several flexible AC transmission system (FACTS) devices and controllable components to accommodate the adversities and guarantee convergence. For that reason, market consolidated power system software s possesses features that handle the characteristics of large-scale networks that are not accounted in tradition power flow formulations. In this context, Brazil is a country of continental dimension and relies on controllable components to handle electrical instability. In Brazilian electric systems representation, it s common to use these components to get more realistic results, handle bad voltage profiles and guarantee convergence. Examples of control actions in Brazilian load flow analysis are: (i) consideration of reactive generation limits and voltage profile limits, (ii) usage of shunt controller and load tap changer transformers (LTC) to handle voltage instability, and (iii) usage phase shifting control to handle active power flow in transmission lines. The present report focusses on two contributions to the power system community: (i) a Julia Language package for parsing the Brazilian power system file (ANAREDE file) into Julia - PWF.jl - and (ii) a Julia package for performing power flow analysis with control action options – ControlPowerFlow.jl, which are both available as open-source projects. The practical validation of those packages is demonstrated in this work through case studies.