Use new "steam_load_lbs" column name from pudl (#437) (#438)

* update steam column name

* update version

Author: grgmiller

CITATION.cff

@@ -22,6 +22,6 @@ authors:
 identifiers:
   - type: doi
     value: 'https://zenodo.org/doi/10.5281/zenodo.7062459'
-version: 0.7.0
+version: 0.7.1
 license: MIT
-date-released: '2025-12-23'
+date-released: '2026-01-23'

LICENSE.md

@@ -1,6 +1,6 @@
 MIT License
 
-Copyright (c) 2022-2025 Singularity Energy
+Copyright (c) 2022-2026 Singularity Energy
 
 Permission is hereby granted, free of charge, to any person obtaining a copy
 of this software and associated documentation files (the "Software"), to deal

notebooks/manual_data/manually_identify_crosswalk_updates.ipynb

@@ -65,7 +65,7 @@
     "        [\n",
     "            \"gross_generation_mwh\",\n",
     "            \"fuel_consumed_mmbtu\",\n",
-    "            \"steam_load_1000_lb\",\n",
+    "            \"steam_load_lbs\",\n",
     "            \"co2_mass_lb\",\n",
     "        ]\n",
     "    ]\n",

notebooks/manual_data/update_steam_units_to_remove.ipynb

@@ -65,7 +65,7 @@
     "        [\n",
     "            \"gross_generation_mwh\",\n",
     "            \"fuel_consumed_mmbtu\",\n",
-    "            \"steam_load_1000_lb\",\n",
+    "            \"steam_load_lbs\",\n",
     "            \"co2_mass_lb\",\n",
     "        ]\n",
     "    ]\n",

notebooks/work_in_progress/GH326_investigate_incomplete_timeseries.ipynb

@@ -283,7 +283,7 @@
     "    \"operating_datetime_utc\",\n",
     "    \"operating_time_hours\",\n",
     "    \"gross_load_mw\",\n",
-    "    \"steam_load_1000_lbs\",\n",
+    "    \"steam_load_lbs\",\n",
     "    \"co2_mass_tons\",\n",
     "    \"co2_mass_measurement_code\",\n",
     "    \"nox_mass_lbs\",\n",

notebooks/work_in_progress/clean_cems_outliers.ipynb

@@ -207,7 +207,7 @@
     "# TODO: here we're assuming that all columns are bad if net gen is bad, and that all bad rows have bad net gen.\n",
     "dat_cols = [\n",
     "    \"gross_generation_mwh\",\n",
-    "    \"steam_load_1000_lb\",\n",
+    "    \"steam_load_lbs\",\n",
     "    \"fuel_consumed_mmbtu\",\n",
     "    \"co2_mass_lb\",\n",
     "    \"nox_mass_lb\",\n",
@@ -320,7 +320,7 @@
     "# Get per-plant, per-variable median and IQR\n",
     "numeric_cols = [\n",
     "    \"gross_generation_mwh\",\n",
-    "    \"steam_load_1000_lb\",\n",
+    "    \"steam_load_lbs\",\n",
     "    \"fuel_consumed_mmbtu\",\n",
     "    \"co2_mass_lb\",\n",
     "    \"ch4_mass_lb\",\n",

pyproject.toml

@@ -5,7 +5,7 @@ build-backend = "hatchling.build"
 
 [project]
 name = "oge"
-version = "0.7.0"
+version = "0.7.1"
 requires-python = ">=3.11,<3.12"
 readme = "README.md"
 authors = [

src/oge/column_checks.py

@@ -77,7 +77,7 @@
         "datetime_utc",
         "operating_time_hours",
         "gross_generation_mwh",
-        "steam_load_1000_lb",
+        "steam_load_lbs",
         "fuel_consumed_mmbtu",
         "co2_mass_lb",
         "nox_mass_lb",
@@ -100,7 +100,7 @@
         "gross_generation_mwh",
         "gtn_method",
         "net_generation_mwh",
-        # "steam_load_1000_lb",
+        # "steam_load_lbs",
         "fuel_consumed_mmbtu",
         "fuel_consumed_for_electricity_mmbtu",
         "co2_mass_lb",
@@ -486,7 +486,7 @@ def get_dtypes() -> dict:
         "so2_removal_efficiency_annual": "float64",
         "so2_removal_efficiency_at_full_load": "float64",
         "start_year": "Int32",
-        "steam_load_1000_lb": "float64",
+        "steam_load_lbs": "float64",
         "subplant_id": "Int32",
         "subplant_primary_fuel": "string",
         "subplant_primary_fuel_from_capacity_mw": "string",
@@ -737,7 +737,7 @@ def get_dtypes() -> dict:
         "startup_source_code_3": "string",
         "startup_source_code_4": "string",
         "state": "string",
-        "steam_load_1000_lbs": "float64",
+        "steam_load_lbs": "float64",
         "steam_plant_type_code": "Int32",
         "stoker_tech": "boolean",
         "street_address": "string",

src/oge/data_cleaning.py

@@ -1020,7 +1020,7 @@ def remove_unmapped_fuel(cems: pd.DataFrame, year: int) -> pd.DataFrame:
     # calculate annual totals by unit
     annual_cems = (
         cems.groupby(["plant_id_eia", "emissions_unit_id_epa"], dropna=False)[
-            ["gross_generation_mwh", "steam_load_1000_lb", "fuel_consumed_mmbtu"]
+            ["gross_generation_mwh", "steam_load_lbs", "fuel_consumed_mmbtu"]
         ]
         .sum()
         .reset_index()
@@ -1064,7 +1064,7 @@ def remove_unmapped_fuel(cems: pd.DataFrame, year: int) -> pd.DataFrame:
     fuel_only_unmapped = annual_cems[
         (annual_cems["generator_id"].isna())
         & (annual_cems["gross_generation_mwh"] == 0)
-        & (annual_cems["steam_load_1000_lb"] == 0)
+        & (annual_cems["steam_load_lbs"] == 0)
         & (annual_cems["fuel_consumed_mmbtu"] > 0)
     ]
     if len(fuel_only_unmapped) > 0:
@@ -1257,7 +1257,7 @@ def identify_and_remove_steam_only_units(cems: pd.DataFrame, year: int) -> pd.Da
     # calculate annual totals by unit
     annual_cems = (
         cems.groupby(["plant_id_eia", "emissions_unit_id_epa"], dropna=False)[
-            ["gross_generation_mwh", "steam_load_1000_lb", "fuel_consumed_mmbtu"]
+            ["gross_generation_mwh", "steam_load_lbs", "fuel_consumed_mmbtu"]
         ]
         .sum()
         .reset_index()
@@ -1284,7 +1284,7 @@ def identify_and_remove_steam_only_units(cems: pd.DataFrame, year: int) -> pd.Da
     steam_only_unmapped = annual_cems[
         (annual_cems["generator_id"].isna())
         & (annual_cems["gross_generation_mwh"] == 0)
-        & (annual_cems["steam_load_1000_lb"] > 0)
+        & (annual_cems["steam_load_lbs"] > 0)
         & (annual_cems["fuel_consumed_mmbtu"] > 0)
     ]
     if len(steam_only_unmapped) > 0:
@@ -1314,7 +1314,7 @@ def identify_and_remove_steam_only_units(cems: pd.DataFrame, year: int) -> pd.Da
     # are still not entirely certain how to interpret steam data in CEMS, so its
     # inclusion could result in potentially anomalous results
     mapped_steam = annual_cems[
-        (annual_cems["steam_load_1000_lb"] > 0) & (~annual_cems["generator_id"].isna())
+        (annual_cems["steam_load_lbs"] > 0) & (~annual_cems["generator_id"].isna())
     ]
     if len(mapped_steam) > 0:
         logger.warning(
@@ -2440,7 +2440,7 @@ def aggregate_cems_to_subplant(cems):
 
     cems_columns_to_aggregate = [
         "gross_generation_mwh",
-        # "steam_load_1000_lb",
+        # "steam_load_lbs",
         "fuel_consumed_mmbtu",
         "co2_mass_lb",
         "ch4_mass_lb",

src/oge/data_pipeline.py

@@ -231,7 +231,7 @@ def main(args):
     # NOTE: steam load columns may be used in the future
     cems.drop(
         columns=[
-            "steam_load_1000_lb",
+            "steam_load_lbs",
             "co2_mass_measurement_code",
             "nox_mass_measurement_code",
             "so2_mass_measurement_code",