| ... | ... | @@ -1497,17 +1497,10 @@ 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 with raw output
|
|
|
|
#Example from CMIP7 experiment with r2b5atm, r2b7oce
|
|
|
|
#Example with raw output 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
|
| ... | ... | @@ -1517,7 +1510,7 @@ 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" />
|
|
|
|
<img src="uploads/d15c9520e8a37b9c94b879786ff45006/image.png" alt="mixed layer depth" width="450" />
|
|
|
|
|
|
|
|
|
|
|
|
```python
|
| ... | ... | |
