PSSE Exporter Part 1: Add Tests and Formalize Exporter Implementation

.gitignore

@@ -1,4 +1,7 @@
 test_simulation_results/*
+# TODO reevaluate next few entries
+test/Raw_Export
+test/test_psse_exports
 #Files generated by invoking Julia with --code-coverage
 *.jl.cov
 *.jl.*.cov

Project.toml

@@ -7,19 +7,23 @@ version = "0.7.0"
 DataFrames = "a93c6f00-e57d-5684-b7b6-d8193f3e46c0"
 DataStructures = "864edb3b-99cc-5e75-8d2d-829cb0a9cfe8"
 Dates = "ade2ca70-3891-5945-98fb-dc099432e06a"
-DelimitedFiles = "8bb1440f-4735-579b-a4ab-409b98df4dab"
 InfrastructureSystems = "2cd47ed4-ca9b-11e9-27f2-ab636a7671f1"
-JSON = "682c06a0-de6a-54ab-a142-c8b1cf79cde6"
+JSON3 = "0f8b85d8-7281-11e9-16c2-39a750bddbf1"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
+Logging = "56ddb016-857b-54e1-b83d-db4d58db5568"
 NLsolve = "2774e3e8-f4cf-5e23-947b-6d7e65073b56"
 PowerNetworkMatrices = "bed98974-b02a-5e2f-9fe0-a103f5c450dd"
 PowerSystems = "bcd98974-b02a-5e2f-9ee0-a103f5c450dd"
 SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 
 [compat]
 DataFrames = "1"
+DataStructures = "0.18"
+Dates = "1"
 InfrastructureSystems = "2"
+JSON3 = "1"
 LinearAlgebra = "1"
+Logging = "1"
 NLsolve = "4"
 PowerNetworkMatrices = "^0.11"
 PowerSystems = "^4"

src/PowerFlows.jl

@@ -8,15 +8,22 @@ export ACPowerFlow
 export PTDFDCPowerFlow
 export vPTDFDCPowerFlow
 export write_results
+export PSSEExporter
+export update_exporter!
+export write_export
+export get_psse_export_paths
 
+import Logging
 import DataFrames
 import PowerSystems
+import PowerSystems: System
 import LinearAlgebra
 import NLsolve
 import SparseArrays
 import InfrastructureSystems
 import PowerNetworkMatrices
 import SparseArrays: SparseMatrixCSC
+import JSON3
 
 const IS = InfrastructureSystems
 const PSY = PowerSystems
@@ -31,15 +38,13 @@ include("ac_power_flow.jl")
 include("ac_power_flow_jacobian.jl")
 include("nlsolve_ac_powerflow.jl")
 include("post_processing.jl")
+include("psse_export.jl")
 
-# PSSE Exporter
-import PowerSystems: System
+# Old PSSE Exporter
 import DataFrames: DataFrame
 import Dates
+import DataStructures
 import DataStructures: OrderedDict
-import DelimitedFiles
-import JSON
 include("psse_exporter/support_tools.jl")
-include("psse_exporter/psse_exporter.jl")
 
 end

src/post_processing.jl

@@ -692,3 +692,32 @@ function write_results(
 
     return Dict("bus_results" => bus_df, "flow_results" => branch_df)
 end
+
+"""
+Modify the values in the given `System` to correspond to the given `PowerFlowData` such that
+if a new `PowerFlowData` is constructed from the resulting system it is the same as `data`.
+See also `write_powerflow_solution!`.
+"""
+function update_system!(sys::PSY.System, data::PowerFlowData)
+    for bus in PSY.get_components(PSY.Bus, sys)
+        if bus.bustype == PSY.ACBusTypes.REF
+            # For REF bus, voltage and angle are fixed; update active and reactive
+            P_gen = data.bus_activepower_injection[data.bus_lookup[PSY.get_number(bus)]]
+            Q_gen = data.bus_reactivepower_injection[data.bus_lookup[PSY.get_number(bus)]]
+            _power_redistribution_ref(sys, P_gen, Q_gen, bus,
+                DEFAULT_MAX_REDISTRIBUTION_ITERATIONS)
+        elseif bus.bustype == PSY.ACBusTypes.PV
+            # For PV bus, active and voltage are fixed; update reactive and angle
+            Q_gen = data.bus_reactivepower_injection[data.bus_lookup[PSY.get_number(bus)]]
+            _reactive_power_redistribution_pv(sys, Q_gen, bus,
+                DEFAULT_MAX_REDISTRIBUTION_ITERATIONS)
+            PSY.set_angle!(bus, data.bus_angles[data.bus_lookup[PSY.get_number(bus)]])
+        elseif bus.bustype == PSY.ACBusTypes.PQ
+            # For PQ bus, active and reactive are fixed; update voltage and angle
+            Vm = data.bus_magnitude[data.bus_lookup[PSY.get_number(bus)]]
+            PSY.set_magnitude!(bus, Vm)
+            PSY.set_angle!(bus, data.bus_angles[data.bus_lookup[PSY.get_number(bus)]])
+        end
+    end
+    # TODO
+end