Difficulty with utilizing PowerModels.update data! function

I think what you are looking for is something like this code below. I would caution, setting up an AC PowerFlow study from a DC dispatch solution has a number of subtle points (e.g. line losses) that should be treated with care.

using PowerModels, Ipopt

case_file = joinpath(dirname(pathof(PowerModels)), "..", "test", "data", "matpower", "case9.m")

network_data = parse_file(case_file)

result = solve_opf(network_data, DCPPowerModel, Ipopt.Optimizer)
print_summary(result["solution"])

for (i,gen) in result["solution"]["gen"]
    network_data["gen"][i]["pg"] = gen["pg"]
end

result2 =solve_pf(network_data, ACPPowerModel, Ipopt.Optimizer)
print_summary(result2["solution"])
1 Like