Add CO2 emissions forcing for esm-hist

[penguian]

Closes #112

Tested using the u-dq819/112-esm-hist@342286 workflow. Workflow output is at gadi:/scratch/tm70/pcl851/cylc-run/u-dq819.112-esm-hist/run1. Ancillary output is at gadi:/g/data/tm70/pcl851/CMIP7/esm1p6_ancil/2025.12.23/modern/historical-emissions/atmosphere/forcing/global.N96/2025.12.23/CO2_fluxes_1849_2023_cmip7.anc. Please pay attention to the calendar and dates in the ancillary file.

[penguian]

Re-tested using the u-dq819/112-esm-hist@342289 workflow, which creates all forcings except ozone. Workflow output is at gadi:/scratch/tm70/pcl851/cylc-run/u-dq819.112-esm-hist/run3. Ancillary output is at gadi:/g/data/tm70/pcl851/CMIP7/esm1p6_ancil/2025.12.23:

/g/data/tm70/pcl851/CMIP7/esm1p6_ancil/2025.12.23
└── modern
    ├── historical
    │   └── atmosphere
    │       ├── aerosol
    │       │   └── global.N96
    │       │       └── 2025.12.23
    │       │           ├── BC_1849_2023_cmip7.anc
    │       │           ├── Bio_1849_2023_cmip7.anc
    │       │           ├── OCFF_1849_2023_cmip7.anc
    │       │           └── scycl_1849_2023_cmip7.anc
    │       ├── forcing
    │       │   └── resolution_independent
    │       │       └── 2025.12.23
    │       │           ├── TSI_CMIP7_ESM
    │       │           └── volcts_cmip7.dat
    │       └── land
    │           └── biogeochemistry
    │               └── global.N96
    │                   └── 2025.12.23
    │                       └── Ndep_1849_2023_cmip7.anc
    ├── historical-emissions
    │   └── atmosphere
    │       └── forcing
    │           └── global.N96
    │               └── 2025.12.23
    │                   └── CO2_fluxes_1849_2023_cmip7.anc
    └── pre-industrial
        └── atmosphere
            ├── aerosol
            │   └── global.N96
            │       └── 2025.12.23
            │           ├── BC_1850_cmip7.anc
            │           ├── Bio_1850_cmip7.anc
            │           ├── OCFF_1850_cmip7.anc
            │           └── scycl_1850_cmip7.anc
            └── land
                └── biogeochemistry
                    └── global.N96
                        └── 2025.12.23
                            └── Ndep_1850_cmip7.anc

27 directories, 13 files

[MartinDix]

Looks like this is missing the aircraft emissions. These need to be summed over levels and added to the surface emissions. Since the supplied data includes 1949 and 2023 might as well use that directly rather than copying 1850 and 2022.

[penguian]

Tested using the u-dq819/112-esm-hist@342289 workflow, which creates all forcings except ozone. Workflow output is at gadi:/scratch/tm70/pcl851/cylc-run/u-dq819.112-esm-hist/run9. Ancillary output is at gadi:/g/data/tm70/pcl851/CMIP7/esm1p6_ancil/2025.12.25:

/g/data/tm70/pcl851/CMIP7/esm1p6_ancil/2025.12.25
└── modern
    ├── historical
    │   └── atmosphere
    │       ├── aerosol
    │       │   └── global.N96
    │       │       └── 2025.12.25
    │       │           ├── BC_1849_2023_cmip7.anc
    │       │           ├── Bio_1849_2023_cmip7.anc
    │       │           ├── OCFF_1849_2023_cmip7.anc
    │       │           └── scycl_1849_2023_cmip7.anc
    │       ├── forcing
    │       │   └── resolution_independent
    │       │       └── 2025.12.25
    │       │           ├── TSI_CMIP7_ESM
    │       │           └── volcts_cmip7.dat
    │       └── land
    │           └── biogeochemistry
    │               └── global.N96
    │                   └── 2025.12.25
    │                       └── Ndep_1849_2023_cmip7.anc
    ├── historical-emissions
    │   └── atmosphere
    │       └── forcing
    │           └── global.N96
    │               └── 2025.12.25
    │                   └── CO2_fluxes_1849_2023_cmip7.anc
    └── pre-industrial
        └── atmosphere
            ├── aerosol
            │   └── global.N96
            │       └── 2025.12.25
            │           ├── BC_1850_cmip7.anc
            │           ├── Bio_1850_cmip7.anc
            │           ├── OCFF_1850_cmip7.anc
            │           └── scycl_1850_cmip7.anc
            └── land
                └── biogeochemistry
                    └── global.N96
                        └── 2025.12.25
                            └── Ndep_1850_cmip7.anc

27 directories, 13 files

[penguian]

The current input path for CO2 emissions is
gadi:/g/data/qv56/replicas/input4MIPs/CMIP7/CMIP/PNNL-JGCRI/CEDS-CMIP-2025-04-18/atmos/mon/CO2_em_AIR_anthro/gn/v20250421/CO2-em-AIR-anthro_input4MIPs_emissions_CMIP_CEDS-CMIP-2025-04-18_gn_{date_range}.nc
for aircraft emissions and
gadi:/g/data/qv56/replicas/input4MIPs/CMIP7/CMIP/PNNL-JGCRI/CEDS-CMIP-2025-04-18/atmos/mon/CO2_em_anthro/gn/v20250421/CO2-em-anthro_input4MIPs_emissions_CMIP_CEDS-CMIP-2025-04-18_gn_{date_range}.nc
for the remaining CEDS emissions.

See https://input4mips-cvs.readthedocs.io/en/latest/dataset-overviews/anthropogenic-slcf-co2-emissions/
and see CMIP7/esm1p6/atmosphere/aerosol/cmip7_aerosol_anthro.py for the file path construction.

[blimlim]

I've attempted to reproduce @MartinDix's calculations with the updated outputs. It looks like the totals now match well: