Extract the dual values for the E8_AnnualEmissionsLimit constraint from the CPLEX solution of an OSeMOSYS model

get_dual_osemosys.py

import pandas as pd
import sys
from logging import getLogger

logger = getLogger(__name__)


def main(file_path: str) -> pd.DataFrame:
    """Read CPLEX solution file and return dual values for E8_AnnualEmissionsLimit constraint.
    """
    df = pd.read_xml(file_path, xpath=".//constraint", parser="etree")
    logger.info(f"Reading CPLEX solution file from {file_path}")
    df[["Variable", "Index"]] = df["name"].str.split("(", expand=True)
    df["Index"] = df["Index"].str.replace(")", "", regex=False)
    df[['Variable', 'Index', 'index', 'status', 'slack', 'dual']][df['Variable'] == 'E8_AnnualEmissionsLimit']
    return df


if __name__ == "__main__":

    args = sys.argv[1:]
    if len(args) != 1:
        logger.error("Please provide the path to the CPLEX solution file e.g. python get_cplex_dual.py my_solution_file.sol")
        sys.exit(1)

    file_path = args[0]

    df = main(file_path)
    logger.info("Writing CPLEX solution file to E8_AnnualEmissionsLimit.csv")
    df.to_csv('E8_AnnualEmissionsLimit.csv', index=False)