Changes
Page history
Update Hands on ICON data analysis for beginners
authored
Jan 29, 2025
by
Dian Putrasahan
Show whitespace changes
Inline
Side-by-side
Hands-on-ICON-data-analysis-for-beginners.md
View page @
f7bfc7b3
...
@@ -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
...
...
...
...