diff --git a/dev/Untitled1.ipynb b/dev/Untitled1.ipynb
index d6e9667b..7b49e66e 100644
--- a/dev/Untitled1.ipynb
+++ b/dev/Untitled1.ipynb
@@ -81,10 +81,10 @@
     "        #{\"model\":\"remind\", \"pathway\":\"SSP2-Base\", \"year\":2010},\n",
     "        #{\"model\":\"image\", \"pathway\":\"SSP2-Base\", \"year\":2020},\n",
     "        #{\"model\":\"remind\", \"pathway\":\"SSP2-Base\", \"year\":2050},\n",
-    "        #{\"model\":\"remind\", \"pathway\":\"SSP2-PkBudg1150\", \"year\":2050},\n",
-    "        {\"model\":\"image\", \"pathway\":\"SSP2-RCP26\", \"year\":2040},\n",
-    "        {\"model\":\"image\", \"pathway\":\"SSP2-RCP26\", \"year\":2060},\n",
-    "        {\"model\":\"image\", \"pathway\":\"SSP2-RCP19\", \"year\":2050},\n",
+    "        {\"model\":\"remind\", \"pathway\":\"SSP2-PkBudg1150\", \"year\":2050},\n",
+    "        #{\"model\":\"image\", \"pathway\":\"SSP2-RCP26\", \"year\":2040},\n",
+    "        #{\"model\":\"image\", \"pathway\":\"SSP2-RCP26\", \"year\":2060},\n",
+    "        #{\"model\":\"image\", \"pathway\":\"SSP2-RCP19\", \"year\":2050},\n",
     "        #{\"model\":\"image\", \"pathway\":\"SSP2-Base\", \"year\":2050},\n",
     "    ],\n",
     "    source_db=\"ecoinvent-3.10-cutoff\", # <-- name of the database in the BW2 project. Must be a string.\n",
@@ -104,7 +104,7 @@
      "name": "stderr",
      "output_type": "stream",
      "text": [
-      "Processing scenarios: 100%|████████████| 3/3 [06:13<00:00, 124.63s/it]"
+      "Processing scenarios: 100%|█████████████| 1/1 [01:38<00:00, 98.83s/it]"
      ]
     },
     {
@@ -192,7 +192,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 6,
+   "execution_count": 5,
    "id": "dff4efb1-69ab-4bd5-8d52-4038df180a85",
    "metadata": {},
    "outputs": [
@@ -211,55 +211,7 @@
      "text": [
       "Writing activities to SQLite3 database:\n",
       "0% [##############################] 100% | ETA: 00:00:00\n",
-      "Total time elapsed: 00:00:41\n"
-     ]
-    },
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "Title: Writing activities to SQLite3 database:\n",
-      "  Started: 07/13/2024 18:54:03\n",
-      "  Finished: 07/13/2024 18:54:45\n",
-      "  Total time elapsed: 00:00:41\n",
-      "  CPU %: 88.20\n",
-      "  Memory %: 12.60\n",
-      "Created database: cement rcp26 2040\n",
-      "Running all checks...\n",
-      "Minor anomalies found: check the change report.\n"
-     ]
-    },
-    {
-     "name": "stderr",
-     "output_type": "stream",
-     "text": [
-      "Writing activities to SQLite3 database:\n",
-      "0% [##############################] 100% | ETA: 00:00:00\n",
-      "Total time elapsed: 00:00:52\n"
-     ]
-    },
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "Title: Writing activities to SQLite3 database:\n",
-      "  Started: 07/13/2024 19:03:00\n",
-      "  Finished: 07/13/2024 19:03:53\n",
-      "  Total time elapsed: 00:00:52\n",
-      "  CPU %: 89.80\n",
-      "  Memory %: 8.77\n",
-      "Created database: cement rcp26 2060\n",
-      "Running all checks...\n",
-      "Minor anomalies found: check the change report.\n"
-     ]
-    },
-    {
-     "name": "stderr",
-     "output_type": "stream",
-     "text": [
-      "Writing activities to SQLite3 database:\n",
-      "0% [##############################] 100% | ETA: 00:00:00\n",
-      "Total time elapsed: 00:00:53\n"
+      "Total time elapsed: 00:00:36\n"
      ]
     },
     {
@@ -267,12 +219,12 @@
      "output_type": "stream",
      "text": [
       "Title: Writing activities to SQLite3 database:\n",
-      "  Started: 07/13/2024 19:12:37\n",
-      "  Finished: 07/13/2024 19:13:30\n",
-      "  Total time elapsed: 00:00:53\n",
-      "  CPU %: 89.00\n",
-      "  Memory %: 7.57\n",
-      "Created database: cement rcp19 2\n",
+      "  Started: 07/13/2024 20:50:02\n",
+      "  Finished: 07/13/2024 20:50:39\n",
+      "  Total time elapsed: 00:00:36\n",
+      "  CPU %: 90.50\n",
+      "  Memory %: 11.02\n",
+      "Created database: ei_cutoff_3.10_remind_SSP2-PkBudg1150_2050 2024-07-13\n",
       "Generate scenario report.\n",
       "Report saved under /Users/romain/GitHub/premise/dev/export/scenario_report.\n",
       "Generate change report.\n",
@@ -281,202 +233,413 @@
     }
    ],
    "source": [
-    "ndb.write_db_to_brightway([\"cement rcp26 2040\", \"cement rcp26 2060\", \"cement rcp19 2\"])"
+    "ndb.write_db_to_brightway()"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 7,
+   "execution_count": 8,
    "id": "a3b56b3f-9583-4c90-ad21-551c5530e27c",
    "metadata": {},
    "outputs": [
     {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "aluminium, primary, liquid\n",
-      "cooling energy, at -25 °c\n",
-      "1,1-dimethylcyclopentane\n",
-      "aluminium, primary, liquid\n",
-      "cooling energy, at -45 °c\n",
-      "cooling energy, at -45 °c\n",
-      "cooling energy, at -25 °c\n",
-      "fraction 8 from naphtha separation\n",
-      "fraction 8 from naphtha separation\n",
-      "fraction 1 from naphtha separation\n",
-      "2,3-dimethylbutane\n",
-      "heptane\n",
-      "fraction 7 from naphtha separation\n",
-      "2,3-dimethylbutane\n",
-      "cooling energy, at -55 °c\n",
-      "heptane\n",
-      "electricity, low voltage\n",
-      "cooling energy, at -55 °c\n",
-      "cooling energy, at -15 °c\n",
-      "2-methylpentane\n",
-      "sawnwood, paraná pine, dried (u=10%)\n",
-      "hexane\n",
-      "hexane\n",
-      "residual wood, dry\n",
-      "formaldehyde\n",
-      "sawlog and veneer log, paraná pine, measured as solid wood under bark\n",
-      "cooling energy, at -25 °c\n",
-      "cooling energy, at -100 °c\n",
-      "ventilation duct, elbow 90°, steel, 100x50 mm\n",
-      "cooling energy, at -160 °c\n",
-      "cooling energy, at -15 °c\n",
-      "cooling energy, at -45 °c\n",
-      "steam, in chemical industry\n",
-      "ventilation duct, elbow 90°, steel, 100x50 mm\n",
-      "cooling energy, at -100 °c\n",
-      "sawnwood, paraná pine, dried (u=10%)\n",
-      "residual wood, dry\n",
-      "aluminium, primary, liquid\n",
-      "sawnwood, paraná pine, dried (u=10%)\n",
-      "sawlog and veneer log, paraná pine, measured as solid wood under bark\n",
-      "fraction 1 from naphtha separation\n",
-      "cooling energy, at -100 °c\n",
-      "sawlog and veneer log, paraná pine, measured as solid wood under bark\n",
-      "aluminium, primary, liquid\n",
-      "steam, in chemical industry\n",
-      "sawnwood, paraná pine, dried (u=10%)\n",
-      "2-methylpentane\n",
-      "sawnwood, paraná pine, dried (u=10%)\n",
-      "cooling energy, at -160 °c\n",
-      "fraction 7 from naphtha separation\n",
-      "ventilation duct, elbow 90°, steel, 100x50 mm\n",
-      "electricity, low voltage\n",
-      "1,1-dimethylcyclopentane\n",
-      "aluminium, primary, liquid\n",
-      "sawlog and veneer log, paraná pine, measured as solid wood under bark\n",
-      "formaldehyde\n",
-      "methylcyclopentane\n",
-      "methylcyclohexane\n",
-      "methylcyclopentane\n",
-      "methylcyclohexane\n",
-      "electricity, high voltage\n",
-      "electricity, high voltage\n",
-      "electricity, high voltage\n",
-      "electricity, high voltage\n",
-      "electricity, high voltage\n",
-      "electricity, high voltage\n",
-      "electricity, high voltage\n",
-      "electricity, high voltage\n",
-      "electricity, high voltage\n",
-      "electricity, high voltage\n",
-      "electricity, high voltage\n",
-      "electricity, high voltage\n",
-      "electricity, high voltage\n",
-      "electricity, high voltage\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "natural gas, high pressure\n",
-      "transport, passenger car\n",
-      "transport, passenger car\n",
-      "transport, passenger car\n",
-      "transport, passenger car\n",
-      "transport, passenger car\n",
-      "transport, passenger car\n",
-      "transport, passenger car\n",
-      "transport, passenger car\n",
-      "transport, passenger car\n",
-      "transport, passenger car\n",
-      "transport, passenger car\n",
-      "transport, passenger car\n",
-      "transport, passenger car\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, freight, lorry\n",
-      "transport, passenger bus\n",
-      "transport, passenger bus\n",
-      "transport, passenger bus\n",
-      "transport, passenger bus\n",
-      "transport, passenger bus\n",
-      "transport, passenger bus\n",
-      "transport, passenger bus\n",
-      "transport, passenger bus\n",
-      "transport, passenger bus\n",
-      "transport, passenger bus\n",
-      "transport, passenger bus\n",
-      "transport, passenger bus\n",
-      "transport, passenger bus\n"
-     ]
+     "data": {
+      "text/html": [
+       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
+       "<defs>\n",
+       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
+       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
+       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
+       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
+       "</symbol>\n",
+       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
+       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
+       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
+       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
+       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
+       "</symbol>\n",
+       "</defs>\n",
+       "</svg>\n",
+       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
+       " *\n",
+       " */\n",
+       "\n",
+       ":root {\n",
+       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
+       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
+       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
+       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
+       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
+       "  --xr-background-color: var(--jp-layout-color0, white);\n",
+       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
+       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
+       "}\n",
+       "\n",
+       "html[theme=dark],\n",
+       "body[data-theme=dark],\n",
+       "body.vscode-dark {\n",
+       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
+       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
+       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
+       "  --xr-border-color: #1F1F1F;\n",
+       "  --xr-disabled-color: #515151;\n",
+       "  --xr-background-color: #111111;\n",
+       "  --xr-background-color-row-even: #111111;\n",
+       "  --xr-background-color-row-odd: #313131;\n",
+       "}\n",
+       "\n",
+       ".xr-wrap {\n",
+       "  display: block !important;\n",
+       "  min-width: 300px;\n",
+       "  max-width: 700px;\n",
+       "}\n",
+       "\n",
+       ".xr-text-repr-fallback {\n",
+       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
+       "  display: none;\n",
+       "}\n",
+       "\n",
+       ".xr-header {\n",
+       "  padding-top: 6px;\n",
+       "  padding-bottom: 6px;\n",
+       "  margin-bottom: 4px;\n",
+       "  border-bottom: solid 1px var(--xr-border-color);\n",
+       "}\n",
+       "\n",
+       ".xr-header > div,\n",
+       ".xr-header > ul {\n",
+       "  display: inline;\n",
+       "  margin-top: 0;\n",
+       "  margin-bottom: 0;\n",
+       "}\n",
+       "\n",
+       ".xr-obj-type,\n",
+       ".xr-array-name {\n",
+       "  margin-left: 2px;\n",
+       "  margin-right: 10px;\n",
+       "}\n",
+       "\n",
+       ".xr-obj-type {\n",
+       "  color: var(--xr-font-color2);\n",
+       "}\n",
+       "\n",
+       ".xr-sections {\n",
+       "  padding-left: 0 !important;\n",
+       "  display: grid;\n",
+       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
+       "}\n",
+       "\n",
+       ".xr-section-item {\n",
+       "  display: contents;\n",
+       "}\n",
+       "\n",
+       ".xr-section-item input {\n",
+       "  display: none;\n",
+       "}\n",
+       "\n",
+       ".xr-section-item input + label {\n",
+       "  color: var(--xr-disabled-color);\n",
+       "}\n",
+       "\n",
+       ".xr-section-item input:enabled + label {\n",
+       "  cursor: pointer;\n",
+       "  color: var(--xr-font-color2);\n",
+       "}\n",
+       "\n",
+       ".xr-section-item input:enabled + label:hover {\n",
+       "  color: var(--xr-font-color0);\n",
+       "}\n",
+       "\n",
+       ".xr-section-summary {\n",
+       "  grid-column: 1;\n",
+       "  color: var(--xr-font-color2);\n",
+       "  font-weight: 500;\n",
+       "}\n",
+       "\n",
+       ".xr-section-summary > span {\n",
+       "  display: inline-block;\n",
+       "  padding-left: 0.5em;\n",
+       "}\n",
+       "\n",
+       ".xr-section-summary-in:disabled + label {\n",
+       "  color: var(--xr-font-color2);\n",
+       "}\n",
+       "\n",
+       ".xr-section-summary-in + label:before {\n",
+       "  display: inline-block;\n",
+       "  content: '►';\n",
+       "  font-size: 11px;\n",
+       "  width: 15px;\n",
+       "  text-align: center;\n",
+       "}\n",
+       "\n",
+       ".xr-section-summary-in:disabled + label:before {\n",
+       "  color: var(--xr-disabled-color);\n",
+       "}\n",
+       "\n",
+       ".xr-section-summary-in:checked + label:before {\n",
+       "  content: '▼';\n",
+       "}\n",
+       "\n",
+       ".xr-section-summary-in:checked + label > span {\n",
+       "  display: none;\n",
+       "}\n",
+       "\n",
+       ".xr-section-summary,\n",
+       ".xr-section-inline-details {\n",
+       "  padding-top: 4px;\n",
+       "  padding-bottom: 4px;\n",
+       "}\n",
+       "\n",
+       ".xr-section-inline-details {\n",
+       "  grid-column: 2 / -1;\n",
+       "}\n",
+       "\n",
+       ".xr-section-details {\n",
+       "  display: none;\n",
+       "  grid-column: 1 / -1;\n",
+       "  margin-bottom: 5px;\n",
+       "}\n",
+       "\n",
+       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
+       "  display: contents;\n",
+       "}\n",
+       "\n",
+       ".xr-array-wrap {\n",
+       "  grid-column: 1 / -1;\n",
+       "  display: grid;\n",
+       "  grid-template-columns: 20px auto;\n",
+       "}\n",
+       "\n",
+       ".xr-array-wrap > label {\n",
+       "  grid-column: 1;\n",
+       "  vertical-align: top;\n",
+       "}\n",
+       "\n",
+       ".xr-preview {\n",
+       "  color: var(--xr-font-color3);\n",
+       "}\n",
+       "\n",
+       ".xr-array-preview,\n",
+       ".xr-array-data {\n",
+       "  padding: 0 5px !important;\n",
+       "  grid-column: 2;\n",
+       "}\n",
+       "\n",
+       ".xr-array-data,\n",
+       ".xr-array-in:checked ~ .xr-array-preview {\n",
+       "  display: none;\n",
+       "}\n",
+       "\n",
+       ".xr-array-in:checked ~ .xr-array-data,\n",
+       ".xr-array-preview {\n",
+       "  display: inline-block;\n",
+       "}\n",
+       "\n",
+       ".xr-dim-list {\n",
+       "  display: inline-block !important;\n",
+       "  list-style: none;\n",
+       "  padding: 0 !important;\n",
+       "  margin: 0;\n",
+       "}\n",
+       "\n",
+       ".xr-dim-list li {\n",
+       "  display: inline-block;\n",
+       "  padding: 0;\n",
+       "  margin: 0;\n",
+       "}\n",
+       "\n",
+       ".xr-dim-list:before {\n",
+       "  content: '(';\n",
+       "}\n",
+       "\n",
+       ".xr-dim-list:after {\n",
+       "  content: ')';\n",
+       "}\n",
+       "\n",
+       ".xr-dim-list li:not(:last-child):after {\n",
+       "  content: ',';\n",
+       "  padding-right: 5px;\n",
+       "}\n",
+       "\n",
+       ".xr-has-index {\n",
+       "  font-weight: bold;\n",
+       "}\n",
+       "\n",
+       ".xr-var-list,\n",
+       ".xr-var-item {\n",
+       "  display: contents;\n",
+       "}\n",
+       "\n",
+       ".xr-var-item > div,\n",
+       ".xr-var-item label,\n",
+       ".xr-var-item > .xr-var-name span {\n",
+       "  background-color: var(--xr-background-color-row-even);\n",
+       "  margin-bottom: 0;\n",
+       "}\n",
+       "\n",
+       ".xr-var-item > .xr-var-name:hover span {\n",
+       "  padding-right: 5px;\n",
+       "}\n",
+       "\n",
+       ".xr-var-list > li:nth-child(odd) > div,\n",
+       ".xr-var-list > li:nth-child(odd) > label,\n",
+       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
+       "  background-color: var(--xr-background-color-row-odd);\n",
+       "}\n",
+       "\n",
+       ".xr-var-name {\n",
+       "  grid-column: 1;\n",
+       "}\n",
+       "\n",
+       ".xr-var-dims {\n",
+       "  grid-column: 2;\n",
+       "}\n",
+       "\n",
+       ".xr-var-dtype {\n",
+       "  grid-column: 3;\n",
+       "  text-align: right;\n",
+       "  color: var(--xr-font-color2);\n",
+       "}\n",
+       "\n",
+       ".xr-var-preview {\n",
+       "  grid-column: 4;\n",
+       "}\n",
+       "\n",
+       ".xr-index-preview {\n",
+       "  grid-column: 2 / 5;\n",
+       "  color: var(--xr-font-color2);\n",
+       "}\n",
+       "\n",
+       ".xr-var-name,\n",
+       ".xr-var-dims,\n",
+       ".xr-var-dtype,\n",
+       ".xr-preview,\n",
+       ".xr-attrs dt {\n",
+       "  white-space: nowrap;\n",
+       "  overflow: hidden;\n",
+       "  text-overflow: ellipsis;\n",
+       "  padding-right: 10px;\n",
+       "}\n",
+       "\n",
+       ".xr-var-name:hover,\n",
+       ".xr-var-dims:hover,\n",
+       ".xr-var-dtype:hover,\n",
+       ".xr-attrs dt:hover {\n",
+       "  overflow: visible;\n",
+       "  width: auto;\n",
+       "  z-index: 1;\n",
+       "}\n",
+       "\n",
+       ".xr-var-attrs,\n",
+       ".xr-var-data,\n",
+       ".xr-index-data {\n",
+       "  display: none;\n",
+       "  background-color: var(--xr-background-color) !important;\n",
+       "  padding-bottom: 5px !important;\n",
+       "}\n",
+       "\n",
+       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
+       ".xr-var-data-in:checked ~ .xr-var-data,\n",
+       ".xr-index-data-in:checked ~ .xr-index-data {\n",
+       "  display: block;\n",
+       "}\n",
+       "\n",
+       ".xr-var-data > table {\n",
+       "  float: right;\n",
+       "}\n",
+       "\n",
+       ".xr-var-name span,\n",
+       ".xr-var-data,\n",
+       ".xr-index-name div,\n",
+       ".xr-index-data,\n",
+       ".xr-attrs {\n",
+       "  padding-left: 25px !important;\n",
+       "}\n",
+       "\n",
+       ".xr-attrs,\n",
+       ".xr-var-attrs,\n",
+       ".xr-var-data,\n",
+       ".xr-index-data {\n",
+       "  grid-column: 1 / -1;\n",
+       "}\n",
+       "\n",
+       "dl.xr-attrs {\n",
+       "  padding: 0;\n",
+       "  margin: 0;\n",
+       "  display: grid;\n",
+       "  grid-template-columns: 125px auto;\n",
+       "}\n",
+       "\n",
+       ".xr-attrs dt,\n",
+       ".xr-attrs dd {\n",
+       "  padding: 0;\n",
+       "  margin: 0;\n",
+       "  float: left;\n",
+       "  padding-right: 10px;\n",
+       "  width: auto;\n",
+       "}\n",
+       "\n",
+       ".xr-attrs dt {\n",
+       "  font-weight: normal;\n",
+       "  grid-column: 1;\n",
+       "}\n",
+       "\n",
+       ".xr-attrs dt:hover span {\n",
+       "  display: inline-block;\n",
+       "  background: var(--xr-background-color);\n",
+       "  padding-right: 10px;\n",
+       "}\n",
+       "\n",
+       ".xr-attrs dd {\n",
+       "  grid-column: 2;\n",
+       "  white-space: pre-wrap;\n",
+       "  word-break: break-all;\n",
+       "}\n",
+       "\n",
+       ".xr-icon-database,\n",
+       ".xr-icon-file-text2,\n",
+       ".xr-no-icon {\n",
+       "  display: inline-block;\n",
+       "  vertical-align: middle;\n",
+       "  width: 1em;\n",
+       "  height: 1.5em !important;\n",
+       "  stroke-width: 0;\n",
+       "  stroke: currentColor;\n",
+       "  fill: currentColor;\n",
+       "}\n",
+       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.DataArray &#x27;value&#x27; (region: 13)&gt; Size: 104B\n",
+       "array([0.24319538, 0.11075869, 0.37354394, 0.1023092 , 0.41335137,\n",
+       "       0.10468423, 0.19437585, 0.19803933, 0.10542634, 0.1051091 ,\n",
+       "       0.09484875, 0.39999575, 0.14865214])\n",
+       "Coordinates:\n",
+       "  * region     (region) object 104B &#x27;CAZ&#x27; &#x27;CHA&#x27; &#x27;EUR&#x27; ... &#x27;SSA&#x27; &#x27;USA&#x27; &#x27;World&#x27;\n",
+       "    year       int64 8B 2050\n",
+       "    variables  &lt;U6 24B &#x27;cement&#x27;</pre><div class='xr-wrap' style='display:none'><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'>'value'</div><ul class='xr-dim-list'><li><span class='xr-has-index'>region</span>: 13</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-f602e43c-d79c-429e-9ec5-379dea27c84e' class='xr-array-in' type='checkbox' checked><label for='section-f602e43c-d79c-429e-9ec5-379dea27c84e' title='Show/hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-array-preview xr-preview'><span>0.2432 0.1108 0.3735 0.1023 0.4134 ... 0.1051 0.09485 0.4 0.1487</span></div><div class='xr-array-data'><pre>array([0.24319538, 0.11075869, 0.37354394, 0.1023092 , 0.41335137,\n",
+       "       0.10468423, 0.19437585, 0.19803933, 0.10542634, 0.1051091 ,\n",
+       "       0.09484875, 0.39999575, 0.14865214])</pre></div></div></li><li class='xr-section-item'><input id='section-76d18dfd-3edb-4c62-9d10-a45e6b21e369' class='xr-section-summary-in' type='checkbox'  checked><label for='section-76d18dfd-3edb-4c62-9d10-a45e6b21e369' class='xr-section-summary' >Coordinates: <span>(3)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>region</span></div><div class='xr-var-dims'>(region)</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>&#x27;CAZ&#x27; &#x27;CHA&#x27; &#x27;EUR&#x27; ... &#x27;USA&#x27; &#x27;World&#x27;</div><input id='attrs-b7437c28-1a31-4bc4-b0d4-7082f5122190' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-b7437c28-1a31-4bc4-b0d4-7082f5122190' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-770b5715-8c37-4964-948c-6051a2cefbd7' class='xr-var-data-in' type='checkbox'><label for='data-770b5715-8c37-4964-948c-6051a2cefbd7' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;CAZ&#x27;, &#x27;CHA&#x27;, &#x27;EUR&#x27;, &#x27;IND&#x27;, &#x27;JPN&#x27;, &#x27;LAM&#x27;, &#x27;MEA&#x27;, &#x27;NEU&#x27;, &#x27;OAS&#x27;, &#x27;REF&#x27;,\n",
+       "       &#x27;SSA&#x27;, &#x27;USA&#x27;, &#x27;World&#x27;], dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>year</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>2050</div><input id='attrs-fdd9de87-bb6c-4cdd-94c8-5834bf5ad0d5' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-fdd9de87-bb6c-4cdd-94c8-5834bf5ad0d5' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-6c1c6d7a-9470-45c9-8589-f91b1778cdc5' class='xr-var-data-in' type='checkbox'><label for='data-6c1c6d7a-9470-45c9-8589-f91b1778cdc5' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array(2050)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>variables</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U6</div><div class='xr-var-preview xr-preview'>&#x27;cement&#x27;</div><input id='attrs-95259748-6c7b-499a-8459-3cab307c1f5d' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-95259748-6c7b-499a-8459-3cab307c1f5d' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-440f804c-f6ef-4671-b92b-9fe6de93e75e' class='xr-var-data-in' type='checkbox'><label for='data-440f804c-f6ef-4671-b92b-9fe6de93e75e' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array(&#x27;cement&#x27;, dtype=&#x27;&lt;U6&#x27;)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-94c57ce9-ac59-4e3a-ade9-e2f859e7cb61' class='xr-section-summary-in' type='checkbox'  ><label for='section-94c57ce9-ac59-4e3a-ade9-e2f859e7cb61' class='xr-section-summary' >Indexes: <span>(1)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-index-name'><div>region</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-8bbd4778-01fe-49ee-84fa-b9fd1d5a0110' class='xr-index-data-in' type='checkbox'/><label for='index-8bbd4778-01fe-49ee-84fa-b9fd1d5a0110' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Index([&#x27;CAZ&#x27;, &#x27;CHA&#x27;, &#x27;EUR&#x27;, &#x27;IND&#x27;, &#x27;JPN&#x27;, &#x27;LAM&#x27;, &#x27;MEA&#x27;, &#x27;NEU&#x27;, &#x27;OAS&#x27;, &#x27;REF&#x27;,\n",
+       "       &#x27;SSA&#x27;, &#x27;USA&#x27;, &#x27;World&#x27;],\n",
+       "      dtype=&#x27;object&#x27;, name=&#x27;region&#x27;))</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-c441cfa1-afb7-48bb-ae55-5f001e6db7fb' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-c441cfa1-afb7-48bb-ae55-5f001e6db7fb' class='xr-section-summary'  title='Expand/collapse section'>Attributes: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'></dl></div></li></ul></div></div>"
+      ],
+      "text/plain": [
+       "<xarray.DataArray 'value' (region: 13)> Size: 104B\n",
+       "array([0.24319538, 0.11075869, 0.37354394, 0.1023092 , 0.41335137,\n",
+       "       0.10468423, 0.19437585, 0.19803933, 0.10542634, 0.1051091 ,\n",
+       "       0.09484875, 0.39999575, 0.14865214])\n",
+       "Coordinates:\n",
+       "  * region     (region) object 104B 'CAZ' 'CHA' 'EUR' ... 'SSA' 'USA' 'World'\n",
+       "    year       int64 8B 2050\n",
+       "    variables  <U6 24B 'cement'"
+      ]
+     },
+     "execution_count": 8,
+     "metadata": {},
+     "output_type": "execute_result"
     }
    ],
    "source": [
-    "for x in ndb.scenarios[1][\"unmatched category flows\"]:\n",
-    "    print(x[1])"
+    "ndb.scenarios[0][\"iam data\"].carbon_capture_rate.sel(variables=\"cement\", year=2050)"
    ]
   },
   {
diff --git a/premise/cement.py b/premise/cement.py
index 1ba94124..34410f39 100644
--- a/premise/cement.py
+++ b/premise/cement.py
@@ -599,25 +599,60 @@ def build_clinker_production_datasets(self) -> list:
                             dataset["exchanges"].append(ccs_exc)
 
                             # Update CO2 exchanges
-                            for exc in dataset["exchanges"]:
-                                if (
-                                    exc["name"].lower().startswith("carbon dioxide")
-                                    and exc["type"] == "biosphere"
-                                ):
+                            for exc in ws.biosphere(
+                                dataset,
+                                ws.contains("name", "Carbon dioxide"),
+                            ):
+                                if exc["name"] == "Carbon dioxide, fossil":
                                     exc["amount"] *= 1 - carbon_capture_rate
 
-                                    if "non-fossil" in exc["name"].lower():
-                                        dataset["log parameters"].update(
-                                            {
-                                                "new biogenic CO2": exc["amount"],
-                                            }
-                                        )
-                                    else:
-                                        dataset["log parameters"].update(
-                                            {
-                                                "new fossil CO2": exc["amount"],
-                                            }
-                                        )
+                                    dataset["log parameters"].update(
+                                        {
+                                            "new fossil CO2": exc["amount"],
+                                        }
+                                    )
+
+                                if exc["name"] == "Carbon dioxide, non-fossil":
+                                    exc["amount"] *= 1 - carbon_capture_rate
+
+                                    dataset["log parameters"].update(
+                                        {
+                                            "new biogenic CO2": exc["amount"],
+                                        }
+                                    )
+
+                                # add a flow of "Carbon dioxide, in air" to reflect
+                                # the permanent storage of biogenic CO2
+                                biogenic_CO2_reduction = (
+                                    dataset["log parameters"]["initial biogenic CO2"]
+                                    - dataset["log parameters"]["new biogenic CO2"]
+                                )
+                                dataset["exchanges"].append(
+                                    {
+                                        "uncertainty type": 0,
+                                        "loc": biogenic_CO2_reduction,
+                                        "amount": biogenic_CO2_reduction,
+                                        "type": "biosphere",
+                                        "name": "Carbon dioxide, in air",
+                                        "unit": "kilogram",
+                                        "categories": (
+                                            "natural resource",
+                                            "in air",
+                                        ),
+                                        "comment": "Permanent storage of biogenic CO2",
+                                        "input": (
+                                            "biosphere3",
+                                            self.biosphere_dict[
+                                                (
+                                                    "Carbon dioxide, in air",
+                                                    "natural resource",
+                                                    "in air",
+                                                    "kilogram",
+                                                )
+                                            ],
+                                        ),
+                                    }
+                                )
 
                         dataset["exchanges"] = [v for v in dataset["exchanges"] if v]