OutOfBoundsError when particle reaches the North Pole (~ 89.5 deg N)

[eavellashaw]

Hi!

I try to advect particles in the North-Atlantic Ocean and follow them to the Arctic, but as soon as a particle tries to cross the last latitude corresponding to the pole (in my case 89.50000080000001), I get an OutOfBoundsError as follows:

Traceback (most recent call last): File "/aos/home/eashaw/Desktop/OneParticle_Adv.py", line 162, in <module> pset.execute(kernels, File "/aos/home/eashaw/miniconda3/envs/py3_parcels/lib/python3.11/site-packages/parcels/particleset/baseparticleset.py", line 450, in execute self.kernel.execute(self, endtime=next_time, dt=dt, recovery=recovery, output_file=output_file, File "/aos/home/eashaw/miniconda3/envs/py3_parcels/lib/python3.11/site-packages/parcels/kernel/kernelsoa.py", line 227, in execute recovery_kernel(p, self.fieldset, p.time) File "/aos/home/eashaw/miniconda3/envs/py3_parcels/lib/python3.11/site-packages/parcels/tools/statuscodes.py", line 202, in recovery_kernel_out_of_bounds raise OutOfBoundsError(particle, fieldset) parcels.tools.statuscodes.OutOfBoundsError: 0 Particle P[0](lon=19.108109, lat=89.499947, depth=214.370804, temperature=0.327767, salinity=34.611294, uvel=-0.000001, vvel=0.000000, wvel=-0.000001, time=136570200.000000) Time: 0185-05-16 04:10:00, timestep dt: 600.000000 Out-of-bounds sampling by particle at (19.108109, 89.499947, 214.370804)

My velocity fields come from the GFDL CM2-1deg model (B grid), but the same error rises with velocity fields obtained from GLORYS with an ORCA12 grid.

Is there a way to avoid this issue? The only solution I found for the moment is to implement a recovery kernel killing the particle as soon as this error rises for the model not to crash.

Thanks in advance,

Esteban Avella Shaw