Hands-on-ICON-data-analysis-for-beginners.md

@@ -1497,12 +1497,25 @@ import matplotlib.pyplot as plt
 import numpy as np
 import pyicon as pyic
 
-#Example from CMIP7 experiment with r2b4atm, r2b6oce
-outdir='/work/mh0287/b364148/Icon/Git_lev/icon.nwp.hdext.intno/experiments/tvp1774_001/'
-dso=xr.open_mfdataset(outdir+'tvp1774_001_oce_def_1[8-9]*.nc')
+##Example from CMIP7 experiment with r2b4atm, r2b6oce
+#outdir='/work/mh0287/b364148/Icon/Git_lev/icon.nwp.hdext.intno/experiments/tvp1774_001/'
+#dso=xr.open_mfdataset(outdir+'tvp1774_001_oce_def_1[8-9]*.nc')
+
+#Example with raw output
+#Example from CMIP7 experiment with r2b5atm, r2b7oce
+expid='slo1826'
+outdir='/work/mh0287/m211032/Icon/Git_lev/icon.XPP.20240717/build.intel-hdext/experiments/'+expid+'/'
+dso=xr.open_mfdataset(outdir+expid+'_oce_def_1[8-9]*.nc')
 
 #Plot mixed layer depth
 dso.mlotst10.isel(time=6).pyic.plot()
+
+#Example with catalogue
+import intake
+cat = intake.open_catalog("https://catalog.gitlab-pages.dkrz.de/ruby/main.yaml")
+ds = cat.dragon.hel24312_r5b7_ctrl.oce_qps.to_dask()
+ds.mlotst.isel(time=6).pyic.plot()
+
 ```
 <img src="https://pad.gwdg.de/uploads/b4238dc7-70de-4c5e-8788-cab007755f8b.png" alt="10m wind speed" width="450" />