ERA5 initial conditions and forcings
Hi,
Quick follow-up on initial conditions:
Are there plans (or could we get) ic_*.nc files for all years 1979–2022, or just more samples (e.g. one per decade or Jan 1 each year)?
Any update on extending forcings and ICs to 2023–2025 with recent ERA5 data?
If extracting from the GCS zarr bucket is the intended path, a small code snippet showing how to pull & save a 6-hourly snapshot as ic_YYYY.nc would save a lot of time.
Also is this the code to generate the forcing and initial conditions? https://github.com/ai2cm/ace/tree/main/scripts/era5/pipeline
I have some issues in run it locally, such as co2 gs is not available.
Thanks for the awesome work!
Hi,
Thanks for letting us know your interest. We don't have immediate plans to make a longer length of forcing available for this checkpoint, or provide additional ICs. As you note you can make ICs from the data available on the GCS bucket documented here. Something like this should work:
import xarray as xr
ds = xr.open_zarr('gs://ai2cm-public-requester-pays/2024-11-13-ai2-climate-emulator-v2-amip/data/era5-1deg-1940-2022.zarr')
ds.sel(time=['1990-01-01T00:00:00', '1990-02-01T00:00:00', #other IC dates here]).to_netcdf('ic_1990.nc')
The data processing script linked in the post is approximately how these files were generated, but is not part of the public API of the fme package (this is true for anything in scripts) and so is best used for guidance in writing your own pipelines rather than executing directly. The specific CO2 data you mention is available via public sources (eg NOAA GML).
Thanks.