pyseobnr ignores lmax

If I run:

import lal
import lalsimulation as lalsim
import numpy as np

import astropy.units as u
import lalsimulation.gwsignal.core.waveform as gwsignal_waveform
from lalsimulation.gwsignal.models import pyseobnr_model

params_dict = {
    'mass1' : 50*u.solMass,
    'mass2' : 50*u.solMass,
    'spin1x' : 0*u.dimensionless_unscaled,
    'spin1y' : 0*u.dimensionless_unscaled,
    'spin1z' : 0*u.dimensionless_unscaled,
    'spin2x' : 0*u.dimensionless_unscaled,
    'spin2y' : 0*u.dimensionless_unscaled,
    'spin2z' : 0*u.dimensionless_unscaled,
    'deltaT' : 1./4096 * u.s,
    'f22_start' : 20*u.Hz,
    'f22_ref': 20*u.Hz,
    'distance' : 1*u.Mpc,
    'eccentricity' : 0*u.dimensionless_unscaled,
    'meanPerAno' : 0*u.rad,
    'lmax': 5,
    'condition' : 0,
}

gen = pyseobnr_model.SEOBNRv5PHM()
hlm = gwsignal_waveform.GenerateTDModes(params_dict, gen)
print (hlm.keys())

I only get

dict_keys([(2, -2), (2, -1), (2, 0), (2, 1), (2, 2), (3, -3), (3, -2), (3, -1), (3, 0), (3, 1), (3, 2), (3, 3), (4, -4), (4, -3), (4, -2), (4, -1), (4, 0), (4, 1), (4, 2), (4, 3), (4, 4)])