Allow bulk writing and reading of time series

src/PowerSystems.jl

@@ -274,6 +274,7 @@ export TimeSeriesCounts
 export get_dynamic_components
 
 export parse_file
+export open_time_series_store!
 export add_time_series!
 export remove_time_series!
 export check_time_series_consistency

src/base.jl

@@ -704,6 +704,34 @@ function add_service!(
     return
 end
 
+"""
+Open the time series store for bulk additions or reads. This is recommended before calling
+add_time_series! many times because of the overhead associated with opening and closing an
+HDF5 file.
+
+This is not necessary for an in-memory time series store.
+
+# Examples
+```julia
+# Assume there is a system with an array of components and SingleTimeSeries
+open_time_series_store!(sys, "r+") do
+    for (component, ts) in zip(components, single_time_series)
+        add_time_series!(sys, component, ts)
+    end
+end
+```
+julia>
+"""
+function open_time_series_store!(
+    func::Function,
+    sys::System,
+    mode = "r",
+    args...;
+    kwargs...,
+)
+    IS.open_time_series_store!(func, sys.data, mode, args...; kwargs...)
+end
+
 """
 Add time series data from a metadata file or metadata descriptors.
 

test/test_system.jl

@@ -248,6 +248,39 @@ end
     end
 end
 
+@testset "Test bulk add of time series" begin
+    sys = System(100.0)
+    bus = ACBus(nothing)
+    bus.bustype = ACBusTypes.REF
+    add_component!(sys, bus)
+    components = []
+    len = 2
+    component = ThermalStandard(nothing)
+    component.name = "gen"
+    component.bus = bus
+    add_component!(sys, component)
+    initial_time = Dates.DateTime("2020-09-01")
+    resolution = Dates.Hour(1)
+    len = 24
+    timestamps = range(initial_time; length = len, step = resolution)
+    arrays = [TimeSeries.TimeArray(timestamps, rand(len)) for _ in 1:5]
+    ts_name = "test"
+
+    open_time_series_store!(sys, "r+") do
+        for (i, ta) in enumerate(arrays)
+            ts = SingleTimeSeries(; data = ta, name = "$(ts_name)_$(i)")
+            add_time_series!(sys, component, ts)
+        end
+    end
+
+    open_time_series_store!(sys, "r") do
+        for (i, expected_array) in enumerate(arrays)
+            ts = IS.get_time_series(IS.SingleTimeSeries, component, "$(ts_name)_$(i)")
+            @test ts.data == expected_array
+        end
+    end
+end
+
 @testset "Test set_name! of system component" begin
     sys = System(100.0)
     bus = ACBus(nothing)