Make an iris cube¶

In [8]:

from iris.coords import DimCoord
from iris.cube import Cube
import numpy as np

Define cube coordinates¶

In [9]:

ysize = 20
xsize = 10

lats = np.linspace(-90, 90, ysize)
lons = np.linspace(0, 360,  xsize)

latitude  = DimCoord(lats, standard_name='latitude',  units='degrees')
longitude = DimCoord(lons, standard_name='longitude', units='degrees')

Build cube¶

In [10]:

long_name = 'something'
units     = 'kg'

cube = Cube(    np.zeros((ysize, xsize), np.float32), 
                dim_coords_and_dims=[(latitude, 0),(longitude, 1)], 
                long_name=long_name, units=units 
           )

In [11]:

print(cube)

something / (kg)                    (latitude: 20; longitude: 10)
     Dimension coordinates:
          latitude                           x              -
          longitude                          -              x

Insert data in cube¶

make sure the shape matches!

In [ ]:

cube.data = my_numpy_array