Polytope vertical profile

The following example first retrieves vertical profiles on pressure levels using polytope, then converts the resulting CoverageJSON data to Xarray.

import earthkit.data as ekd

request = {
    "class": "od",
    "stream": "enfo",
    "type": "pf",
    "date": -1,
    "time": "0000",
    "levtype": "pl",
    "expver": 1,
    "domain": "g",
    "param": "203/133",
    "number": "1",
    "step": "0",
    "levelist": "1/to/1000",
    "feature": {
        "type": "verticalprofile",
        "points": [[38.9, -9.1]],
    },
}

ds = ekd.from_source("polytope", "ecmwf-mars", request, stream=False, address='polytope.ecmwf.int')
ds.to_xarray()

2024-12-18 09:10:29 - INFO - Key read from /Users/cgr/.polytopeapirc
2024-12-18 09:10:29 - INFO - Sending request...
{'request': 'class: od\n'
            'date: -1\n'
            'domain: g\n'
            'expver: 1\n'
            'feature:\n'
            '  points:\n'
            '  - - 38.9\n'
            '    - -9.1\n'
            '  type: verticalprofile\n'
            'levelist: 1/to/1000\n'
            'levtype: pl\n'
            "number: '1'\n"
            'param: 203/133\n'
            "step: '0'\n"
            'stream: enfo\n'
            "time: '0000'\n"
            'type: pf\n',
 'verb': 'retrieve'}
2024-12-18 09:10:29 - INFO - Polytope user key found in session cache for user cgr
2024-12-18 09:10:30 - INFO - Request accepted. Please poll ./76f67d4d-f386-43ff-b11c-39e6ea9c404f for status
2024-12-18 09:10:30 - INFO - Polytope user key found in session cache for user cgr
2024-12-18 09:10:30 - INFO - Checking request status (./76f67d4d-f386-43ff-b11c-39e6ea9c404f)...
2024-12-18 09:10:30 - INFO - The current status of the request is 'queued'
2024-12-18 09:10:30 - INFO - The current status of the request is 'processing'
2024-12-18 09:10:33 - INFO - The current status of the request is 'processed'

<xarray.Dataset> Size: 640B
Dimensions:   (x: 1, y: 1, number: 1, datetime: 1, t: 1, z: 22)
Coordinates:
  * x         (x) float64 8B 38.91
  * y         (y) float64 8B 350.9
  * number    (number) int64 8B 1
  * datetime  (datetime) <U20 80B '2024-12-17T00:00:00Z'
  * t         (t) int64 8B 0
  * z         (z) int64 176B 1 2 3 5 7 10 20 30 ... 400 500 600 700 850 925 1000
Data variables:
    q         (x, y, number, datetime, t, z) float64 176B 4.254e-06 ... 0.005184
    o3        (x, y, number, datetime, t, z) float64 176B 4.228e-06 ... 6.039...
Attributes:
    class:          od
    Forecast date:  2024-12-17T00:00:00Z
    domain:         g
    expver:         0001
    levelist:       1
    levtype:        pl
    number:         1
    stream:         enfo
    type:           pf

xarray.Dataset

• Dimensions:
  • x: 1
  • y: 1
  • number: 1
  • datetime: 1
  • t: 1
  • z: 22
• Coordinates: (6)
  • x
    (x)
    float64
    38.91

    array([38.910368])

  • y
    (y)
    float64
    350.9

    array([350.889193])

  • number
    (number)
    int64
    1

    array([1])

  • datetime
    (datetime)
    <U20
    '2024-12-17T00:00:00Z'

    array(['2024-12-17T00:00:00Z'], dtype='<U20')

  • t
    (t)
    int64
    0

    array([0])

  • z
    (z)
    int64
    1 2 3 5 7 ... 600 700 850 925 1000

    array([   1,    2,    3,    5,    7,   10,   20,   30,   50,   70,  100,  150,
            200,  250,  300,  400,  500,  600,  700,  850,  925, 1000])

• Data variables: (2)
  • q
    (x, y, number, datetime, t, z)
    float64
    4.254e-06 4.205e-06 ... 0.005184

    type :

    Parameter

    units :

    kg kg**-1

    long_name :

    q

    array([[[[[[4.25423787e-06, 4.20500510e-06, 4.04536877e-06,
                3.70361522e-06, 3.59550177e-06, 3.53840187e-06,
                3.36708308e-06, 3.05463527e-06, 2.89008494e-06,
                2.92804998e-06, 3.00258398e-06, 4.47779894e-06,
                1.73449516e-05, 3.82959843e-05, 5.53429127e-05,
                6.47187233e-04, 2.22325325e-03, 3.93819809e-03,
                5.02681732e-03, 7.14731216e-03, 3.43990326e-03,
                5.18369675e-03]]]]]])

  • o3
    (x, y, number, datetime, t, z)
    float64
    4.228e-06 6.801e-06 ... 6.039e-08

    type :

    Parameter

    units :

    kg kg**-1

    long_name :

    o3

    array([[[[[[4.22841276e-06, 6.80094581e-06, 9.99311214e-06,
                1.27267540e-05, 1.20746226e-05, 1.08049953e-05,
                8.61148692e-06, 7.40709220e-06, 3.24864413e-06,
                2.65042974e-06, 1.42876524e-06, 7.21018296e-07,
                1.00844773e-07, 8.01956048e-08, 8.11123755e-08,
                7.92133505e-08, 7.48914317e-08, 7.68486643e-08,
                7.79750914e-08, 7.98123949e-08, 7.46204520e-08,
                6.03884125e-08]]]]]])

• Indexes: (6)
  • x
    PandasIndex

    PandasIndex(Index([38.910368186756], dtype='float64', name='x'))

  • y
    PandasIndex

    PandasIndex(Index([350.889192886457], dtype='float64', name='y'))

  • number
    PandasIndex

    PandasIndex(Index([1], dtype='int64', name='number'))

  • datetime
    PandasIndex

    PandasIndex(Index(['2024-12-17T00:00:00Z'], dtype='object', name='datetime'))

  • t
    PandasIndex

    PandasIndex(Index([0], dtype='int64', name='t'))

  • z
    PandasIndex

    PandasIndex(Index([   1,    2,    3,    5,    7,   10,   20,   30,   50,   70,  100,  150,
            200,  250,  300,  400,  500,  600,  700,  850,  925, 1000],
          dtype='int64', name='z'))

• Attributes: (9)

  class :

  od

  Forecast date :

  2024-12-17T00:00:00Z

  domain :

  g

  expver :

  0001

  levelist :

  1

  levtype :

  pl

  number :

  1

  stream :

  enfo

  type :

  pf