Looping over a subset of times in the `FieldSet`

[gmacgilchrist]

I wish to evolve particles in a velocity field that loops multiple times over a specific year (for example, a year with certain anomalous conditions). I thought I would be able to do this by setting time_periodic to timedelta(365) and setting the initialization time to be the start of the year that I wish to loop. However, this fails with:

ValueError: Time series provided is longer than the time_periodic parameter

The implication is that I can only loop when I wish to do so across the time boundaries of the FieldSet.

Is there a specific reason for this constraint? Would it be possible to remove it, and thus have the potential to loop over a subset of the FieldSet?

I am running with parcels version 2.3.0.

[erikvansebille]

Thanks for bringing this up, @gmacgilchrist. I agree that this is not very clear in the documentation.

The design philosophy of the time_periodic functionality is that the looping will be over the full time-range of the FieldSet. If a user wants to run over a subset in time, they will have to make sure that the FieldSet is shortened (e.g., by providing fewer input netcdf files).

The reason why Parcels still supports a timedelta object in time_periodic is that Parcels doesn't know what the time interval is between the last snapshot (t_n) and the first snapshot (t₀). Parcels could assume that it's the same as between e.g. t₀ and t₁, but we want to keep the flexibility to support FieldSets with variable time intervals.

So, this is not a bug but a design choice. I will transfer this Issue to a Discussion