Remove FieldSet.add_field()

src/parcels/_core/fieldset.py

@@ -139,29 +139,6 @@ def time_interval(self):
 
         return overlap
 
-    def add_field(self, field: Field, name: str | None = None):
-        """Add a :class:`parcels.field.Field` object to the FieldSet.
-
-        Parameters
-        ----------
-        field : parcels.field.Field
-            Field object to be added
-        name : str
-            Name of the :class:`parcels.field.Field` object to be added. Defaults
-            to name in Field object.
-        """
-        if not isinstance(field, (Field, VectorField)):
-            raise ValueError(f"Expected `field` to be a Field or VectorField object. Got {type(field)}")
-        assert_compatible_calendars((*self.fields.values(), field))
-
-        name = field.name if name is None else name
-
-        if name in self.fields:
-            raise ValueError(f"FieldSet already has a Field with name '{name}'")
-
-        self.fields[name] = field
-        _warn_if_fields_use_different_meshes(self.fields.values())
-
     def to_windowed_arrays(self, *, max_levels: int | None = None):
         """Wrap dask-backed field data in rolling time-window caches.
 

tests/test_compat_v3.py

@@ -49,7 +49,7 @@ def assert_valid_v3_particlefile_structure(ds: xr.Dataset):
     assert ds["lat"].attrs["axis"] == "Y"  # attrs are copied accross correctly
 
 
-@settings(deadline=timedelta(seconds=0.3))
+@settings(deadline=timedelta(seconds=1))
 @example(buf=example_particlefile())
 @given(buf=pst.particlefile_output())
 def test_particlefile_to_v3_zarr(buf):

tests/test_fieldset.py

@@ -303,11 +303,6 @@ def test_fieldset_grid_deduplication():
     ...
 
 
-def test_fieldset_add_field_after_pset():
-    # ? Should it be allowed to add fields (normal or vector) after a ParticleSet has been initialized?
-    ...
-
-
 def test_fieldset_from_icon():
     ds = convert.icon_to_ugrid(datasets_unstructured["icon_square_delaunay_uniform_z_coordinate"])
     fieldset = FieldSet.from_ugrid_conventions(ds)