MIR: regridding Xarray data¶

This example shows how to interpolate Xarray data with the mir backend.

To make this notebook work earthkit-data and earthkit-plots have to be installed.

[1]:

import earthkit.data as ekd
import earthkit.plots as ekp

import earthkit.geo as ekg

Please note the Xarray support in regridding is experimental. The grid spec of the input data can only inferred in some special cases, otherwise it has to be specified manually.

Xarray generated from GRIB with earthkit-data¶

In the first example we fetch GRIB data on 0.25°x0.25° grid on a subarea and convert to Xarray with earthkit-data.

[2]:

ds = ekd.from_source("sample", "era5_tquv_pl_subarea.grib").to_xarray()
ds

[2]:

<xarray.Dataset> Size: 69MB
Dimensions:                  (forecast_reference_time: 4, level: 3,
                              latitude: 301, longitude: 601)
Coordinates:
  * forecast_reference_time  (forecast_reference_time) datetime64[ns] 32B 201...
  * level                    (level) int64 24B 700 850 925
  * latitude                 (latitude) float64 2kB 90.0 89.75 ... 15.25 15.0
  * longitude                (longitude) float64 5kB -100.0 -99.75 ... 50.0
Data variables:
    q                        (forecast_reference_time, level, latitude, longitude) float64 17MB ...
    t                        (forecast_reference_time, level, latitude, longitude) float64 17MB ...
    u                        (forecast_reference_time, level, latitude, longitude) float64 17MB ...
    v                        (forecast_reference_time, level, latitude, longitude) float64 17MB ...
Attributes:
    Conventions:  CF-1.8
    institution:  ECMWF

xarray.Dataset

Dimensions:
- forecast_reference_time: 4
- level: 3
- latitude: 301
- longitude: 601
Coordinates: (4)
- forecast_reference_time
  (forecast_reference_time)
  datetime64[ns]
  2016-09-26 ... 2016-09-26T18:00:00
  standard_name :
  forecast_reference_time
  long_name :
  initial time of forecast
```
array(['2016-09-26T00:00:00.000000000', '2016-09-26T06:00:00.000000000',
       '2016-09-26T12:00:00.000000000', '2016-09-26T18:00:00.000000000'],
      dtype='datetime64[ns]')
```
- level
  (level)
  int64
  700 850 925
  standard_name :
  air_pressure
  long_name :
  pressure
  units :
  hectopascal
  positive :
  down
```
array([700, 850, 925])
```
- latitude
  (latitude)
  float64
  90.0 89.75 89.5 ... 15.5 15.25 15.0
  units :
  degrees_north
  standard_name :
  latitude
  long_name :
  latitude
```
array([90.  , 89.75, 89.5 , ..., 15.5 , 15.25, 15.  ], shape=(301,))
```
- longitude
  (longitude)
  float64
  -100.0 -99.75 -99.5 ... 49.75 50.0
  units :
  degrees_east
  standard_name :
  longitude
  long_name :
  longitude
```
array([-100.  ,  -99.75,  -99.5 , ...,   49.5 ,   49.75,   50.  ], shape=(601,))
```
Data variables: (4)
- q
  (forecast_reference_time, level, latitude, longitude)
  float64
  ...
  standard_name :
  specific_humidity
  long_name :
  Specific humidity
  units :
  dimensionless
  level_type :
  pressure
  _earthkit :
  {'message': b'GRIB\x00\x00l\x01\x00\x004\x80b\x91\xff\x80\x85d\x02\xbc\x10\t\x1a\x00\x00\x01\x00\x00\x00\x00\x00\x00\x15\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x01\x17\x02\x04\x010001\x00\x00\x00\x00\x00 \x00\xff\x00\x02Y\x01-\x01_\x90\x81\x86\xa0\x80\x00:\x98\x00\xc3P\x00\xfa\x00\xfa\x00\x00\x00\x00\x00\x00\x00\x0c\x00\x80\x02D\xb9}n\x00\x007777', 'bitsPerValue': 16, 'grid_spec': '{"area": [90, -100, 15, 50], "grid": [0.25, 0.25]}'}
```
[2170812 values with dtype=float64]
```
- t
  (forecast_reference_time, level, latitude, longitude)
  float64
  ...
  standard_name :
  air_temperature
  long_name :
  Temperature
  units :
  kelvin
  level_type :
  pressure
  _earthkit :
  {'message': b'GRIB\x00\x00l\x01\x00\x004\x80b\x91\xff\x80\x82d\x02\xbc\x10\t\x1a\x00\x00\x01\x00\x00\x00\x00\x00\x00\x15\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x01\x17\x02\x04\x010001\x00\x00\x00\x00\x00 \x00\xff\x00\x02Y\x01-\x01_\x90\x81\x86\xa0\x80\x00:\x98\x00\xc3P\x00\xfa\x00\xfa\x00\x00\x00\x00\x00\x00\x00\x0c\x00\x80\x02D\xb9}n\x00\x007777', 'bitsPerValue': 16, 'grid_spec': '{"area": [90, -100, 15, 50], "grid": [0.25, 0.25]}'}
```
[2170812 values with dtype=float64]
```
- u
  (forecast_reference_time, level, latitude, longitude)
  float64
  ...
  standard_name :
  eastward_wind
  long_name :
  U component of wind
  units :
  meter / second
  level_type :
  pressure
  _earthkit :
  {'message': b'GRIB\x00\x00l\x01\x00\x004\x80b\x91\xff\x80\x83d\x02\xbc\x10\t\x1a\x00\x00\x01\x00\x00\x00\x00\x00\x00\x15\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x01\x17\x02\x04\x010001\x00\x00\x00\x00\x00 \x00\xff\x00\x02Y\x01-\x01_\x90\x81\x86\xa0\x80\x00:\x98\x00\xc3P\x00\xfa\x00\xfa\x00\x00\x00\x00\x00\x00\x00\x0c\x00\x80\x02D\xb9}n\x00\x007777', 'bitsPerValue': 16, 'grid_spec': '{"area": [90, -100, 15, 50], "grid": [0.25, 0.25]}'}
```
[2170812 values with dtype=float64]
```
- v
  (forecast_reference_time, level, latitude, longitude)
  float64
  ...
  standard_name :
  northward_wind
  long_name :
  V component of wind
  units :
  meter / second
  level_type :
  pressure
  _earthkit :
  {'message': b'GRIB\x00\x00l\x01\x00\x004\x80b\x91\xff\x80\x84d\x02\xbc\x10\t\x1a\x00\x00\x01\x00\x00\x00\x00\x00\x00\x15\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x01\x17\x02\x04\x010001\x00\x00\x00\x00\x00 \x00\xff\x00\x02Y\x01-\x01_\x90\x81\x86\xa0\x80\x00:\x98\x00\xc3P\x00\xfa\x00\xfa\x00\x00\x00\x00\x00\x00\x00\x0c\x00\x80\x02D\xb9}n\x00\x007777', 'bitsPerValue': 16, 'grid_spec': '{"area": [90, -100, 15, 50], "grid": [0.25, 0.25]}'}
```
[2170812 values with dtype=float64]
```
Attributes: (2)
Conventions :
CF-1.8
institution :
ECMWF

The data generated in this way has the grid spec that can be accessed via the earthkit accessor.

[3]:

ds.earthkit.grid_spec

[3]:

{'area': [90, -100, 15, 50], 'grid': [0.25, 0.25]}

Now, we interpolate the data to a 5°x5° regular latitude-longitude grid. Since the grid spec can be inferred from the dataset we only need to specify the out_grid.

[4]:

# the target grid is 5x5 degree regular latitude-longitude grid
out_grid = {"grid": [5, 5]}

# perform interpolation
ds_r = ekg.regrid(ds, out_grid=out_grid, interpolation="linear")
ds_r

[4]:

<xarray.Dataset> Size: 191kB
Dimensions:                  (forecast_reference_time: 4, level: 3,
                              latitude: 16, longitude: 31)
Coordinates:
  * forecast_reference_time  (forecast_reference_time) datetime64[ns] 32B 201...
  * level                    (level) int64 24B 700 850 925
  * latitude                 (latitude) float64 128B 90.0 85.0 ... 20.0 15.0
  * longitude                (longitude) float64 248B -100.0 -95.0 ... 45.0 50.0
Data variables:
    q                        (forecast_reference_time, level, latitude, longitude) float64 48kB ...
    t                        (forecast_reference_time, level, latitude, longitude) float64 48kB ...
    u                        (forecast_reference_time, level, latitude, longitude) float64 48kB ...
    v                        (forecast_reference_time, level, latitude, longitude) float64 48kB ...

xarray.Dataset

Dimensions:
- forecast_reference_time: 4
- level: 3
- latitude: 16
- longitude: 31

Coordinates: (4)

forecast_reference_time

(forecast_reference_time)

datetime64[ns]

2016-09-26 ... 2016-09-26T18:00:00

standard_name :: forecast_reference_time
long_name :: initial time of forecast

array(['2016-09-26T00:00:00.000000000', '2016-09-26T06:00:00.000000000',
       '2016-09-26T12:00:00.000000000', '2016-09-26T18:00:00.000000000'],
      dtype='datetime64[ns]')

level
(level)
int64
700 850 925
standard_name :
air_pressure
long_name :
pressure
units :
hectopascal
positive :
down
```
array([700, 850, 925])
```

latitude

(latitude)

float64

90.0 85.0 80.0 ... 25.0 20.0 15.0

array([90., 85., 80., 75., 70., 65., 60., 55., 50., 45., 40., 35., 30., 25.,
       20., 15.])

longitude

(longitude)

float64

-100.0 -95.0 -90.0 ... 45.0 50.0

array([-100.,  -95.,  -90.,  -85.,  -80.,  -75.,  -70.,  -65.,  -60.,  -55.,
        -50.,  -45.,  -40.,  -35.,  -30.,  -25.,  -20.,  -15.,  -10.,   -5.,
          0.,    5.,   10.,   15.,   20.,   25.,   30.,   35.,   40.,   45.,
         50.])

Data variables: (4)

(forecast_reference_time, level, latitude, longitude)

float64

0.001419 0.001419 ... 0.01652

standard_name :: specific_humidity
long_name :: Specific humidity
units :: dimensionless
level_type :: pressure
_earthkit :: {'message': b'GRIB\x00\x00l\x01\x00\x004\x80b\x91\xff\x80\x85d\x02\xbc\x10\t\x1a\x00\x00\x01\x00\x00\x00\x00\x00\x00\x15\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x01\x17\x02\x04\x010001\x00\x00\x00\x00\x00 \x00\xff\x00\x00\x1f\x00\x10\x01_\x90\x81\x86\xa0\x80\x00:\x98\x00\xc3P\x13\x88\x13\x88\x00\x00\x00\x00\x00\x00\x00\x0c\x08\x00\x00D\xb9}n\x00\x007777', 'bitsPerValue': 16, 'grid_spec': '{"area": [90, -100, 15, 50], "grid": [5, 5]}'}

array([[[[0.00141868, 0.00141868, 0.00141868, ..., 0.00141868,
          0.00141868, 0.00141868],
         [0.00160369, 0.00221714, 0.00290093, ..., 0.0014883 ,
          0.0010179 , 0.00080809],
         [0.00240263, 0.00299176, 0.0029026 , ..., 0.00040969,
          0.00030002, 0.0006102 ],
         ...,
         [0.00822291, 0.0071629 , 0.00646505, ..., 0.00387391,
          0.00369391, 0.00369176],
         [0.00992569, 0.00833186, 0.00830874, ..., 0.00578937,
          0.00371226, 0.00358185],
         [0.00789627, 0.00765714, 0.00932607, ..., 0.00655111,
          0.00445708, 0.00628623]],

        [[0.00278865, 0.00278865, 0.00278865, ..., 0.00278865,
          0.00278865, 0.00278865],
         [0.00323926, 0.00343333, 0.00348912, ..., 0.00332652,
          0.00324498, 0.00311671],
         [0.00353108, 0.00337802, 0.00405847, ..., 0.00079547,
          0.00083886, 0.00097524],
...
         [0.01244843, 0.01115573, 0.0091306 , ..., 0.00449288,
          0.00361073, 0.00485766],
         [0.0095197 , 0.0127579 , 0.01093018, ..., 0.0069958 ,
          0.0044781 , 0.00433505],
         [0.01305926, 0.01232684, 0.01200307, ..., 0.0069176 ,
          0.00496161, 0.00657571]],

        [[0.00274271, 0.00274271, 0.00274271, ..., 0.00274271,
          0.00274271, 0.00274271],
         [0.00219292, 0.00226731, 0.00246376, ..., 0.00206465,
          0.00294012, 0.00316138],
         [0.00174851, 0.00198836, 0.00259346, ..., 0.00390715,
          0.00390572, 0.00262303],
         ...,
         [0.01386971, 0.01455206, 0.01341433, ..., 0.00455231,
          0.00389714, 0.00399727],
         [0.00951952, 0.01555867, 0.0159697 , ..., 0.00955672,
          0.00462622, 0.00438303],
         [0.01565069, 0.01568121, 0.01256413, ..., 0.00855727,
          0.00496144, 0.01651806]]]], shape=(4, 3, 16, 31))

(forecast_reference_time, level, latitude, longitude)

float64

263.4 263.4 263.4 ... 305.6 298.7

standard_name :: air_temperature
long_name :: Temperature
units :: kelvin
level_type :: pressure
_earthkit :: {'message': b'GRIB\x00\x00l\x01\x00\x004\x80b\x91\xff\x80\x82d\x02\xbc\x10\t\x1a\x00\x00\x01\x00\x00\x00\x00\x00\x00\x15\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x01\x17\x02\x04\x010001\x00\x00\x00\x00\x00 \x00\xff\x00\x00\x1f\x00\x10\x01_\x90\x81\x86\xa0\x80\x00:\x98\x00\xc3P\x13\x88\x13\x88\x00\x00\x00\x00\x00\x00\x00\x0c\x08\x00\x00D\xb9}n\x00\x007777', 'bitsPerValue': 16, 'grid_spec': '{"area": [90, -100, 15, 50], "grid": [5, 5]}'}

array([[[[263.43481445, 263.43481445, 263.43481445, ..., 263.43481445,
          263.43481445, 263.43481445],
         [264.56762695, 264.91333008, 265.9855957 , ..., 266.49926758,
          266.34301758, 266.11645508],
         [263.83911133, 265.7824707 , 265.96801758, ..., 265.59301758,
          265.16333008, 264.05004883],
         ...,
         [282.5793457 , 282.3137207 , 281.65356445, ..., 286.5090332 ,
          285.17114258, 283.65356445],
         [284.50512695, 283.71801758, 282.14575195, ..., 285.8293457 ,
          284.29418945, 284.85473633],
         [283.86450195, 282.99731445, 282.7824707 , ..., 285.09301758,
          285.05981445, 283.98168945]],

        [[268.47851562, 268.47851562, 268.47851562, ..., 268.47851562,
          268.47851562, 268.47851562],
         [270.27539062, 271.03515625, 270.22265625, ..., 272.47070312,
          270.78710938, 269.484375  ],
         [270.6875    , 269.85742188, 272.44921875, ..., 271.25976562,
          270.64453125, 269.64257812],
...
         [290.80004883, 290.69262695, 290.9387207 , ..., 301.86450195,
          302.64379883, 298.48754883],
         [293.27270508, 290.6496582 , 292.09692383, ..., 300.28051758,
          300.20239258, 299.3293457 ],
         [291.18481445, 290.95825195, 294.9543457 , ..., 299.12231445,
          300.65551758, 297.29614258]],

        [[268.93676758, 268.93676758, 268.93676758, ..., 268.93676758,
          268.93676758, 268.93676758],
         [265.62817383, 266.03637695, 267.01879883, ..., 272.14770508,
          271.26489258, 270.99926758],
         [264.04223633, 264.63012695, 268.02661133, ..., 273.11645508,
          272.91137695, 272.58129883],
         ...,
         [295.8918457 , 294.6496582 , 295.04418945, ..., 307.4230957 ,
          309.13598633, 303.80004883],
         [295.6965332 , 294.20629883, 296.63598633, ..., 304.75317383,
          304.77856445, 306.47387695],
         [294.39770508, 294.49926758, 300.10473633, ..., 303.40356445,
          305.59887695, 298.6965332 ]]]], shape=(4, 3, 16, 31))

(forecast_reference_time, level, latitude, longitude)

float64

-0.0001373 -0.0001373 ... 2.664

standard_name :: eastward_wind
long_name :: U component of wind
units :: meter / second
level_type :: pressure
_earthkit :: {'message': b'GRIB\x00\x00l\x01\x00\x004\x80b\x91\xff\x80\x83d\x02\xbc\x10\t\x1a\x00\x00\x01\x00\x00\x00\x00\x00\x00\x15\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x01\x17\x02\x04\x010001\x00\x00\x00\x00\x00 \x00\xff\x00\x00\x1f\x00\x10\x01_\x90\x81\x86\xa0\x80\x00:\x98\x00\xc3P\x13\x88\x13\x88\x00\x00\x00\x00\x00\x00\x00\x0c\x08\x00\x00D\xb9}n\x00\x007777', 'bitsPerValue': 16, 'grid_spec': '{"area": [90, -100, 15, 50], "grid": [5, 5]}'}

array([[[[-1.37329102e-04, -1.37329102e-04, -1.37329102e-04, ...,
          -1.37329102e-04, -1.37329102e-04, -1.37329102e-04],
         [ 3.36900330e+00,  2.78501892e+00,  2.81822205e+00, ...,
           1.11346283e+01,  1.10565033e+01,  1.05721283e+01],
         [-6.17324829e-01,  1.67759705e+00,  2.03501892e+00, ...,
           8.58580017e+00,  6.99595642e+00,  5.80064392e+00],
         ...,
         [-3.10365295e+00, -3.46498108e+00, -5.89981079e-01, ...,
           3.83384705e+00, -2.13294983e+00,  1.38658142e+00],
         [-2.97865295e+00, -4.11732483e+00, -3.78529358e+00, ...,
          -2.47279358e+00, -2.65762329e-01, -2.81387329e-01],
         [-2.54043579e-01, -2.48646545e+00, -1.41810608e+00, ...,
          -4.84974670e+00, -1.03380280e+01, -4.69935608e+00]],

        [[-8.54492188e-04, -8.54492188e-04, -8.54492188e-04, ...,
          -8.54492188e-04, -8.54492188e-04, -8.54492188e-04],
         [ 3.80578613e+00,  4.12219238e+00,  4.64562988e+00, ...,
           1.25968018e+01,  1.10850830e+01,  9.61242676e+00],
         [ 2.69641113e+00, -1.46765137e+00,  1.82727051e+00, ...,
           6.41320801e+00,  6.63781738e+00,  5.98352051e+00],
...
           5.23757935e+00, -4.34295654e-01,  2.77664185e+00],
         [-2.09249878e+00,  1.89578247e+00, -3.77413940e+00, ...,
          -1.50265503e+00, -6.53195190e+00, -3.74093628e+00],
         [-8.63983154e-01, -3.99093628e+00, -1.93429565e+00, ...,
          -1.79367065e+00,  3.44543457e-02,  4.68048096e-01]],

        [[-3.50952148e-04, -3.50952148e-04, -3.50952148e-04, ...,
          -3.50952148e-04, -3.50952148e-04, -3.50952148e-04],
         [ 7.12269592e+00,  6.34535217e+00,  5.58363342e+00, ...,
           5.85511780e+00,  6.53285217e+00,  8.56996155e+00],
         [ 4.30238342e+00, -1.01165771e-02,  2.22425842e+00, ...,
           7.90785217e+00,  1.01246490e+01,  7.91371155e+00],
         ...,
         [-9.76913452e-01, -4.87339783e+00, -1.76206970e+00, ...,
           2.45082092e+00,  1.98867798e-01,  6.40274048e-01],
         [-2.08042908e+00,  1.72816467e+00, -2.34996033e+00, ...,
          -5.27694702e-01, -3.46910095e+00, -5.68980408e+00],
         [-9.67147827e-01, -3.63121033e+00, -1.14878845e+00, ...,
           3.19300842e+00, -2.24960327e-01,  2.66371155e+00]]]],
      shape=(4, 3, 16, 31))

(forecast_reference_time, level, latitude, longitude)

float64

0.0004425 0.0004425 ... -0.142

standard_name :: northward_wind
long_name :: V component of wind
units :: meter / second
level_type :: pressure
_earthkit :: {'message': b'GRIB\x00\x00l\x01\x00\x004\x80b\x91\xff\x80\x84d\x02\xbc\x10\t\x1a\x00\x00\x01\x00\x00\x00\x00\x00\x00\x15\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x01\x17\x02\x04\x010001\x00\x00\x00\x00\x00 \x00\xff\x00\x00\x1f\x00\x10\x01_\x90\x81\x86\xa0\x80\x00:\x98\x00\xc3P\x13\x88\x13\x88\x00\x00\x00\x00\x00\x00\x00\x0c\x08\x00\x00D\xb9}n\x00\x007777', 'bitsPerValue': 16, 'grid_spec': '{"area": [90, -100, 15, 50], "grid": [5, 5]}'}

array([[[[ 4.42504883e-04,  4.42504883e-04,  4.42504883e-04, ...,
           4.42504883e-04,  4.42504883e-04,  4.42504883e-04],
         [ 1.10844269e+01,  9.42231750e+00,  7.25630188e+00, ...,
          -1.16440887e+01, -1.31616669e+01, -1.43003387e+01],
         [ 5.06880188e+00,  1.08109894e+01,  6.35200500e+00, ...,
          -6.08744812e+00, -1.05698700e+01, -1.26401825e+01],
         ...,
         [ 3.43794250e+00,  2.40278625e+00,  1.36372375e+00, ...,
           7.32661438e+00,  2.82075500e+00, -1.49292450e+01],
         [-4.72612000e+00, -1.61666870e-01,  3.96919250e+00, ...,
          -1.42729187e+00, -2.75346375e+00,  2.54348755e-01],
         [-2.07572937e+00,  4.77005005e-01,  2.26802063e+00, ...,
          -1.75151062e+00, -7.15190125e+00, -4.50932312e+00]],

        [[-5.95092773e-04, -5.95092773e-04, -5.95092773e-04, ...,
          -5.95092773e-04, -5.95092773e-04, -5.95092773e-04],
         [ 7.69081116e+00,  5.66151428e+00,  8.03456116e+00, ...,
          -1.01412201e+01, -1.01900482e+01, -9.92832947e+00],
         [-4.84970093e-01,  9.44471741e+00,  2.17909241e+00, ...,
          -5.09629822e+00, -5.92051697e+00, -6.77793884e+00],
...
          -7.10606384e+00, -1.73692322e+00, -1.30220795e+01],
         [-1.98497009e+00,  1.50331116e+00,  5.54042053e+00, ...,
           1.35803223e-03, -9.18563843e-01, -5.18614197e+00],
         [-2.59043884e+00,  6.50526428e+00, -3.37982178e-02, ...,
           5.13076782e-01, -1.59434509e+00,  2.46034241e+00]],

        [[ 6.25610352e-04,  6.25610352e-04,  6.25610352e-04, ...,
           6.25610352e-04,  6.25610352e-04,  6.25610352e-04],
         [ 8.26568604e-02, -3.54843140e-01, -1.89056396e-02, ...,
           2.50648499e+00,  1.25453186e+00, -2.01304626e+00],
         [-9.61851501e+00, -7.53843689e+00, -3.43882751e+00, ...,
          -3.78257751e+00, -9.88218689e+00, -1.29310150e+01],
         ...,
         [-9.79843140e-01, -6.77337646e-02,  2.52015686e+00, ...,
          -4.12437439e+00, -1.08335876e+00, -1.00599213e+01],
         [-2.00914001e+00,  2.70765686e+00,  9.63516235e-01, ...,
          -2.34898376e+00, -5.69686890e-01, -4.55210876e+00],
         [-8.70468140e-01,  2.86195374e+00, -7.94525146e-02, ...,
           2.38906860e-01, -1.83140564e+00, -1.41952515e-01]]]],
      shape=(4, 3, 16, 31))

The resulting dataset still has the earthkit accessor with the correct grid spec set on it.

[5]:

ds_r.earthkit.grid_spec

[5]:

{'area': [90, -100, 15, 50], 'grid': [5, 5]}

We check the result visually by plotting the same data from the input and output dataset.

[6]:

ekp.quickplot(ds.t.isel(forecast_reference_time=0, level=0))
ekp.quickplot(ds_r.t.isel(forecast_reference_time=0, level=0))

[6]:

<earthkit.plots.components.maps.Map at 0x17ec7d310>

../_images/tutorials_mir_regrid_xarray_14_1.png

../_images/tutorials_mir_regrid_xarray_14_2.png

Xarray generated from NetCDF¶

In this example the data is available on the same grid as in the previous example but the format is NetCDF. We fetch the file and load it into Xarray.

[7]:

ds = ekd.from_source("sample", "era5_tquv_pl_subarea.nc").to_xarray()
ds

[7]:

<xarray.Dataset> Size: 35MB
Dimensions:         (valid_time: 4, pressure_level: 3, latitude: 301,
                     longitude: 601)
Coordinates:
  * valid_time      (valid_time) datetime64[ns] 32B 2016-09-26 ... 2016-09-26...
  * pressure_level  (pressure_level) float64 24B 925.0 850.0 700.0
  * latitude        (latitude) float64 2kB 90.0 89.75 89.5 ... 15.5 15.25 15.0
  * longitude       (longitude) float64 5kB -100.0 -99.75 -99.5 ... 49.75 50.0
Data variables:
    q               (valid_time, pressure_level, latitude, longitude) float32 9MB ...
    t               (valid_time, pressure_level, latitude, longitude) float32 9MB ...
    u               (valid_time, pressure_level, latitude, longitude) float32 9MB ...
    v               (valid_time, pressure_level, latitude, longitude) float32 9MB ...
Attributes:
    GRIB_centre:             ecmf
    GRIB_centreDescription:  European Centre for Medium-Range Weather Forecasts
    GRIB_subCentre:          0
    Conventions:             CF-1.7
    institution:             European Centre for Medium-Range Weather Forecasts
    history:                 2026-03-12T10:30 GRIB to CDM+CF via cfgrib-0.9.1...

xarray.Dataset

Dimensions:
- valid_time: 4
- pressure_level: 3
- latitude: 301
- longitude: 601
Coordinates: (4)
- valid_time
  (valid_time)
  datetime64[ns]
  2016-09-26 ... 2016-09-26T18:00:00
  long_name :
  time
  standard_name :
  time
```
array(['2016-09-26T00:00:00.000000000', '2016-09-26T06:00:00.000000000',
       '2016-09-26T12:00:00.000000000', '2016-09-26T18:00:00.000000000'],
      dtype='datetime64[ns]')
```
- pressure_level
  (pressure_level)
  float64
  925.0 850.0 700.0
  long_name :
  pressure
  units :
  hPa
  positive :
  down
  stored_direction :
  decreasing
  standard_name :
  air_pressure
```
array([925., 850., 700.])
```
- latitude
  (latitude)
  float64
  90.0 89.75 89.5 ... 15.5 15.25 15.0
  units :
  degrees_north
  standard_name :
  latitude
  long_name :
  latitude
  stored_direction :
  decreasing
```
array([90.  , 89.75, 89.5 , ..., 15.5 , 15.25, 15.  ], shape=(301,))
```
- longitude
  (longitude)
  float64
  -100.0 -99.75 -99.5 ... 49.75 50.0
  units :
  degrees_east
  standard_name :
  longitude
  long_name :
  longitude
```
array([-100.  ,  -99.75,  -99.5 , ...,   49.5 ,   49.75,   50.  ], shape=(601,))
```
Data variables: (4)
- q
  (valid_time, pressure_level, latitude, longitude)
  float32
  ...
  GRIB_paramId :
  133
  GRIB_dataType :
  an
  GRIB_numberOfPoints :
  180901
  GRIB_typeOfLevel :
  isobaricInhPa
  GRIB_stepUnits :
  1
  GRIB_stepType :
  instant
  GRIB_gridType :
  regular_ll
  GRIB_uvRelativeToGrid :
  0
  GRIB_NV :
  0
  GRIB_Nx :
  601
  GRIB_Ny :
  301
  GRIB_cfName :
  specific_humidity
  GRIB_cfVarName :
  q
  GRIB_gridDefinitionDescription :
  Latitude/Longitude Grid
  GRIB_iDirectionIncrementInDegrees :
  0.25
  GRIB_iScansNegatively :
  0
  GRIB_jDirectionIncrementInDegrees :
  0.25
  GRIB_jPointsAreConsecutive :
  0
  GRIB_jScansPositively :
  0
  GRIB_latitudeOfFirstGridPointInDegrees :
  90.0
  GRIB_latitudeOfLastGridPointInDegrees :
  15.0
  GRIB_longitudeOfFirstGridPointInDegrees :
  -100.0
  GRIB_longitudeOfLastGridPointInDegrees :
  50.0
  GRIB_missingValue :
  3.4028234663852886e+38
  GRIB_name :
  Specific humidity
  GRIB_shortName :
  q
  GRIB_totalNumber :
  0
  GRIB_units :
  kg kg**-1
  long_name :
  Specific humidity
  units :
  kg kg**-1
  standard_name :
  specific_humidity
```
[2170812 values with dtype=float32]
```
- t
  (valid_time, pressure_level, latitude, longitude)
  float32
  ...
  GRIB_paramId :
  130
  GRIB_dataType :
  an
  GRIB_numberOfPoints :
  180901
  GRIB_typeOfLevel :
  isobaricInhPa
  GRIB_stepUnits :
  1
  GRIB_stepType :
  instant
  GRIB_gridType :
  regular_ll
  GRIB_uvRelativeToGrid :
  0
  GRIB_NV :
  0
  GRIB_Nx :
  601
  GRIB_Ny :
  301
  GRIB_cfName :
  air_temperature
  GRIB_cfVarName :
  t
  GRIB_gridDefinitionDescription :
  Latitude/Longitude Grid
  GRIB_iDirectionIncrementInDegrees :
  0.25
  GRIB_iScansNegatively :
  0
  GRIB_jDirectionIncrementInDegrees :
  0.25
  GRIB_jPointsAreConsecutive :
  0
  GRIB_jScansPositively :
  0
  GRIB_latitudeOfFirstGridPointInDegrees :
  90.0
  GRIB_latitudeOfLastGridPointInDegrees :
  15.0
  GRIB_longitudeOfFirstGridPointInDegrees :
  -100.0
  GRIB_longitudeOfLastGridPointInDegrees :
  50.0
  GRIB_missingValue :
  3.4028234663852886e+38
  GRIB_name :
  Temperature
  GRIB_shortName :
  t
  GRIB_totalNumber :
  0
  GRIB_units :
  K
  long_name :
  Temperature
  units :
  K
  standard_name :
  air_temperature
```
[2170812 values with dtype=float32]
```
- u
  (valid_time, pressure_level, latitude, longitude)
  float32
  ...
  GRIB_paramId :
  131
  GRIB_dataType :
  an
  GRIB_numberOfPoints :
  180901
  GRIB_typeOfLevel :
  isobaricInhPa
  GRIB_stepUnits :
  1
  GRIB_stepType :
  instant
  GRIB_gridType :
  regular_ll
  GRIB_uvRelativeToGrid :
  0
  GRIB_NV :
  0
  GRIB_Nx :
  601
  GRIB_Ny :
  301
  GRIB_cfName :
  eastward_wind
  GRIB_cfVarName :
  u
  GRIB_gridDefinitionDescription :
  Latitude/Longitude Grid
  GRIB_iDirectionIncrementInDegrees :
  0.25
  GRIB_iScansNegatively :
  0
  GRIB_jDirectionIncrementInDegrees :
  0.25
  GRIB_jPointsAreConsecutive :
  0
  GRIB_jScansPositively :
  0
  GRIB_latitudeOfFirstGridPointInDegrees :
  90.0
  GRIB_latitudeOfLastGridPointInDegrees :
  15.0
  GRIB_longitudeOfFirstGridPointInDegrees :
  -100.0
  GRIB_longitudeOfLastGridPointInDegrees :
  50.0
  GRIB_missingValue :
  3.4028234663852886e+38
  GRIB_name :
  U component of wind
  GRIB_shortName :
  u
  GRIB_totalNumber :
  0
  GRIB_units :
  m s**-1
  long_name :
  U component of wind
  units :
  m s**-1
  standard_name :
  eastward_wind
```
[2170812 values with dtype=float32]
```
- v
  (valid_time, pressure_level, latitude, longitude)
  float32
  ...
  GRIB_paramId :
  132
  GRIB_dataType :
  an
  GRIB_numberOfPoints :
  180901
  GRIB_typeOfLevel :
  isobaricInhPa
  GRIB_stepUnits :
  1
  GRIB_stepType :
  instant
  GRIB_gridType :
  regular_ll
  GRIB_uvRelativeToGrid :
  0
  GRIB_NV :
  0
  GRIB_Nx :
  601
  GRIB_Ny :
  301
  GRIB_cfName :
  northward_wind
  GRIB_cfVarName :
  v
  GRIB_gridDefinitionDescription :
  Latitude/Longitude Grid
  GRIB_iDirectionIncrementInDegrees :
  0.25
  GRIB_iScansNegatively :
  0
  GRIB_jDirectionIncrementInDegrees :
  0.25
  GRIB_jPointsAreConsecutive :
  0
  GRIB_jScansPositively :
  0
  GRIB_latitudeOfFirstGridPointInDegrees :
  90.0
  GRIB_latitudeOfLastGridPointInDegrees :
  15.0
  GRIB_longitudeOfFirstGridPointInDegrees :
  -100.0
  GRIB_longitudeOfLastGridPointInDegrees :
  50.0
  GRIB_missingValue :
  3.4028234663852886e+38
  GRIB_name :
  V component of wind
  GRIB_shortName :
  v
  GRIB_totalNumber :
  0
  GRIB_units :
  m s**-1
  long_name :
  V component of wind
  units :
  m s**-1
  standard_name :
  northward_wind
```
[2170812 values with dtype=float32]
```
Attributes: (6)
GRIB_centre :
ecmf
GRIB_centreDescription :
European Centre for Medium-Range Weather Forecasts
GRIB_subCentre :
0
Conventions :
CF-1.7
institution :
European Centre for Medium-Range Weather Forecasts
history :
2026-03-12T10:30 GRIB to CDM+CF via cfgrib-0.9.15.1/ecCodes-2.42.0 with {"source": "tmpfsay9t6q/data.grib", "filter_by_keys": {"stream": ["oper"], "stepType": ["instant"]}, "encode_cf": ["parameter", "time", "geography", "vertical"]}

Although the dataset has the earthkit accessor (it is always available when earthkit-data is installed) the grid spec is not available on it.

[8]:

ds.earthkit.grid_spec  # returns None

Currently, earthkit-geo is not able infer the grid spec if it is not available on the earthkit accessor. As a result, in order to make the interpolation work we need to provide the grid spec via the in_grid kwarg to regrid().

[9]:

in_grid = {"grid": [0.25, 0.25], "area": [90.0, -100.0, 15.0, 50.0]}

# the target grid is 5x5 degree regular latitude-longitude grid
out_grid = {"grid": [5, 5]}

ds_r = ekg.regrid(ds, in_grid=in_grid, out_grid=out_grid, interpolation="linear")
ds_r

[9]:

<xarray.Dataset> Size: 96kB
Dimensions:         (valid_time: 4, pressure_level: 3, latitude: 16,
                     longitude: 31)
Coordinates:
  * valid_time      (valid_time) datetime64[ns] 32B 2016-09-26 ... 2016-09-26...
  * pressure_level  (pressure_level) float64 24B 925.0 850.0 700.0
  * latitude        (latitude) float64 128B 90.0 85.0 80.0 ... 25.0 20.0 15.0
  * longitude       (longitude) float64 248B -100.0 -95.0 -90.0 ... 45.0 50.0
Data variables:
    q               (valid_time, pressure_level, latitude, longitude) float32 24kB ...
    t               (valid_time, pressure_level, latitude, longitude) float32 24kB ...
    u               (valid_time, pressure_level, latitude, longitude) float32 24kB ...
    v               (valid_time, pressure_level, latitude, longitude) float32 24kB ...

xarray.Dataset

Dimensions:
- valid_time: 4
- pressure_level: 3
- latitude: 16
- longitude: 31

Coordinates: (4)

valid_time

(valid_time)

datetime64[ns]

2016-09-26 ... 2016-09-26T18:00:00

long_name :: time
standard_name :: time

array(['2016-09-26T00:00:00.000000000', '2016-09-26T06:00:00.000000000',
       '2016-09-26T12:00:00.000000000', '2016-09-26T18:00:00.000000000'],
      dtype='datetime64[ns]')

pressure_level
(pressure_level)
float64
925.0 850.0 700.0
long_name :
pressure
units :
hPa
positive :
down
stored_direction :
decreasing
standard_name :
air_pressure
```
array([925., 850., 700.])
```

latitude

(latitude)

float64

90.0 85.0 80.0 ... 25.0 20.0 15.0

array([90., 85., 80., 75., 70., 65., 60., 55., 50., 45., 40., 35., 30., 25.,
       20., 15.])

longitude

(longitude)

float64

-100.0 -95.0 -90.0 ... 45.0 50.0

array([-100.,  -95.,  -90.,  -85.,  -80.,  -75.,  -70.,  -65.,  -60.,  -55.,
        -50.,  -45.,  -40.,  -35.,  -30.,  -25.,  -20.,  -15.,  -10.,   -5.,
          0.,    5.,   10.,   15.,   20.,   25.,   30.,   35.,   40.,   45.,
         50.])

Data variables: (4)

(valid_time, pressure_level, latitude, longitude)

float32

0.003074 0.003074 ... 0.006416

GRIB_paramId :: 133
GRIB_dataType :: an
GRIB_numberOfPoints :: 180901
GRIB_typeOfLevel :: isobaricInhPa
GRIB_stepUnits :: 1
GRIB_stepType :: instant
GRIB_gridType :: regular_ll
GRIB_uvRelativeToGrid :: 0
GRIB_NV :: 0
GRIB_Nx :: 601
GRIB_Ny :: 301
GRIB_cfName :: specific_humidity
GRIB_cfVarName :: q
GRIB_gridDefinitionDescription :: Latitude/Longitude Grid
GRIB_iDirectionIncrementInDegrees :: 0.25
GRIB_iScansNegatively :: 0
GRIB_jDirectionIncrementInDegrees :: 0.25
GRIB_jPointsAreConsecutive :: 0
GRIB_jScansPositively :: 0
GRIB_latitudeOfFirstGridPointInDegrees :: 90.0
GRIB_latitudeOfLastGridPointInDegrees :: 15.0
GRIB_longitudeOfFirstGridPointInDegrees :: -100.0
GRIB_longitudeOfLastGridPointInDegrees :: 50.0
GRIB_missingValue :: 3.4028234663852886e+38
GRIB_name :: Specific humidity
GRIB_shortName :: q
GRIB_totalNumber :: 0
GRIB_units :: kg kg**-1
long_name :: Specific humidity
units :: kg kg**-1
standard_name :: specific_humidity

array([[[[0.00307444, 0.00307444, 0.00307444, ..., 0.00307444,
          0.00307444, 0.00307444],
         [0.00281265, 0.00304201, 0.00326326, ..., 0.0032666 ,
          0.00317791, 0.00295713],
         [0.00299242, 0.00369766, 0.00389269, ..., 0.00129106,
          0.00148275, 0.00135448],
         ...,
         [0.01354387, 0.01518467, 0.01445034, ..., 0.00574377,
          0.00430801, 0.00583246],
         [0.01278332, 0.01550939, 0.0140498 , ..., 0.01537207,
          0.00397995, 0.00607326],
         [0.01478508, 0.01405456, 0.01484707, ..., 0.00938299,
          0.00475624, 0.01519516]],

        [[0.00278865, 0.00278865, 0.00278865, ..., 0.00278865,
          0.00278865, 0.00278865],
         [0.00323926, 0.00343333, 0.00348912, ..., 0.00332652,
          0.00324498, 0.00311671],
         [0.00353108, 0.00337802, 0.00405847, ..., 0.00079547,
          0.00083886, 0.00097524],
...
         [0.01244843, 0.01115573, 0.0091306 , ..., 0.00449288,
          0.00361073, 0.00485766],
         [0.0095197 , 0.0127579 , 0.01093018, ..., 0.0069958 ,
          0.0044781 , 0.00433505],
         [0.01305926, 0.01232684, 0.01200307, ..., 0.0069176 ,
          0.00496161, 0.00657571]],

        [[0.0016844 , 0.0016844 , 0.0016844 , ..., 0.0016844 ,
          0.0016844 , 0.0016844 ],
         [0.00084445, 0.00161026, 0.00209234, ..., 0.00050709,
          0.00039337, 0.00061343],
         [0.00026033, 0.00150106, 0.00177405, ..., 0.00136135,
          0.00121424, 0.00047538],
         ...,
         [0.00857995, 0.00706933, 0.00701783, ..., 0.00467727,
          0.00360415, 0.00098059],
         [0.00879571, 0.0097079 , 0.0060787 , ..., 0.00614355,
          0.00439856, 0.00365446],
         [0.0083668 , 0.00853298, 0.00892375, ..., 0.00507019,
          0.00490782, 0.00641558]]]], shape=(4, 3, 16, 31), dtype=float32)

(valid_time, pressure_level, latitude, longitude)

float32

269.7 269.7 269.7 ... 285.6 284.2

GRIB_paramId :: 130
GRIB_dataType :: an
GRIB_numberOfPoints :: 180901
GRIB_typeOfLevel :: isobaricInhPa
GRIB_stepUnits :: 1
GRIB_stepType :: instant
GRIB_gridType :: regular_ll
GRIB_uvRelativeToGrid :: 0
GRIB_NV :: 0
GRIB_Nx :: 601
GRIB_Ny :: 301
GRIB_cfName :: air_temperature
GRIB_cfVarName :: t
GRIB_gridDefinitionDescription :: Latitude/Longitude Grid
GRIB_iDirectionIncrementInDegrees :: 0.25
GRIB_iScansNegatively :: 0
GRIB_jDirectionIncrementInDegrees :: 0.25
GRIB_jPointsAreConsecutive :: 0
GRIB_jScansPositively :: 0
GRIB_latitudeOfFirstGridPointInDegrees :: 90.0
GRIB_latitudeOfLastGridPointInDegrees :: 15.0
GRIB_longitudeOfFirstGridPointInDegrees :: -100.0
GRIB_longitudeOfLastGridPointInDegrees :: 50.0
GRIB_missingValue :: 3.4028234663852886e+38
GRIB_name :: Temperature
GRIB_shortName :: t
GRIB_totalNumber :: 0
GRIB_units :: K
long_name :: Temperature
units :: K
standard_name :: air_temperature

array([[[[269.66333, 269.66333, 269.66333, ..., 269.66333, 269.66333,
          269.66333],
         [268.58325, 269.6067 , 270.53833, ..., 270.4817 , 270.30786,
          269.9563 ],
         [269.4309 , 271.77856, 274.40747, ..., 271.45044, 270.63013,
          270.927  ],
         ...,
         [299.13013, 295.98755, 295.2102 , ..., 303.42114, 303.94458,
          303.01294],
         [295.0735 , 294.9641 , 296.67896, ..., 302.7688 , 301.92505,
          304.22583],
         [295.15942, 295.00317, 297.1438 , ..., 301.99146, 303.36255,
          300.14966]],

        [[268.47852, 268.47852, 268.47852, ..., 268.47852, 268.47852,
          268.47852],
         [270.2754 , 271.03516, 270.22266, ..., 272.4707 , 270.7871 ,
          269.48438],
         [270.6875 , 269.85742, 272.44922, ..., 271.25977, 270.64453,
          269.64258],
...
         [290.80005, 290.69263, 290.93872, ..., 301.8645 , 302.6438 ,
          298.48755],
         [293.2727 , 290.64966, 292.09692, ..., 300.28052, 300.2024 ,
          299.32935],
         [291.1848 , 290.95825, 294.95435, ..., 299.1223 , 300.65552,
          297.29614]],

        [[263.2678 , 263.2678 , 263.2678 , ..., 263.2678 , 263.2678 ,
          263.2678 ],
         [261.09396, 261.35178, 261.78928, ..., 267.51974, 267.2131 ,
          266.5471 ],
         [258.9006 , 259.97287, 260.4553 , ..., 266.9006 , 267.22678,
          267.39865],
         ...,
         [281.58615, 282.47092, 281.66232, ..., 286.72873, 286.24435,
          285.1877 ],
         [284.5217 , 282.67795, 282.41818, ..., 286.04318, 283.93967,
          284.6428 ],
         [282.88693, 282.32053, 283.20334, ..., 285.467  , 285.64865,
          284.20724]]]], shape=(4, 3, 16, 31), dtype=float32)

(valid_time, pressure_level, latitude, longitude)

float32

7.629e-05 7.629e-05 ... -4.351

GRIB_paramId :: 131
GRIB_dataType :: an
GRIB_numberOfPoints :: 180901
GRIB_typeOfLevel :: isobaricInhPa
GRIB_stepUnits :: 1
GRIB_stepType :: instant
GRIB_gridType :: regular_ll
GRIB_uvRelativeToGrid :: 0
GRIB_NV :: 0
GRIB_Nx :: 601
GRIB_Ny :: 301
GRIB_cfName :: eastward_wind
GRIB_cfVarName :: u
GRIB_gridDefinitionDescription :: Latitude/Longitude Grid
GRIB_iDirectionIncrementInDegrees :: 0.25
GRIB_iScansNegatively :: 0
GRIB_jDirectionIncrementInDegrees :: 0.25
GRIB_jPointsAreConsecutive :: 0
GRIB_jScansPositively :: 0
GRIB_latitudeOfFirstGridPointInDegrees :: 90.0
GRIB_latitudeOfLastGridPointInDegrees :: 15.0
GRIB_longitudeOfFirstGridPointInDegrees :: -100.0
GRIB_longitudeOfLastGridPointInDegrees :: 50.0
GRIB_missingValue :: 3.4028234663852886e+38
GRIB_name :: U component of wind
GRIB_shortName :: u
GRIB_totalNumber :: 0
GRIB_units :: m s**-1
long_name :: U component of wind
units :: m s**-1
standard_name :: eastward_wind

array([[[[ 7.62939453e-05,  7.62939453e-05,  7.62939453e-05, ...,
           7.62939453e-05,  7.62939453e-05,  7.62939453e-05],
         [ 2.65437317e+00,  3.41218567e+00,  4.39656067e+00, ...,
           7.12702942e+00,  6.54499817e+00,  6.22077942e+00],
         [ 2.25398254e+00, -7.21893311e-02,  9.49295044e-01, ...,
           6.43367004e+00,  7.09577942e+00,  5.29695129e+00],
         ...,
         [-9.90158081e-01, -6.96672058e+00, -6.50773621e+00, ...,
           2.81648254e+00, -1.63078308e+00,  3.29695129e+00],
         [-3.59298706e-01, -3.45500183e+00, -7.86320496e+00, ...,
          -5.09689331e-01, -2.35539246e+00, -8.81632996e+00],
         [ 4.69512939e-02, -4.45304871e+00, -8.00018311e-02, ...,
           2.97468567e+00,  6.77810669e-01, -1.42501831e-01]],

        [[-8.54492188e-04, -8.54492188e-04, -8.54492188e-04, ...,
          -8.54492188e-04, -8.54492188e-04, -8.54492188e-04],
         [ 3.80578613e+00,  4.12219238e+00,  4.64562988e+00, ...,
           1.25968018e+01,  1.10850830e+01,  9.61242676e+00],
         [ 2.69641113e+00, -1.46765137e+00,  1.82727051e+00, ...,
           6.41320801e+00,  6.63781738e+00,  5.98352051e+00],
...
           5.23757935e+00, -4.34295654e-01,  2.77664185e+00],
         [-2.09249878e+00,  1.89578247e+00, -3.77413940e+00, ...,
          -1.50265503e+00, -6.53195190e+00, -3.74093628e+00],
         [-8.63983154e-01, -3.99093628e+00, -1.93429565e+00, ...,
          -1.79367065e+00,  3.44543457e-02,  4.68048096e-01]],

        [[ 1.67846680e-04,  1.67846680e-04,  1.67846680e-04, ...,
           1.67846680e-04,  1.67846680e-04,  1.67846680e-04],
         [ 3.57829285e+00,  3.00798035e+00,  2.24040222e+00, ...,
           1.96329803e+01,  1.78517303e+01,  1.52521210e+01],
         [ 8.93376160e+00,  4.48649597e+00,  3.14274597e+00, ...,
           7.46696472e+00,  7.56852722e+00,  6.33610535e+00],
         ...,
         [-1.83186340e+00, -5.51545715e+00, -1.99592590e+00, ...,
           1.04118347e+00, -1.46858215e+00, -3.16194153e+00],
         [-3.66194153e+00,  1.39860535e+00, -4.20491028e+00, ...,
          -6.50373840e+00, -5.66584778e+00,  2.31462097e+00],
         [-1.17366028e+00, -2.03694153e+00, -2.82209778e+00, ...,
          -7.27326965e+00, -7.96467590e+00, -4.35139465e+00]]]],
      shape=(4, 3, 16, 31), dtype=float32)

(valid_time, pressure_level, latitude, longitude)

float32

0.0007172 0.0007172 ... -4.286

GRIB_paramId :: 132
GRIB_dataType :: an
GRIB_numberOfPoints :: 180901
GRIB_typeOfLevel :: isobaricInhPa
GRIB_stepUnits :: 1
GRIB_stepType :: instant
GRIB_gridType :: regular_ll
GRIB_uvRelativeToGrid :: 0
GRIB_NV :: 0
GRIB_Nx :: 601
GRIB_Ny :: 301
GRIB_cfName :: northward_wind
GRIB_cfVarName :: v
GRIB_gridDefinitionDescription :: Latitude/Longitude Grid
GRIB_iDirectionIncrementInDegrees :: 0.25
GRIB_iScansNegatively :: 0
GRIB_jDirectionIncrementInDegrees :: 0.25
GRIB_jPointsAreConsecutive :: 0
GRIB_jScansPositively :: 0
GRIB_latitudeOfFirstGridPointInDegrees :: 90.0
GRIB_latitudeOfLastGridPointInDegrees :: 15.0
GRIB_longitudeOfFirstGridPointInDegrees :: -100.0
GRIB_longitudeOfLastGridPointInDegrees :: 50.0
GRIB_missingValue :: 3.4028234663852886e+38
GRIB_name :: V component of wind
GRIB_shortName :: v
GRIB_totalNumber :: 0
GRIB_units :: m s**-1
long_name :: V component of wind
units :: m s**-1
standard_name :: northward_wind

array([[[[ 7.17163086e-04,  7.17163086e-04,  7.17163086e-04, ...,
           7.17163086e-04,  7.17163086e-04,  7.17163086e-04],
         [ 5.28392029e+00,  6.75657654e+00,  8.26048279e+00, ...,
          -8.82154846e+00, -9.05787659e+00, -8.81764221e+00],
         [-2.06373596e+00,  4.51438904e+00,  8.47015381e-02, ...,
          -5.37818909e+00, -5.02662659e+00, -5.24342346e+00],
         ...,
         [-2.84439087e-01,  3.22337341e+00,  1.38548279e+00, ...,
          -3.52857971e+00,  3.85618591e+00, -8.35279846e+00],
         [-5.75454712e-01, -1.04751587e-01,  3.46360779e+00, ...,
           5.30014038e-01,  1.13548279e+00, -1.08840485e+01],
         [ 9.44671631e-02,  4.97727966e+00, -2.08267212e-01, ...,
          -2.26686096e+00, -1.39186096e+00, -4.66079712e-01]],

        [[-5.95092773e-04, -5.95092773e-04, -5.95092773e-04, ...,
          -5.95092773e-04, -5.95092773e-04, -5.95092773e-04],
         [ 7.69081116e+00,  5.66151428e+00,  8.03456116e+00, ...,
          -1.01412201e+01, -1.01900482e+01, -9.92832947e+00],
         [-4.84970093e-01,  9.44471741e+00,  2.17909241e+00, ...,
          -5.09629822e+00, -5.92051697e+00, -6.77793884e+00],
...
          -7.10606384e+00, -1.73692322e+00, -1.30220795e+01],
         [-1.98497009e+00,  1.50331116e+00,  5.54042053e+00, ...,
           1.35803223e-03, -9.18563843e-01, -5.18614197e+00],
         [-2.59043884e+00,  6.50526428e+00, -3.37982178e-02, ...,
           5.13076782e-01, -1.59434509e+00,  2.46034241e+00]],

        [[ 6.10351562e-04,  6.10351562e-04,  6.10351562e-04, ...,
           6.10351562e-04,  6.10351562e-04,  6.10351562e-04],
         [ 3.69396973e+00,  3.12756348e+00,  4.63342285e+00, ...,
          -7.41345215e+00, -9.87438965e+00, -1.11712646e+01],
         [-6.65368652e+00, -6.27087402e+00,  9.84985352e-01, ...,
          -6.43688965e+00, -1.08802490e+01, -1.47591553e+01],
         ...,
         [ 8.91235352e-01,  3.18969727e-01,  7.79907227e-01, ...,
           1.89318848e+00, -8.90014648e-01, -1.54993896e+01],
         [-3.14978027e+00,  1.19982910e+00,  6.86389160e+00, ...,
          -2.25134277e+00, -1.18103027e+00, -7.00329590e+00],
         [-4.26696777e+00,  6.48498535e+00,  1.56115723e+00, ...,
          -5.03295898e-01, -5.82751465e+00, -4.28649902e+00]]]],
      shape=(4, 3, 16, 31), dtype=float32)

We check the result visually by plotting the same data from the input and output dataset.

[10]:

ekp.quickplot(ds.t.isel(valid_time=0, pressure_level=0))
ekp.quickplot(ds_r.t.isel(valid_time=0, pressure_level=0))

[10]:

<earthkit.plots.components.maps.Map at 0x17f021e00>

../_images/tutorials_mir_regrid_xarray_23_1.png

../_images/tutorials_mir_regrid_xarray_23_2.png

[ ]: