Plot installed capacity

In [1]:

Copied!





"""
Installed  Capacity
===================================

This is an example using the core model.
"""
"""
Installed  Capacity
===================================

This is an example using the core model.
"""

Out[1]:

'\nInstalled  Capacity\n===================================\n\nThis is an example using the core model.\n'

In [2]:

Copied!

import matplotlib.pyplot as plt
from energyscope.energyscope import Energyscope
from energyscope.models import core

import matplotlib.pyplot as plt
from energyscope.energyscope import Energyscope
from energyscope.models import core

In [3]:

Copied!

thresh = 1
figsize = (10, 6)
thresh = 1
figsize = (10, 6)

In [4]:

Copied!

# Create a model and calculate the results
es_core = Energyscope(model=core)
results = es_core.calc()
# Create a model and calculate the results
es_core = Energyscope(model=core)
results = es_core.calc()

[INFO] Activating AMPL license with UUID

---------------------------------------------------------------------------
SystemError                               Traceback (most recent call last)
File /builds/energyscope/energyscope/src/energyscope/energyscope.py:41, in Energyscope.es_model(self)
     40             modules.activate(ampl_uuid)
---> 41         self.__es_model = AMPL()
     42 except SystemError:
     43     # Try to create the object a second time to prevent errors when starting `ampl_lic`

File ~/.local/share/hatch/env/virtual/energyscope/4BIeM0-F/docs/lib/python3.10/site-packages/amplpy/ampl.py:105, in AMPL.__init__(self, environment)
    104 if environment is None:
--> 105     self._impl = amplpython.AMPL()
    106 else:

SystemError: License not valid.
Message: 
License file /root/.local/share/hatch/env/virtual/energyscope/4BIeM0-F/docs/lib/python3.10/site-packages/ampl_module_base/bin/ampl.lic:
Today = 20260223; found license file "/root/.local/share/hatch/env/virtual/energyscope/4BIeM0-F/docs/lib/python3.10/site-packages/ampl_module_base/bin/ampl.lic":
# Bundle #7221.7759 expiring 20260215: ME-409 Energy conversion and renewable energy, Prof. Fran?ois Mar?chal, EPFL
699c7aa2 de55a760 d0520895 371c367c 70249cf6 78eeb7b2 660707cd 82640652
f3e6a306 69c75deb 028e8b63 230c0451 26e61c83 960cc9d8 d790f152 32f080fe
172152d5 58d5ba02 f35b656f c7f96298 663afdc7 5c107dce b1d043a1 71a6b0b5
1145f193 67f0b132 75bbb785 516ea73a 66fc7b36 b4d169ec a5793166 a7ca55bd
e4f419e7 8f286a77 24e4ee80 95369081 f642d101 dff63ab9 a3134e9c b5313012
b1ef53d9 1845eebf e0dc2f6f 7719c263 f360ef2d 988352f3 65e8dca1 a26287c5
e3566e89 725fc641 246cefd9 b7e917fb 575d3f0c 08e7b718 f54e3b9f 33957656
7384f7b0 92058c20 57b29e35 8397911c d665f1df bc449a42 a73eb3c3 d5108181
062aa5d8 2b2241c5 06aa15a8 414ba20e 6591b8ce 96087be8 50ea1ef0 f5f8e772
d59fe7e0 b434bc37 377db5e8 6134f09f 462ffa28 85769a8b 272648dc 8699728b
964cf499 705750de 8449d3e6 275e70f0 8018c075 559bfb3b 55a902ad f3f1a5d1
c4637893 4282727d 74

Computed 1-1-0-eb4da5d1
runner--pdbxrphm-project-62195384-concurrent-0
32:6C:CD:60:38:18

HOSTINFO.NPROC.4
HOSTINFO.CORES.4
HOSTINFO.PLATFORM.linux64

No ampl license for this machine.



During handling of the above exception, another exception occurred:

SystemError                               Traceback (most recent call last)
Cell In[4], line 3
      1 # Create a model and calculate the results
      2 es_core = Energyscope(model=core)
----> 3 results = es_core.calc()

File /builds/energyscope/energyscope/src/energyscope/energyscope.py:83, in Energyscope.calc(self, ds, parser)
     79 def calc(self, ds: Dataset = None, parser: Callable[[AMPL], Result] = parse_result) -> Result:
     80     """
     81     Calls AMPL with `df` as .dat and returns the parsed result.
     82     """
---> 83     if self.es_model.getSets().__len__() == 0:  # Check if AMPL instance is empty
     84         self._initial_run(ds=ds)
     86     # Solve the model

File /builds/energyscope/energyscope/src/energyscope/energyscope.py:44, in Energyscope.es_model(self)
     41                 self.__es_model = AMPL()
     42         except SystemError:
     43             # Try to create the object a second time to prevent errors when starting `ampl_lic`
---> 44             self.__es_model = AMPL()
     45 return self.__es_model

File ~/.local/share/hatch/env/virtual/energyscope/4BIeM0-F/docs/lib/python3.10/site-packages/amplpy/ampl.py:105, in AMPL.__init__(self, environment)
    103     environment = Environment()
    104 if environment is None:
--> 105     self._impl = amplpython.AMPL()
    106 else:
    107     self._impl = amplpython.AMPL(environment._impl)

SystemError: License not valid.
Message: 
License file /root/.local/share/hatch/env/virtual/energyscope/4BIeM0-F/docs/lib/python3.10/site-packages/ampl_module_base/bin/ampl.lic:
Today = 20260223; found license file "/root/.local/share/hatch/env/virtual/energyscope/4BIeM0-F/docs/lib/python3.10/site-packages/ampl_module_base/bin/ampl.lic":
# Bundle #7221.7759 expiring 20260215: ME-409 Energy conversion and renewable energy, Prof. Fran?ois Mar?chal, EPFL
699c7aa2 de55a760 d0520895 371c367c 70249cf6 78eeb7b2 660707cd 82640652
f3e6a306 69c75deb 028e8b63 230c0451 26e61c83 960cc9d8 d790f152 32f080fe
172152d5 58d5ba02 f35b656f c7f96298 663afdc7 5c107dce b1d043a1 71a6b0b5
1145f193 67f0b132 75bbb785 516ea73a 66fc7b36 b4d169ec a5793166 a7ca55bd
e4f419e7 8f286a77 24e4ee80 95369081 f642d101 dff63ab9 a3134e9c b5313012
b1ef53d9 1845eebf e0dc2f6f 7719c263 f360ef2d 988352f3 65e8dca1 a26287c5
e3566e89 725fc641 246cefd9 b7e917fb 575d3f0c 08e7b718 f54e3b9f 33957656
7384f7b0 92058c20 57b29e35 8397911c d665f1df bc449a42 a73eb3c3 d5108181
062aa5d8 2b2241c5 06aa15a8 414ba20e 6591b8ce 96087be8 50ea1ef0 f5f8e772
d59fe7e0 b434bc37 377db5e8 6134f09f 462ffa28 85769a8b 272648dc 8699728b
964cf499 705750de 8449d3e6 275e70f0 8018c075 559bfb3b 55a902ad f3f1a5d1
c4637893 4282727d 74

Computed 1-1-0-eb4da5d1
runner--pdbxrphm-project-62195384-concurrent-0
32:6C:CD:60:38:18

HOSTINFO.NPROC.4
HOSTINFO.CORES.4
HOSTINFO.PLATFORM.linux64

No ampl license for this machine.

In [5]:

Copied!





F = results.variables["F"]
storages = results.sets['STORAGE_TECH']["STORAGE_TECH"].values

F_storage = F[F.index.get_level_values(0).isin(storages)]
F_tech = F[~F.index.get_level_values(0).isin(storages)]

# Sort by size
F_storage = F_storage.sort_values(by=F_storage.columns[0], ascending=False)
F_storage = F_storage[F_storage.values[:,0]>thresh]
F_tech = F_tech[F_tech.values[:,0]>thresh]
F_tech = F_tech.sort_values(by=F_tech.columns[0], ascending=False)


# Calculate subplot widths proportional to number of bars
n_tech = len(F_tech)
n_storage = len(F_storage)
total = n_tech + n_storage
widths = [n_tech, n_storage] if total > 0 else [1, 1]

fig, axes = plt.subplots(1, 2, figsize=figsize, gridspec_kw={'width_ratios': widths})

# Technologies subplot (left)
axes[0].bar(x=range(n_tech), height=F_tech.values[:, 0])
axes[0].set_xticks(range(n_tech))
axes[0].set_xticklabels(F_tech.index, rotation=90)
axes[0].set_ylabel('Installed Capacity of Conversion Technologies [GW]')
axes[0].set_xlim(-1, n_tech)
axes[0].set_title('Conversion Technologies', fontsize=20)

# Storages subplot (right)
axes[1].bar(x=range(n_storage), height=F_storage.values[:, 0], color='orange')
axes[1].set_xticks(range(n_storage))
axes[1].set_xticklabels(F_storage.index, rotation=90)
axes[1].set_ylabel('Installed Capacity of Storages [GWh]')
axes[1].set_xlim(-1, n_storage)
axes[1].set_title('Storages', fontsize=20)

plt.tight_layout()
plt.show()
F = results.variables["F"]
storages = results.sets['STORAGE_TECH']["STORAGE_TECH"].values

F_storage = F[F.index.get_level_values(0).isin(storages)]
F_tech = F[~F.index.get_level_values(0).isin(storages)]

# Sort by size
F_storage = F_storage.sort_values(by=F_storage.columns[0], ascending=False)
F_storage = F_storage[F_storage.values[:,0]>thresh]
F_tech = F_tech[F_tech.values[:,0]>thresh]
F_tech = F_tech.sort_values(by=F_tech.columns[0], ascending=False)


# Calculate subplot widths proportional to number of bars
n_tech = len(F_tech)
n_storage = len(F_storage)
total = n_tech + n_storage
widths = [n_tech, n_storage] if total > 0 else [1, 1]

fig, axes = plt.subplots(1, 2, figsize=figsize, gridspec_kw={'width_ratios': widths})

# Technologies subplot (left)
axes[0].bar(x=range(n_tech), height=F_tech.values[:, 0])
axes[0].set_xticks(range(n_tech))
axes[0].set_xticklabels(F_tech.index, rotation=90)
axes[0].set_ylabel('Installed Capacity of Conversion Technologies [GW]')
axes[0].set_xlim(-1, n_tech)
axes[0].set_title('Conversion Technologies', fontsize=20)

# Storages subplot (right)
axes[1].bar(x=range(n_storage), height=F_storage.values[:, 0], color='orange')
axes[1].set_xticks(range(n_storage))
axes[1].set_xticklabels(F_storage.index, rotation=90)
axes[1].set_ylabel('Installed Capacity of Storages [GWh]')
axes[1].set_xlim(-1, n_storage)
axes[1].set_title('Storages', fontsize=20)

plt.tight_layout()
plt.show()

---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
Cell In[5], line 1
----> 1 F = results.variables["F"]
      2 storages = results.sets['STORAGE_TECH']["STORAGE_TECH"].values
      4 F_storage = F[F.index.get_level_values(0).isin(storages)]

NameError: name 'results' is not defined