Getting Started

Creating a PSRIO Script

In this guide, we will show you how to create PSRIO recipes from scratch using output data from example cases inside the SDDP installation folder. For the purpose of this tutorial, we will use example case 20, which consists of a single system with three thermal plants supplying energy to consumers over 12 stages, with deterministic execution.

Example

To calculate the average demand marginal cost over the study horizon from the cmgdem.csv file, we first load the system collection with the following command:

system = System();

Then, we load the marginal cost data from the file:

cmgdem = system:load("cmgdem");

The data can be seen in the image below:

The processing is accomplished using the method aggregate_stages, with the parameter BY_AVERAGE():

cmgdem_avg = cmgdem:aggregate_stages(BY_AVERAGE());

We save the expression with the following method, and the PSRIO will output the cmgdem_avg.bin and cmgdem_avg.hdr.

cmgdem_avg:save("cmgdem_avg");

If we want to save a CSV file, we should use the additional parameter:

cmgdem_avg:save("cmgdem_avg", { csv = true });

Let's see what the final result looks like:

As we can see, the stages collapsed into one single value representing the average value of the demand marginal costs over the study horizon.

Example

Suppose we want to extract the total amount of energy generated by the plants based on the information received from SDDP about the generations in each stage of each plant, which we found in the file gerter.csv. First, we load the thermal collection with the following command:

thermal = Thermal();

Then, we load the thermal generation data from the file:

gerter = thermal:load("gerter");

Check the data contained in the file:

Notice that we managed to load the generation data by informing the output filename. Now, we can make operations with the thermal generation data. We want to sum the values of each plant regarding one stage so that there is only one value per stage representing the total generation of the system.

total_gerter = gerter:aggregate_agents(BY_SUM(), "Total thermal");

total_gerter:save("total_gerter", { csv = true });

The first instruction is to aggregate agents, summing their respective generations for each stage. In the subsequent command, we are saving the data stored in the variable total_gerter to the file total_gerter.csv. The processed data looks like this:

"The most effective approach to visualize the output is through data plotting. To start, let's define a chart object."

chart = Chart("My first chart");

We indicate to the chart object that we want to plot the data in the file variable total_gerter, which we have just created.

chart:add_line(total_gerter);

After that, we will define a tab object to hold chart objects information and then add the chart and save it to the file thermal_report.html.

tab = Tab("My first tab");
tab:push(chart);

dashboard = Dashboard();
dashboard:push(tab);
dashboard:save("thermal_report");

By clicking in the file thermal_report.html, that's what we get:

In short, the recipe for this work is as follows:

-- Load thermal collections and thermal generation data
thermal = Thermal();
gerter = thermal:load("gerter");

-- Aggregate agents generations by sum and save to output file
total_gerter = gerter:aggregate_agents(BY_SUM(), "Total thermal");
total_gerter:save("total_gerter", { csv = true });

-- Create chart object
chart = Chart("My first chart");

-- Add total generation info
chart:add_line(total_gerter);

-- Create dashboad and tab object
dashboard = Dashboard();
tab = Tab("My first tab");

-- Add chart to tab and tab to dashboard
tab:push(chart);
dashboard:push(tab);

-- Save dashboard to html file
tab:save("thermal_report");

There are two ways to run the script above-mentioned. In the following sections, we will explain each one of them.

Running a script

Graph 4

Graph 4 is a powerful graphing tool that allows users to make PSRIO scripts within an editor. Once the SDDP execution is complete, you can open Graph 4 by following these steps:

Select the PSRIO editor option from the menu on your screen's left side:

You can create a new script, write down the code, save it, and you will be ready to go:

SDDP

Another way to run a PSRIO script is by using SDDP. You can automate running a PSRIO script by saving an sddp.lua file in the case folder. The SDDP model will run this script and produce the outputs and tab defined on it.

Following either method, you can quickly run the script and get the desired outputs.

Command Prompt

To run PSRIO from the command prompt (CMD), you can use the PSRIO.exe executable. Typically, you can find it in the path "C:\PSR\GraphModule\Oper\psrplot\psrio." To execute PSRIO.exe, you need to provide the following parameters:

<executable_path>\PSRIO.exe -r <script_path>\<script_name>.lua <sddp_case_path>

Here's an example command:

"C:\PSR\GraphModule\Oper\psrplot\psrio\PSRIO.exe" -r "D:\scripts\psrio.lua" "D:\case"

If you want to specify an output path, you can include it as a parameter in the psrio.exe call:

<executable_path>\PSRIO.exe -r <script_path>\<script_name>.lua -o <output_path> <sddp_case_path>

Here's an example command with an output path:

"C:\PSR\GraphModule\Oper\psrplot\psrio\PSRIO.exe" -r "D:\scripts\psrio.lua" -o "D:\outputs" "D:\case"

Changelog

1.2.1 (xxxx-xx-xx)

https://psrio.psr-inc.com/1.2.1/index.html

Added

Added (Collection)

Added (Input Data)

Syntax Changes

Changed

Fixed

Removed

Submodules

1.2.0 (2024-10-21)

https://psrio.psr-inc.com/1.2.0/index.html

Added

• Interactive charts with animation controls
• Spreadsheet() and spreadsheet:add(exp, string)
• exp:rainflow_count(exp1, range_min, range_max)
• exp:initial_date()
• exp:final_date()
• exp:data_info()
• exp:select_blocks({ i, j, k }) and exp:select_blocks(i, k)
• circuit:bus_to(f)
• circuit:bus_from(f)
• chart:enable_controls()
• chart:invert_axes()
• chart:xlabel_datetime_format(string)
• chart:add_line(exp, { data_format = '...' }), ...
• chart:add_box_plot_categories(exp1, exp2, exp3, exp4, exp5, label, options)
• chart:add_heatmap_agents(exp)
• chart:enable_vertical_zoom()
• chart:horizontal_legend
• BY_AVERAGE_EXCLUDING(double) aggregation
• BY_SUM_EXCLUDING(double) aggregation
• BY_PERCENTILE_EXCLUDING(int, double) aggregation
• PSR.console_verbose_level(i)
• PSR.console_verbose_level()
• --skip_typical_days_validation argument

Added (Collection)

• ElectrificationFixedConverter
• ElectrificationFixedCommodity

Added (Input Data)

• circuit.decommissioned
• hydro.ramp_up
• hydro.ramp_down
• hydro.bid_price
• expansion_project.integration_cost
• expansion_project.invest_cost_unit
• expansion_project.invest_cost
• expansion_project.annualized_invest_cost
• thermal.commitment_type
• battery.charge_ramp
• battery.discharge_ramp
• battery.max_reserve
• battery.bid_price
• renewable.max_reserve
• renewable.bid_price

Submodules

• lua 5.4.6 -> 5.4.7

Changelog

1.1.0 (2024-05-22)

https://psrio.psr-inc.com/1.1.0/index.html

Added

• exp:stage_profile_day(f)
• exp:cumsum_stages()
• exp:select_scenarios_kth_largest(i)
• exp:select_scenarios_kth_smallest(i)
• exp:convert()
• exp:aggregate_stages(..., Profile.PER_SEMESTER)
• exp:aggregate_stages(..., Profile.SEMESTER)
• exp:reshape_stages(Profile.PER_SEMESTER)
• chart:add_line_categories(...)
• chart:add_line_stacking_categories(...)
• chart:add_spline_categories(...)
• chart:add_column_stacking_categories(...)
• chart:add_column_percent_categories(...)
• --load_format argument

Added (Input Data)

• Battery().max_storage
• Battery().min_storage
• Battery().om_cost

Syntax Changes

• chart:add_categories(...) to chart:add_column_categories(...)

Removed

• --binf_rule argument

Submodules

• PSRClasses 6a8cd0db -> c4e86064

Changelog

1.0.0 (2024-03-17)

https://psrio.psr-inc.com/1.0.0/index.html

Added

• responsive containers to side by side charts
• marketplace documentation
• exp:save(..., { index = false })
• toml:get_keys()
• toml:get_table(key)
• toml:get_table(key, fallback)
• toml:get_table_array(key)
• --binf_rule argument
• hydro.discharge_ramp_up
• hydro.discharge_ramp_down
• chart:add_waterfall(...)
• collection:force_load(file_name)
• concentrated_solar_power.system

Added (Input Data)

• demand.is_enabled

Changed

• major performance improvements in exp:select_agents(query) method
• major performance improvements in exp:remove_zeros() method
• minor performance improvements in exp:convert(...) method
• toml:get_string(key, fallback) where fallback is dynamic
• toml:get_string_array(key, fallback) where fallback is dynamic
• toml:get_integer(key, fallback) where fallback is dynamic
• toml:get_integer_array(key, fallback) where fallback is dynamic
• toml:get_real(key, fallback) where fallback is dynamic
• toml:get_real_array(key, fallback) where fallback is dynamic
• toml:get_boolean(key, fallback) where fallback is dynamic

Fixed

• GasNode collection
• collection:print_elements() index
• error in exp:aggregate_blocks() when exp is null

Removed

• --load_only_csv argument

Submodules

• bootstrap 5.3.2 -> 5.3.3
• bootstrap-icons 1.11.1 -> 1.11.2
• github-markdown-css 5.4.0 -> 5.5.1
• json 3.11.2 -> 3.11.3
• psrplot 0.11.0 -> 0.12.0

Changelog

0.27.0 (2023-12-12)

https://psrio.psr-inc.com/0.27.0/index.html

Added

• datetime library
• table extension
• exp:collection()
• exp:study_index()
• exp:reshape_stages(Profile.QUARTER)
• exp:to_block() and exp:to_hour() with automatic unit solver
• chart:add_thermal_merit_order_curve(...)
• collection:ids()
• collection:get_study_index()
• study:get_relationship_array(...)
• toml:get_real_array(key)
• BY_RAINFLOW_COUNT() aggregation
• Reader class
• debug(...) and trace(...)
• base/sddp/deficit_cost_segment.lua
• --modeloutput, --relative_tolerance, --absolute_tolerance arguments
• PSR.interpolate_colors(color1, color2, color3, n)

Added (Collection)

• OptPriceAgent collection
• OptPriceContract collection
• OptPriceLoad collection
• OptPricePlant collection
• OptPriceSystem collection

Added (Input Data)

• battery.system
• hydro.initial_condition
• hydro.initial_condition_type
• powerinjection.system
• system.risk_aversion_curve

Syntax Changes

• study:get_relationship(...) to study:get_relationship_map(...)
• powerinjection.hour_capacity -> powerinjection.capacity
• powerinjection.hour_price -> powerinjection.price

Changed

• internal changes in reshape_stages

Fixed

• GasEmission collection
• error when loading an csv file inside a subdirectory
• error in concatenate(...)
• error saving an expression that has stage type = quarter

Removed

• writer:open()

Submodules

• github-markdown-css 5.3.0 -> 5.4.0
• PSRClasses 5aa8e761 -> 57b055d9
• zlib v1.2.13 -> v1.3

Changelog

0.26.0 (2023-10-18)

https://psrio.psr-inc.com/0.26.0/index.html

Added

• exp:shift_values(shift)
• exp:aggregate_stages_weighted(by, weights, profile)
• exp:assignments_mask(k)
• safe_divide(exp1, exp2)
• chart:add_line_exclude_zeros(...)
• chart:add_heatmap_series_exclude_zeros(...)
• study:final_year_without_buffer_years()
• study:get_parameter(model, key, fallback)
• toml:get_string_array(...)
• GenericConstraintInterpolation collection
• $/kW, m3 as favorite units
• base/sddp/unitary_variable_cost.lua

Changed

• charts front-end improvements
• psrio.pmk to PSRIO.pmk
• psrio.pmd to PSRIO.pmd
• psrio.ver to PSRIO.ver
• psrio.exe to PSRIO.exe
• psrio to PSRIO (linux)

Fixed

• AWS S3 header reading when there are duplicated agents
• unit when exp:force_unit(...) is the last operation
• optgen study horizon
• binary expression backwards relationship
• base/sddp/reservoir_stored_energy.lua

Submodules

• PSRClasses 2cfa7e4a -> 5aa8e761
• bootstrap 5.3.1 -> 5.3.2
• bootstrap icons 1.11.0 -> 1.11.1
• github-markdown-css 5.2.0 -> 5.3.0
• KaTeX 0.16.8 -> 0.16.9
• highlightjs 11.8.0 -> 11.9.0

Changelog

0.25.0 (2023-08-11)

https://psrio.psr-inc.com/0.25.0/index.html

Added

• PSRIO Marketplace
• PSRIO Extensions
• collections relationship rework with an internal graph search
• collection:file_exists(path)
• PSR.source_line(level)
• expansion_project.capacity
• renewable.system
• chart:add_heatmap_series(...)
• chart:add_cumulative_distribution_function(...)
• chart:push_annotation(x, y, text)
• xLabel and yLabel for the chart options
• #tab length operator
• tab:push_table(...)
• tab:push_divider()
• exp:aggregate_agents_by_label(...)
• exp:reorder_agents(...)
• exp:bus_from() and exp:bus_to()
• exp:select_and_rename_agent(agent, label)
• exp:set(stage, scenario, block, agent, data)
• exp:add_agent_right(label)
• exp:add_agent_left(label)
• study:get_relationship(collection_from, collection_to)
• study.deficit_segment_1, ..., study.deficit_segment_4
• study.deficit_cost_1, ..., study.deficit_cost_4
• study:initial_hydrology_stage() and study:initial_hydrology_year()
• hydro_gauging_station:final_hydrology_year()

Fixed

• histogram chart
• errors handling huge graf files

Removed

• GaugingStation collection (users should use HydroGaugingStation collection)
• exp:save(..., { remove_zeros = true }), the users should use the exp:remove_zeros() method

Submodules

• bootstrap 5.2.3 -> 5.3.1
• fast_float 4.0.0 -> 5.2.0
• GitHub Flavored Markdown 0.29.0.gfm.11 -> 0.29.0.gfm.13
• highlightjs 11.7.0 -> 11.8.0
• KaTeX 0.16.4 -> 0.16.8
• Lua 5.4.4 -> 5.4.6
• PSRClasses 5ce7fe6e -> 2cfa7e4a

Changelog

0.24.0 (2023-05-12)

https://psrio.psr-inc.com/0.24.0/index.html

Added

• improved NCP integration
• chart:add_box_plot(...)
• dashboard:set_refresh_rate(...)
• toml:get_integer_array(key)
• exp:select_agents(collection) when exp is from a generic collection
• exp:uncouple_blocks(clear_first_data)
• exp:uncouple_stages_blocks()
• exp:uncouple_stages_blocks(clear_first_data)
• exp:spread_hours()
• exp:aggregate_broadcast_blocks(f, profile)
• collection length operator: #collection
• --upload_dashboards argument (require AWS_ACCESS_KEY_ID and AWS_SECRET_ACCESS_KEY)
• --incremental_progress_bar argument
• circuit_flow.lua, thermal_operation.lua and thermal_dispatch_factor.lua
• hydro.min_total_outflow_hourly
• hydro.max_total_outflow_hourly
• hydro.min_operative_storage_hourly
• hydro.max_operative_storage_hourly
• hydro.flood_control_hourly
• hydro.alert_storage_hourly
• hydro.min_spillage_hourly
• hydro.max_spillage_hourly
• hydro.min_bio_spillage_hourly
• hydro.target_storage_hourly
• hydro.max_turbining_hourly

Syntax Changes

• toml:get_array_size(...) to toml:get_table_array_size(...)
• toml:get_array(...) to toml:get_table_array(...)
• thermal_generation_running_hours.lua to thermal_plant_running_hours.lua
• concentratedsolarpower.hour_generation to concentratedsolarpower.hour_scenarios

Changed

• if the binary expression relationship can not be found, try the other way around
• collection:load_table and collection:load_table_without_header behavior when the file does not exist; now return an empty table

Fixed

• binary expression relationship issue that occurred when the collection consisted of only one agent
• missing trim in collection:load_table(...)
• chart:add_probability_of_exceedance(...) and chart:add_probability_of_exceedance(...) last point (100%)
• collection:load_table_without_header loading (#133)
• issue that occurred when attempting to save a dashboard or report within a nested folder
• error that occurred while running a script with a report in dependency mode

Removed

• Dashboard("title")
• dashboard + dashboard operator
• dashboard:push(dashboard)
• dashboard:push(chart)
• dashboard:push(markdown)
• area.imported
• area.exported
• circuit.international_cost_right
• circuit.international_cost_left
• dclink.capacity_right
• dclink.capacity_left
• fuelcontract.max_offtake
• interconnection.capacity_right
• interconnection.capacity_left
• interconnection.cost_right
• interconnection.cost_left
• renewable.hour_generation
• renewable.block_generation
• renewable_gauging_station.hour_generation
• renewable_gauging_station.block_generation
• renewable_gauging_station.hour_historical_generation

Submodules

• Bootstrap Icons 1.10.2 -> 1.10.4
• fast_float 3.10.0 -> 4.0.0
• GitHub Flavored Markdown 0.29.0.gfm.9 -> 0.29.0.gfm.11
• PSRClasses 62b06eab -> 5ce7fe6e

Changelog

0.23.0 (2023-04-03)

https://psrio.psr-inc.com/0.23.0/index.html

Added

• dashboard loading overlay
• HydroGaugingStation collection (alias to GaugingStation, to be deprecated in the future)
• collection:avids()
• relational:add(header, vector)
• exp:spread_stages(int)
• exp:replace_scenarios(src, {int, int, ...})
• study:stages_without_buffer_years() and study:buffer_stages()
• hydro.storage_x_elevation__storage_1, ..., hydro.storage_x_elevation__storage_5
• hydro.storage_x_elevation__elevation_1, ..., hydro.storage_x_elevation__elevation_5

Changed

• BY_PERCENTILE(p) to be the same as Excel PERCENTILE.EXC

Fixed

• fuel.emission_factor
• collection:labels()
• --magent intermittent error

Removed

• writer:open(), now the open function is called in the constructor

Submodules

• PSRClasses c72d7788 -> 62b06eab

Changelog

0.22.0 (2023-03-09)

https://psrio.psr-inc.com/0.22.0/index.html

Added

• new AWS S3 saving implementation
• exp:aggregate_scenarios_blocks(f)
• exp:aggregate_topology(f, topology, min_level, max_level)
• chart:add_heatmap(exp)
• chart:add_column_range(exp, exp)
• chart:add_error_bar(exp, exp)
• chart:set_digits(int) and timeseries:set_digits(int)
• chart:add_probability_of_nonexceedance(...)
• Topology.SPILLED_FROM
• --load_only_csv argument

Syntax Changes

• Collection.DCLINK to Collection.DC_LINK

Changed

• improved heatmap rendering performance
• improved the significant digits in the fast_csv
• collection:labels(), collection:codes(), and collection:inputs() now it returns a table

Fixed

• error in collection:get_files(...) and collection:get_directories(...) when the directory does not exist
• error when trying to save data in a subfolder (#93)
• error loading netplan bases
• error when converting units with negative stages (#96)

Removed

• chart:add_heatmap_hourly(...)
• awsintegration

Submodules

• fast_float 3.8.1 -> 3.10.0
• GitHub Flavored Markdown 0.29.0.gfm.7 -> 0.29.0.gfm.9
• Github Markdown CSS 5.1.0 -> 5.2.0
• PSRClasses 1ce193dd -> 66c9973f
• PSRPlot 0.4 -> 0.7

Changelog

0.21.0 (2023-01-16)

https://psrio.psr-inc.com/0.21.0/index.html

Added

• NCP support
• exp:aggregate_blocks() with automatic aggregation function selection
• GenericConstraint collection
• exp:map_scenarios(filename)
• --load_from_output_path reading in collection:load_toml(...)
• dashboard:set_title(string)
• generic:create_writer(filename)
• writer:open()
• writer:is_open()
• writer:write(string)
• writer:write_line(string)
• writer:close()
• --log_append argument

Syntax Changes

• --logname to --log_name argument
• toml:get_bool(...) to toml:get_boolean(...)
• toml:get_int(...) to toml:get_integer(...)
• toml:get_double(...) to toml:get_real(...)

Fixed

• intermittent error in Load collection
• relational:add(string, study, from, to) when the collection is empty
• error in aggregate_stages when first_stage > 1` (#86)
• error using aggregate_stages(..., Profile.PER_WEEK) when exp has weekly stages (#89)

Removed

• some trace messages
• loadfile and dofile functions

Submodules

• PSRClasses 07a23f47 -> 1ce193dd

Changelog

0.20.0 (2022-12-23)

https://psrio.psr-inc.com/0.20.0/index.html

Added

• study:has_flexible_demand()
• study:has_hourly_renewable_scenarios()

Syntax Changes

• all base functions now receive the study index, ex: sddpgrxxd(i, suffix)
• battery.existing to battery.state
• fuel_reservoir.maxinjection to fuel_reservoir.max_injection
• fuel_reservoir.maxinjection_chronological to fuel_reservoir.max_injection_constraint
• hydro.min_total_outflow_penalty to hydro.min_total_outflow_unit_violation_cost
• hydro.max_total_outflow_penalty to hydro.max_total_outflow_unit_violation_cost
• hydro.min_operative_storage_penalty to hydro.min_operative_storage_unit_violation_cost
• hydro.max_operative_storage_penalty to hydro.max_operative_storage_unit_violation_cost
• hydro.alert_storage_penalty to hydro.alert_storage_unit_violation_cost
• hydro.min_spillage_penalty to hydro.min_spillage_unit_violation_cost
• hydro.max_spillage_penalty to hydro.max_spillage_unit_violation_cost
• hydro.min_bio_spillage_penalty to hydro.min_bio_spillage_unit_violation_cost
• hydro.min_turbining_penalty to hydro.min_turbining_unit_violation_cost
• hydro.max_turbining_penalty to hydro.max_turbining_unit_violation_cost
• hydro.min_total_outflow_penalty_type to hydro.min_total_outflow_violation_type
• hydro.max_total_outflow_penalty_type to hydro.max_total_outflow_violation_type
• hydro.min_operative_storage_penalty_type to hydro.min_operative_storage_violation_type
• hydro.max_operative_storage_penalty_type to hydro.max_operative_storage_violation_type
• hydro.alert_storage_penalty_type to hydro.alert_storage_violation_type
• hydro.min_spillage_penalty_type to hydro.min_spillage_violation_type
• hydro.max_spillage_penalty_type to hydro.max_spillage_violation_type
• hydro.min_bio_spillage_penalty_type to hydro.min_bio_spillage_violation_type
• hydro.min_turbining_penalty_type to hydro.min_turbining_violation_type
• hydro.max_turbining_penalty_type to hydro.max_turbining_violation_type
• renewable.hour_generation to renewable.hour_scenarios
• renewable.block_generation to renewable.block_scenarios
• renewable_gauging_station.hour_generation to renewable_gauging_station.hour_scenarios
• renewable_gauging_station.block_generation to renewable_gauging_station.block_scenarios
• renewable_gauging_station.hour_historical_generation to renewable_gauging_station.hour_historical_scenarios

Changed

• disabled tab behavior when there are sub tabs

Fixed

• exp:to_hour(BY_AVERAGE())
• hblock loading when --model none and there is no duraci
• error adding the same tab to a dashboard more the once
• error converting to csv when the binary has negative stages
• intermittent error loading an output on linux
• system.load_level_length floating point approximation

Submodules

• PSRClasses updated
• Highlight.js 11.6.0 -> 11.7.0
• KaTeX 0.16.3 -> 0.16.4

Changelog

0.19.0 (2022-12-15)

https://psrio.psr-inc.com/0.19.0/index.html

Added

• new AWS S3 downloader
• exp:select_scenarios(int, int)
• tab:set_disabled()
• tab:set_collapsed(boolean)
• study:has_hour_block_map()
• study:get_parameter(string, string)
• concatenate_blocks(...)
• info(std::vector<std::string>)
• ConcentratedSolarPower collection (variables: hour_generation)
• ElectrificationDemand collection
• ElectrificationDemandSegment collection (variables: hour, block, cost, hour_price)
• ElectrificationNetwork collection
• ElectrificationNode collection
• ElectrificationProcess collection
• ElectrificationProducer collection
• ElectrificationStorage collection
• ElectrificationTransport collection
• GasEmission collection (variables: code, cost)
• GasNode collection (variables: code, production_cost_constraint)
• Load collection (variables: code, bus, hour_bus)
• ReservoirSet collection (variables: code, security_energy, flood_control_energy, alert_energy)
• dclink.wheeling_cost_from and dclink.wheeling_cost_to
• fuel.min_consumption, fuel.max_consumption and fuel.availability
• system.carbon_credit_cost
• thermal.startup_cost_constraint, thermal.spinning_reserve, thermal.max_reserve, and thermal.bid_price
• balancing_area.regulation_up and balancing_area.regulation_down
• demand.variable_block_duration
• demand_segment.hour_price and demand_segment.block_price
• error check in exp:to_block() and exp:to_hour() when the case has no hour block map

Syntax Changes

• study:is_hourly_load() to study:has_hourly_load()
• circuit.existing to circuit.state
• dclink.existing to dclink.state
• generation_constraint.existing to generation_constraint.state
• generation_constraint.capacity to generation_constraint.data
• hydro.existing to hydro.state
• renewable.existing to renewable.state
• thermal.existing to thermal.state

Fixed

• docs publish (tag)
• errors in concatenate_stages and concatenate_scenarios when first_stage > 1
• add_index_translation
• active tab when the first ones are disabled
• error when creating the output path

Removed

• thermal.forced_generation
• include folder (julia users should use the PSRIO.jl package)

Submodules

• AWS 1.9.379 added
• zlib 1.2.13 added
• AWS Integration updated
• Bootstrap 5.2.2 -> 5.2.3
• Bootstrap Icons 1.9.1 -> 1.10.2
• S3Integrations removed

Changelog

0.18.0 (2022-11-21)

https://psrio.psr-inc.com/0.18.0/index.html

Added

• new documentation (https://psrio.psr-inc.com/)
• BY_ORDER(i)
• collection:inputs()
• study:stage_type(), study:blocks(stage), and study:initial_stage()
• --load_from_output_path reading in collection:load_table(...)
• exp:moving(f, window)
• Profile.PER_QUARTER and Profile.QUARTER
• hydro.alert_storage_penalty and hydro.alert_storage_penalty_type
• hydro.min_total_outflow_penalty and hydro.min_total_outflow_penalty_type
• hydro.max_total_outflow_penalty and hydro.max_total_outflow_penalty_type
• hydro.min_spillage_penalty and hydro.min_spillage_penalty_type
• hydro.max_spillage_penalty and hydro.max_spillage_penalty_type
• hydro.max_turbining_penalty and hydro.max_turbining_penalty_type
• hydro.min_bio_spillage_penalty and hydro.min_bio_spillage_penalty_type
• hydro.min_turbining_penalty and hydro.min_turbining_penalty_type
• hydro.min_operative_storage_penalty and hydro.min_operative_storage_penalty_type
• hydro.max_operative_storage_penalty and hydro.max_operative_storage_penalty_type

Syntax Changes

• study.blocks to study.blocks_per_stage
• hydro.capacity to hydro.max_generation
• hydro.capacity_maintenance to hydro.max_generation_available
• hydro.vmin to hydro.min_storage
• hydro.vmax to hydro.max_storage
• hydro.qmin to hydro.min_turbining_outflow
• hydro.qmax to hydro.max_turbining_outflow
• hydro.vmin_chronological_historical_scenarios_nodata to hydro.min_operative_storage_historical_scenarios_nodata
• hydro.vmin_chronological_historical_scenarios to hydro.min_operative_storage_historical_scenarios
• hydro.vmin_chronological to hydro.min_operative_storage
• hydro.vmax_chronological_historical_scenarios_nodata to hydro.max_operative_storage_historical_scenarios_nodata
• hydro.vmax_chronological_historical_scenarios to hydro.max_operative_storage_historical_scenarios
• hydro.vmax_chronological to hydro.max_operative_storage

Changed

• internal log structure
• internal compiler structure

Fixed

• temporary psrclasses files leftovers error
• relationship error
• --model none parameters check
• exp:to_hour(f, "hblock file") and exp:to_block(f, "hblock file") error when the file does not exist

Submodules

• KaTeX 0.16.2 -> 0.16.3

Changelog

0.17.0 (2022-10-19)

Added

• exp:sort_agents_descending() and exp:sort_agents_ascending()
• exp:uncouple_stages(clear_first_stage)
• ExpansionConstraint and ExpansionDecision collections
• base/sddp/vere15.lua
• hydro.max_turbining, hydro.spinning_reserve, and hydro.max_reserve
• thermal.alternative_fuel
• $ to favorite units
• CI/CD: psr-update-modules

Syntax Changes

• chart:add_heatmap(exp) to chart:add_heatmap_hourly(exp)
• generic:create(values, unit) to generic:create(label, unit, values)
• exp:sort_agents to exp:sort_agents_labels
• demand.maximum_increase to demand.max_increase
• demand.maximum_decrease to demand.max_decrease
• demand.maximum_curtailment to demand.max_curtailment
• hydro.omcost to hydro.om_cost
• hydro.FOR to hydro.forced_outage_rate
• hydro.COR to hydro.historical_outage_factor
• expansion_project.omcost to expansion_project.om_cost
• renewable.omcost to renewable.om_cost
• system.duraci to system.load_level_length
• system.hblock to system.hour_block_map
• thermal.omcost to thermal.om_cost
• thermal.germin to thermal.min_generation
• thermal.germin_maintenance to thermal.min_generation_available
• thermal.capacity to thermal.max_generation
• thermal.capacity_maintenance to thermal.max_generation_available
• thermal.FOR to thermal.forced_outage_rate
• thermal.COR to thermal.historical_outage_factor
• thermal.cesp1 to thermal.specific_consumption_segment_1
• thermal.cesp2 to thermal.specific_consumption_segment_2
• thermal.cesp3 to thermal.specific_consumption_segment_3
• thermal.transport_cost to thermal.fuel_transportation_cost
• thermal.must_run to thermal.operation_mode
• thermal.minimum_uptime to thermal.min_uptime
• thermal.minimum_downtime to thermal.min_downtime
• thermal.maximum_startups to thermal.max_startups
• thermal.maximum_shutdowns to thermal.max_shutdowns

Fixed

• output selection
• aggregate_stages when the first_stage > 1
• missing columns when the series are added individually
• error when adding a chart more than once in the dashboard

Submodules

• Bootstrap 5.2.1 -> 5.2.2

Changelog

0.16.0 (2022-09-21)

Added

• Tab class
• dashboard submenu support
• dashboard:set_icon(url) and dashboard:set_logo(url)
• exp:spread_blocks(blocks)
• exp:select_agents_by_regex(regex) and exp:remove_agents_by_regex(regex)

Fixed

• data literal when using a --model none
• exp:to_block(f, "hblock file")
• exp:spread_stages()

Removed

• Dashboard concatenation operator (!)

Submodules

• Bootstrap 5.2.1
• GitHub Flavored Markdown 0.29.0.gfm.6

Changelog

0.15.0 (2022-09-02)

Added

• improved CI/CD
• PSR.assert_version(version), ex: PSR.assert_version(=0.15.0), PSR.assert_version(>0.15.0)
• collection:codes()
• exp:fill(value)
• BY_KTH_LARGEST(p) and BY_KTH_LARGEST_INDEX(p)
• BY_KTH_SMALLEST(p) and BY_KTH_SMALLEST_INDEX(p)
• BY_MAX_INDEX() and BY_MIN_INDEX()
• exp:reshape_stages(Profile.MONTH)

Syntax Changes

• BY_NTH_ELEMENT(p) to BY_PERCENTILE_INDEX(p)

Changed

• BY_NPV(pu) to BY_NPV(percent), ex: BY_NPV(0.2) to BY_NPV(20)

Fixed

• exp:select_agents({...}) with duplicate agents
• exp:reshape_stages(Profile.DAY) when exp has weekly stages

Submodules

• GitHub Flavored Markdown 0.29.0.gfm.5
• KaTeX 0.16.2
• JSON 3.11.2

Changelog

0.14.3 (2022-08-19)

Added

• CI/CD: notify on slack

Fixed

• CI/CD: publish on master
• hourly scenarios data dimensions check
• exp:aggregate_stages(f, profile) initial stage

0.14.0 (2022-08-17)

Added

• multi-language code syntax highlighting
• CI/CD via github actions
• require(...) access to .lua files next to the recipes
• exp:circuit_from() and exp:circuit_to() to map buses to circuits
• exp:dclink_from() and exp:dclink_to() to map buses to dclinks
• relationship from buses to dclinks in exp:aggregate_agents()
• exp:uncouple_blocks()
• exp:aggregate_slices(f, size)

Fixed

• units check in report
• data reading with "0/1" unit

Changelog

0.13.0 (2022-07-28)

Added

• subhourly support (1 min, 5 min, 10 min, 15 min, 20 min, 30 min)
• TimeSeries and Relational classes
• FlowController collection
• chart:add_line_stacking(exp)
• #chart length
• study:initial_year() and study:final_year()
• collection.code input data
• exp:to_list(stage, scenario, block) and exp:to_int_list(stage, scenario, block)
• exp:remove_scenarios({int, ...})
• exp:spread_stages()
• compare:set_relative_tolerance(double) and compare:set_absolute_tolerance(double)

Syntax Changes

• Compare(title, rtol = 1e-4, atol = 0) to Compare(title)

Changed

• lua 5.4.2 to 5.4.4
• sol 3.2.3 to 3.3.0
• feather-icons replaced by lucide.dev
• collection:load_table(...) and collection:load_table_without_header(...) from now on it considers all data as string
• update to 0.29.0.gfm.4

Fixed

• utf8 encoding
• error in concatenate_scenarios(...)
• ternary operator when there were negative stages

Removed

• psrclasses epoch check
• study:save_relationship(...)
• --return_errors argument

Changelog

0.12.0 (2022-05-29)

Added

• raw HTML parsing in markdown
• chart:add_categories(exp, label) and chart:add_categories(exp1, exp2, label)
• chart:add_line(exp, {color = {"#9b5de5", "#f15bb5", ...}})
• study:save_relationship("filename", from, to)
• Topology.TURBINED_FROM
• concatenate_scenarios(exp1, exp2, ...) and clamp(exp, low, hi)
• exp:repeat_on_buffer_years()
• exp:to_hour(f, "hblock file") and exp:to_block(f, "hblock file")
• exp:force_blocks()
• exp:set_initial_stage(int)
• exp:select_agents(collection, {label1, label2, ...})
• exp:aggregate_stages(f, {stage1, stage2, ...}) and exp:aggregate_scenarios_stages(f, {stage1, stage2, ...})
• exp:cumsum_agents()
• exp:fill_agents({int or string, int or string, ...}, value)
• collection:load_table_without_header("filename")
• psrclasses-netplan.dat loading
• --load_from_output_path argument

Syntax Changes

• PSR.colors({"#ff0029", "#377eb8", ...}) to PSR.set_global_colors({"#ff0029", "#377eb8", ...})

Changed

• improvements in dashboard frontend
• incremental progress bar in psrcloud

Fixed

• error in exp:select_agents(query) when the attributes were different
• error in exp:select_agents(collection, {label}) when there were multiple agents with the same label

Changelog

0.11.0 (2022-04-11)

Added

• dashboard frontend rework
• new unit matching algorithm
• input data dimensions check rework
• --label name, --label code, and --label avid arguments
• --ignore_hrbmap argument
• exp:save(filename, {label="name"}), exp:save(filename, {label="code"}), and exp:save(filename, {label="avid"})
• exp:first_stage(), exp:last_stage(), and exp:code(index)
• Chart(title, subtitle)
• chart:add_spline(exp)
• chart:add_area_spline(exp), chart:add_area_spline_stacking(exp), chart:add_area_spline_percent(exp), and chart:add_area_spline_range(exp1, exp2)
• exp:remove_agent(index), exp:remove_agent(label), and exp:rename_agent(label)
• exp:select_agent_by_code(code) and exp:select_agents_by_code({code1, code2, ...})
• exp:remove_agent_by_code(code) and exp:remove_agents_by_code({code1, code2, ...})
• exp:rename_agents_with_codes()
• exp:select_agents(collection, label)
• exp:spread_scenarios(scenarios)
• collection:language()
• study:get_parameter(key, double)
• PSR.colors({"#ff0029", "#377eb8", ...}) and PSR.interpolate_colors("#ff0000", "#00ff00", 4)
• toml:get_array_size(key) and toml:get_array(key, i)
• sfd, R$ units

Syntax Changes

• exp:select_stages(stage) to exp:select_stage(stage)
• exp:select_stages(first_stage, -1) to exp:select_first_stage(first_stage)
• exp:select_stages(-1, last_stage) to exp:select_last_stage(last_stage)
• study:macro_agents(...) to collection:load_tag(...)
• study:load_table(filename) to collection:load_table(filename)
• study:load_toml(filename) to collection:load_toml(filename)
• study:get_parameter(key) to study:get_parameter(key, int)
• Compare(title, rtol = 1e-4) to Compare(title, rtol = 1e-4, atol = 0)

Fixed

• fast_csv warning when the csv was locked
• saving in a directory that does not exist

Removed

• --avid argument
• chart:add_line_stacking(exp) and chart:add_line_percent(exp)

Changelog

0.10.1 (2022-02-14)

Fixed

• intermittent error in dependencies mode

0.10.0 (2022-02-07)

Added

• improved error messages
• linear domain axis chart (stage_type = 0)
• chart:add_heatmap(exp) and chart:save(filename, title)
• dashboard:set_icon(string) based on feathericons
• collection:cloudname()
• collection:get_directories(regex) and collection:get_directories(subpath, regex)
• exp:select_agent(index) and exp:select_agent(label)
• exp:set_stage_type(int)
• exp:save(filename, {fast_csv = true}) and exp:save(filename, {bin = true})
• exp:add_prefix(string)
• exp:rename_agent(string)
• exp:has_agent(label)
• exp:stages_to_agents() and exp:stages_to_agents(selected)
• PSR.sleep(seconds) PSR.version()
• info(table)
• all aggregation functions in exp:to_block(function)

Syntax Changes

• study:load_table(filename) no longer assume the extension is a .csv
• exp:rename_agents_with_suffix(string) to exp:add_suffix(string)

Fixed

• exp:select_scenarios({int, ...}) when the first stage was different than 1
• `chart:add_area_range(exp, exp) with data considering leap years

Removed

• exp:save(filename, {tmp = true})

Changelog

0.9.0 (2022-01-03)

Added

• guards in the loading method to protect from psrclasses unknown error
• report:add_header(string)
• generic:create({value1, value2, ...}) and generic:create({value1, value2, ...}, "unit")
• study:load_table(filename) and collection:get_files(regex)
• PSR.studies()
• negative epoch
• relationship thermal/hydro/renewable plant to network areas
• CO, tCO units

Syntax Changes

• add_index_translation to PSR.add_index_translation

Changed

• improvements in security
• improvements in report classes
• improvements in units classes
• improvements in makefile and psrcloud deploy
• argument --recipes can be a directory, and the psrio will search for all *.lua files

Fixed

• CVAR aggregation when the values were equal
• error in dependencies mode when saving a not selected output
• error loading an output with upper case filename in linux
• memory constraint on loading large hourly scenarios
• error in BY_SUM(), BY_MIN(), and BY_MAX() when there were no values to aggregate
• rollback to lua 5.4.2 (https://github.com/ThePhD/sol2/issues/1266)

Changelog

0.8.0 (2021-11-18)

Added

• LaTeX math support in markdown ($$ ... $$)
• exp:select_agents(std::vector<std::string>)
• exp:replace(exp)
• int = study:get_parameter(string)
• exp:select_largest_agents(k) and exp:select_smallest_agents(k)
• study:openings()
• exp:sort_agents()
• Compare(title, epsilon = 1e-4)
• 100 max errors in Compare
• sddpconfig.dat reading
• exp:force_collection(collection)
• parm reference relationship
• --return_errors argument

Syntax Changes

• study.stages to study:stages()
• study.stages_per_year to study:stages_per_year()
• study.scenarios to study:scenarios()
• study:parent_path() to study:dirname()
• collection:load(filename, true) to collection:load_or_error(filename)

Changed

• performance improvement in concatenate(exp1, exp2, ...)
• performance improvement in exp:select_agents(query)

Fixed

• stderr when the case does not exists
• unit in power binary expression (^)
• study.discount_rate unit
• error using remove_zeros and horizon

Changelog

0.7.0 (2021-09-24)

Added

• new dashboard structure with markdown support
• multiple cases support
• output comparison with report files
• lua helpers, e.g. require("lua/spairs")
• exp:rename_agents_with_suffix(string)
• exp:select_block(block)
• collection:load(filename, true) to raise error if the file is not found
• exp:force_unit(unit)
• exp = exp:save_cache()
• exp:aggregate_scenarios_stages(f)
• exp:remove_zeros()
• exp:aggregate_agents(f, group_size)
• exp:uncouple_stages()
• study:path() and study:parent_path()
• exp:set_initial_year(year)
• do_string(code) (meta programming)
• shift_stages(stages)
• topology enumerate and exp:aggregate_topology(f, topology)
• info(msg), warning(msg), error(msg)
• forward and backward reading

Changed

• exp:agents() return a vector of lua objects
• unit conversion in logical binary expressions

Fixed

• exp:aggregate_stages() in monthly-hourly cases
• missing blocks in binary expression
• eps in logical binary operations (lq, le, gt, ge)
• exp:to_hour() and exp:to_block() when the input has no blocks
• stderr not showing in the log in release mode
• optgen study loading

Changelog

0.6.0 (2021-06-04)

Added

• interactive mode (REPL)
• psrio-base and require("<base>")
• include("<file>") with relative directory
• runtime input data loading: collection:load_vector(attribute, unit) and collection:load_parameter(attribute, unit)
• exp:select_stages(stage), exp:select_stages_by_year(initial_year, final_year), exp:reset_stages()
• exp:select_scenarios(scenario) and exp:select_scenarios({scenario}) [per stage]
• exp:select_agents(query)
• exp:round(digits)
• exp:force_hourly()
• exp:save("<file>", { horizon = true })
• attributes getters:
  • exp:loaded(), exp:unit(), exp:agents(), exp:agents_size(),exp:agent(int i), exp:scenarios()
  • exp:stage_type(), exp:initial_year(), exp:initial_stage(), exp:stages(), exp:week(stage), exp:month(stage), exp:year(stage)
  • exp:blocks(int), exp:has_blocks(), exp:is_hourly(), exp:to_list(), exp:to_int_list()
• study:is_hourly() and study:is_genesys() methods
• pu unit
• add_index_translation(...) and macro_agents() functions
• single binary support
• relationship from buses to circuits in exp:aggregate_agents()

Changed

• added files in indice.grf even if it is not in indexdat.fmt
• deleted bin/hdr files in CSV mode

Fixed

• several fixes of methods with first_stage > 1
• error handling when loading an output with the wrong collection
• error in concatenate when expression agents = 0
• input hourly data loading when there is no hour-block mapping file
• exp:to_block() and historical scenarios performance

Removed

• variable_by_scenario metadata (variable_by_scenario = false is now scenarios = 1)
• --skip argument, now it is automatic

Changelog

0.5.0 (2021-01-27)

Added

• Lua
• dashboard support
• --horizon <date> to change to initial date
• --model none to load an empty study

Syntax Changes

• several breaking changes due to Lua introduction

Removed

• bison scanner and parser
• verbose 1f, 2f, and 3f (use 1, 2, or 3)

Changelog

0.4.1 (2021-01-12)

Added

• first and last stage attributes to all variables and operators
• select_stages(exp, first_stage, last_stage)
• rename_agents(exp, [string, string, ...])
• argument --logname <string> to rename the log files (psrio.log and psrio-psrclasses.log to <string>.log and <string>-psrclasses.log)

Changed

• crop_stages(exp) is now select_stages(exp)

Fixed

• indexdat.fmt concurrent access

0.4.0 (2020-12-01)

Changed

• agents selection by order or name: select_agents(exp, {int or string, int or string, ...})
• scenarios selection parameter in aggregation: aggregate_scenarios(type, exp, {int, int, ...})

Removed

• argument --append

Examples

Thermal Generation per Fuel

thermal = Thermal();
gerter = thermal:load("gerter");
gerter:aggregate_agents(BY_SUM(), Collection.FUEL):save("gerter_per_fuel");

Select System and Technology

renewable = Renewable();
renewable_generation = renewable:load("gergnd");

solar = renewable.tech_type:eq(2);
nordeste = renewable.system:eq(4);

renewable_generation:select_agents(solar & nordeste):save("renewable_generation_solar_nordeste");

Renewable Generation per Technology

renewable = Renewable();
renewable_generation = renewable:load("gergnd");

general = renewable.tech_type:eq(0);
wind = renewable.tech_type:eq(1);
solar = renewable.tech_type:eq(2);
biomass = renewable.tech_type:eq(3);
small_hydro = renewable.tech_type:eq(4);
concentrated_solar_power = renewable.tech_type:eq(5);

concatenate(
    renewable_generation:select_agents(general):aggregate_agents(BY_SUM(), "General"),
    renewable_generation:select_agents(wind):aggregate_agents(BY_SUM(), "Wind"),
    renewable_generation:select_agents(solar):aggregate_agents(BY_SUM(), "Solar"),
    renewable_generation:select_agents(biomass):aggregate_agents(BY_SUM(), "Biomass"),
    renewable_generation:select_agents(small_hydro):aggregate_agents(BY_SUM(), "Small Hydro"),
    renewable_generation:select_agents(concentrated_solar_power):aggregate_agents(BY_SUM(), "CSP")
):save("renewable_generation_per_technology");

Select All December Stages

generic = Generic();
output = generic:load("cmgdem");
first_stage = output:first_stage();
last_stage = output:last_stage();

selected_stages = {};
for stage = first_stage, last_stage do
    if output:month(stage) == 12 then
        table.insert(selected_stages, stage);
    end
end

output:stages_to_agents(selected_stages):save("example");

Custom Chart 1

local system = System();
local cmgdem = system:load("cmgdem"):aggregate_blocks(BY_AVERAGE());
local p10 = cmgdem:aggregate_scenarios(BY_PERCENTILE(10)):aggregate_agents(BY_SUM(), "p10");
local p25 = cmgdem:aggregate_scenarios(BY_PERCENTILE(25)):aggregate_agents(BY_SUM(), "p25");
local p50 = cmgdem:aggregate_scenarios(BY_PERCENTILE(50)):aggregate_agents(BY_SUM(), "p50");
local p75 = cmgdem:aggregate_scenarios(BY_PERCENTILE(75)):aggregate_agents(BY_SUM(), "p75");
local p90 = cmgdem:aggregate_scenarios(BY_PERCENTILE(90)):aggregate_agents(BY_SUM(), "p90");
local chart = Chart("Custom Chart 1");

chart:add_area_range(p10, p90, { color = "blue", lineWidth = 0, fillOpacity = 0.2 });
chart:add_area_range(p25, p75, { color = "blue", lineWidth = 0, fillOpacity = 0.2 });
chart:add_line(p50, { color = "red", lineWidth = 2, dashStyle = "dash" });

Custom Chart 2

local system = System();
local cmgdem = system:load("cmgdem"):aggregate_blocks(BY_AVERAGE());
local p10 = cmgdem:aggregate_scenarios(BY_PERCENTILE(10)):aggregate_agents(BY_SUM(), "p10");
local p90 = cmgdem:aggregate_scenarios(BY_PERCENTILE(90)):aggregate_agents(BY_SUM(), "p90");
local avg = cmgdem:aggregate_scenarios(BY_AVERAGE()):aggregate_agents(BY_SUM(), "Average");
local chart = Chart("Custom Chart 2");

local options = {
    xGridLineWidth = 1
};
chart:add_error_bar(p10, p90, options);

local options = {
    color = "green",
    lineWidth = 0,
    marker = { enabled = true, symbol = "rect" }
};
chart:add_line(p10, options);

local options = {
    color = "red",
    lineWidth = 0,
    marker = { enabled = true, symbol = "triangle" }
};
chart:add_line(p90, options);

local options = {
    color = "blue",
    lineWidth = 0,
    marker = { enabled = true, symbol = "circle" },
    dataLabels = { enabled = true, format = "{point.y:,.0f}" }
};
chart:add_line(avg, options);

Custom Chart 3

local system = System();
local cmgdem = system:load("cmgdem"):aggregate_blocks(BY_AVERAGE());
local p10 = cmgdem:aggregate_scenarios(BY_PERCENTILE(10)):aggregate_agents(BY_SUM(), "p10");
local p50 = cmgdem:aggregate_scenarios(BY_PERCENTILE(50)):aggregate_agents(BY_SUM(), "p50");
local p90 = cmgdem:aggregate_scenarios(BY_PERCENTILE(90)):aggregate_agents(BY_SUM(), "p90");
local chart = Chart("Custom Chart 3");

chart:add_error_bar(p10, p90, { color = "black" });
chart:add_line(p50, { color = "grey", marker = { enabled = true, symbol = "circle" } });

Collection Syntax

A collection refers to a set of characteristics associated with a specific group of data. It stores information such as agent names, stage resolutions, and the number of scenarios within the group. For instance, the collection Hydro() contains information about hydro plants in the study.

Collection	Constructor	Enumerate
Area	Area()	Collection.AREA
Balancing Area	BalancingArea()	Collection.BALANCING_AREA
Balancing Area Hydro	BalancingAreaHydro()	Collection.BALANCING_AREA_HYDRO
Balancing Area Thermal	BalancingAreaThermal()	Collection.BALANCING_AREA_THERMAL
Battery	Battery()	Collection.BATTERY
Bus	Bus()	Collection.BUS
Circuit	Circuit()	Collection.CIRCUIT
Circuits Sum	CircuitsSum()	Collection.CIRCUITS_SUM
Concentrated Solar Power	ConcentratedSolarPower()	Collection.CONCENTRATED_SOLAR_POWER
DC Link	DCLink()	Collection.DC_LINK
Demand	Demand()	Collection.DEMAND
Demand Segment	DemandSegment()	Collection.DEMAND_SEGMENT
Expansion Capacity	ExpansionCapacity()	Collection.EXPANSION_CAPACITY
Expansion Project	ExpansionProject()	Collection.EXPANSION_PROJECT
Fuel	Fuel()	Collection.FUEL
Fuel Consumption	FuelConsumption()	Collection.FUEL_CONSUMPTION
Fuel Contract	FuelContract()	Collection.FUEL_CONTRACT
Fuel Reservoir	FuelReservoir()	Collection.FUEL_RESERVOIR
Generator	Generator()	Collection.GENERATOR
Generation Constraint	GenerationConstraint()	Collection.GENERATION_CONSTRAINT
Generic	Generic()	---
Hydro	Hydro()	Collection.HYDRO
Hydro Gauging Station	HydroGaugingStation()	Collection.HYDRO_GAUGING_STATION
Interconnection	Interconnection()	Collection.INTERCONNECTION
Interconnection Sum	InterconnectionSum()	Collection.INTERCONNECTION_SUM
Power Injection	PowerInjection()	Collection.POWER_INJECTION
Renewable	Renewable()	Collection.RENEWABLE
Renewable Gauging Station	RenewableGaugingStation()	Collection.RENEWABLE_GAUGING_STATION
Reserve Generation Constraint	ReserveGenerationConstraint()	Collection.RESERVE_GENERATION_CONSTRAINT
Reservoir Set	ReservoirSet()	Collection.RESERVOIR_SET
Study	Study()	---
System	System()	Collection.SYSTEM
Thermal	Thermal()	Collection.THERMAL

Loading an Output

PSRIO can load output data from the PSR models, particularly those in the graph format that can be loaded by PSRIO using the load method. Please note that data should be loaded using the specific collection associated with them.

Operator	Syntax
Load method	output = Collection():load("filename")

The generic collection can load any output with any agent type.

Example 1

The following example loads two outputs, gerhid and fprodt, considering the agents as hydro plants collection:

hydro = Hydro();
gerhid = hydro:load("gerhid");
fprodt = hydro:load("fprodt");

Example 2

The following example loads two outputs, cmgdem and demand, considering the agents as system collection:

system = System();
cmgdem = system:load("cmgdem");
demand = system:load("demand");

Example 3

The following example loads two outputs, gerter and coster, considering the agents as thermal plants collection:

thermal = Thermal();
gerter = thermal:load("gerter");
coster = thermal:load("coster");

Example 4

The following example loads two outputs, objcop and outdfact, considering the agents as generic:

generic = Generic();
objcop = generic:load("objcop");
outdfact = generic:load("outdfact");

Loading a Generic CSV ( with header )

Operator	Syntax
Load table method	table = Study():load_table("filename.csv")

The following example loads a generic csv table and iterates through the rows:

Example

file.csv:

Name         , Node, Interval, Startup, Participation
H_1_1        , 1   , 1       , 1      , 0.79
R_1_3        , 1   , 1       , 1      , 0.20
CONTRATO_FLAT, 1   , 1       , 1      , 1

script:

study = Study();
table = study:load_table("file.csv");

info("Name,Node,Interval,Startup,Participation");
for i=1,#table do
    info(
        table[i]["Name"] .. "," .. 
        table[i]["Node"] .. "," .. 
        table[i]["Interval"] .. "," .. 
        table[i]["Startup"] .. "," .. 
        table[i]["Participation"]
    );
end

Loading a Generic CSV ( without header )

Operator	Syntax
Load table method	table = Study():load_table_without_header("filename.csv")

The following example loads a generic csv table and iterates through the rows:

Example 1

file.csv:

Name         , Node, Interval, Startup, Participation
H_1_1        , 1   , 1       , 1      , 0.79
R_1_3        , 1   , 1       , 1      , 0.20
CONTRATO_FLAT, 1   , 1       , 1      , 1

script:

study = Study();
table = study:load_table("file.csv");

for i=1,#table do
    info(
        table[i][1] .. "," .. 
        table[i][2] .. "," .. 
        table[i][3] .. "," .. 
        table[i][4] .. "," .. 
        table[i][5]
    );
end

Loading an Input

• Area
• Battery
• Bus
• Circuit
• Concentrated Solar Power
• DC Link
• Demand
• Demand Segment
• Fuel
• Hydro
• Hydro (Tables)
• Hydro Gauging Station
• Interconnection
• Power Injection
• Renewable
• Renewable Gauging Station
• Reservoir Set
• Study
• System
• Thermal

PSRIO is capable of loading a majority of input data from the PSR model. Input data should be loaded using the specific collection associated with each data type. The ´.´ operator must be utilized to perform this loading, as illustrated below.

Area

Data	Unit
Area().code	---
Area().import_limit	MW
Area().export_limit	MW

Battery

Data	Unit
Battery().code	---
Battery().state	---
Battery().capacity	MW

Bus

Data	Unit
Bus().code	---
Bus().voltage_level	kV

Circuit

Data	Unit
Circuit().code	---
Circuit().state	---
Circuit().resistance	%
Circuit().reactance	%
Circuit().capacity	MW
Circuit().emergency_capacity	MW
Circuit().DLR_factor	pu
Circuit().monitored	---
Circuit().monitored_contingencies	---
Circuit().is_dc	---
Circuit().international_cost_from	$/MWh
Circuit().international_cost_to	$/MWh

Concentrated Solar Power

Data	Unit
ConcentratedSolarPower().hour_scenarios	pu

DC Link

Data	Unit
DCLink().code	---
DCLink().state	---
DCLink().capacity_from	MW
DCLink().capacity_to	MW
DCLink().wheeling_cost_from	$/MWh
DCLink().wheeling_cost_to	$/MWh

Demand

Data	Unit
Demand().code	---
Demand().is_elastic	---
Demand().inelastic_hour	MW
Demand().inelastic_block	GWh
Demand().is_flexible	---
Demand().max_increase	pu
Demand().max_decrease	pu
Demand().curtailment_cost	$/MWh
Demand().max_curtailment	pu
Demand().variable_block_duration	h

Demand Segment

Data	Unit
DemandSegment().hour	MW
DemandSegment().block	GWh
DemandSegment().hour_price	$/MWh
DemandSegment().block_price	$/MWh

Fuel

Data	Unit
Fuel().code	---
Fuel().cost	$/gal
Fuel().emission_factor	tCO2/MWh
Fuel().min_consumption	kgal
Fuel().max_consumption	kgal
Fuel().availability	kgal

Hydro

Data	Unit
Hydro().code	---
Hydro().state	---
Hydro().units	---
Hydro().system	---
Hydro().max_generation	MW
Hydro().max_generation_available	MW
Hydro().forced_outage_rate	%
Hydro().historical_outage_factor	%
Hydro().min_storage	hm3
Hydro().max_storage	hm3
Hydro().min_turbining_outflow	m3/s
Hydro().max_turbining_outflow	m3/s
Hydro().om_cost	$/MWh
Hydro().irrigation	m3/s
Hydro().min_total_outflow_modification	m3/s
Hydro().target_storage_tolerance	%
Hydro().disconsider_in_stored_and_inflow_energy	---
Hydro().mean_production_coefficient	MW/(m3/s)
Hydro().loss_factor	pu
Hydro().min_total_outflow	m3/s
Hydro().min_total_outflow_unit_violation_cost	k$/m3/s
Hydro().min_total_outflow_violation_type	---
Hydro().max_total_outflow	m3/s
Hydro().max_total_outflow_unit_violation_cost	k$/m3/s
Hydro().max_total_outflow_violation_type	---
Hydro().min_operative_storage	hm3
Hydro().min_operative_storage_unit_violation_cost	k$/hm3
Hydro().min_operative_storage_violation_type	---
Hydro().max_operative_storage	hm3
Hydro().max_operative_storage_unit_violation_cost	k$/hm3
Hydro().max_operative_storage_violation_type	---
Hydro().flood_control	hm3
Hydro().alert_storage	hm3
Hydro().alert_storage_unit_violation_cost	k$/hm3
Hydro().alert_storage_violation_type	---
Hydro().min_spillage	m3/s
Hydro().min_spillage_unit_violation_cost	k$/m3/s
Hydro().min_spillage_violation_type	---
Hydro().max_spillage	m3/s
Hydro().max_spillage_unit_violation_cost	k$/m3/s
Hydro().max_spillage_violation_type	---
Hydro().min_bio_spillage	%
Hydro().min_bio_spillage_unit_violation_cost	k$/hm3
Hydro().min_bio_spillage_violation_type	---
Hydro().target_storage	hm3
Hydro().max_turbining	m3/s
Hydro().max_turbining_unit_violation_cost	k$/m3/s
Hydro().max_turbining_violation_type	---
Hydro().min_turbining_unit_violation_cost	k$/m3/s
Hydro().min_turbining_violation_type	---
Hydro().spinning_reserve	%
Hydro().max_reserve	MW

Hydro (tables)

Data	Unit
Hydro().storage_x_elevation__storage_1	hm3
Hydro().storage_x_elevation__storage_2	hm3
Hydro().storage_x_elevation__storage_3	hm3
Hydro().storage_x_elevation__storage_4	hm3
Hydro().storage_x_elevation__storage_5	hm3
Hydro().storage_x_elevation__elevation_1	masl
Hydro().storage_x_elevation__elevation_2	masl
Hydro().storage_x_elevation__elevation_3	masl
Hydro().storage_x_elevation__elevation_4	masl
Hydro().storage_x_elevation__elevation_5	masl

Hydro Gauging Station

Data	Unit
HydroGaugingStation().code	---
HydroGaugingStation().inflow	m3/s
HydroGaugingStation().forward	m3/s
HydroGaugingStation().backward	m3/s
HydroGaugingStation().hour_inflow_historical_scenarios_nodata	---
HydroGaugingStation().hour_inflow_historical_scenarios	m3/s
HydroGaugingStation().hour_inflow	m3/s

Interconnection

Data	Unit
Interconnection().code	---
Interconnection().state	---
Interconnection().capacity_from	MW
Interconnection().capacity_to	MW
Interconnection().cost_from	$/MWh
Interconnection().cost_to	$/MWh

Power Injection

Data	Unit
PowerInjection().code	---
PowerInjection().hour_capacity	MW
PowerInjection().hour_price	$/MWh

Renewable

Data	Unit
Renewable().code	---
Renewable().state	---
Renewable().units	---
Renewable().tech_type	---
Renewable().capacity	MW
Renewable().om_cost	$/MWh
Renewable().hour_scenarios	pu
Renewable().block_scenarios	pu
Renewable().operation_factor	pu

Renewable Gauging Station

Data	Unit
RenewableGaugingStation().code	---
RenewableGaugingStation().hour_scenarios	pu
RenewableGaugingStation().block_scenarios	pu
RenewableGaugingStation().hour_historical_scenarios	pu

Reservoir Set

Data	Unit
ReservoirSet().code	---
ReservoirSet().security_energy	MWh
ReservoirSet().flood_control_energy	MWh
ReservoirSet().alert_energy	MWh

Study

Data	Unit
Study().deficit_segment_1	%
Study().deficit_segment_2	%
Study().deficit_segment_3	%
Study().deficit_segment_4	%
Study().deficit_cost_1	$/MWh
Study().deficit_cost_2	$/MWh
Study().deficit_cost_3	$/MWh
Study().deficit_cost_4	$/MWh

System

Data	Unit
System().code	---
System().load_level_length	---
System().hour_block_map	---
System().sensitivity	---
System().carbon_credit_cost	$/tCO2

Thermal

Data	Unit
Thermal().code	---
Thermal().state	---
Thermal().units	---
Thermal().system	---
Thermal().min_generation	MW
Thermal().min_generation_available	MW
Thermal().min_generation_constraint	MW
Thermal().max_generation	MW
Thermal().max_generation_available	MW
Thermal().forced_outage_rate	%
Thermal().historical_outage_factor	%
Thermal().startup_cost	k$
Thermal().startup_cost_constraint	k$
Thermal().om_cost	$/MWh
Thermal().specific_consumption_segment_1	gal/MWh
Thermal().specific_consumption_segment_2	gal/MWh
Thermal().specific_consumption_segment_3	gal/MWh
Thermal().fuel_transportation_cost	$/gal
Thermal().operation_mode	---
Thermal().emission_coefficient	pu
Thermal().ramp_up	MW/min
Thermal().ramp_down	MW/min
Thermal().min_uptime	hour
Thermal().min_downtime	hour
Thermal().max_startups	---
Thermal().max_shutdowns	---
Thermal().shutdown_cost	k$
Thermal().alternative_fuel	---
Thermal().spinning_reserve	%
Thermal().max_reserve	MW
Thermal().bid_price	$/MWh

Study Attributes

An important collection that will be available in all cases is the Study() collection. This collection contains essential information about the case, such as the number of stages and scenarios, the horizon, and the resolution of blocks/hours. Below are the methods related to the Study() collection.

Stages

Method	Return Type
Study():stage_type()	number
Study():initial_stage()	number
Study():initial_year()	number
Study():final_year()	number
Study():stages()	number
Study():stages_per_year()	number
Study():stages_without_buffer_years()	number

Blocks

Method	Return Type
Study():blocks(stage)	number
Study():is_hourly()	boolean
Study():has_hourly_load()	boolean
Study():has_hour_block_map()	boolean

Scenarios

Method	Return Type
Study():scenarios()	number
Study():openings()	number

Others

Method	Return Type
Study():get_parameter(key, fallback_integer)	integer
Study():get_parameter(key, fallback_double)	double

Directories

Method	Return Type
Study():dirname()	string
Study():cloudname()	string

Expressions Attributes

Attributes are properties that characterize an expression. For example, expressions related to the results of a PSR model, whose data can change over stages, scenarios, and blocks, have a set of attributes in PSRIO that can be useful in some tasks. In the example below, we load the thermal generation:

thermal = Thermal();
gerter = thermal:load("gerter");

And the log shows some expression's attributes information:

[info] Loading gerter (stages: 12 [1:12] [month] [10/2016], blocks: hour, scenarios: 1200, unit: GWh, agents: 3 [thermal])

From the log, it is possible to see that gerter has:

• 12 stages in a month-level resolution, where the first stage is 1 and the last is 12
• Initial stage is 10 (october) and the initial year is 2016
• Hourly discretization
• 1200 scenarios
• Generation data given in GWh
• 3 agents (thermal plants)

Expression Attributes

All attributes related to stages, block, scenarios, agents and units of an expression and their respective methods are presented in the tables below.

Stages

Method	Return Type
exp:first_stage()	number
exp:last_stage()	number
exp:stage_type()	number
exp:initial_stage()	number
exp:initial_year()	number
exp:final_year()	number
exp:week(stage)	number
exp:month(stage)	number
exp:year(stage)	number
exp:stage(index)	number
exp:stages()	number

Example

stages = gerter:stages() -- 12

stage1 = gerter:stage(1) -- 10
stage2 = gerter:stage(2) -- 11
stage3 = gerter:stage(3) -- 12
stage4 = gerter:stage(4) -- 13

stage_type = gerter:stage_type() -- 2 (monthly)

initial_stage = gerter:initial_stage() -- 10
initial_year = gerter:initial_year() -- 2016
final_year = gerter:final_year() -- 2017

month1 = gerter:month(1) -- 10
month2 = gerter:month(2) -- 11
month3 = gerter:month(3) -- 12
month4 = gerter:month(4) --  1

year1 = gerter:year(1) -- 2016
year4 = gerter:year(4) -- 2017

Blocks

Method	Return Type
exp:blocks(stage)	number
exp:has_blocks()	boolean
exp:is_hourly()	boolean
exp:hour_discretization()	number

Example

blocks1 = gerter:blocks(1) -- 744
blocks2 = gerter:blocks(2) -- 720

has_blocks = gerter:has_blocks() -- true

is_hourly = gerter:is_hourly() -- true

Scenarios

Method	Return Type
exp:scenarios()	number

Example

scenarios = gerter:scenarios() -- 1200

Agents

Method	Return Type
exp:agents_size()	number
exp:agent(index)	string
exp:code(index)	number
exp:agents()	table of strings
exp:codes()	table of numbers
exp:has_agent(string)	boolean

Example

agents = gerter:agents() -- { "Thermal 1", "Thermal 2", "Thermal 3" }
for i, agent in ipairs(agents) do
    info(agent .. " at index " .. i);
end

size = gerter:agents_size() -- 3
for i = 1,size do
    agent = gerter:agent(i);
    info(agent .. " at index " .. i);
end

[info] Thermal 1 at index 1
[info] Thermal 2 at index 2
[info] Thermal 3 at index 3

Unit

Method	Return Type
exp:unit()	string

Example

unit = gerter:unit() -- "GWh"

Unary Expressions

Minus

$$\operatorname{exp}=-\operatorname{exp1}$$

Absolute Value

$$\operatorname{exp}=\operatorname{exp1:abs}()$$

Example

circuit = Circuit();
cirflw = circuit:load("cirflw");
abs_cirflw = cirflw:abs();

Round

$$\operatorname{exp}=\operatorname{exp1:round}(\operatorname{int})$$

Fill

$$\operatorname{exp}=\operatorname{exp1:fill}(\operatorname{double})$$

Unit Conversion

$$\operatorname{exp}=\operatorname{exp1:convert}(\operatorname{string})$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
gerterMW = gerter:convert("MW");

Force Unit

$$\operatorname{force\_unit}=\operatorname{exp1:convert}(\operatorname{string})$$

Example

hydro = Hydro();
final_cota = hydro:load("cotfin);
final_cota_in_changed_unit = final_cota:force_unit("msmn");

Binary Expressions

Addition

$$\operatorname{exp}=\operatorname{exp1}+\operatorname{exp2}$$

Example

hydro = Hydro();
spilled_outflow = hydro:load("qverti");
turbined_outflow = hydro:load("qturbi");
total_outflow = spilled_outflow + turbined_outflow;

Subtraction

$$\operatorname{exp}=\operatorname{exp1}-\operatorname{exp2}$$

Example

hydro = Hydro();
useful_storage = hydro.max_storage - hydro.min_storage;

Multiplication

$$\operatorname{exp}=\operatorname{exp1}*\operatorname{exp2}$$

Example

renewable = Renewable();
renewable_generation = renewable:load("gergnd");
renewable_om_cost = renewable_generation * renewable.om_cost;

Division

$$\operatorname{exp}=\operatorname{exp1}/\operatorname{exp2}$$

Example

circuit = Circuit();
circuit_flow = circuit:load("cirflw");
circuit_loading = circuit_flow / circuit.capacity;

Safe Division

$$\operatorname{exp}=\operatorname{safe\_divide}(\operatorname{exp1}, \operatorname{exp2})$$

Example

renewable = Renewable();

generation = generic:load("gergnd");
spillage = generic:load("vergnd");
spillage_proportion = safe_divide(spillage / (spillage + generation));

Power

$$\operatorname{exp}=\operatorname{exp1}\char`^\operatorname{exp2}$$

Equal to

$$\operatorname{exp}=\operatorname{exp1:eq}(\operatorname{exp2})$$

Example

circuit = Circuit();
circuit_flow = circuit:load("cirflw");
circuit_in_max_charge = circuit_flow:abs():eq(circuit.capacity);

Not Equal to

$$\operatorname{exp}=\operatorname{exp1:ne}(\operatorname{exp2})$$

Example

circuit = Circuit();
circuit_flow = circuit:load("cirflw");
circuit_not_in_max_charge = circuit_flow:abs():ne(circuit.capacity);

Less-than

$$\operatorname{exp}=\operatorname{exp1:lt}(\operatorname{exp2})$$

Example

circuit = Circuit();
circuit_flow = circuit:load("cirflw");
circuit_less_than_the_max_charge = circuit_flow:abs():lt(circuit.capacity);

Less-than-or-equals to

$$\operatorname{exp}=\operatorname{exp1:le}(\operatorname{exp2})$$

Example

circuit = Circuit();
circuit_flow = circuit:load("cirflw");
circuit_less_or_equal_the_max_charge = circuit_flow:abs():le(circuit.capacity);

Greater-than

$$\operatorname{exp}=\operatorname{exp1:gt}(\operatorname{exp2})$$

Example

circuit = Circuit();
circuit_flow = circuit:load("cirflw");
circuit_greater_than_the_max_charge = circuit_flow:abs():gt(circuit.capacity);

Greater-than-or-equals to

$$\operatorname{exp}=\operatorname{exp1:ge}(\operatorname{exp2})$$

Example

circuit = Circuit();
circuit_flow = circuit:load("cirflw");
circuit_greater_than_or_equal_the_max_charge = circuit_flow:abs():ge(circuit.capacity);

And

$$\operatorname{exp}=\operatorname{exp1}\&\operatorname{exp2}$$

Or

$$\operatorname{exp}=\operatorname{exp1}|\operatorname{exp2}$$

Maximum

$$\operatorname{exp}=\operatorname{max}(\operatorname{exp1}, \operatorname{exp2})$$

Example

hydro = Hydro();
turbined_outflow = hydro:load("qturbi");
minimum_turbined_outflow = max(hydro.min_turbining_outflow - turbined_outflow, 0);

Minimum

$$\operatorname{exp}=\operatorname{min}(\operatorname{exp1}, \operatorname{exp2})$$

Example

hydro = Hydro();
production_factor = hydro:load("fprodt");
available_hydro_capacity = min(hydro.max_turbining_outflow * production_factor, hydro.max_generation_available);

Ternary Expressions

Conditional

$$\operatorname{exp}=\operatorname{ifelse}(\operatorname{exp1}, \operatorname{exp2}, \operatorname{exp3})$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
gerter_gt_zero = ifelse(gerter:gt(0.0), 1, 0);

Aggregate Functions

Syntax
BY_SUM()
BY_MULTIPLICATION()
BY_AVERAGE()
BY_CVAR_L(number)
BY_CVAR_R(number)
BY_STDDEV()
BY_REPEATING()
BY_FIRST_VALUE()
BY_ORDER(number)
BY_LAST_VALUE()
BY_MAX()
BY_MAX_INDEX()
BY_MIN()
BY_MIN_INDEX()
BY_KTH_LARGEST(number)
BY_KTH_LARGEST_INDEX(number)
BY_KTH_SMALLEST(number)
BY_KTH_SMALLEST_INDEX(number)
BY_PERCENTILE(number)
BY_PERCENTILE_INDEX(number)

Scenarios

This section will cover how to apply the methods found in the aggregate_functions tab to data scenarios.

Aggregate Scenarios

$$\operatorname{exp}=\operatorname{exp1:aggregate\_scenarios}(\operatorname{f})$$

Example

system = System();
cmgdem = system:load("cmgdem");
cmgdem_avg = cmgdem:aggregate_scenarios(BY_AVERAGE());
cmgdem_p90 = cmgdem:aggregate_scenarios(BY_PERCENTILE(90));

Aggregate Selected Scenarios

$$\operatorname{exp}=\operatorname{exp1:aggregate\_scenarios}(\operatorname{f}, \{\operatorname{int}, \operatorname{int}, ...\})$$

Example

system = System();
cmgdem = system:load("cmgdem");
cmgdem_max = cmgdem:aggregate_scenarios(BY_MAX(), {1, 2, 3, 4, 5});

Select One Scenario

$$\operatorname{exp}=\operatorname{exp1:select\_scenario}(\operatorname{int})$$

Example

system = System();
cmgdem = system:load("cmgdem");
cmgdem_scenario32 = cmgdem:select_scenario(32);

Select Multiple Scenarios

$$\operatorname{exp}=\operatorname{exp1:select\_scenarios}(\{\operatorname{int}, \operatorname{int}, ...\})$$

Select Scenarios Range

$$\operatorname{exp}=\operatorname{exp1:select\_scenarios}(\operatorname{int}, \operatorname{int})$$

Remove Multiple Scenarios

$$\operatorname{exp}=\operatorname{exp1:remove\_scenarios}(\{\operatorname{int}, \operatorname{int}, ...\})$$

Concatenate Scenarios

$$\operatorname{exp}=\operatorname{concatenate\_scenarios}(\{\operatorname{exp1}, \operatorname{exp2}, ...\})$$

Example

hydro = Hydro();

gerhid = hydro:load("gerhid");
gerhid_scenario_5 = gerhid:select_scenario(5);
gerhid_scenario_15 = gerhid:select_scenario(15);
gerhid_scenario_25 = gerhid:select_scenario(25);
gerhid_scenarios = concatenate_scenarios(
    gerhid_scenario_5,
    gerhid_scenario_15,
    gerhid_scenario_25
);

Blocks/Hours

This section will cover how to apply the methods found in the aggregate_functions tab to data bolcks\hours.

Aggregate Blocks/Hours

$$\operatorname{exp}=\operatorname{exp1:aggregate\_blocks}(\operatorname{f})$$

Example

system = System();
cmgdem = system:load("cmgdem");
cmgdem_agg = cmgdem:aggregate_blocks(BY_AVERAGE());

renewable = Renewable();
gergnd = renewable:load("gergnd");
gergnd_agg = gergnd:aggregate_blocks(BY_SUM());

Select One Block/Hour

$$\operatorname{exp}=\operatorname{exp1:select\_block}(\operatorname{int})$$

Example

system = System();
cmgdem = system:load("cmgdem");
cmgdem_block21 = cmgdem:select_block(21);

Map Blocks into Hours

$$\operatorname{exp}=\operatorname{exp1:to\_hour}(\operatorname{f})$$

Where f is BY_AVERAGE() or BY_REPEATING().

Example

system = System();
cmgdem_block = system:load("cmgdem");
cmgdem_hourly = cmgdem_block:to_hour(BY_REPEATING());

Map Hours into Blocks

$$\operatorname{exp}=\operatorname{exp1:to\_block}(\operatorname{f})$$

Where f is BY_AVERAGE() or BY_SUM().

Example

thermal = Thermal();
gerter_hourly = thermal:load("gerter");
gerter_block = gerter_hourly:to_block(BY_SUM());

Stages

This section will cover how to apply the methods found in the aggregate_functions tab to data stages.

Aggregate Stages

$$\operatorname{exp}=\operatorname{exp1:aggregate\_stages}(\operatorname{f})$$

Aggregate Stages into a Profile

$$\operatorname{exp}=\operatorname{exp1:aggregate\_stages}(\operatorname{f}, \operatorname{profile})$$

Where profile is the following enumerate:

Profiles
Profile.STAGE
Profile.WEEK
Profile.MONTH
Profile.QUARTER
Profile.YEAR
Profile.PER_WEEK
Profile.PER_MONTH
Profile.PER_QUARTER
Profile.PER_YEAR

Profile.STAGE

The Profile.STAGE is the default value to characterize the aggregation if the user does not inform any profile. The data associated with each stage of the study horizon is aggregated.

exp1	exp (Profile.STAGE)
n (daily)	1 (daily)
n (weekly)	1 (weekly)
n (monthly)	1 (monthly)
n (yearly)	1 (yearly)

Profile.WEEK and Profile.PER_WEEK

When the data has a daily resolution and the aggregation profile is Profile.WEEK, PSRIO will aggregate the data for each day of the weeks in the study period, i.e., that data regarding all Mondays in the data set will be aggregated into one value and the same for Tuesday, Wednesday, and so on. With a daily resolution and the aggregation Profile.PER_WEEK, PSRIO will aggregate the data related to each week of the study.

When the data is weekly and we request is Profile.WEEK, the data associated with each week is aggregated in one week. If the request is the aggregation Profile.PER_WEEK, PSRIO will do nothing and the data will remains the same.

These aggregation profiles are not defined for monthly and yearly resolution data.

exp1	exp (Profile.WEEK)	exp (Profile.PER_WEEK)
n (daily)	7 (daily)	n/7 (weekly)
n (weekly)	1 (weekly)	n (weekly)
n (monthly)	❌	❌
n (yearly)	❌	❌

Example

system = System();
cmgdem = system:load("cmgdem");
cmgdem_agg = cmgdem:aggregate_stages(BY_AVERAGE(), Profile.WEEK);

Example

system = System();
cmgdem = system:load("cmgdem");
cmgdem_agg = cmgdem:aggregate_stages(BY_AVERAGE(), Profile.PER_WEEK);

Profile.MONTH and Profile.PER_MONTH

Similar to Profile.WEEK and Profile.PER_WEEK, when the data has daily discretization and we request Profile.MONTH, PSRIO will aggregate the data for each day of the months in the study period. For example, PSRIO will aggregate all 1^st days of each month and the same for the others months. If we request the Profile.PER_MONTH, PSRIO will aggregate the data related to each month.

When the data has month-level discretization and we request Profile.MONTH, PSRIO will aggregate the data associated with each month. If we request Profile.PER_MONTH, nothing is done to the data.

These aggregation profiles are not defined for weekly and yearly resolution data.

exp1	exp (Profile.MONTH)	exp (Profile.PER_MONTH)
n (daily)	31 (daily)	n/~30 (monthly)
n (weekly)	❌	❌
n (monthly)	1 (monthly)	n (monthly)
n (yearly)	❌	❌

Example

system = System();
cmgdem = system:load("cmgdem");
cmgdem_agg = cmgdem:aggregate_stages(BY_AVERAGE(), Profile.MONTH);

Example

system = System();
cmgdem = system:load("cmgdem");
cmgdem_agg = cmgdem:aggregate_stages(BY_AVERAGE(), Profile.PER_MONTH);

Profile.YEAR and Profile.PER_YEAR

When the data has a daily resolution and we request a Profile.YEAR, the data related to each day of years is aggregated. For example, PSRIO will aggregate all January 1^st days in the study period, which is done for the other days of the year. If we request a Profile.PER_YEAR, the data associated with each year is aggregated. The same happens when the data has week, month, or year-level resolution if Profile.PER_YEAR is selected.

If the data has a year resolution and profile.YEAR is selected. the data related to the same year is aggregated.

The profile.YEAR profile is not defined for weekly and monthly resolution data.

exp1	exp (Profile.YEAR)	exp (Profile.PER_YEAR)
n (daily)	365 (daily)	n/365 (yearly)
n (weekly)	52 (weekly)	n/52 (yearly)
n (monthly)	12 (monthly)	n/12 (yearly)
n (yearly)	1 (yearly)	1 (yearly)

Example

system = System();
defcit = system:load("defcit");
defcit_per_year = defcit:aggregate_stages(BY_SUM(), Profile.YEAR);

Example

system = System()
defcit = system:load("defcit")
defcit_per_year = defcit:aggregate_stages(BY_SUM(), Profile.PER_YEAR);

Select One Stage

$$\operatorname{exp}=\operatorname{exp1:select\_stage}(\operatorname{int})$$

generic = Generic();
objcop = generic:load("objcop");
objcop_1st_stage = objcop:select_stage(1);

Select the First Stage

$$\operatorname{exp}=\operatorname{exp1:select\_first\_stage}(\operatorname{int})$$

generic = Generic();
objcop = generic:load("objcop");
objcop_2nd_to_n_stage = objcop:select_first_stage(2);

Select the Last Stage

$$\operatorname{exp}=\operatorname{exp1:select\_last\_stage}(\operatorname{int})$$

generic = Generic();
objcop = generic:load("objcop");
objcop_1st_to_2nd_stage = objcop:select_last_stage(2);

Select the First and the Last Stages

$$\operatorname{exp}=\operatorname{exp1:select\_stages}(\operatorname{int}, \operatorname{int})$$

generic = Generic();
objcop = generic:load("objcop");
objcop_2nd_to_4th_stage = objcop:select_stages(2, 4);

Select Stages by a First and a Last Years

$$\operatorname{exp}=\operatorname{exp1:select\_stages\_by\_year}(\operatorname{int}, \operatorname{int})$$

generic = Generic();
objcop = generic:load("objcop");
objcop_2030_2040 = objcop:select_stages_by_year(2030, 2040);

Select Stages by a Year

$$\operatorname{exp}=\operatorname{exp1:select\_stages\_by\_year}(\operatorname{int})$$

generic = Generic();
objcop = generic:load("objcop");
objcop_2035 = objcop:select_stages_by_year(2035);

Reshape Stages

$$\operatorname{exp}=\operatorname{exp1:reshape\_stages}(\operatorname{Profile.DAILY})$$

For an hourly represented data, the stage resolution can changed from a week, month or year level to a daily one using that method

exp1	exp
n (daily-hourly)	n (daily-hourly)
n (weekly-hourly)	7n (daily-hourly)
n (monthly-hourly)	~30n (daily-hourly)
n (yearly-hourly)	365n (daily-hourly)

Example

thermal = Thermal();
gerter = thermal:load("gerter");
gerter_daily = gerter:reshape_stages(Profile.DAILY);

Reset Stages

$$\operatorname{exp}=\operatorname{exp1:reset\_stages}()$$

The reset method sets the initial stage to 1.

Example

thermal = Thermal();
gerter = thermal:load("gerter");
gerter_cut = gerter:select_stages(10,30);
gerter_reset = gerter_cut:reset_stages();

Concatenate Stages

$$\operatorname{exp}=\operatorname{concatenate\_stages}(\{\operatorname{exp1}, \operatorname{exp2}, ...\})$$

Example

hydro = Hydro();

gerhid = hydro:load("gerhid");
gerhid_stage_5 = gerhid:select_stage(5);
gerhid_stage_15 = gerhid:select_stage(15);
gerhid_stage_25 = gerhid:select_stage(25);
gerhid_stages = concatenate_stages(
    gerhid_stage_5,
    gerhid_stage_15,
    gerhid_stage_25
);

Set Initial Stage

$$\operatorname{exp} = \operatorname{set\_initial\_stage}(\operatorname{n})$$

Where n is the desired intial stage.

Example

hydro = Hydro();

gerhid = hydro:load("gerhid");
gerhid_initial_stage = gerhid:set_initial_stage(5);

Set Initial Year

$$\operatorname{exp} = \operatorname{exp1:set\_initial\_year}(\operatorname{n})$$

Where n is the desired intial year.

Example

hydro = Hydro();

gerhid = hydro:load("gerhid");
gerhid_ini_yea = gerhid:set_initial_stage(2021);

Uncouple Stages

$$\operatorname{exp}=\operatorname{exp1:uncouple\_stages}()$$

Example

hydro = Hydro();
gerhid = hydro:load("gerhid");
gerhid_unc = gerhid:uncouple_stages();

Uncouple Stages and Blocks

$$\operatorname{exp}=\operatorname{exp1:uncouple\_stages\_blocks}()$$

Example

hydro = Hydro();
gerhid = hydro:load("gerhid");
gerhid_unc = gerhid:uncouple_stages_blocks();

Agents

This section will cover how to apply the methods found in the aggregate_functions tab to data agents.

Aggregate All Agents

$$\operatorname{exp}=\operatorname{exp1:aggregate\_agents}(\operatorname{f}, \operatorname{string})$$

Example

hydro = Hydro();
gerhid = hydro:load("gerhid");
gerhid_sum = gerhid:aggregate_agents(BY_SUM(), "Total Hydro");

Aggregate Agents into Collection

$$\operatorname{exp}=\operatorname{exp1:aggregate\_agents}(\operatorname{f}, \operatorname{collection})$$

Where collection is Collection Enumerate.

Example

hydro = Hydro();
gerhid = hydro:load("gerhid");
gerhid_systems = gerhid:aggregate_agents(BY_SUM(), Collection.SYSTEM);
gerhid_buses = gerhid:aggregate_agents(BY_SUM(), Collection.BUSES);

Select One Agent by Name or Index

$$\operatorname{exp}=\operatorname{exp1:select\_agent}(\text{string or int})$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
gerter_t1 = gerter:select_agent("Thermal 1");
gerter_t2 = gerter:select_agent(2);

Select Mulitple Agents by Names or Indices

$$\operatorname{exp}=\operatorname{exp1:select\_agents}(\{\text{string or int}, \text{string or int}, ...\})$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
gerter_t1_and_t2 = gerter:select_agents({"Thermal 1", 2});

Select Agents within a Collection

$$\operatorname{exp}=\operatorname{exp1:select\_agents}(\operatorname{collection})$$

Where collection is Collection Enumerate.

Example

expansion_project = ExpansionProject()
outidec = expansion_project:load("outidec");
outidec_dclinks = outidec:select_agents(Collection.DCLINK);

Select Agents within a Collection Element

$$\operatorname{exp}=\operatorname{exp1:select\_agents}(\operatorname{collection},\operatorname{string})$$

Where collection is Collection Enumerate and string the name of a element that belongs to Collection.

Example

expansion_project = ExpansionProject()
outidec = expansion_project:load("outidec");
outidec_from_S1_system = outidec:select_agents(Collection.SYSTEM, "S1");

Select Agents with a Query

$$\operatorname{exp}=\operatorname{exp1:select\_agents}(\operatorname{string})$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
non_zero_gerter = gerter:select_agents(gerter:ne(0));

Select Agents by Regex

$$\operatorname{exp}=\operatorname{exp1:select\_agents\_by\_regex}(\operatorname{string})$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
UFV_agents = gerter:select_agents_by_regex("(UFV_)(.*)");

Select Agent by Code

$$\operatorname{exp}=\operatorname{exp1:select\_agent\_by\_code}(\operatorname{string})$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
agent = gerter:select_agent_by_code("1");

Select Agents by Code

$$\operatorname{exp}=\operatorname{exp1:select\_agents\_by\_code}(\{\operatorname{string}\})$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
agents = gerter:select_agent_by_code({"1","2"});

Remove One Agent by Name or Index

$$\operatorname{exp}=\operatorname{exp1:remove\_agent}(\text{string or int})$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
gerter_t2_and_t3 = gerter:remove_agent("Thermal 1");
gerter_t1_and_t2 = gerter:remove_agent(3);

Remove Mulitple Agents by Names or Indices

$$\operatorname{exp}=\operatorname{exp1:remove\_agents}(\{\text{string or int}, \text{string or int}, ...\})$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
gerter_t1 = gerter:remove_agents({"Thermal 2", 3});

Rename One Agent

$$\operatorname{exp}=\operatorname{exp1:rename\_agent}(\text{string})$$

Rename Mulitple Agents with One Name

$$\operatorname{exp}=\operatorname{exp1:rename\_agents}(\text{string})$$

Rename Mulitple Agents with Mutiple Names

$$\operatorname{exp}=\operatorname{exp1:rename\_agents}(\{\text{string}, \text{string}, ...\})$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
gerter_renamed = gerter:rename_agents({"T1", "T2", "T3"});

Rename Mulitple Agents by Adding a Suffix

$$\operatorname{exp}=\operatorname{exp1:add\_suffix}(\text{string})$$

Rename Mulitple Agents by Adding a Prefix

$$\operatorname{exp}=\operatorname{exp1:add\_prefix}(\text{string})$$

Concatenate Agents

$$\operatorname{exp}=\operatorname{concatenate}(\{\operatorname{exp1}, \operatorname{exp2}, ...\})$$

Example 1

hydro = Hydro();
thermal = Thermal();
renewable = Renewable();

gerhid = hydro:load("gerhid");
gerter = thermal:load("gerter");
gergnd = renewable:load("gergnd");
generation = concatenate(gerhid, gerter, gergnd);

Aggregate Topology

$$\operatorname{exp}=\operatorname{exp1:aggregate\_topology}(\operatorname{f}, \operatorname{topology})$$

Where topology is the following enumerate:

Topologies
Topology.CONTROLLED_BY
Topology.TURBINED_TO
Topology.TURBINED_FROM
Topology.SPILLED_TO
Topology.SPILLED_FROM
Topology.FILTRATION_TO
Topology.ASSOCIATED_RESERVOIR
Topology.NEUTRAL
Topology.STORED_ENERGY_TO

Example

hydro = Hydro();
production_factor = hydro:load("fprodt");
production_factor_accumulated = production_factor:aggregate_topology(BY_SUM(), Topology.TURBINED_TO);

Aggregate Topology with Min and Max Levels

$$\operatorname{exp}=\operatorname{exp1:aggregate\_topology}(\operatorname{f}, \operatorname{topology}, \operatorname{min}, \operatorname{max})$$

Example

hydro = Hydro();
spillage = hydro:load("qverti");
spillage_parents = spillage:aggregate_topology(BY_SUM(), Topology.TURBINED_FROM, 1, 1);

Replace

$$\operatorname{exp}=\operatorname{exp1:replace}(\operatorname{exp2})$$

The agents data from exp1 will be replace by exp2 data with the same agents name.

Example

thermal = Thermal();
potter_agent_1 = thermal:load("potter"):select_agent(1);
gerter = thermal:load("gerter");
gerter_replaced = gerter:replace(potter_agent_1);

Select Smallest Agents

$$\operatorname{exp}=\operatorname{exp1:select\_smallest\_agents}(\operatorname{n})$$

Where n is the number of agents that will be select.

Example

thermal = Thermal();
gerter = thermal:load("gerter");
gerter_smallest = gerter:select_smallest_agents(5);

Select Largest Agents

$$\operatorname{exp}=\operatorname{exp1:select\_largest\_agents}(\operatorname{n})$$

Where n is the number of agents that will be select.

Example

thermal = Thermal();
gerter = thermal:load("gerter");
gerter_largest = gerter:select_largest_agents(5);

Cumulative Sum Agents

$$\operatorname{exp}=\operatorname{exp1:cumsum\_agents}()$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
gerter_sum = gerter:cumsum_agents();

Remove Zeros

$$\operatorname{exp}=\operatorname{exp1:remove\_zeros}()$$

Remove agents with all data equal to zero.

Example

thermal = Thermal();
gerter = thermal:load("gerter");
gerter_without_zeros = gerter:remove_zeros();

Others Dimensions

This section will introduce other types of aggregation that can be applied to stages data.

Moving Function

$$\operatorname{exp}=\operatorname{exp1:moving}(\operatorname{f}, \operatorname{int})$$

This operator can be used to calculate the moving average, as demonstrated in the following example.

Example

system = System();
cmgdem = system:load("cmgdem");
cmgdem_12months = cmgdem:moving(BY_AVERAGE(), 12);

Saving

Files in PSRIO can be saved in various formats, including:

-Binary Format: This format stores two files, .bin and .hdr, representing the binary data optimized for computer processing.
-Single Binary Format: This format is similar to binary format; however, it stores the data in only one file, .dat.
-Comma-Separated Values Format (CSV): This format stores the data in a .csv file.

Save Options

Description	Syntax	Default
Save output as BIN/HDR	bin	true
Save output as CSV	csv	false
Save output as DAT	dat	false

Save

$$\operatorname{exp1:save}(\operatorname{string})$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
gerter:save("example");

Save with Options

$$\operatorname{exp1:save}(\operatorname{string}, \{\operatorname{options}\})$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
gerter:save("example", { csv = true });

Save and Load into an Expression

$$\operatorname{exp}=\operatorname{exp1:save\_and\_load}(\operatorname{string})$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
example = gerter:save_and_load("example");

Save with Options and Load into an Expression

$$\operatorname{exp}=\operatorname{exp1:save\_and\_load}(\operatorname{string}, \{\operatorname{options}\})$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
example = gerter:save_and_load("example", { csv = true });

Save Cache

$$\operatorname{exp}=\operatorname{exp1:save\_cache}()$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
cache = gerter:save_cache();

Add Index Translation

$$\operatorname{PSR.add\_index\_translation}(\operatorname{id}, \operatorname{name}, \operatorname{description})$$

Example

thermal = Thermal();
gerter = thermal:load("gerter");
gerter:aggregate_agents(BY_SUM(), Collection.BUS):save("gerter_per_bus");

PSR.add_index_translation("gerter_per_bus", "Thermal generation per bus", "Thermal plants generation per bus")

Dashboard

PSRIO provides a set of methods that allow users to analyze their results simply and easily.

Create a Tab

$$\operatorname{tab}=\operatorname{Tab}(\text{string})$$

Example

tab = Tab("Tab title");

Set the Tab Icon

$$\operatorname{tab:set\_icon}(\text{string})$$

https://lucide.dev/

Example

tab = Tab("Tab title");
tab:set_icon("home");

Disable the Tab

$$\operatorname{tab:set\_disabled}()$$

Example

tab = Tab("Tab title");
tab:set_disabled();

Set the Collapse Flag

$$\operatorname{tab:set\_collapsed}(\text{boolean})$$

Example

tab = Tab("Tab title");
tab:set_collapsed(true);

Push a Chart to a Tab

$$\operatorname{tab:push}(\text{chart})$$

Example

tab = Tab("Tab title");
chart = Chart("Chart title");
tab:push(chart);

Push a Markdown to a Tab

$$\operatorname{tab:push}(\text{string})$$

Example

tab = Tab("Tab title");
tab:push("# Heading 1");

Push a Tab to a Tab

$$\operatorname{tab:push}(\text{tab})$$

Example

tab = Tab("Tab title");
subtab = Tab("Sub Tab title");
tab:push(subtab);

Create a Dashboard

$$\operatorname{dashboard}=\operatorname{Dashboard}()$$

Example

dashboard = Dashboard();

Push a Tab to a Dashboard

$$\operatorname{dashboard:push}(\operatorname{tab})$$

Example

tab = Tab("Tab title");
dashboard = Dashboard();
dashboard:push(tab);

Save a Dashboard

$$\operatorname{dashboard:save}(\text{string})$$

Example

dashboard = Dashboard();
dashboard:save("example");

Charts

Charts are the essential elements of dashboards. They provide a visual way of analyzing the data generated by PSR models. To create a chart object, we should use one of the following methods:

Method	Syntax
Create a Chart	chart = Chart()
Create a Chart with a title	chart = Chart("title")
Create a Chart with a title and a subtitle	chart = Chart("title", "subtitle")

After creating the chart object, we can start to push the data in it. PSRIO gives multiple options for chart types, which can be seen below:

Method	Syntax
Line	chart:add_line(exp)
Line Categories	chart:add_line_categories(exp)
Spline	chart:add_spline(exp)
Column	chart:add_column(exp)
Column Categories	chart:add_column_categories(exp, "label")
Column Stacking	chart:add_column_stacking(exp)
Column Stacking Categories	chart:add_column_stacking_categories(exp, "label")
Column Percent	chart:add_column_percent(exp)
Column Percent Categories	chart:add_column_percent_categories(exp)
Column Range	chart:add_column_range(exp1, exp2)
Column Range Categories	chart:add_column_range_categories(exp1, exp2)
Area	chart:add_area(exp)
Area Stacking	chart:add_area_stacking(exp)
Area Percent	chart:add_area_percent(exp)
Area Range	chart:add_area_range(exp1, exp2)
Area Spline	chart:add_area_spline(exp)
Area Spline Stacking	chart:add_area_spline_stacking(exp)
Area Spline Percent	chart:add_area_spline_percent(exp)
Area Spline Range	chart:add_area_spline_range(exp1, exp2)
Error Bar	chart:add_error_bar(exp1, exp2)
Pie	chart:add_pie(exp)
Histogram	chart:add_histogram(exp)
Heatmap	chart:add_heatmap(exp)
Heatmap Series	chart:add_heatmap_series(exp)
Scatter	chart:add_scatter(exp1, exp2, "label")
Probability of Exceedance	chart:add_probability_of_exceedance(exp)
Probability of Non Exceedance	chart:add_probability_of_nonexceedance(exp)
Cumulative Distribution Function	chart:add_cumulative_distribution_function(exp)
Box Plot	chart:add_box_plot(exp1, exp2, exp3, exp4, exp5)

Examples

local thermal = Thermal();
local gerter = thermal:load("gerter");

local hydro = Hydro();
local gerhid = hydro:load("gerhid");

local renewable = Renewable();
local gergnd = renewable:load("gergnd");

local tab = Tab("Tutorial");

-- push charts to tab here --

local dashboard = Dashboard();
dashboard:push(tab);
dashboard:save("tutorial");

Example: Line

local chart = Chart("Line");
chart:add_line(gerter:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):select_largest_agents(5));
tab:push(chart);

Example: Spline

local chart = Chart("Spline");
chart:add_spline(gerter:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):select_largest_agents(5));
tab:push(chart);

Example: Column

local chart = Chart("Column");
chart:add_column(gerter:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):select_largest_agents(5));
tab:push(chart);

Example: Column Categories

local chart = Chart("Column Categories");
chart:add_column_categories(gerter:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):aggregate_agents(BY_SUM(), Collection.SYSTEM), "Thermal");
chart:add_column_categories(gerhid:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):aggregate_agents(BY_SUM(), Collection.SYSTEM), "Hydro");
chart:add_column_categories(gergnd:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):aggregate_agents(BY_SUM(), Collection.SYSTEM), "Renewable");
tab:push(chart);

Example: Column Stacking

local chart = Chart("Column Stacking");
chart:add_column_stacking(gerter:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):select_largest_agents(5));
tab:push(chart);

Example: Column Stacking Categories

local chart = Chart("Column Stacking Categories");
chart:add_column_stacking_categories(gerter:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):aggregate_agents(BY_SUM(), Collection.SYSTEM), "Thermal");
chart:add_column_stacking_categories(gerhid:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):aggregate_agents(BY_SUM(), Collection.SYSTEM), "Hydro");
chart:add_column_stacking_categories(gergnd:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):aggregate_agents(BY_SUM(), Collection.SYSTEM), "Renewable");
tab:push(chart);

Example: Column Percent

local chart = Chart("Column Percent");
chart:add_column_percent(gerter:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):select_largest_agents(5));
tab:push(chart);

Example: Column Range

local chart = Chart("Column Range");
chart:add_column_range(
    gerter:aggregate_blocks(BY_SUM()):aggregate_agents(BY_SUM(), "p10"):aggregate_scenarios(BY_PERCENTILE(10)),
    gerter:aggregate_blocks(BY_SUM()):aggregate_agents(BY_SUM(), "p90"):aggregate_scenarios(BY_PERCENTILE(90))
);
tab:push(chart);

Example: Area

local chart = Chart("Area");
chart:add_area(gerter:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):select_largest_agents(5));
tab:push(chart);

Example: Area Stacking

local chart = Chart("Area Stacking");
chart:add_area_stacking(gerter:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):select_largest_agents(5));
tab:push(chart);

Example: Area Percent

local chart = Chart("Area Percent");
chart:add_area_percent(gerter:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):select_largest_agents(5));
tab:push(chart);

Example: Area Range

local chart = Chart("Area Range");
chart:add_area_range(
    gerter:aggregate_blocks(BY_SUM()):aggregate_agents(BY_SUM(), "p10"):aggregate_scenarios(BY_PERCENTILE(10)),
    gerter:aggregate_blocks(BY_SUM()):aggregate_agents(BY_SUM(), "p90"):aggregate_scenarios(BY_PERCENTILE(90))
);
tab:push(chart);

Example: Area Spline

local chart = Chart("Area Spline");
chart:add_area_spline(gerter:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):select_largest_agents(5));
tab:push(chart);

Example: Area Spline Stacking

local chart = Chart("Area Spline Stacking");
chart:add_area_spline_stacking(gerter:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):select_largest_agents(5));
tab:push(chart);

Example: Area Spline Percent

local chart = Chart("Area Spline Percent");
chart:add_area_spline_percent(gerter:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):select_largest_agents(5));
tab:push(chart);

Example: Area Spline Range

local chart = Chart("Area Spline Range");
chart:add_area_spline_range(
    gerter:aggregate_blocks(BY_SUM()):aggregate_agents(BY_SUM(), "p10"):aggregate_scenarios(BY_PERCENTILE(10)),
    gerter:aggregate_blocks(BY_SUM()):aggregate_agents(BY_SUM(), "p90"):aggregate_scenarios(BY_PERCENTILE(90))
);
tab:push(chart);

Example: Error Bar

local chart = Chart("Error Bar");
chart:add_error_bar(
    gerter:aggregate_blocks(BY_SUM()):aggregate_agents(BY_SUM(), "p10"):aggregate_scenarios(BY_PERCENTILE(10)),
    gerter:aggregate_blocks(BY_SUM()):aggregate_agents(BY_SUM(), "p90"):aggregate_scenarios(BY_PERCENTILE(90))
);
tab:push(chart);

Example: Pie

local chart = Chart("Pie");
chart:add_pie(gerter:aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):aggregate_stages(BY_SUM()));
tab:push(chart);

Example: Histogram

local chart = Chart("Histogram");
chart:add_histogram(gerter:aggregate_blocks(BY_SUM()):aggregate_agents(BY_SUM(), "Total"));
tab:push(chart);

Example: Heatmap (Hourly)

local chart = Chart("Heatmap (Hourly)");
chart:add_heatmap(gerter:aggregate_scenarios(BY_AVERAGE()):aggregate_agents(BY_SUM(), "total"));
tab:push(chart);

Example: Heatmap (Daily)

local chart = Chart("Heatmap (Daily)");
chart:add_heatmap(gerter:reshape_stages(Profile.DAILY):aggregate_blocks(BY_SUM()):aggregate_scenarios(BY_AVERAGE()):aggregate_agents(BY_SUM(), "total"));
tab:push(chart);