Exctract a level from a 3D field

# for figures in notebook
%matplotlib inline
# import & initialize epygram
import epygram
epygram.init_env()
# path to file
import os
INPUTS_DIR = os.path.join('..', 'inputs')

# open resource as a meta-CL-resource, CL = Combine Levels
r = epygram.resources.meta_resource(os.path.join(INPUTS_DIR, 'ICMSHAROM+0022'), 'r', 'CL')
# in this case, field need to be defined with GRIB2 nomenclature

f = r.readfield({'discipline':0, 'parameterCategory':0, 'parameterNumber':0, 'typeOfFirstFixedSurface':119})

f.sp2gp()

f18 = f.getlevel(18)

print(type(f))
print(type(f18))
print(f18.geometry.vcoordinate.levels)

<class 'epygram.fields.D3Field.D3Field'>
<class 'epygram.fields.H2DField.H2DField'>
[18]