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:40, in Energyscope.es_model(self)
     39             modules.activate(ampl_uuid)
---> 40         self.__es_model = AMPL()
     41 except SystemError:
     42     # 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 = 20260325; 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
69c3f78b 1dc11d18 d62162c9 14b14c35 0b418d08 25d60012 109f1051 34cbe8c2
b2b238a0 7e7e2a0c 24bb31bf 03e30000 95774846 e2e0cb23 c31ce5c0 7fde72be
3fcea157 b63ec984 afd17bc4 a0415fe5 c47e7030 17f03473 70c96730 17c3be80
6aaac0d0 540aa293 50554f11 9d1b6641 bca973e3 b516b670 ef005e43 22e22327
d7644845 f5a97855 51107c33 f0d7fa9f c300c8c5 624ca5c6 bcb217b7 ccdfc395
a5c8cf15 80b8d144 5fcda354 c434de26 ab3a2353 b6bcbf01 26db94a4 f494113a
78a40913 e40c4907 ee17f051 4265cccb d2818757 43f0b4e3 5304bac3 345ceb18
022df8c4 502724f7 83873517 2549a305 f95881a7 87fbe651 88347da3 4648f731
742be6c7 518ac423 44a71683 3ffabd4b 51d34103 56f3de77 f04e9ba7 edac4e15
c4b1aea4 b622f7b3 a0dfafb7 b540c741 e4a44fd4 a52e46e4 71e668d4 1f0cc914
85cebf60 52b613a6 200c0040 2feeda30 84ddc616 9596d246 1717cfc0 1211a016
d7f91e10 b267cf9d 85

Computed 1-1-0-10220046
runner--pdbxrphm-project-62195384-concurrent-0
FA:61:DB:37:C3:4C

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:82, in Energyscope.calc(self, ds, parser)
     78 def calc(self, ds: Dataset = None, parser: Callable[[AMPL], Result] = parse_result) -> Result:
     79     """
     80     Calls AMPL with `df` as .dat and returns the parsed result.
     81     """
---> 82     if self.es_model.getSets().__len__() == 0:  # Check if AMPL instance is empty
     83         self._initial_run(ds=ds)
     85     # Solve the model

File /builds/energyscope/energyscope/src/energyscope/energyscope.py:43, in Energyscope.es_model(self)
     40                 self.__es_model = AMPL()
     41         except SystemError:
     42             # Try to create the object a second time to prevent errors when starting `ampl_lic`
---> 43             self.__es_model = AMPL()
     44 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 = 20260325; 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
69c3f78b 1dc11d18 d62162c9 14b14c35 0b418d08 25d60012 109f1051 34cbe8c2
b2b238a0 7e7e2a0c 24bb31bf 03e30000 95774846 e2e0cb23 c31ce5c0 7fde72be
3fcea157 b63ec984 afd17bc4 a0415fe5 c47e7030 17f03473 70c96730 17c3be80
6aaac0d0 540aa293 50554f11 9d1b6641 bca973e3 b516b670 ef005e43 22e22327
d7644845 f5a97855 51107c33 f0d7fa9f c300c8c5 624ca5c6 bcb217b7 ccdfc395
a5c8cf15 80b8d144 5fcda354 c434de26 ab3a2353 b6bcbf01 26db94a4 f494113a
78a40913 e40c4907 ee17f051 4265cccb d2818757 43f0b4e3 5304bac3 345ceb18
022df8c4 502724f7 83873517 2549a305 f95881a7 87fbe651 88347da3 4648f731
742be6c7 518ac423 44a71683 3ffabd4b 51d34103 56f3de77 f04e9ba7 edac4e15
c4b1aea4 b622f7b3 a0dfafb7 b540c741 e4a44fd4 a52e46e4 71e668d4 1f0cc914
85cebf60 52b613a6 200c0040 2feeda30 84ddc616 9596d246 1717cfc0 1211a016
d7f91e10 b267cf9d 85

Computed 1-1-0-10220046
runner--pdbxrphm-project-62195384-concurrent-0
FA:61:DB:37:C3:4C

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