implement this

function run_powerflow!(inputs...)
    logger = PSY.configure_logging(
        console_level = PSY_CONSOLE_LEVEL,
        file_level = PSY_FILE_LEVEL,
    )
    try
        Logging.with_logger(logger) do
            PowerFlows.run_powerflow!(inputs...)
            return
        end
    finally
        close(logger)
        PSY.configure_logging(
            console_level = NODE_CONSOLE_LEVEL,
            file_level = NODE_FILE_LEVEL,
        )
    end
end

Implement modeling of more complex load models.

Hi, there are two bugs when reporting Pload and Qload of PQ buses:

Line 143 should have a minus before Vm_squared, as follows:

reactive_power += - Vm_squared * imag(PSY.get_Y(l))

Finally, lines 419 and 420 give results in pu, they need to be converted to MW and MVAr:

P_load_vect[ix] += P_admittance .* sys_basepower
Q_load_vect[ix] += Q_admittance .* sys_basepower

I was trying to run muti-period DC power flow on a System with StandardLoads and it fails because PSY.get_active_power(load::PSY.StandardLoad) isn't defined.

This is why it fails :

Here's the stacktrace:
julia> data_1 = PowerFlowData(DCPowerFlow(), pr_sys, timesteps);
ERROR: MethodError: no method matching get_active_power(::StandardLoad)
Closest candidates are:
get_active_power(::HydroDispatch) at ~/.conda-envs/julia17-sdhulipa/share/julia/packages/PowerSystems/9jm8v/src/models/generated/HydroDispatch.jl:120
get_active_power(::PowerSystems.ExponentialLoad) at ~/.conda-envs/julia17-sdhulipa/share/julia/packages/PowerSystems/9jm8v/src/models/generated/ExponentialLoad.jl:105
get_active_power(::InterruptiblePowerLoad) at ~/.conda-envs/julia17-sdhulipa/share/julia/packages/PowerSystems/9jm8v/src/models/generated/InterruptiblePowerLoad.jl:100
...
Stacktrace:
[1] _get_withdrawals!(bus_activepower_withdrawals::Vector{Float64}, bus_reactivepower_withdrawals::Vector{Float64}, bus_lookup::Dict{Int64, Int64}, sys::System)
@ PowerFlows ~/.conda-envs/julia17-sdhulipa/share/julia/packages/PowerFlows/GpXz5/src/common.jl:57
[2] PowerFlowData(::DCPowerFlow, sys::System, timesteps::Int64, timestep_names::Vector{String})
@ PowerFlows ~/.conda-envs/julia17-sdhulipa/share/julia/packages/PowerFlows/GpXz5/src/PowerFlowData.jl:278
[3] PowerFlowData(::DCPowerFlow, sys::System, timesteps::Int64)
@ PowerFlows ~/.conda-envs/julia17-sdhulipa/share/julia/packages/PowerFlows/GpXz5/src/PowerFlowData.jl:229
[4] top-level scope
@ REPL[7398]:1

Hi, I was wondering if you have experience with running power flow on power systems with some thousands buses and branches? I am working with a TSO and they seemed quite interested in having their system running on an open source software.

As far as I can tell the system does not contain any component types not supported by the library. However, even though I export the raw data with the solution obtained by PSS/E it does not converge. Have you had any similar experiences and have any tips on what I should look for?

Another thing I noticed is that several of the buses have a angles outside the +- 90 degrees, so I need to run the import runchecks=False to be able to import the system.

I guess this is more of a question than an issue, but I hope it is a question that is relevant to more than me and that you have some tips.

Best regards Sigurd

This issue is used to trigger TagBot; feel free to unsubscribe.

If you haven't already, you should update your TagBot.yml to include issue comment triggers.
Please see this post on Discourse for instructions and more details.

If you'd like for me to do this for you, comment TagBot fix on this issue.
I'll open a PR within a few hours, please be patient!