Climate dt example

Climate DT Full Fields Extraction Example¶

This example demonstrates how to extract a Full Field from the DestinE Climate DT.

In [ ]:

Copied!





import earthkit.data

request = {
    'activity': 'ScenarioMIP',
    'class': 'd1',
    'dataset': 'climate-dt',
    'date': '20200102',
    'experiment': 'SSP3-7.0',
    'expver': '0001',
    'generation': '1',
    'levtype': 'sfc',
    'model': 'IFS-NEMO',
    'param': '134/165/166',
    'realization': '1',
    'resolution': 'standard',
    'stream': 'clte',
    'time': '0100', 
    'type': 'fc'
}


ds = earthkit.data.from_source("polytope", "destination-earth", request, address="polytope.lumi.apps.dte.destination-earth.eu", stream=False)
import earthkit.data

request = {
    'activity': 'ScenarioMIP',
    'class': 'd1',
    'dataset': 'climate-dt',
    'date': '20200102',
    'experiment': 'SSP3-7.0',
    'expver': '0001',
    'generation': '1',
    'levtype': 'sfc',
    'model': 'IFS-NEMO',
    'param': '134/165/166',
    'realization': '1',
    'resolution': 'standard',
    'stream': 'clte',
    'time': '0100', 
    'type': 'fc'
}


ds = earthkit.data.from_source("polytope", "destination-earth", request, address="polytope.lumi.apps.dte.destination-earth.eu", stream=False)

In [3]:

Copied!

ds.ls()
ds.ls()

Out[3]:

	centre	shortName	typeOfLevel	level	dataDate	dataTime	dataType	number	gridType
0	ecmf	sp	surface	0	20200102	100	fc	None	healpix
1	ecmf	10u	heightAboveGround	10	20200102	100	fc	None	healpix
2	ecmf	10v	heightAboveGround	10	20200102	100	fc	None	healpix

Visualise¶

The following visualisation uses earthkit-plots.

In [4]:

Copied!

import earthkit.plots

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

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

ce05761e06ddc639c16cad8e47ca64ded634218b6c6bca2bc72c14e472d38825.npz:   0%|          | 0.00/84.5M [00:00<?, ?B…

No description has been provided for this image

In [ ]: