Skip to content

Commit f2519fe

Browse files
Fixing RK45 kernel
1 parent 2c3f924 commit f2519fe

File tree

1 file changed

+2
-2
lines changed

1 file changed

+2
-2
lines changed

src/parcels/kernels/advection.py

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -29,7 +29,7 @@ def _constrain_dt_to_within_time_interval(time_interval, time, dt):
2929
if time_interval:
3030
dt = np.where(time + dt <= time_interval.right, dt, time_interval.right - time)
3131
dt = np.where(time + dt >= time_interval.left, dt, time - time_interval.left)
32-
return dt
32+
return dt.astype("timedelta64[s]")
3333

3434

3535
def AdvectionRK2(particles, fieldset): # pragma: no cover
@@ -156,7 +156,7 @@ def AdvectionRK45(particles, fieldset): # pragma: no cover
156156
"""
157157
dt = _constrain_dt_to_within_time_interval(fieldset.time_interval, particles.time, particles.dt)
158158
dt_flt = dt / np.timedelta64(1, "s") # TODO: improve API for converting dt to seconds
159-
sign_dt = np.sign(dt)
159+
sign_dt = np.sign(dt_flt)
160160

161161
c = [1.0 / 4.0, 3.0 / 8.0, 12.0 / 13.0, 1.0, 1.0 / 2.0]
162162
A = [

0 commit comments

Comments
 (0)