Open CommonClimate opened 2 years ago
Update: Dan Lunt called the AR6 Fig 3.2 to my attention:
. This would be a great figure to reproduce in this workflow (in general terms). To simplify, we could focus on LMR2.1 as the background , and put up whatever PMIP4 GMST we have (e.g. MPI past2k). An interesting application would be to show how to plot not just GMST, but e.g. JJA temperature at PAGES 2k sites, showcasing xarray selection semantics (e.g. isel
).
This should be straightforward to do once LMR2.1 is cloud-hosted. In the interim, we can use the netCDF file for PAGES 2k (2019) that @fzhu2e generated.
For this plot, use gmt_MCruns_ensemble_full_LMRv2.1.nc
For the gridded fields, you may re-use the VICS dashboard code and (if it's easy) enable sliders for year selection
To easily plot the LMR/PAGES2k ensemble: https://fzhu2e.github.io/pens/ or pyleoclim.
@CommonClimate new notebook (with wqrts and all): https://github.com/LinkedEarth/PaleoBooks/blob/LMR-CMIP6-eVolv2k/notebooks/CMIP6_LMR.ipynb
Since it's in the proposal and we have all the pieces in hand, let's have a notebook (or section of a notebook) that facilitates a comparison of simulated and reconstructed temperature fields in the spirit of [IPCC AR5 Fig 5.9]( .
Master word: KEEP IT SIMPLE. Just use the simulations you have vs LMR, and show how to map time slice averages, like Little Ice Age - Medieval Climate Anomaly (so 1550:1850 - 950:1250 AD; definitions vary and we could not care less, the point is to show how to easily plot time slice comparisons of PMIP4 output vs reconstructions.) That's it. No fancy science, as the science behind this particular comparison is quite lame, but the general illustration of using xarray's isel + plot capabilities for mapping will be useful enough to get re-used a bunch, I predict.
@jordanplanders perhaps this would be a good time to draw up a checklist of what notebooks you have in this general area, and their % completion?
I say we can close this once the CMIP6-LMR repo is published as a JupyterBook
Hi @jordanplanders ! You will have the honor of inaugurating this repo with a brand new tutorial: Model-data comparison over the Common Era (see the "Climate of the Past Millennium" section of the proposal), leveraging recent additions to the Cloud.
Here is what @khider and I thought would be good to try:
MPI-ESM1-2-LR.past2k.r1i1p1f1
past2k experiment surface air temperature data, now available in the cloud (thanks to our awesome pangeo-forge colleagues). It is available here, and I have personally tested the following code chunk:dataset_public_url = 'https://ncsa.osn.xsede.org/Pangeo/pangeo-forge-test/prod/recipe-run-820/pangeo-forge/cmip6-feedstock/CMIP6.PMIP.MPI-M.MPI-ESM1-2-LR.past2k.r1i1p1f1.Amon.tas.gn.dummy.zarr' mapper = fsspec.get_mapper(dataset_public_url) ds = xr.open_zarr(mapper, consolidated=True) ds ds['tas'].mean(axis=0).plot()