This is a little example on how to work with CMIP6 data and plot it the data was first prepared by using CDOs: ONI_all.sh this script reads in SST-data (tos*.nc), selects the NINO3.4 region, calculates anomalies to a centered base-period, calculates a 3month-running mean and averages it over the region the ONI is the outcome (this does not work in the example directory though, since not all the data is there) #---------------------------------------------------------------------------- following are two python scripts as examples on how to read in and plot the data these should also work within this example directory ONI_plot.py reads in a list of files containing ONI (these files were created by a ONI_all.sh, which calculates the ONI from SST-data (tos-files) with cdo) plots time series of ONI for each file (model) in: ./exampledata/*.nc out: ./exampledata/plots/.pdf #-------------------------------------------------------------------------------- ONI_trend.py reads in a list of files containing ONI (these files were created by a ONI_all.sh, which calculates the ONI from SST-data (tos-files) with cdo) plots ONI and 30year moving average for each model (window is easily adjustable) plots moving average and its mean for whole ensemble in: ./exampledata/*.nc out ./exampledata/plots/Trend*.pdf