| ... | ... | @@ -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" />
|
|
|
|
|
| ... | ... | |
