Update Hands on ICON data analysis for beginners authored by Dian Putrasahan's avatar Dian Putrasahan
...@@ -1497,17 +1497,10 @@ import matplotlib.pyplot as plt ...@@ -1497,17 +1497,10 @@ import matplotlib.pyplot as plt
import numpy as np import numpy as np
import pyicon as pyic import pyicon as pyic
##Example from CMIP7 experiment with r2b4atm, r2b6oce #Example with raw output from CMIP7 experiment with r2b5atm, r2b7oce
#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' expid='slo1826'
outdir='/work/mh0287/m211032/Icon/Git_lev/icon.XPP.20240717/build.intel-hdext/experiments/'+expid+'/' 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') dso=xr.open_mfdataset(outdir+expid+'_oce_def_1[8-9]*.nc')
#Plot mixed layer depth
dso.mlotst10.isel(time=6).pyic.plot() dso.mlotst10.isel(time=6).pyic.plot()
#Example with catalogue #Example with catalogue
...@@ -1517,7 +1510,7 @@ ds = cat.dragon.hel24312_r5b7_ctrl.oce_qps.to_dask() ...@@ -1517,7 +1510,7 @@ ds = cat.dragon.hel24312_r5b7_ctrl.oce_qps.to_dask()
ds.mlotst.isel(time=6).pyic.plot() 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" /> <img src="uploads/d15c9520e8a37b9c94b879786ff45006/image.png" alt="mixed layer depth" width="450" />
```python ```python
... ...
......