Operational Full Fields Extraction Example¶
This example demonstrates how to extract a Full Field from the ECMWF Operational Forecast.
In [ ]:
Copied!
import earthkit.data
request = {
"class": "od",
"stream" : "enfo",
"type" : "pf",
"date" : -1, # Note: date must be within the last two days
"time" : "0000",
"levtype" : "sfc",
"expver" : "0001",
"domain" : "g",
"param" : "164/167/169",
"number" : "1",
"step": "0",
}
ds = earthkit.data.from_source(
"polytope",
"ecmwf-mars",
request,
stream=False,
address='polytope.ecmwf.int',
)
import earthkit.data
request = {
"class": "od",
"stream" : "enfo",
"type" : "pf",
"date" : -1, # Note: date must be within the last two days
"time" : "0000",
"levtype" : "sfc",
"expver" : "0001",
"domain" : "g",
"param" : "164/167/169",
"number" : "1",
"step": "0",
}
ds = earthkit.data.from_source(
"polytope",
"ecmwf-mars",
request,
stream=False,
address='polytope.ecmwf.int',
)
In [2]:
Copied!
ds.ls()
ds.ls()
Out[2]:
| centre | shortName | typeOfLevel | level | dataDate | dataTime | stepRange | dataType | number | gridType | |
|---|---|---|---|---|---|---|---|---|---|---|
| 0 | ecmf | tcc | surface | 0 | 20251204 | 0 | 0 | pf | 1 | reduced_gg |
| 1 | ecmf | 2t | surface | 0 | 20251204 | 0 | 0 | pf | 1 | reduced_gg |
| 2 | ecmf | ssrd | surface | 0 | 20251204 | 0 | 0 | pf | 1 | reduced_gg |
Visualise¶
The following visualisation uses earthkit-plots.
In [ ]:
Copied!
import earthkit.plots
earthkit.plots.quickplot(ds[1]).show()
import earthkit.plots
earthkit.plots.quickplot(ds[1]).show()
867262f33dd699f318e74cfb5938737afb74a5311b086c2621ba48f4626637eb.npz: 0%| | 0.00/115M [00:00<?, ?B/…
In [ ]:
Copied!