In #556 Layered HFB packages have been added. A user can define multiple linestring and per linestring a layer can be specified, However the layer specified for a linestring is applied to all linestrings

barrier_gdf = gpd.GeoDataFrame(
  geometry=
    [
      linestring1,
      linestring2
    ],
  data=
    {
      "hydraulic_characteristic": [1e-3, 1e-3,],
      "layer": [1, 2,]
    },
)


hfb = imod.mf6.LayeredHorizontalFlowBarrierHydraulicCharacteristic(barrier_gdf)

In the snippet above I would expect that linestring1 is applied to layer 1 and that linestring2 is applied to layer 2. However linestring1 is applied to layers 1 and 2 and the same applies for linestrin2