Skip to content

External forcings file

The external forcing .ext file contains the forcing data for a D-Flow FM model. This includes open boundaries, lateral discharges, source/sink terms, bubble screens, and meteorological forcings. The documentation below only concerns the 'new' format (ExtForceFileNew in the MDU file).

Class hierarchy

The [Meteo] and [Spatial] blocks (this file) and the [Initial] / [Parameter] blocks (the initial field file) are spatial-field blocks that share their field validation. The shared logic lives in two plain validator mixins, while INIBasedModel is mixed in only at the two concrete block bases.

Hold "Ctrl" to enable pan & zoom
classDiagram
    direction TB

    class INIBasedModel
    class OperandInterpolationValidators {
        <<mixin>>
    }
    class LocationTypeDataFileTypeValidators {
        <<mixin>>
    }
    class SpatialForcingBase {
        <<abstract>>
    }
    class AbstractSpatialField {
        <<abstract>>
    }
    class Meteo
    class Spatial
    class InitialField
    class ParameterField

    OperandInterpolationValidators <|-- SpatialForcingBase
    OperandInterpolationValidators <|-- AbstractSpatialField
    INIBasedModel <|-- SpatialForcingBase
    LocationTypeDataFileTypeValidators <|-- AbstractSpatialField
    INIBasedModel <|-- AbstractSpatialField
    SpatialForcingBase <|-- Meteo
    SpatialForcingBase <|-- Spatial
    LocationTypeDataFileTypeValidators <|-- Spatial
    AbstractSpatialField <|-- InitialField
    AbstractSpatialField <|-- ParameterField

    note for OperandInterpolationValidators "operand + interpolationMethod validators (all four blocks)"
    note for LocationTypeDataFileTypeValidators "locationType + dataFileType validators (not Meteo)"
    note for SpatialForcingBase "file-model resolution, shared comments, unknown-keyword policy"
    note for AbstractSpatialField "averagingType + datafile validators"

Reading the diagram:

  • OperandInterpolationValidators and LocationTypeDataFileTypeValidators are independent plain mixins (not models, and neither inherits the other). The former holds the operand / interpolationMethod validators shared by all four blocks; the latter holds the locationType / dataFileType validators shared by Spatial, InitialField and ParameterField (not Meteo, which has neither field).
  • INIBasedModel is inherited only by the two concrete block bases: SpatialForcingBase (Meteo / Spatial) and AbstractSpatialField (InitialField / ParameterField).
  • Each block base inherits exactly the validator mixins it needs: SpatialForcingBase inherits OperandInterpolationValidators; AbstractSpatialField inherits both mixins; and Spatial adds LocationTypeDataFileTypeValidators on top of SpatialForcingBase.
  • averagingType is validated per-class (Spatial, AbstractSpatialField) rather than in a shared mixin, because Meteo reaches OperandInterpolationValidators too and stores averagingType as a raw integer.

Model

Models for the external forcings file (new format) of D-Flow FM.

Boundary

Bases: INIBasedModel

A [Boundary] block for use inside an external forcings file.

I.e., a ExtModel.

All lowercased attributes match with the boundary input as described in UM Sec.C.5.2.1.

Source code in hydrolib/core/dflowfm/ext/models.py
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
class Boundary(INIBasedModel):
    """A `[Boundary]` block for use inside an external forcings file.

    I.e., a [ExtModel][hydrolib.core.dflowfm.ext.models.ExtModel].

    All lowercased attributes match with the boundary input as described in
    [UM Sec.C.5.2.1](https://content.oss.deltares.nl/delft3dfm1d2d/D-Flow_FM_User_Manual_1D2D.pdf#subsection.C.5.2.1).
    """

    _header: Literal["Boundary"] = "Boundary"
    quantity: str = Field(alias="quantity")
    nodeid: str | None = Field(None, alias="nodeId")
    locationfile: Annotated[
        DiskOnlyFileModel, BeforeValidator(set_default_disk_only_file_model)
    ] = Field(default_factory=lambda: DiskOnlyFileModel(None), alias="locationFile")
    forcingfile: ForcingModel | DiskOnlyFileModel = Field(alias="forcingFile")
    bndwidth1d: float | None = Field(None, alias="bndWidth1D")
    bndbldepth: float | None = Field(None, alias="bndBlDepth")
    returntime: float | None = Field(None, alias="returnTime")
    operand: Operand | None = Field(None, alias="operand")

    def is_intermediate_link(self) -> bool:
        return True

    @field_validator("forcingfile", mode="before")
    @classmethod
    def validate_forcingfile(cls, data: Any) -> Any:
        if isinstance(data, (str, Path)):
            data = resolve_file_model(data, ForcingModel)
        elif not isinstance(data, (ForcingModel, DiskOnlyFileModel)):
            raise TypeError(
                "Forcing file must be a ForcingModel or a path to a forcing file."
            )
        return data

    @classmethod
    def _is_valid_locationfile_data(
        cls, elem: Union[None, str, Path, DiskOnlyFileModel]
    ) -> bool:
        return isinstance(elem, Path) or (
            isinstance(elem, DiskOnlyFileModel) and elem.filepath is not None
        )

    @classmethod
    def _exclude_from_validation(cls, input_data: Optional[dict] = None) -> Set:
        unknown_keywords = ["return_time"]
        return set(unknown_keywords)

    @model_validator(mode="before")
    @classmethod
    def rename_return_time_field(cls, values: Dict) -> Dict:
        """Renames the deprecated return_time field to returnTime.

        Args:
            values (Dict): Dictionary with raw, unvalidated input values.

        Returns:
            Dict: Validated dictionary of values for Boundary.
        """
        if "return_time" in values:
            values["returnTime"] = values.pop("return_time")
        return values

    @model_validator(mode="before")
    @classmethod
    def check_nodeid_or_locationfile_present(cls, values: Dict) -> Dict:
        """Verifies that either nodeid or locationfile properties have been set.

        Args:
            values (Dict): Dictionary with values already validated.

        Raises:
            ValueError: When none of the values are present.

        Returns:
            Dict: Validated dictionary of values for Boundary.
        """
        node_id = values.get("nodeid")
        location_file = values.get("locationfile")
        if str_is_empty_or_none(node_id) and not cls._is_valid_locationfile_data(
            location_file
        ):
            raise ValueError(
                "Either nodeId or locationFile fields should be specified."
            )
        return values

    def _get_identifier(self, data: dict) -> Optional[str]:
        """
        Retrieves the identifier for a boundary, which is the nodeid.

        Args:
            data (dict): Dictionary of values for this boundary.

        Returns:
            str: The nodeid value or None if not found.
        """
        return data.get("nodeid")

    @property
    def forcing(self) -> Union[ForcingBase, None]:
        """Retrieves the corresponding forcing data for this boundary.

        Returns:
            ForcingBase: The corresponding forcing data, or None when no matching forcing block is
            found or when the forcing file has not been parsed (e.g. loaded with `recurse=False`,
            in which case it is a `DiskOnlyFileModel` placeholder).
        """
        result = None
        if isinstance(self.forcingfile, ForcingModel):
            for forcing in self.forcingfile.forcing:
                if self.nodeid == forcing.name and any(
                    quantity.quantity.startswith(self.quantity)
                    for quantity in forcing.quantityunitpair
                ):
                    result = forcing
                    break

        return result

    @model_validator(mode="before")
    @classmethod
    def validate_locationfile(cls, data: Any) -> Any:
        file_location = data.get("locationfile") or data.get("locationFile")
        data.pop("locationFile", None)  # Remove alias if present

        # Convert string to DiskOnlyFileModel if needed
        if isinstance(file_location, (str, Path)):
            data["locationfile"] = DiskOnlyFileModel(file_location)
        return data

    @field_validator("operand", mode="before")
    @classmethod
    def validate_operand(cls, v: Any):
        return enum_value_parser(v, Operand, Operand.legacy_alternatives())

forcing property

Retrieves the corresponding forcing data for this boundary.

Returns:

Name Type Description
ForcingBase Union[ForcingBase, None]

The corresponding forcing data, or None when no matching forcing block is

Union[ForcingBase, None]

found or when the forcing file has not been parsed (e.g. loaded with recurse=False,

Union[ForcingBase, None]

in which case it is a DiskOnlyFileModel placeholder).

check_nodeid_or_locationfile_present(values) classmethod

Verifies that either nodeid or locationfile properties have been set.

Parameters:

Name Type Description Default
values Dict

Dictionary with values already validated.

required

Raises:

Type Description
ValueError

When none of the values are present.

Returns:

Name Type Description
Dict Dict

Validated dictionary of values for Boundary.

Source code in hydrolib/core/dflowfm/ext/models.py
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
@model_validator(mode="before")
@classmethod
def check_nodeid_or_locationfile_present(cls, values: Dict) -> Dict:
    """Verifies that either nodeid or locationfile properties have been set.

    Args:
        values (Dict): Dictionary with values already validated.

    Raises:
        ValueError: When none of the values are present.

    Returns:
        Dict: Validated dictionary of values for Boundary.
    """
    node_id = values.get("nodeid")
    location_file = values.get("locationfile")
    if str_is_empty_or_none(node_id) and not cls._is_valid_locationfile_data(
        location_file
    ):
        raise ValueError(
            "Either nodeId or locationFile fields should be specified."
        )
    return values

rename_return_time_field(values) classmethod

Renames the deprecated return_time field to returnTime.

Parameters:

Name Type Description Default
values Dict

Dictionary with raw, unvalidated input values.

required

Returns:

Name Type Description
Dict Dict

Validated dictionary of values for Boundary.

Source code in hydrolib/core/dflowfm/ext/models.py
239
240
241
242
243
244
245
246
247
248
249
250
251
252
@model_validator(mode="before")
@classmethod
def rename_return_time_field(cls, values: Dict) -> Dict:
    """Renames the deprecated return_time field to returnTime.

    Args:
        values (Dict): Dictionary with raw, unvalidated input values.

    Returns:
        Dict: Validated dictionary of values for Boundary.
    """
    if "return_time" in values:
        values["returnTime"] = values.pop("return_time")
    return values

BoundaryError

Bases: Exception

BoundaryError.

Source code in hydrolib/core/dflowfm/ext/models.py
1390
1391
1392
1393
1394
1395
class BoundaryError(Exception):
    """BoundaryError."""

    def __init__(self, error_message: str):
        """Initialize with an error message."""
        super().__init__(error_message)

__init__(error_message)

Initialize with an error message.

Source code in hydrolib/core/dflowfm/ext/models.py
1393
1394
1395
def __init__(self, error_message: str):
    """Initialize with an error message."""
    super().__init__(error_message)

BubbleScreen

Bases: CoordinateValidator, INIBasedModel

A [BubbleScreen] block for use inside an external forcings file.

Represents a bubble screen (air curtain): a row of nozzles on the waterway bed that releases air bubbles to create a hydraulic barrier. Location is specified either via a locationFile (a plain .pli) or via inline numCoordinates/xCoordinates/yCoordinates. Unlike SourceSink, BubbleScreen polylines carry no z information; vertical placement is given exclusively by the required zLevel field.

Attributes:

Name Type Description
id str

Unique identifier for the block within the ext file.

name str

Human-readable label. Defaults to an empty string.

locationfile Annotated[DiskOnlyFileModel, BeforeValidator(set_default_disk_only_file_model)]

Reference to an external .pli polyline file that defines the screen location. Mutually exclusive with inline coordinates.

numcoordinates int | None

Number of vertices in the inline polyline.

xcoordinates list[float] | None

X-coordinates of the inline polyline vertices.

ycoordinates list[float] | None

Y-coordinates of the inline polyline vertices.

zlevel float

Depth of the nozzle row in model vertical coordinate units.

discharge ForcingData

Volumetric air-flow rate. Accepts a scalar float, the string "realtime", or a path to a .bc forcing file.

Raises:

Type Description
ValidationError

If neither locationFile nor the full set of inline coordinate fields (numCoordinates, xCoordinates, yCoordinates) is provided.

ValidationError

If numCoordinates does not match the length of xCoordinates or yCoordinates.

Examples:

  • Build a block with inline coordinates and inspect the geometry:
    >>> from hydrolib.core.dflowfm.ext.models import BubbleScreen
    >>> screen = BubbleScreen(
    ...     id="air_curtain",
    ...     numcoordinates=2,
    ...     xcoordinates=[100.0, 200.0],
    ...     ycoordinates=[50.0, 50.0],
    ...     zlevel=-3.5,
    ...     discharge=1.2,
    ... )
    >>> screen.id
    'air_curtain'
    >>> screen.zlevel
    -3.5
    >>> screen.xcoordinates
    [100.0, 200.0]
    
  • Build a block that references an external polyline file:
    >>> from pathlib import Path
    >>> from hydrolib.core.base.models import DiskOnlyFileModel
    >>> from hydrolib.core.dflowfm.ext.models import BubbleScreen
    >>> screen = BubbleScreen(
    ...     id="canal_lock",
    ...     locationfile=DiskOnlyFileModel(filepath=Path("lock_screen.pli")),
    ...     zlevel=-6.0,
    ...     discharge=0.8,
    ... )
    >>> screen.id
    'canal_lock'
    >>> str(screen.locationfile.filepath)
    'lock_screen.pli'
    
  • Omitting both location forms raises ValidationError:
    >>> from pydantic import ValidationError
    >>> from hydrolib.core.dflowfm.ext.models import BubbleScreen
    >>> try:
    ...     BubbleScreen(id="bad", zlevel=-5.0, discharge=1.0)
    ... except ValidationError:
    ...     print("location is required")
    location is required
    
See Also

SourceSink: Sister block for water and solute injection or extraction. ExtModel: The new-format external forcings file model that hosts [BubbleScreen] blocks under its bubblescreen list.

Source code in hydrolib/core/dflowfm/ext/models.py
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
class BubbleScreen(CoordinateValidator, INIBasedModel):
    """A `[BubbleScreen]` block for use inside an external forcings file.

    Represents a bubble screen (air curtain): a row of nozzles on the waterway
    bed that releases air bubbles to create a hydraulic barrier. Location is
    specified either via a `locationFile` (a plain `.pli`) or via inline
    `numCoordinates`/`xCoordinates`/`yCoordinates`. Unlike `SourceSink`,
    BubbleScreen polylines carry no z information; vertical placement is given
    exclusively by the required `zLevel` field.

    Attributes:
        id: Unique identifier for the block within the ext file.
        name: Human-readable label. Defaults to an empty string.
        locationfile: Reference to an external `.pli` polyline file that
            defines the screen location. Mutually exclusive with inline
            coordinates.
        numcoordinates: Number of vertices in the inline polyline.
        xcoordinates: X-coordinates of the inline polyline vertices.
        ycoordinates: Y-coordinates of the inline polyline vertices.
        zlevel: Depth of the nozzle row in model vertical coordinate units.
        discharge: Volumetric air-flow rate. Accepts a scalar `float`,
            the string `"realtime"`, or a path to a `.bc` forcing file.

    Raises:
        ValidationError: If neither `locationFile` nor the full set of inline
            coordinate fields (`numCoordinates`, `xCoordinates`,
            `yCoordinates`) is provided.
        ValidationError: If `numCoordinates` does not match the length of
            `xCoordinates` or `yCoordinates`.

    Examples:
        - Build a block with inline coordinates and inspect the geometry:
            ```python
            >>> from hydrolib.core.dflowfm.ext.models import BubbleScreen
            >>> screen = BubbleScreen(
            ...     id="air_curtain",
            ...     numcoordinates=2,
            ...     xcoordinates=[100.0, 200.0],
            ...     ycoordinates=[50.0, 50.0],
            ...     zlevel=-3.5,
            ...     discharge=1.2,
            ... )
            >>> screen.id
            'air_curtain'
            >>> screen.zlevel
            -3.5
            >>> screen.xcoordinates
            [100.0, 200.0]

            ```
        - Build a block that references an external polyline file:
            ```python
            >>> from pathlib import Path
            >>> from hydrolib.core.base.models import DiskOnlyFileModel
            >>> from hydrolib.core.dflowfm.ext.models import BubbleScreen
            >>> screen = BubbleScreen(
            ...     id="canal_lock",
            ...     locationfile=DiskOnlyFileModel(filepath=Path("lock_screen.pli")),
            ...     zlevel=-6.0,
            ...     discharge=0.8,
            ... )
            >>> screen.id
            'canal_lock'
            >>> str(screen.locationfile.filepath)
            'lock_screen.pli'

            ```
        - Omitting both location forms raises `ValidationError`:
            ```python
            >>> from pydantic import ValidationError
            >>> from hydrolib.core.dflowfm.ext.models import BubbleScreen
            >>> try:
            ...     BubbleScreen(id="bad", zlevel=-5.0, discharge=1.0)
            ... except ValidationError:
            ...     print("location is required")
            location is required

            ```

    See Also:
        SourceSink: Sister block for water and solute injection or extraction.
        ExtModel: The new-format external forcings file model that hosts
            `[BubbleScreen]` blocks under its `bubblescreen` list.

    """

    _header: Literal["BubbleScreen"] = "BubbleScreen"
    id: str = Field(alias="id")
    name: str = Field("", alias="name")
    locationfile: Annotated[
        DiskOnlyFileModel, BeforeValidator(set_default_disk_only_file_model)
    ] = Field(default_factory=lambda: DiskOnlyFileModel(None), alias="locationFile")
    numcoordinates: int | None = Field(None, alias="numCoordinates")
    xcoordinates: list[float] | None = Field(None, alias="xCoordinates")
    ycoordinates: list[float] | None = Field(None, alias="yCoordinates")
    zlevel: float = Field(alias="zLevel")
    discharge: ForcingData = Field(alias="discharge")

    def is_intermediate_link(self) -> bool:
        return True

    @field_validator("discharge", mode="before")
    @classmethod
    def resolve_forcing_reference(cls, v):
        """Resolve string/Path inputs to scalar, RealTime enum, or ForcingModel."""
        return _resolve_forcing_data(v)

    @model_validator(mode="before")
    @classmethod
    def validate_location_specification(cls, values):
        locationfile = values.get("locationfile", values.get("locationFile"))
        numcoordinates = values.get("numcoordinates", values.get("numCoordinates"))
        xcoordinates = values.get("xcoordinates", values.get("xCoordinates"))
        ycoordinates = values.get("ycoordinates", values.get("yCoordinates"))

        has_locationfile = locationfile is not None
        if hasattr(locationfile, "filepath"):
            has_locationfile = locationfile.filepath is not None

        all_inline_present = (
            numcoordinates is not None
            and xcoordinates is not None
            and ycoordinates is not None
        )
        lengths_match = all_inline_present and (
            _coordinate_length(xcoordinates)
            == _coordinate_length(ycoordinates)
            == int(numcoordinates)
        )

        if all_inline_present and not lengths_match:
            raise ValueError(
                f"`numCoordinates` ({int(numcoordinates)}) does not match the length of "
                f"`xCoordinates` ({_coordinate_length(xcoordinates)}) and `yCoordinates` "
                f"({_coordinate_length(ycoordinates)}) for the BubbleScreen block "
                f"`{values.get('id')}`."
            )

        if not (has_locationfile or lengths_match):
            raise ValueError(
                "Either `locationFile` or the combination of `numCoordinates`, "
                "`xCoordinates`, and `yCoordinates` must be provided "
                f"for the BubbleScreen block `{values.get('id')}`."
            )

        return values

    @model_validator(mode="before")
    @classmethod
    def validate_locationfile(cls, data: Any) -> Any:
        file_location = data.get("locationfile") or data.get("locationFile")
        data.pop("locationFile", None)

        if isinstance(file_location, (str, Path)):
            data["locationfile"] = DiskOnlyFileModel(file_location)
        else:
            data["locationfile"] = file_location
        return data

resolve_forcing_reference(v) classmethod

Resolve string/Path inputs to scalar, RealTime enum, or ForcingModel.

Source code in hydrolib/core/dflowfm/ext/models.py
742
743
744
745
746
@field_validator("discharge", mode="before")
@classmethod
def resolve_forcing_reference(cls, v):
    """Resolve string/Path inputs to scalar, RealTime enum, or ForcingModel."""
    return _resolve_forcing_data(v)

ExtGeneral

Bases: INIGeneral

The external forcing file's [General] section with file meta-data.

Source code in hydrolib/core/dflowfm/ext/models.py
1280
1281
1282
1283
1284
1285
class ExtGeneral(INIGeneral):
    """The external forcing file's `[General]` section with file meta-data."""

    _header: Literal["General"] = "General"
    fileversion: str = Field("3.00", alias="fileVersion")
    filetype: Literal["extForce"] = Field("extForce", alias="fileType")

ExtModel

Bases: INIModel

The overall external forcings model that contains the contents of one external forcings file (new format).

This model is typically referenced under a FMModel.external_forcing.extforcefilenew.

Attributes:

Name Type Description
general ExtGeneral

[General] block with file metadata.

boundary List[Boundary]

List of [Boundary] blocks for all boundary conditions.

lateral List[Lateral]

List of [Lateral] blocks for all lateral discharges.

sourcesink List[SourceSink]

List of [SourceSink] blocks for all source/sink terms.

bubblescreen List[BubbleScreen]

List of [BubbleScreen] blocks for all bubble screens (air curtains).

meteo List[Meteo]

List of [Meteo] blocks for legacy meteorological forcings. Deprecated: use spatial instead.

spatial List[Spatial]

List of [Spatial] blocks for spatial forcings (meteo, initial conditions, and spatial parameters).

Source code in hydrolib/core/dflowfm/ext/models.py
1288
1289
1290
1291
1292
1293
1294
1295
1296
1297
1298
1299
1300
1301
1302
1303
1304
1305
1306
1307
1308
1309
1310
1311
1312
1313
1314
1315
1316
1317
1318
1319
1320
1321
1322
1323
1324
1325
1326
1327
1328
1329
1330
1331
1332
1333
1334
1335
1336
1337
1338
1339
1340
1341
1342
1343
1344
1345
1346
1347
1348
1349
1350
1351
1352
1353
1354
1355
1356
1357
1358
1359
1360
1361
1362
1363
1364
1365
1366
1367
1368
1369
1370
1371
class ExtModel(INIModel):
    """
    The overall external forcings model that contains the contents of one external forcings file (new format).

    This model is typically referenced under a [FMModel][hydrolib.core.dflowfm.mdu.models.FMModel]`.external_forcing.extforcefilenew`.

    Attributes:
        general (ExtGeneral): `[General]` block with file metadata.
        boundary (List[Boundary]): List of `[Boundary]` blocks for all boundary conditions.
        lateral (List[Lateral]): List of `[Lateral]` blocks for all lateral discharges.
        sourcesink (List[SourceSink]): List of `[SourceSink]` blocks for all source/sink terms.
        bubblescreen (List[BubbleScreen]): List of `[BubbleScreen]` blocks for all bubble
            screens (air curtains).
        meteo (List[Meteo]): List of `[Meteo]` blocks for legacy meteorological forcings.
            Deprecated: use `spatial` instead.
        spatial (List[Spatial]): List of `[Spatial]` blocks for spatial forcings (meteo,
            initial conditions, and spatial parameters).
    """

    general: ExtGeneral = ExtGeneral()
    boundary: Annotated[List[Boundary], BeforeValidator(make_list)] = Field(
        default_factory=list
    )
    lateral: Annotated[List[Lateral], BeforeValidator(make_list)] = Field(
        default_factory=list
    )
    sourcesink: Annotated[List[SourceSink], BeforeValidator(make_list)] = Field(
        default_factory=list
    )
    bubblescreen: Annotated[List[BubbleScreen], BeforeValidator(make_list)] = Field(
        default_factory=list
    )
    meteo: Annotated[List[Meteo], BeforeValidator(make_list)] = Field(
        default_factory=list
    )
    spatial: Annotated[List[Spatial], BeforeValidator(make_list)] = Field(
        default_factory=list
    )
    serializer_config: INISerializerConfig = INISerializerConfig(
        section_indent=0, property_indent=0
    )

    @model_validator(mode="after")
    def _warn_on_meteo(self) -> "ExtModel":
        """Emit a DeprecationWarning when [Meteo] blocks are present.

        The `[Meteo]` block is superseded by the `[Spatial]` block. New models
        should use `ExtModel.spatial` instead of `ExtModel.meteo`.
        """
        if self.meteo:
            warnings.warn(
                "`ExtModel.meteo` is deprecated; use `ExtModel.spatial` instead. "
                "`[Meteo]` blocks should be replaced by `[Spatial]` blocks.",
                DeprecationWarning,
                stacklevel=2,
            )
        return self

    @property
    def n_forcing_blocks(self) -> int:
        """Total number of forcing blocks held across all block types.

        Sums every `[Boundary]`, `[Lateral]`, `[SourceSink]`, `[BubbleScreen]`,
        `[Meteo]` and `[Spatial]` block, whether produced by conversion or loaded
        from an existing file. Use this to decide whether the model has any content
        worth writing; counting the individual lists by hand is error-prone and has
        silently dropped block types before.
        """
        return (
            len(self.boundary)
            + len(self.lateral)
            + len(self.sourcesink)
            + len(self.bubblescreen)
            + len(self.meteo)
            + len(self.spatial)
        )

    @classmethod
    def _ext(cls) -> str:
        return ".ext"

    @classmethod
    def _filename(cls) -> str:
        return "bnd"

n_forcing_blocks property

Total number of forcing blocks held across all block types.

Sums every [Boundary], [Lateral], [SourceSink], [BubbleScreen], [Meteo] and [Spatial] block, whether produced by conversion or loaded from an existing file. Use this to decide whether the model has any content worth writing; counting the individual lists by hand is error-prone and has silently dropped block types before.

Lateral

Bases: CoordinateValidator, INIBasedModel

A [Lateral] block for use inside an external forcings file.

I.e., a ExtModel.

All lowercased attributes match with the lateral input as described in UM Sec.C.5.2.2.

Source code in hydrolib/core/dflowfm/ext/models.py
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
class Lateral(CoordinateValidator, INIBasedModel):
    """A `[Lateral]` block for use inside an external forcings file.

    I.e., a [ExtModel][hydrolib.core.dflowfm.ext.models.ExtModel].

    All lowercased attributes match with the lateral input as described in
    [UM Sec.C.5.2.2](https://content.oss.deltares.nl/delft3dfm1d2d/D-Flow_FM_User_Manual_1D2D.pdf#subsection.C.6.3.2).
    """

    _header: Literal["Lateral"] = "Lateral"
    id: str = Field(alias="id")
    name: str = Field("", alias="name")
    locationtype: str | None = Field(None, alias="locationType")
    nodeid: str | None = Field(None, alias="nodeId")
    branchid: str | None = Field(None, alias="branchId")
    chainage: float | None = Field(None, alias="chainage")
    numcoordinates: int | None = Field(None, alias="numCoordinates")
    xcoordinates: list[float] | None = Field(None, alias="xCoordinates")
    ycoordinates: list[float] | None = Field(None, alias="yCoordinates")
    locationfile: Annotated[
        DiskOnlyFileModel, BeforeValidator(set_default_disk_only_file_model)
    ] | None = Field(None, alias="locationFile")
    applytransport: int | None = Field(None, alias="applyTransport")
    discharge: ForcingData = Field(alias="discharge")

    def is_intermediate_link(self) -> bool:
        return True

    @field_validator("discharge", mode="before")
    @classmethod
    def validate_discharge(cls, v):
        return _resolve_forcing_data(v)

    @model_validator(mode="before")
    def validate_that_location_specification_is_correct(cls, values: Dict) -> Dict:
        """Validate that one complete location specification is provided for `[Lateral]`.

        Validation is delegated to :class:`LocationValidatorUtils` with
        `minimum_num_coordinates=1` and `validate_location_file=True`.
        With that configuration, the input is accepted when it matches one of
        these alternatives:

        - `locationFile`
        - `nodeId`
        - `branchId` together with `chainage`
        - `numCoordinates` + `xCoordinates` + `yCoordinates`

        The helper applies the detailed consistency checks for each alternative
        (including coordinate-length checks and `locationType` validation).
        """
        location_validator = LocationValidatorUtils(
            values,
            config=LocationValidationConfiguration(
                minimum_num_coordinates=1,
                validate_location_file=True
            ),
        )
        values = location_validator.validate()
        return values

    def _get_identifier(self, data: dict) -> Optional[str]:
        return data.get("id") or data.get("name")

    @field_validator("applytransport", mode="after")
    @classmethod
    def validate_applytransport(cls, v: Any) -> int | None:
        if v not in (0, 1, None):
            raise ValueError(f"applyTransport must be 0 or 1, got '{v}'.")
        return v

    @field_validator("locationtype", mode="before")
    @classmethod
    def validate_location_type(cls, v: Any) -> LocationType:
        return enum_value_parser(v, LocationType)

validate_that_location_specification_is_correct(values)

Validate that one complete location specification is provided for [Lateral].

Validation is delegated to :class:LocationValidatorUtils with minimum_num_coordinates=1 and validate_location_file=True. With that configuration, the input is accepted when it matches one of these alternatives:

  • locationFile
  • nodeId
  • branchId together with chainage
  • numCoordinates + xCoordinates + yCoordinates

The helper applies the detailed consistency checks for each alternative (including coordinate-length checks and locationType validation).

Source code in hydrolib/core/dflowfm/ext/models.py
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
@model_validator(mode="before")
def validate_that_location_specification_is_correct(cls, values: Dict) -> Dict:
    """Validate that one complete location specification is provided for `[Lateral]`.

    Validation is delegated to :class:`LocationValidatorUtils` with
    `minimum_num_coordinates=1` and `validate_location_file=True`.
    With that configuration, the input is accepted when it matches one of
    these alternatives:

    - `locationFile`
    - `nodeId`
    - `branchId` together with `chainage`
    - `numCoordinates` + `xCoordinates` + `yCoordinates`

    The helper applies the detailed consistency checks for each alternative
    (including coordinate-length checks and `locationType` validation).
    """
    location_validator = LocationValidatorUtils(
        values,
        config=LocationValidationConfiguration(
            minimum_num_coordinates=1,
            validate_location_file=True
        ),
    )
    values = location_validator.validate()
    return values

LateralError

Bases: Exception

LateralError.

Source code in hydrolib/core/dflowfm/ext/models.py
1406
1407
1408
1409
1410
1411
class LateralError(Exception):
    """LateralError."""

    def __init__(self, error_message: str):
        """Initialize with an error message."""
        super().__init__(error_message)

__init__(error_message)

Initialize with an error message.

Source code in hydrolib/core/dflowfm/ext/models.py
1409
1410
1411
def __init__(self, error_message: str):
    """Initialize with an error message."""
    super().__init__(error_message)

Meteo

Bases: SpatialForcingBase

A [Meteo] block for use inside an external forcings file.

I.e., a ExtModel.

All lowercased attributes match with the meteo input as described in UM Sec.C.5.2.3.

Source code in hydrolib/core/dflowfm/ext/models.py
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
class Meteo(SpatialForcingBase):
    """A `[Meteo]` block for use inside an external forcings file.

    I.e., a [ExtModel][hydrolib.core.dflowfm.ext.models.ExtModel].

    All lowercased attributes match with the meteo input as described in
    [UM Sec.C.5.2.3](https://content.oss.deltares.nl/delft3dfm1d2d/D-Flow_FM_User_Manual_1D2D.pdf#subsection.C.5.2.3).
    """

    class Comments(SpatialForcingBase.Comments):
        """Comments for the Meteo block fields.

        Inherits the shared descriptions from `SpatialForcingBase.Comments`; only
        the `forcing*` file keywords and the `extrapolationSearchRadius` wording are
        specific to this block.
        """

        forcingfile: Optional[str] = Field(
            "Name of file containing the forcing for this meteo quantity.",
            alias="forcingFile",
        )
        forcingfiletype: Optional[str] = Field(
            "Type of forcingFile.", alias="forcingFileType"
        )
        forcingvariablename: Optional[str] = Field(
            "Variable name used in forcingfile associated with this forcing. See UM Section C.5.3",
            alias="forcingVariableName",
        )
        targetmaskfile: Optional[str] = Field(
            "Name of <*.pol> file to be used as mask. Grid parts inside any polygon will receive the meteo forcing.",
            alias="targetMaskFile",
        )
        extrapolationsearchradius: Optional[str] = Field(
            "Maximum search radius for nearest neighbor extrapolation in space.",
            alias="extrapolationSearchRadius",
        )

    comments: Comments = Comments()

    _header: Literal["Meteo"] = "Meteo"
    quantity: str = Field(alias="quantity")
    forcingfile: Union[TimModel, ForcingModel, DiskOnlyFileModel, PolyFile] = Field(
        alias="forcingFile"
    )
    forcingvariablename: str | None = Field(None, alias="forcingVariableName")
    forcingfiletype: MeteoForcingFileType = Field(alias="forcingFileType")
    targetmaskfile: PolyFile | None = Field(None, alias="targetMaskFile")
    targetmaskinvert: bool | None = Field(None, alias="targetMaskInvert")
    interpolationmethod: MeteoInterpolationMethod | None = Field(
        None, alias="interpolationMethod"
    )
    operand: Operand | None = Field(Operand.override.value, alias="operand")
    extrapolationallowed: bool | None = Field(None, alias="extrapolationAllowed")
    extrapolationsearchradius: float | None = Field(
        None, alias="extrapolationSearchRadius"
    )
    averagingtype: int | None = Field(None, alias="averagingType")
    averagingnummin: PositiveInt | None = Field(None, alias="averagingNumMin")
    averagingpercentile: float | None = Field(None, alias="averagingPercentile")

    # Deprecated camelCase aliases — intentional case clash with the fields above; remove in 2.0.0 (docs/migration.md).
    forcingVariableName = DeprecatedAttributeAlias(  # NOSONAR S1845
        "forcingvariablename", removed_in="2.0.0", since="1.1.0"
    )
    extrapolationAllowed = DeprecatedAttributeAlias(  # NOSONAR S1845
        "extrapolationallowed", removed_in="2.0.0", since="1.1.0"
    )
    extrapolationSearchRadius = DeprecatedAttributeAlias(  # NOSONAR S1845
        "extrapolationsearchradius", removed_in="2.0.0", since="1.1.0"
    )
    averagingType = DeprecatedAttributeAlias(  # NOSONAR S1845
        "averagingtype", removed_in="2.0.0", since="1.1.0"
    )
    averagingNumMin = DeprecatedAttributeAlias(  # NOSONAR S1845
        "averagingnummin", removed_in="2.0.0", since="1.1.0"
    )
    averagingPercentile = DeprecatedAttributeAlias(  # NOSONAR S1845
        "averagingpercentile", removed_in="2.0.0", since="1.1.0"
    )

    @model_validator(mode="before")
    @classmethod
    def choose_file_model(cls, values: Dict[str, Any]) -> Dict[str, Any]:
        """Select the right class for the forcingFile parameter based on forcingFileType.

        Uses the shared `SpatialForcingBase._resolve_file_models` helper against
        this block's `forcingFile`/`forcingFileType` keywords.
        """
        return cls._resolve_file_models(
            values,
            ("forcingfile", "forcingFile"),
            ("forcingfiletype", "forcingFileType"),
        )

    def is_intermediate_link(self) -> bool:
        return True

    @field_validator("forcingfiletype", mode="before")
    @classmethod
    def forcingfiletype_validator(cls, v):
        return enum_value_parser(v, MeteoForcingFileType)

Comments

Bases: Comments

Comments for the Meteo block fields.

Inherits the shared descriptions from SpatialForcingBase.Comments; only the forcing* file keywords and the extrapolationSearchRadius wording are specific to this block.

Source code in hydrolib/core/dflowfm/ext/models.py
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
class Comments(SpatialForcingBase.Comments):
    """Comments for the Meteo block fields.

    Inherits the shared descriptions from `SpatialForcingBase.Comments`; only
    the `forcing*` file keywords and the `extrapolationSearchRadius` wording are
    specific to this block.
    """

    forcingfile: Optional[str] = Field(
        "Name of file containing the forcing for this meteo quantity.",
        alias="forcingFile",
    )
    forcingfiletype: Optional[str] = Field(
        "Type of forcingFile.", alias="forcingFileType"
    )
    forcingvariablename: Optional[str] = Field(
        "Variable name used in forcingfile associated with this forcing. See UM Section C.5.3",
        alias="forcingVariableName",
    )
    targetmaskfile: Optional[str] = Field(
        "Name of <*.pol> file to be used as mask. Grid parts inside any polygon will receive the meteo forcing.",
        alias="targetMaskFile",
    )
    extrapolationsearchradius: Optional[str] = Field(
        "Maximum search radius for nearest neighbor extrapolation in space.",
        alias="extrapolationSearchRadius",
    )

choose_file_model(values) classmethod

Select the right class for the forcingFile parameter based on forcingFileType.

Uses the shared SpatialForcingBase._resolve_file_models helper against this block's forcingFile/forcingFileType keywords.

Source code in hydrolib/core/dflowfm/ext/models.py
973
974
975
976
977
978
979
980
981
982
983
984
985
@model_validator(mode="before")
@classmethod
def choose_file_model(cls, values: Dict[str, Any]) -> Dict[str, Any]:
    """Select the right class for the forcingFile parameter based on forcingFileType.

    Uses the shared `SpatialForcingBase._resolve_file_models` helper against
    this block's `forcingFile`/`forcingFileType` keywords.
    """
    return cls._resolve_file_models(
        values,
        ("forcingfile", "forcingFile"),
        ("forcingfiletype", "forcingFileType"),
    )

ParameterFieldError

Bases: Exception

ParameterFieldError.

Source code in hydrolib/core/dflowfm/ext/models.py
1398
1399
1400
1401
1402
1403
class ParameterFieldError(Exception):
    """ParameterFieldError."""

    def __init__(self, error_message: str):
        """Initialize with an error message."""
        super().__init__(error_message)

__init__(error_message)

Initialize with an error message.

Source code in hydrolib/core/dflowfm/ext/models.py
1401
1402
1403
def __init__(self, error_message: str):
    """Initialize with an error message."""
    super().__init__(error_message)

SourceSink

Bases: CoordinateValidator, INIBasedModel

A [SourceSink] block for use inside an external forcings file.

I.e., a ExtModel.

All lowercased attributes match with the source-sink input as described in UM Sec.C.5.2.4.

Source code in hydrolib/core/dflowfm/ext/models.py
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
class SourceSink(CoordinateValidator, INIBasedModel):
    """A `[SourceSink]` block for use inside an external forcings file.

    I.e., a [ExtModel][hydrolib.core.dflowfm.ext.models.SourceSink].

    All lowercased attributes match with the source-sink input as described in
    [UM Sec.C.5.2.4](https://content.oss.deltares.nl/delft3dfm1d2d/D-Flow_FM_User_Manual_1D2D.pdf#subsection.C.5.2.4).
    """
    model_config = ConfigDict(extra="allow")

    _header: Literal["SourceSink"] = "SourceSink"
    id: str = Field(alias="id")
    name: str = Field("", alias="name")
    locationfile: Optional[DiskOnlyFileModel] = Field(
        default_factory=lambda: DiskOnlyFileModel(None), alias="locationFile"
    )

    numcoordinates: Optional[int] = Field(None, alias="numCoordinates")
    xcoordinates: Optional[List[float]] = Field(None, alias="xCoordinates")
    ycoordinates: Optional[List[float]] = Field(None, alias="yCoordinates")

    zsource: Optional[Union[float, List[float]]] = Field(None, alias="zSource")
    zsink: Optional[Union[float, List[float]]] = Field(None, alias="zSink")
    area: Optional[float] = Field(None, alias="Area")

    discharge: ForcingData = Field(alias="discharge")
    salinity: Optional[ForcingData] = Field(None, alias="salinity")
    temperature: Optional[ForcingData] = Field(None, alias="temperature")

    def is_intermediate_link(self) -> bool:
        return True

    @field_validator("zsource", "zsink", mode="before")
    @classmethod
    def split_z_values(cls, v):
        """Parse `zSource`/`zSink` from a space-separated string.

        Range sources are written as e.g. `zSource = -7.5 -3.01` in the ext
        file; the parser delivers this as a single string that Pydantic
        cannot coerce to `float` or `List[float]`. An empty or whitespace-only
        string (e.g. `zSource =`) is treated as "not provided" (`None`),
        consistent with `str_is_empty_or_none`.
        """
        result = v
        if isinstance(v, str):
            if str_is_empty_or_none(v):
                result = None
            else:
                parts = v.split()
                if len(parts) == 1:
                    result = float(parts[0])
                else:
                    result = [float(p) for p in parts]
        return result

    @field_validator(
        "discharge", "salinity", "temperature", mode="before"
    )
    @classmethod
    def validate_forcing_data(cls, v):
        return _resolve_forcing_data(v)

    @model_validator(mode="before")
    @classmethod
    def _resolve_dynamic_forcing_deltas(cls, values: Any) -> Any:
        """Apply `_resolve_forcing_data` to dynamic `tracer<...>Delta`/`sedFrac<...>Delta` keys.

        Also renames legacy `salinitydelta`/`temperaturedelta` keys (produced by
        the INI parser from old-format `salinityDelta`/`temperatureDelta`) to the
        current field names `salinity`/`temperature` for backward compatibility.

        Per D-Flow FM User Manual Table C.8 (§C.6.2.4), `tracer<name>Delta` and
        `sedFrac<name>Delta` accept a scalar Double or the name of a `.bc`
        time-series file. The first-class `discharge`/`salinity`/
        `temperature` fields are already handled by `validate_forcing_data`;
        this validator extends the same coercion to the dynamic Delta-suffix
        fields that arrive via `extra="allow"`.

        Legacy dynamic fields (`initialtracer_*`, `tracerbnd*`, `sedfracbnd_*`,
        `initialsedfrac_*`) do not end with `delta` and are left untouched.
        """
        if isinstance(values, dict):
            # Migrate legacy salinityDelta/temperatureDelta keys (any casing) from old-format ext files.
            lowercase_map = {k.lower(): k for k in values}
            if "salinitydelta" in lowercase_map and "salinity" not in values:
                values["salinity"] = values.pop(lowercase_map["salinitydelta"])
            if "temperaturedelta" in lowercase_map and "temperature" not in values:
                values["temperature"] = values.pop(lowercase_map["temperaturedelta"])
            for key in values:
                if _is_dynamic_forcing_delta_key(key):
                    values[key] = _resolve_forcing_data(
                        values[key], allow_realtime=False
                    )
        return values

    @classmethod
    def _exclude_from_validation(cls, input_data: dict | None = None) -> Set:
        input_data = input_data or {}
        fields = cls.model_fields
        dynamic_fields = input_data.get(_DYNAMIC_FIELDS_KEY) or []
        unknown_keywords = [
            key
            for key in input_data.keys()
            if key not in fields
            and (
                key.startswith(SOURCE_SINKS_QUANTITIES_VALID_PREFIXES)
                or _is_dynamic_forcing_delta_key(key)
                or key in dynamic_fields
                or key == _DYNAMIC_FIELDS_KEY
            )
        ]
        return set(unknown_keywords)

    def __init__(self, dynamic_fields: Optional[List[str]] = None, **data):
        """Initialize SourceSink and set dynamic fields as instance attributes.

        When `dynamic_fields` is provided, exactly those names (whose values are passed
        in `data`) are attached onto this instance. When it is omitted, the legacy
        behaviour applies: every key in `data` that starts with one of
        `SOURCE_SINKS_QUANTITIES_VALID_PREFIXES` is attached. In both cases the values
        are stored as-is on the instance (in `model_extra`); no coercion is applied yet.
        Args:
            dynamic_fields: Names of extra fields (whose values are passed in `data`)
                to attach onto this instance. If `None`, the prefix-based detection is
                used instead.
            **data: The regular SourceSink field values, plus the values for any dynamic
                fields.
        """
        # Thread the dynamic field names through validation so that
        # `_exclude_from_validation` can whitelist them as known keywords.
        if dynamic_fields is not None:
            data[_DYNAMIC_FIELDS_KEY] = list(dynamic_fields)

        super().__init__(**data)

        # Drop the reserved key so it does not linger as an instance attribute.
        if self.__pydantic_extra__ is not None:
            self.__pydantic_extra__.pop(_DYNAMIC_FIELDS_KEY, None)

    @model_validator(mode="before")
    def validate_location_specification(cls, values):
        """
        Ensures that either `locationfile` or a valid set of coordinates is provided.

         This validation enforces that at least one of the following conditions is met:
         1. `locationfile` is provided.
         2. The combination of `numcoordinates`, `xcoordinates`, and `ycoordinates` is valid:
             - `xcoordinates` and `ycoordinates` must be lists of equal length.
             - The length of `xcoordinates` and `ycoordinates` must match `numcoordinates`.

         Raises:
             ValueError: If neither `locationfile` nor a valid coordinate set is provided.

         Returns:
             Dict: The validated input values.
        """
        locationfile = values.get("locationfile", values.get("locationFile"))

        numcoordinates = values.get("numcoordinates", values.get("numCoordinates"))
        xcoordinates = values.get("xcoordinates", values.get("xCoordinates"))
        ycoordinates = values.get("ycoordinates", values.get("yCoordinates"))

        has_locationfile = locationfile is not None
        has_coordinates = (
            numcoordinates is not None
            and xcoordinates is not None
            and ycoordinates is not None
            and _coordinate_length(xcoordinates)
            == _coordinate_length(ycoordinates)
            == int(numcoordinates)
        )

        if not (has_locationfile or has_coordinates):
            raise ValueError(
                "Either `locationFile` or the combination of `numCoordinates`, `xCoordinates`, and `yCoordinates` "
                f"must be provided. for the SourceSink block `{values.get('id')}`."
            )

        return values

    @model_validator(mode="before")
    @classmethod
    def validate_locationfile_z_conflict(cls, values):
        """Reject `.pliz` locationFile combined with explicit `zSource`/`zSink`.

        A `.pliz` polyline already encodes vertical placement via its column
        count (3 or 5 cols). Combining it with explicit `zSource`/`zSink` is
        ambiguous and is rejected by the kernel. A plain `.pli` polyline
        carries no z, so combining it with explicit `zSource`/`zSink` is
        allowed (and required for a coupled source-sink specified via
        `locationFile`).
        """
        locationfile = values.get("locationfile", values.get("locationFile"))
        zsource = values.get("zsource", values.get("zSource"))
        zsink = values.get("zsink", values.get("zSink"))

        locationfile_path = None
        if hasattr(locationfile, "filepath"):
            locationfile_path = locationfile.filepath
        elif isinstance(locationfile, (str, Path)):
            locationfile_path = Path(locationfile)

        file_has_z = locationfile_path is not None and str(
            locationfile_path
        ).lower().endswith(".pliz")

        conflict_field = None
        if file_has_z and zsource is not None:
            conflict_field = "zSource"
        elif file_has_z and zsink is not None:
            conflict_field = "zSink"

        if conflict_field is not None:
            raise ValueError(
                f"`locationFile` (a `.pliz` polyline) cannot be combined with explicit "
                f"`{conflict_field}` for the SourceSink block `{values.get('id')}`. "
                f"The `.pliz` file already encodes vertical placement. "
                f"Either remove `{conflict_field}` or use a `.pli` locationFile "
                f"(no z column) so the explicit z-values are unambiguous."
            )

        return values

    @model_validator(mode="before")
    @classmethod
    def validate_locationfile(cls, data: Any) -> Any:
        file_location = data.get("locationfile") or data.get("locationFile")
        data.pop("locationFile", None)  # Remove alias if present

        # Convert string to DiskOnlyFileModel if needed
        if isinstance(file_location, (str, Path)):
            data["locationfile"] = DiskOnlyFileModel(file_location)
        else:
            data["locationfile"] = file_location
        return data

__init__(dynamic_fields=None, **data)

Initialize SourceSink and set dynamic fields as instance attributes.

When dynamic_fields is provided, exactly those names (whose values are passed in data) are attached onto this instance. When it is omitted, the legacy behaviour applies: every key in data that starts with one of SOURCE_SINKS_QUANTITIES_VALID_PREFIXES is attached. In both cases the values are stored as-is on the instance (in model_extra); no coercion is applied yet. Args: dynamic_fields: Names of extra fields (whose values are passed in data) to attach onto this instance. If None, the prefix-based detection is used instead. **data: The regular SourceSink field values, plus the values for any dynamic fields.

Source code in hydrolib/core/dflowfm/ext/models.py
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
def __init__(self, dynamic_fields: Optional[List[str]] = None, **data):
    """Initialize SourceSink and set dynamic fields as instance attributes.

    When `dynamic_fields` is provided, exactly those names (whose values are passed
    in `data`) are attached onto this instance. When it is omitted, the legacy
    behaviour applies: every key in `data` that starts with one of
    `SOURCE_SINKS_QUANTITIES_VALID_PREFIXES` is attached. In both cases the values
    are stored as-is on the instance (in `model_extra`); no coercion is applied yet.
    Args:
        dynamic_fields: Names of extra fields (whose values are passed in `data`)
            to attach onto this instance. If `None`, the prefix-based detection is
            used instead.
        **data: The regular SourceSink field values, plus the values for any dynamic
            fields.
    """
    # Thread the dynamic field names through validation so that
    # `_exclude_from_validation` can whitelist them as known keywords.
    if dynamic_fields is not None:
        data[_DYNAMIC_FIELDS_KEY] = list(dynamic_fields)

    super().__init__(**data)

    # Drop the reserved key so it does not linger as an instance attribute.
    if self.__pydantic_extra__ is not None:
        self.__pydantic_extra__.pop(_DYNAMIC_FIELDS_KEY, None)

split_z_values(v) classmethod

Parse zSource/zSink from a space-separated string.

Range sources are written as e.g. zSource = -7.5 -3.01 in the ext file; the parser delivers this as a single string that Pydantic cannot coerce to float or List[float]. An empty or whitespace-only string (e.g. zSource =) is treated as "not provided" (None), consistent with str_is_empty_or_none.

Source code in hydrolib/core/dflowfm/ext/models.py
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
@field_validator("zsource", "zsink", mode="before")
@classmethod
def split_z_values(cls, v):
    """Parse `zSource`/`zSink` from a space-separated string.

    Range sources are written as e.g. `zSource = -7.5 -3.01` in the ext
    file; the parser delivers this as a single string that Pydantic
    cannot coerce to `float` or `List[float]`. An empty or whitespace-only
    string (e.g. `zSource =`) is treated as "not provided" (`None`),
    consistent with `str_is_empty_or_none`.
    """
    result = v
    if isinstance(v, str):
        if str_is_empty_or_none(v):
            result = None
        else:
            parts = v.split()
            if len(parts) == 1:
                result = float(parts[0])
            else:
                result = [float(p) for p in parts]
    return result

validate_location_specification(values)

Ensures that either locationfile or a valid set of coordinates is provided.

This validation enforces that at least one of the following conditions is met: 1. locationfile is provided. 2. The combination of numcoordinates, xcoordinates, and ycoordinates is valid: - xcoordinates and ycoordinates must be lists of equal length. - The length of xcoordinates and ycoordinates must match numcoordinates.

Raises: ValueError: If neither locationfile nor a valid coordinate set is provided.

Returns: Dict: The validated input values.

Source code in hydrolib/core/dflowfm/ext/models.py
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
@model_validator(mode="before")
def validate_location_specification(cls, values):
    """
    Ensures that either `locationfile` or a valid set of coordinates is provided.

     This validation enforces that at least one of the following conditions is met:
     1. `locationfile` is provided.
     2. The combination of `numcoordinates`, `xcoordinates`, and `ycoordinates` is valid:
         - `xcoordinates` and `ycoordinates` must be lists of equal length.
         - The length of `xcoordinates` and `ycoordinates` must match `numcoordinates`.

     Raises:
         ValueError: If neither `locationfile` nor a valid coordinate set is provided.

     Returns:
         Dict: The validated input values.
    """
    locationfile = values.get("locationfile", values.get("locationFile"))

    numcoordinates = values.get("numcoordinates", values.get("numCoordinates"))
    xcoordinates = values.get("xcoordinates", values.get("xCoordinates"))
    ycoordinates = values.get("ycoordinates", values.get("yCoordinates"))

    has_locationfile = locationfile is not None
    has_coordinates = (
        numcoordinates is not None
        and xcoordinates is not None
        and ycoordinates is not None
        and _coordinate_length(xcoordinates)
        == _coordinate_length(ycoordinates)
        == int(numcoordinates)
    )

    if not (has_locationfile or has_coordinates):
        raise ValueError(
            "Either `locationFile` or the combination of `numCoordinates`, `xCoordinates`, and `yCoordinates` "
            f"must be provided. for the SourceSink block `{values.get('id')}`."
        )

    return values

validate_locationfile_z_conflict(values) classmethod

Reject .pliz locationFile combined with explicit zSource/zSink.

A .pliz polyline already encodes vertical placement via its column count (3 or 5 cols). Combining it with explicit zSource/zSink is ambiguous and is rejected by the kernel. A plain .pli polyline carries no z, so combining it with explicit zSource/zSink is allowed (and required for a coupled source-sink specified via locationFile).

Source code in hydrolib/core/dflowfm/ext/models.py
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
@model_validator(mode="before")
@classmethod
def validate_locationfile_z_conflict(cls, values):
    """Reject `.pliz` locationFile combined with explicit `zSource`/`zSink`.

    A `.pliz` polyline already encodes vertical placement via its column
    count (3 or 5 cols). Combining it with explicit `zSource`/`zSink` is
    ambiguous and is rejected by the kernel. A plain `.pli` polyline
    carries no z, so combining it with explicit `zSource`/`zSink` is
    allowed (and required for a coupled source-sink specified via
    `locationFile`).
    """
    locationfile = values.get("locationfile", values.get("locationFile"))
    zsource = values.get("zsource", values.get("zSource"))
    zsink = values.get("zsink", values.get("zSink"))

    locationfile_path = None
    if hasattr(locationfile, "filepath"):
        locationfile_path = locationfile.filepath
    elif isinstance(locationfile, (str, Path)):
        locationfile_path = Path(locationfile)

    file_has_z = locationfile_path is not None and str(
        locationfile_path
    ).lower().endswith(".pliz")

    conflict_field = None
    if file_has_z and zsource is not None:
        conflict_field = "zSource"
    elif file_has_z and zsink is not None:
        conflict_field = "zSink"

    if conflict_field is not None:
        raise ValueError(
            f"`locationFile` (a `.pliz` polyline) cannot be combined with explicit "
            f"`{conflict_field}` for the SourceSink block `{values.get('id')}`. "
            f"The `.pliz` file already encodes vertical placement. "
            f"Either remove `{conflict_field}` or use a `.pli` locationFile "
            f"(no z column) so the explicit z-values are unambiguous."
        )

    return values

SourceSinkError

Bases: Exception

SourceSinkError.

Source code in hydrolib/core/dflowfm/ext/models.py
1374
1375
1376
1377
1378
1379
class SourceSinkError(Exception):
    """SourceSinkError."""

    def __init__(self, error_message: str):
        """Initialize with an error message."""
        super().__init__(error_message)

__init__(error_message)

Initialize with an error message.

Source code in hydrolib/core/dflowfm/ext/models.py
1377
1378
1379
def __init__(self, error_message: str):
    """Initialize with an error message."""
    super().__init__(error_message)

Spatial

Bases: SpatialForcingBase, LocationTypeDataFileTypeValidators

A [Spatial] block for use inside an external forcings file.

I.e., a ExtModel.

This block replaces both the legacy [Meteo] block (for meteorological forcings) and the [Initial] / [Parameter] blocks in inifield files (for initial conditions and spatial parameters).

All lowercased attributes match with the spatial input as described in UM Sec.C.5.2.3.

Source code in hydrolib/core/dflowfm/ext/models.py
 996
 997
 998
 999
1000
1001
1002
1003
1004
1005
1006
1007
1008
1009
1010
1011
1012
1013
1014
1015
1016
1017
1018
1019
1020
1021
1022
1023
1024
1025
1026
1027
1028
1029
1030
1031
1032
1033
1034
1035
1036
1037
1038
1039
1040
1041
1042
1043
1044
1045
1046
1047
1048
1049
1050
1051
1052
1053
1054
1055
1056
1057
1058
1059
1060
1061
1062
1063
1064
1065
1066
1067
1068
1069
1070
1071
1072
1073
1074
1075
1076
1077
1078
1079
1080
1081
1082
1083
1084
1085
1086
1087
1088
1089
1090
1091
1092
1093
1094
1095
1096
1097
1098
1099
1100
1101
1102
1103
1104
1105
1106
1107
1108
1109
1110
1111
1112
1113
1114
1115
1116
1117
1118
1119
1120
1121
1122
1123
1124
1125
1126
1127
1128
1129
1130
1131
1132
1133
1134
1135
1136
1137
1138
1139
1140
1141
1142
1143
1144
1145
1146
1147
1148
1149
1150
1151
1152
1153
1154
1155
1156
1157
1158
1159
1160
1161
1162
1163
1164
1165
1166
1167
1168
1169
1170
1171
1172
1173
1174
1175
1176
1177
1178
1179
1180
1181
1182
1183
1184
1185
1186
1187
1188
1189
1190
1191
1192
1193
1194
1195
1196
1197
1198
1199
1200
1201
1202
1203
1204
1205
1206
1207
1208
1209
1210
1211
1212
1213
1214
1215
1216
1217
1218
1219
1220
1221
1222
1223
1224
1225
1226
1227
1228
1229
1230
1231
1232
1233
1234
1235
1236
1237
1238
1239
1240
1241
1242
1243
1244
1245
1246
1247
1248
1249
1250
1251
1252
1253
1254
1255
1256
1257
1258
1259
1260
1261
1262
1263
1264
1265
1266
1267
1268
1269
1270
1271
1272
1273
1274
1275
1276
1277
class Spatial(SpatialForcingBase, LocationTypeDataFileTypeValidators):
    """A `[Spatial]` block for use inside an external forcings file.

    I.e., a [ExtModel][hydrolib.core.dflowfm.ext.models.ExtModel].

    This block replaces both the legacy `[Meteo]` block (for meteorological
    forcings) and the `[Initial]` / `[Parameter]` blocks in inifield files
    (for initial conditions and spatial parameters).

    All lowercased attributes match with the spatial input as described in
    [UM Sec.C.5.2.3](https://content.oss.deltares.nl/delft3dfm1d2d/D-Flow_FM_User_Manual_1D2D.pdf#subsection.C.5.2.3).
    """

    class Comments(SpatialForcingBase.Comments):
        """Comments for the Spatial block fields.

        Inherits the shared descriptions from `SpatialForcingBase.Comments`; only
        the `data*` file keywords, the `extrapolationSearchRadius` wording, and the
        fields unique to the Spatial block are declared here.
        """

        datafile: str | None = Field(
            "Name of file containing the data for this spatial quantity.",
            alias="dataFile",
        )
        datafiletype: str | None = Field(
            "Type of dataFile.", alias="dataFileType"
        )
        datavariablename: str | None = Field(
            "Variable name used in dataFile associated with this quantity.",
            alias="dataVariableName",
        )
        targetmaskfile: str | None = Field(
            "Name of <*.pol> file to be used as mask. Grid parts inside any polygon will receive the spatial forcing.",
            alias="targetMaskFile",
        )
        extrapolationsearchradius: str | None = Field(
            "Maximum search radius for nearest neighbour extrapolation in space.",
            alias="extrapolationSearchRadius",
        )
        averagingtype: str | None = Field(
            "Type of averaging, if interpolationMethod=averaging.",
            alias="averagingType",
        )
        averagingrelsize: str | None = Field(
            "Relative search cell size for averaging.", alias="averagingRelSize"
        )
        averagingnummin: str | None = Field(
            "Minimum number of points in averaging. Must be ≥ 1.",
            alias="averagingNumMin",
        )
        averagingpercentile: str | None = Field(
            "Percentile value for which data values to include in averaging. 0.0 means off.",
            alias="averagingPercentile",
        )
        locationtype: str | None = Field(
            "Target location of interpolation.", alias="locationType"
        )
        datavalue: str | None = Field(
            "Constant value to be set inside all model points inside the polygon, "
            "used when no dataFile/dataFileType is specified. "
            "Requires targetMaskFile=*.pol and interpolationMethod=constant.",
            alias="dataValue",
        )
        frictiontype: str | None = Field(
            "Only for quantity=frictionCoefficient. The friction type.", alias="frictionType",
        )
        tracerfallvelocity: str | None = Field(
            "Only for initialtracer<tracername>. Fall velocity of the tracer.",
            alias="tracerFallVelocity",
        )
        tracerdecaytime: str | None = Field(
            "Only for initialtracer<tracername>. Decay time of the tracer.",
            alias="tracerDecayTime",
        )
        targetlayer: str | None = Field(
            "Target layer for the data: bottom, all, or a positive layer number.",
            alias="targetLayer",
        )

    comments: Comments = Comments()

    _header: Literal["Spatial"] = "Spatial"
    quantity: str = Field(alias="quantity")
    datafile: TimModel | ForcingModel | DiskOnlyFileModel | PolyFile | None = Field(
        None, alias="dataFile"
    )
    datafiletype: DataFileType | None = Field(None, alias="dataFileType")
    datavariablename: str | None = Field(None, alias="dataVariableName")
    targetmaskfile: PolyFile | DiskOnlyFileModel | None = Field(None, alias="targetMaskFile")
    targetmaskinvert: bool | None = Field(None, alias="targetMaskInvert")
    interpolationmethod: InterpolationMethod | None = Field(
        None, alias="interpolationMethod"
    )
    operand: Operand | None = Field(Operand.override.value, alias="operand")
    extrapolationallowed: bool | None = Field(False, alias="extrapolationAllowed")
    extrapolationsearchradius: float | None = Field(
        None, alias="extrapolationSearchRadius"
    )
    averagingtype: AveragingType | None = Field(None, alias="averagingType")
    averagingrelsize: NonNegativeFloat | None = Field(None, alias="averagingRelSize")
    averagingnummin: PositiveInt | None = Field(None, alias="averagingNumMin")
    averagingpercentile: NonNegativeFloat | None = Field(None, alias="averagingPercentile")
    locationtype: LocationType | None = Field(
        LocationType.all.value, alias="locationType"
    )
    datavalue: float | None = Field(None, alias="dataValue")
    frictiontype: str | None = Field(None, alias="frictionType")
    tracerfallvelocity: float | None = Field(None, alias="tracerFallVelocity")
    tracerdecaytime: float | None = Field(None, alias="tracerDecayTime")
    targetlayer: TargetLayer | int | None = Field(None, alias="targetLayer")

    @classmethod
    def _normalize_spatial_keys(cls, values: Dict) -> Dict:
        """Normalize camelCase aliases and coerce any unresolved datafile to DiskOnlyFileModel.

        `dataFile` is normally resolved to its concrete file model earlier in the
        validation flow (see `_resolve_file_models`). This fallback only wraps a
        path that is still raw (e.g. no `dataFileType` was supplied) so that the
        field always holds a file model rather than a bare string.
        """
        data_file = values.get("datafile") or values.get("dataFile")
        if isinstance(data_file, (str, Path)):
            data_file = DiskOnlyFileModel(data_file)
            values.pop("dataFile", None)
            values["datafile"] = data_file

        if "dataValue" in values and "datavalue" not in values:
            values["datavalue"] = values.pop("dataValue")
        if "targetMaskFile" in values and "targetmaskfile" not in values:
            values["targetmaskfile"] = values.pop("targetMaskFile")
        return values

    @classmethod
    def _validate_datavalue_path(
        cls, values: Dict, has_datafile: bool, has_datafiletype: bool
    ) -> None:
        """Validate the `dataValue` usage path (constant value inside polygon)."""
        if has_datafile or has_datafiletype:
            raise ValueError(
                "When 'dataValue' is provided, 'dataFile' and 'dataFileType' must not be specified."
            )
        interp = values.get("interpolationmethod") or values.get("interpolationMethod")
        if interp is None:
            values["interpolationmethod"] = InterpolationMethod.constant
        elif str(interp).lower() != str(InterpolationMethod.constant).lower():
            raise ValueError(
                f"When 'dataValue' is provided, 'interpolationMethod' must be "
                f"'{InterpolationMethod.constant}', got '{interp}'."
            )

    @classmethod
    def _validate_datafile_path(
        cls, values: Dict, has_datafile: bool, has_datafiletype: bool
    ) -> None:
        """Validate the `dataFile` usage path and emit deprecation warning when needed."""
        if not has_datafile:
            raise ValueError("'dataFile' is required when 'dataValue' is not specified.")
        if not has_datafiletype:
            raise ValueError("'dataFileType' is required when 'dataValue' is not specified.")

        raw_filetype = values.get("datafiletype") or values.get("dataFileType")
        quantity = values.get("quantity") or ""
        if (
            raw_filetype is not None
            and str(raw_filetype).lower() == DataFileType.polygon
            and not str(quantity).startswith("initialvertical")
        ):
            warnings.warn(
                "Using dataFileType=polygon for 'inside polygon' data is deprecated. "
                "Use dataValue + targetMaskFile=<*.pol> + interpolationMethod=constant instead. "
                "The polygon dataFileType remains supported only for initialvertical* quantities "
                "(e.g. initialverticalsalinityprofile).",
                DeprecationWarning,
                stacklevel=2,
            )

    @classmethod
    def _process_section_values(cls, values):
        """Flatten a Section object into a dictionary of raw values.

        The raw `dataFile` value is left as a path/string so that the subsequent
        `_resolve_file_models` step can select the concrete file model from
        `dataFileType` (rather than being forced to `DiskOnlyFileModel` here).

        Args:
            values: The values to process, which may be a Section object or a dictionary.

        Returns:
            A dictionary containing the processed values.
        """
        return cls._convert_section_to_dict(values)

    @model_validator(mode="before")
    @classmethod
    def validate_datavalue_or_datafile(cls, values: Dict) -> Dict:
        """Validates the two mutually exclusive usage paths of a Spatial block.

        When `dataValue` is provided the block describes a constant value applied
        inside a polygon mask.  In this mode:
        - `dataFile` and `dataFileType` must **not** be specified.
        - `targetMaskFile` (a `.pol` file) is optional but often used.
        - `interpolationMethod` must be `constant` (set automatically when omitted).

        When `dataValue` is absent, `dataFile` and `dataFileType` are both
        required.

        Note: using `dataFileType=polygon` for "inside polygon" initial-condition
        data is **deprecated**.  Use `dataValue` + `targetMaskFile=*.pol` +
        `interpolationMethod=constant` instead.  The `polygon` dataFileType
        remains supported for quantities such as `initialvertical*` (e.g.
        `initialverticalsalinityprofile`) that use polygon files for a different
        purpose and have no new alternative yet.

        The data file is resolved to its concrete file model only on the `dataFile`
        path (`dataValue` absent). This prevents an invalid combination
        (`dataValue` together with `dataFile`) from parsing a file before the
        mutual-exclusion check rejects it, so callers get the exclusion error rather
        than a file-parse error.
        """
        values = cls._process_section_values(values)

        on_datavalue_path = (
            values.get("datavalue") is not None or values.get("dataValue") is not None
        )
        if not on_datavalue_path:
            values = cls._resolve_file_models(
                values,
                ("datafile", "dataFile"),
                ("datafiletype", "dataFileType"),
            )

        values = cls._normalize_spatial_keys(values)

        datavalue = values.get("datavalue")
        has_datafile = (values.get("datafile") or values.get("dataFile")) is not None
        has_datafiletype = (values.get("datafiletype") or values.get("dataFileType")) is not None

        if datavalue is not None:
            cls._validate_datavalue_path(values, has_datafile, has_datafiletype)
        else:
            cls._validate_datafile_path(values, has_datafile, has_datafiletype)

        return values

    def is_intermediate_link(self) -> bool:
        return True

    @field_validator("targetmaskfile", mode="before")
    @classmethod
    def validate_targetmaskfile(cls, v: Any) -> Any:
        if isinstance(v, (str, Path)):
            return resolve_file_model(v, PolyFile)
        return v

    @field_validator("averagingtype", mode="before")
    @classmethod
    def validate_average_type(cls, v):
        return enum_value_parser(v, AveragingType)

    @field_validator("targetlayer", mode="before")
    @classmethod
    def validate_targetlayer(cls, v):
        """Coerce targetLayer to a TargetLayer member or a positive integer layer number.

        Accepts `bottom` and `all` (case-insensitive) or a positive integer. The
        old external-forcings `LAYER` values `-1` and `0` are represented by
        `bottom` and `all` respectively and are not accepted as integers here.
        """
        result = v
        if v is not None and not isinstance(v, TargetLayer):
            text = str(v).strip()
            if text.lower() in (TargetLayer.bottom.value, TargetLayer.all.value):
                result = TargetLayer(text.lower())
            elif text.lstrip("+").isdigit() and int(text) > 0:
                result = int(text)
            else:
                raise ValueError(
                    "targetLayer must be 'bottom', 'all', or a positive integer, "
                    f"got '{v}'."
                )
        return result

Comments

Bases: Comments

Comments for the Spatial block fields.

Inherits the shared descriptions from SpatialForcingBase.Comments; only the data* file keywords, the extrapolationSearchRadius wording, and the fields unique to the Spatial block are declared here.

Source code in hydrolib/core/dflowfm/ext/models.py
1009
1010
1011
1012
1013
1014
1015
1016
1017
1018
1019
1020
1021
1022
1023
1024
1025
1026
1027
1028
1029
1030
1031
1032
1033
1034
1035
1036
1037
1038
1039
1040
1041
1042
1043
1044
1045
1046
1047
1048
1049
1050
1051
1052
1053
1054
1055
1056
1057
1058
1059
1060
1061
1062
1063
1064
1065
1066
1067
1068
1069
1070
1071
1072
1073
1074
class Comments(SpatialForcingBase.Comments):
    """Comments for the Spatial block fields.

    Inherits the shared descriptions from `SpatialForcingBase.Comments`; only
    the `data*` file keywords, the `extrapolationSearchRadius` wording, and the
    fields unique to the Spatial block are declared here.
    """

    datafile: str | None = Field(
        "Name of file containing the data for this spatial quantity.",
        alias="dataFile",
    )
    datafiletype: str | None = Field(
        "Type of dataFile.", alias="dataFileType"
    )
    datavariablename: str | None = Field(
        "Variable name used in dataFile associated with this quantity.",
        alias="dataVariableName",
    )
    targetmaskfile: str | None = Field(
        "Name of <*.pol> file to be used as mask. Grid parts inside any polygon will receive the spatial forcing.",
        alias="targetMaskFile",
    )
    extrapolationsearchradius: str | None = Field(
        "Maximum search radius for nearest neighbour extrapolation in space.",
        alias="extrapolationSearchRadius",
    )
    averagingtype: str | None = Field(
        "Type of averaging, if interpolationMethod=averaging.",
        alias="averagingType",
    )
    averagingrelsize: str | None = Field(
        "Relative search cell size for averaging.", alias="averagingRelSize"
    )
    averagingnummin: str | None = Field(
        "Minimum number of points in averaging. Must be ≥ 1.",
        alias="averagingNumMin",
    )
    averagingpercentile: str | None = Field(
        "Percentile value for which data values to include in averaging. 0.0 means off.",
        alias="averagingPercentile",
    )
    locationtype: str | None = Field(
        "Target location of interpolation.", alias="locationType"
    )
    datavalue: str | None = Field(
        "Constant value to be set inside all model points inside the polygon, "
        "used when no dataFile/dataFileType is specified. "
        "Requires targetMaskFile=*.pol and interpolationMethod=constant.",
        alias="dataValue",
    )
    frictiontype: str | None = Field(
        "Only for quantity=frictionCoefficient. The friction type.", alias="frictionType",
    )
    tracerfallvelocity: str | None = Field(
        "Only for initialtracer<tracername>. Fall velocity of the tracer.",
        alias="tracerFallVelocity",
    )
    tracerdecaytime: str | None = Field(
        "Only for initialtracer<tracername>. Decay time of the tracer.",
        alias="tracerDecayTime",
    )
    targetlayer: str | None = Field(
        "Target layer for the data: bottom, all, or a positive layer number.",
        alias="targetLayer",
    )

validate_datavalue_or_datafile(values) classmethod

Validates the two mutually exclusive usage paths of a Spatial block.

When dataValue is provided the block describes a constant value applied inside a polygon mask. In this mode: - dataFile and dataFileType must not be specified. - targetMaskFile (a .pol file) is optional but often used. - interpolationMethod must be constant (set automatically when omitted).

When dataValue is absent, dataFile and dataFileType are both required.

Note: using dataFileType=polygon for "inside polygon" initial-condition data is deprecated. Use dataValue + targetMaskFile=*.pol + interpolationMethod=constant instead. The polygon dataFileType remains supported for quantities such as initialvertical* (e.g. initialverticalsalinityprofile) that use polygon files for a different purpose and have no new alternative yet.

The data file is resolved to its concrete file model only on the dataFile path (dataValue absent). This prevents an invalid combination (dataValue together with dataFile) from parsing a file before the mutual-exclusion check rejects it, so callers get the exclusion error rather than a file-parse error.

Source code in hydrolib/core/dflowfm/ext/models.py
1189
1190
1191
1192
1193
1194
1195
1196
1197
1198
1199
1200
1201
1202
1203
1204
1205
1206
1207
1208
1209
1210
1211
1212
1213
1214
1215
1216
1217
1218
1219
1220
1221
1222
1223
1224
1225
1226
1227
1228
1229
1230
1231
1232
1233
1234
1235
1236
1237
1238
1239
@model_validator(mode="before")
@classmethod
def validate_datavalue_or_datafile(cls, values: Dict) -> Dict:
    """Validates the two mutually exclusive usage paths of a Spatial block.

    When `dataValue` is provided the block describes a constant value applied
    inside a polygon mask.  In this mode:
    - `dataFile` and `dataFileType` must **not** be specified.
    - `targetMaskFile` (a `.pol` file) is optional but often used.
    - `interpolationMethod` must be `constant` (set automatically when omitted).

    When `dataValue` is absent, `dataFile` and `dataFileType` are both
    required.

    Note: using `dataFileType=polygon` for "inside polygon" initial-condition
    data is **deprecated**.  Use `dataValue` + `targetMaskFile=*.pol` +
    `interpolationMethod=constant` instead.  The `polygon` dataFileType
    remains supported for quantities such as `initialvertical*` (e.g.
    `initialverticalsalinityprofile`) that use polygon files for a different
    purpose and have no new alternative yet.

    The data file is resolved to its concrete file model only on the `dataFile`
    path (`dataValue` absent). This prevents an invalid combination
    (`dataValue` together with `dataFile`) from parsing a file before the
    mutual-exclusion check rejects it, so callers get the exclusion error rather
    than a file-parse error.
    """
    values = cls._process_section_values(values)

    on_datavalue_path = (
        values.get("datavalue") is not None or values.get("dataValue") is not None
    )
    if not on_datavalue_path:
        values = cls._resolve_file_models(
            values,
            ("datafile", "dataFile"),
            ("datafiletype", "dataFileType"),
        )

    values = cls._normalize_spatial_keys(values)

    datavalue = values.get("datavalue")
    has_datafile = (values.get("datafile") or values.get("dataFile")) is not None
    has_datafiletype = (values.get("datafiletype") or values.get("dataFileType")) is not None

    if datavalue is not None:
        cls._validate_datavalue_path(values, has_datafile, has_datafiletype)
    else:
        cls._validate_datafile_path(values, has_datafile, has_datafiletype)

    return values

validate_targetlayer(v) classmethod

Coerce targetLayer to a TargetLayer member or a positive integer layer number.

Accepts bottom and all (case-insensitive) or a positive integer. The old external-forcings LAYER values -1 and 0 are represented by bottom and all respectively and are not accepted as integers here.

Source code in hydrolib/core/dflowfm/ext/models.py
1256
1257
1258
1259
1260
1261
1262
1263
1264
1265
1266
1267
1268
1269
1270
1271
1272
1273
1274
1275
1276
1277
@field_validator("targetlayer", mode="before")
@classmethod
def validate_targetlayer(cls, v):
    """Coerce targetLayer to a TargetLayer member or a positive integer layer number.

    Accepts `bottom` and `all` (case-insensitive) or a positive integer. The
    old external-forcings `LAYER` values `-1` and `0` are represented by
    `bottom` and `all` respectively and are not accepted as integers here.
    """
    result = v
    if v is not None and not isinstance(v, TargetLayer):
        text = str(v).strip()
        if text.lower() in (TargetLayer.bottom.value, TargetLayer.all.value):
            result = TargetLayer(text.lower())
        elif text.lstrip("+").isdigit() and int(text) > 0:
            result = int(text)
        else:
            raise ValueError(
                "targetLayer must be 'bottom', 'all', or a positive integer, "
                f"got '{v}'."
            )
    return result

SpatialError

Bases: Exception

SpatialError.

Source code in hydrolib/core/dflowfm/ext/models.py
1382
1383
1384
1385
1386
1387
class SpatialError(Exception):
    """SpatialError."""

    def __init__(self, error_message: str):
        """Initialize with an error message."""
        super().__init__(error_message)

__init__(error_message)

Initialize with an error message.

Source code in hydrolib/core/dflowfm/ext/models.py
1385
1386
1387
def __init__(self, error_message: str):
    """Initialize with an error message."""
    super().__init__(error_message)

SpatialForcingBase

Bases: OperandInterpolationValidators, INIBasedModel, ABC

Shared behaviour for the [Meteo] and [Spatial] external-forcing blocks.

Meteo (legacy) and Spatial (its successor) share the same data-file model resolution logic and unknown-keyword handling. This abstract base holds that common behaviour so a single fix applies to both.

The operand / interpolationMethod validators are inherited from OperandInterpolationValidators (common to all four spatial-field blocks). Spatial additionally inherits LocationTypeDataFileTypeValidators for the locationType / dataFileType validators it shares with the inifield blocks; Meteo does not, as it has neither field.

Field declarations remain on the concrete subclasses: their keyword names differ (forcing* versus data*) and their serialization order must be preserved, so only behaviour (not fields) is hoisted here.

Source code in hydrolib/core/dflowfm/ext/models.py
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
class SpatialForcingBase(OperandInterpolationValidators, INIBasedModel, ABC):
    """Shared behaviour for the `[Meteo]` and `[Spatial]` external-forcing blocks.

    `Meteo` (legacy) and `Spatial` (its successor) share the same data-file model
    resolution logic and unknown-keyword handling. This abstract base holds that
    common behaviour so a single fix applies to both.

    The `operand` / `interpolationMethod` validators are inherited from
    `OperandInterpolationValidators` (common to all four spatial-field blocks). `Spatial`
    additionally inherits `LocationTypeDataFileTypeValidators` for the `locationType` /
    `dataFileType` validators it shares with the inifield blocks; `Meteo` does not,
    as it has neither field.

    Field declarations remain on the concrete subclasses: their keyword names
    differ (`forcing*` versus `data*`) and their serialization order must be
    preserved, so only behaviour (not fields) is hoisted here.
    """

    class Comments(INIBasedModel.Comments):
        """Comments shared by the `[Meteo]` and `[Spatial]` block fields.

        Only the descriptions that are identical in both blocks live here. The
        file-specific ones (`forcingFile`/`dataFile`, `extrapolationSearchRadius`)
        stay on the subclasses because their wording differs.
        """

        quantity: Optional[str] = Field(
            "Name of the quantity. See UM Section C.5.3", alias="quantity"
        )
        targetmaskinvert: Optional[str] = Field(
            "Flag indicating whether the target mask should be inverted, i.e., outside of all polygons: no or yes.",
            alias="targetMaskInvert",
        )
        interpolationmethod: Optional[str] = Field(
            "Type of (spatial) interpolation.", alias="interpolationMethod"
        )
        operand: Optional[str] = Field(
            "How this data is combined with previous data for the same quantity (if any).",
            alias="operand",
        )
        extrapolationallowed: Optional[str] = Field(
            "Optionally allow nearest neighbour extrapolation in space (0: no, 1: yes). Default off.",
            alias="extrapolationAllowed",
        )

    @classmethod
    def _get_unknown_keyword_error_manager(cls) -> Optional[UnknownKeywordErrorManager]:
        """Neither block currently raises an error on unknown keywords."""
        return None

    @staticmethod
    def _resolve_file_models(
        values: Dict[str, Any], file_keys: tuple, type_keys: tuple
    ) -> Dict[str, Any]:
        """Select the concrete file model for the data/forcing file from its type.

        Mirrors the historical per-class `choose_file_model` bodies: when both a
        file keyword and its type keyword are present and the file is still a raw
        path, the path is resolved into the model class from
        `FILETYPE_FILEMODEL_MAPPING` (types not in the mapping fall back to
        `DiskOnlyFileModel`).

        Args:
            values: Raw, unvalidated input values for the block.
            file_keys: The lowercase and camelCase names of the file keyword,
                e.g. `("datafile", "dataFile")`.
            type_keys: The lowercase and camelCase names of the file-type keyword,
                e.g. `("datafiletype", "dataFileType")`.

        Returns:
            Dict[str, Any]: The (possibly updated) values dictionary.
        """
        if any(key in values for key in type_keys) and any(
            key in values for key in file_keys
        ):
            type_key = type_keys[0] if type_keys[0] in values else type_keys[1]
            file_key = file_keys[0] if file_keys[0] in values else file_keys[1]

            file_type = values.get(type_key)
            file_type = str(file_type).lower() if file_type is not None else None

            raw_path = values.get(file_key)
            if isinstance(raw_path, (Path, str)):
                model = FILETYPE_FILEMODEL_MAPPING.get(file_type)
                if model is None:
                    values[file_key] = DiskOnlyFileModel(raw_path)
                else:
                    values[file_key] = resolve_file_model(raw_path, model)

        return values

Comments

Bases: Comments

Comments shared by the [Meteo] and [Spatial] block fields.

Only the descriptions that are identical in both blocks live here. The file-specific ones (forcingFile/dataFile, extrapolationSearchRadius) stay on the subclasses because their wording differs.

Source code in hydrolib/core/dflowfm/ext/models.py
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
class Comments(INIBasedModel.Comments):
    """Comments shared by the `[Meteo]` and `[Spatial]` block fields.

    Only the descriptions that are identical in both blocks live here. The
    file-specific ones (`forcingFile`/`dataFile`, `extrapolationSearchRadius`)
    stay on the subclasses because their wording differs.
    """

    quantity: Optional[str] = Field(
        "Name of the quantity. See UM Section C.5.3", alias="quantity"
    )
    targetmaskinvert: Optional[str] = Field(
        "Flag indicating whether the target mask should be inverted, i.e., outside of all polygons: no or yes.",
        alias="targetMaskInvert",
    )
    interpolationmethod: Optional[str] = Field(
        "Type of (spatial) interpolation.", alias="interpolationMethod"
    )
    operand: Optional[str] = Field(
        "How this data is combined with previous data for the same quantity (if any).",
        alias="operand",
    )
    extrapolationallowed: Optional[str] = Field(
        "Optionally allow nearest neighbour extrapolation in space (0: no, 1: yes). Default off.",
        alias="extrapolationAllowed",
    )

TargetLayer

Bases: StrEnum

Valid non-numeric values for the targetLayer attribute of a [Spatial] block.

Corresponds to the LAYER value in the old external forcings file: bottom (old -1) and all (old 0). A positive integer layer number is also accepted; see Spatial.targetlayer.

Source code in hydrolib/core/dflowfm/ext/models.py
74
75
76
77
78
79
80
81
82
83
class TargetLayer(StrEnum):
    """Valid non-numeric values for the `targetLayer` attribute of a `[Spatial]` block.

    Corresponds to the `LAYER` value in the old external forcings file: `bottom`
    (old `-1`) and `all` (old `0`). A positive integer layer number is also
    accepted; see `Spatial.targetlayer`.
    """

    bottom = "bottom"
    all = "all"