import earthkit.data
import earthkit.maps
data = earthkit.data.from_source(
'cds',
'reanalysis-era5-single-levels',
{
'product_type': 'reanalysis',
'variable': [
'2m_temperature',
'mean_sea_level_pressure',
],
'year': '1985',
'month': '12',
'day': '25',
'time': '12:00',
},
)
chart = earthkit.maps.Chart(domain="North Atlantic")
chart.plot(data[0], units="celsius")
chart.plot(data[1], units="hPa")
chart.coastlines(resolution="medium")
chart.borders()
chart.gridlines()
chart.subplot_titles()
chart.legend()
chart.show()