Process CO2 emission forcings for ScenarioMIP

[paul-leopardi]

Closes #140

For each scenario, the code in CMIP7/esm1p6/atmosphere/co2/cmip7_ES_CO2_interpolate.py reads both the ScenarioMIP CO2 emission forcing data for CO2 and the ScenarioMIP aircraft CO2 emission forcing data, adds the two, then interpolates the time series and regrids the data to produce an ancillary file.

It turns out that the equivalent code in
CMIP7/esm1p6/atmosphere/co2/cmip7_EH_CO2_interpolate.py for the historical forcings interpolated monthly data from the cubes corresponding the the historical forcing data before collapsing the CO2 emission data by sector and the CO2 aircraft emission data by altitude. When this was tried for ScenarioMIP the task ran out of memory when interpolating the time series. The solution applied to both the EH and ES code is to collapse the respective cubes before interpolating the time series. That is why four files are changed.

The debugging work and a small amount of coding was assisted by Google Antigravity using Gemini 3.6 Flash (high).

[read-the-docs-community]

Documentation build overview

📚 CMIP7-Input | 🛠️ Build #33954817 | 📁 Comparing db9711c against latest (6217a09)

  🔍 Preview build  

1 file changed

± index.html

[paul-leopardi]

I created and ran the u-dq819.scenarios-co2 workflow. Logs are in /scratch/tm70/pcl851/cylc-run/u-dq819.scenarios-co2/run9/log

Output is in /g/data/tm70/pcl851/CMIP7/esm1p6_ancil/2026.08.07
In particular, the esm-scen7-h to esm-scen7-vl output is in

/g/data/tm70/pcl851/CMIP7/esm1p6_ancil/2026.08.07/esm-scenarios/h/atmosphere/forcing/global.N96/2026.08.07/CO2_fluxes_h_2022_2100_cmip7.anc
/g/data/tm70/pcl851/CMIP7/esm1p6_ancil/2026.08.07/esm-scenarios/hl/atmosphere/forcing/global.N96/2026.08.07/CO2_fluxes_hl_2022_2100_cmip7.anc
/g/data/tm70/pcl851/CMIP7/esm1p6_ancil/2026.08.07/esm-scenarios/m/atmosphere/forcing/global.N96/2026.08.07/CO2_fluxes_m_2022_2100_cmip7.anc
/g/data/tm70/pcl851/CMIP7/esm1p6_ancil/2026.08.07/esm-scenarios/vl/atmosphere/forcing/global.N96/2026.08.07/CO2_fluxes_vl_2022_2100_cmip7.anc

The directory name esm-scenarios is arbitrary and can be changed,

[paul-leopardi]

See the comment in #140 as well as https://input4mips-cvs.readthedocs.io/en/latest/dataset-overviews/anthropogenic-slcf-co2-emissions/#scenariomip_2. The workflow needs to be re-run using 1-1-2 rather than 1-1-1 datasets.

[paul-leopardi]

The 1-1-2 datasets are not yet available at /g/data/qv56/replicas/input4MIPs/CMIP7/ScenarioMIP/IIASA-IAMC
I have submitted https://track.nci.org.au/servicedesk/customer/portal/5/HELP-210008 to ask for these to be loaded.