Opendata example

Open Data Full Fields Extraction Example¶

This example demonstrates how to extract a Full Field from the ECMWF Open Data.

In [ ]:

Copied!





import earthkit.data

request = {
    "class": "ai",
    "stream" : "oper",
    "type" : "fc",
    "date" : -1,
    "time" : "0000",
    "levtype" : "sfc",
    "expver" : "0001", 
    "model": "aifs-single",
    "domain" : "g",
    "param" : "166/167/169", 
}

ds = earthkit.data.from_source(
    "polytope",
    "ecmwf-mars",
    request,
    stream=False,
    address='polytope.ecmwf.int',
)
import earthkit.data

request = {
    "class": "ai",
    "stream" : "oper",
    "type" : "fc",
    "date" : -1,
    "time" : "0000",
    "levtype" : "sfc",
    "expver" : "0001", 
    "model": "aifs-single",
    "domain" : "g",
    "param" : "166/167/169", 
}

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	dataType	number	gridType
0	ecmf	10v	heightAboveGround	10	20251204	fc	None	reduced_gg
1	ecmf	2t	heightAboveGround	2	20251204	fc	None	reduced_gg
2	ecmf	ssrd	surface	0	20251204	fc	None	reduced_gg

Visualise¶

The following visualisation uses earthkit-plots.

In [3]:

Copied!

import earthkit.plots

earthkit.plots.quickplot(ds[1]).show()
import earthkit.plots

earthkit.plots.quickplot(ds[1]).show()

8510b03559104792b1fb8a92ec2bdba1e061d1bebc0ed840c6cb3e0babf41015.npz:   0%|          | 0.00/76.2M [00:00<?, ?B…

No description has been provided for this image

In [ ]: