Corrected indexing

Author: sambuddhac

example_systems/1_three_zones_asymmetric_bidirectional_Actual_Symmetric 2/README.md

@@ -0,0 +1,32 @@
+# Three Zones
+
+**Three Zones**, a one-year example with hourly resolution, contains zones representing Massachusetts, Connecticut, and Maine. The ten represented resources include natural gas, solar PV, wind, and lithium-ion battery storage.
+
+To run the model, first navigate to the example directory:
+
+- Using a Julia REPL:
+
+```bash
+$ julia
+julia> cd("example_systems/1_three_zones/")
+```
+
+- Using a terminal or command prompt:
+```bash
+$ cd example_systems/1_three_zones/
+``` 
+   
+Next, ensure that your settings in `settings/genx_settings.yml` are correct. The default settings use the solver `HiGHS`, time domain reduced input data (`TimeDomainReduction: 1`) and minimum capacity requirement policy (`MinCapReq: 1`) as specified in the `policies/Minimum_capacity_requirement.csv` file. Other optional policies include a capacity reserve margin, an energy share requirement (such as renewable portfolio standard (RPS) or clean electricity standard (CES) policies), a CO2 emissions cap, and a maximum capacity requirement policy (see the documentation for more details). For this example, a rate-based carbon cap of 50 gCO<sub>2</sub> per kWh is specified in the `policies/CO2_cap.csv` input file.
+
+Once the settings are confirmed, run the model with the `Run.jl` script in the example directory:
+
+- Using a Julia REPL (recommended)
+```julia
+julia> include("Run.jl")
+```
+- Using a terminal or command prompt:
+```bash
+$ julia Run.jl
+```
+
+Once the model has completed, results will write to the `results` directory.

example_systems/1_three_zones_asymmetric_bidirectional_Actual_Symmetric 2/Run.jl

@@ -0,0 +1,3 @@
+using GenX
+
+run_genx_case!(dirname(@__FILE__))

example_systems/1_three_zones_asymmetric_bidirectional_Actual_Symmetric 2/policies/CO2_cap.csv

@@ -0,0 +1,4 @@
+,Network_zones,CO_2_Cap_Zone_1,CO_2_Cap_Zone_2,CO_2_Cap_Zone_3,CO_2_Max_tons_MWh_1,CO_2_Max_tons_MWh_2,CO_2_Max_tons_MWh_3,CO_2_Max_Mtons_1,CO_2_Max_Mtons_2,CO_2_Max_Mtons_3
+MA,z1,1,0,0,0.05,0,0,0.018,0,0
+CT,z2,0,1,0,0,0.05,0,0,0.025,0
+ME,z3,0,0,1,0,0,0.05,0,0,0.025

example_systems/1_three_zones_asymmetric_bidirectional_Actual_Symmetric 2/policies/Minimum_capacity_requirement.csv

@@ -0,0 +1,4 @@
+MinCapReqConstraint,ConstraintDescription,Min_MW
+1,MA_PV,5000
+2,CT_Wind,10000
+3,All_Batteries,6000

example_systems/1_three_zones_asymmetric_bidirectional_Actual_Symmetric 2/resources/Storage.csv

@@ -0,0 +1,4 @@
+Resource,Zone,Model,New_Build,Can_Retire,Existing_Cap_MW,Existing_Cap_MWh,Max_Cap_MW,Max_Cap_MWh,Min_Cap_MW,Min_Cap_MWh,Inv_Cost_per_MWyr,Inv_Cost_per_MWhyr,Fixed_OM_Cost_per_MWyr,Fixed_OM_Cost_per_MWhyr,Var_OM_Cost_per_MWh,Var_OM_Cost_per_MWh_In,Self_Disch,Eff_Up,Eff_Down,Min_Duration,Max_Duration,Reg_Max,Rsv_Max,Reg_Cost,Rsv_Cost,region,cluster
+MA_battery,1,1,1,0,0,0,-1,-1,0,0,19584,22494,4895,5622,0.15,0.15,0,0.92,0.92,1,10,0,0,0,0,MA,0
+CT_battery,2,1,1,0,0,0,-1,-1,0,0,19584,22494,4895,5622,0.15,0.15,0,0.92,0.92,1,10,0,0,0,0,CT,0
+ME_battery,3,1,1,0,0,0,-1,-1,0,0,19584,22494,4895,5622,0.15,0.15,0,0.92,0.92,1,10,0,0,0,0,ME,0

example_systems/1_three_zones_asymmetric_bidirectional_Actual_Symmetric 2/resources/Thermal.csv

@@ -0,0 +1,4 @@
+Resource,Zone,Model,New_Build,Can_Retire,Existing_Cap_MW,Max_Cap_MW,Min_Cap_MW,Inv_Cost_per_MWyr,Fixed_OM_Cost_per_MWyr,Var_OM_Cost_per_MWh,Heat_Rate_MMBTU_per_MWh,Fuel,Cap_Size,Start_Cost_per_MW,Start_Fuel_MMBTU_per_MW,Up_Time,Down_Time,Ramp_Up_Percentage,Ramp_Dn_Percentage,Min_Power,Reg_Max,Rsv_Max,Reg_Cost,Rsv_Cost,region,cluster
+MA_natural_gas_combined_cycle,1,1,1,0,0,-1,0,65400,10287,3.55,7.43,MA_NG,250,91,2,6,6,0.64,0.64,0.468,0.25,0.5,0,0,MA,1
+CT_natural_gas_combined_cycle,2,1,1,0,0,-1,0,65400,9698,3.57,7.12,CT_NG,250,91,2,6,6,0.64,0.64,0.338,0.133332722,0.266665444,0,0,CT,1
+ME_natural_gas_combined_cycle,3,1,1,0,0,-1,0,65400,16291,4.5,12.62,ME_NG,250,91,2,6,6,0.64,0.64,0.474,0.033333333,0.066666667,0,0,ME,1

example_systems/1_three_zones_asymmetric_bidirectional_Actual_Symmetric 2/resources/Vre.csv

@@ -0,0 +1,5 @@
+Resource,Zone,Num_VRE_Bins,New_Build,Can_Retire,Existing_Cap_MW,Max_Cap_MW,Min_Cap_MW,Inv_Cost_per_MWyr,Fixed_OM_Cost_per_MWyr,Var_OM_Cost_per_MWh,Reg_Max,Rsv_Max,Reg_Cost,Rsv_Cost,region,cluster
+MA_solar_pv,1,1,1,0,0,-1,0,85300,18760,0,0,0,0,0,MA,1
+CT_onshore_wind,2,1,1,0,0,-1,0,1,1,0.1,0,0,0,0,CT,1
+CT_solar_pv,2,1,1,0,0,-1,0,1,1,0,0,0,0,0,CT,1
+ME_onshore_wind,3,1,1,0,0,-1,0,97200,43205,0.1,0,0,0,0,ME,1

example_systems/1_three_zones_asymmetric_bidirectional_Actual_Symmetric 2/resources/policy_assignments/Resource_minimum_capacity_requirement.csv

@@ -0,0 +1,6 @@
+Resource,Min_Cap_1,Min_Cap_2,Min_Cap_3
+MA_solar_pv,1,0,0
+CT_onshore_wind,0,1,0
+MA_battery,0,0,1
+CT_battery,0,0,1
+ME_battery,0,0,1

example_systems/1_three_zones_asymmetric_bidirectional_Actual_Symmetric 2/settings/clp_settings.yml

@@ -0,0 +1,14 @@
+# Clp Solver parameters https://github.com/jump-dev/Clp.jl
+# Common solver settings
+Feasib_Tol: 1e-5              # Primal/Dual feasibility tolerance
+TimeLimit: -1.0               # Terminate after this many seconds have passed. A negative value means no time limit
+Pre_Solve: 0                  # Set to 1 to disable presolve
+Method: 5                     # Solution method: dual simplex (0), primal simplex (1), sprint (2), barrier with crossover (3), barrier without crossover (4), automatic (5)
+
+#Clp-specific solver settings
+DualObjectiveLimit: 1e308     # When using dual simplex (where the objective is monotonically changing), terminate when the objective exceeds this limit
+MaximumIterations: 2147483647 # Terminate after performing this number of simplex iterations
+LogLevel: 1                   # Set to 1, 2, 3, or 4 for increasing output. Set to 0 to disable output
+InfeasibleReturn: 0           # Set to 1 to return as soon as the problem is found to be infeasible (by default, an infeasibility proof is computed as well)
+Scaling: 3                    # 0 -off, 1 equilibrium, 2 geometric, 3 auto, 4 dynamic(later)
+Perturbation: 100             # switch on perturbation (50), automatic (100), don't try perturbing (102)

example_systems/1_three_zones_asymmetric_bidirectional_Actual_Symmetric 2/settings/cplex_settings.yml

@@ -0,0 +1,10 @@
+# CPLEX Solver Parameters
+Feasib_Tol: 1.0e-05 # Constraint (primal) feasibility tolerances.
+Optimal_Tol: 1e-5 # Dual feasibility tolerances.
+Pre_Solve: 1 # Controls presolve level.
+TimeLimit: 110000 # Limits total time solver.
+MIPGap: 1e-3 # Relative (p.u. of optimal) mixed integer optimality tolerance for MIP problems (ignored otherwise).
+Method: 2 # Algorithm used to solve continuous models (including MIP root relaxation).
+BarConvTol: 1.0e-08 # Barrier convergence tolerance (determines when barrier terminates).
+NumericFocus: 0 # Numerical precision emphasis.
+SolutionType: 2 # Solution type for LP or QP.