(xi, yi)
| 419 | speed = np.ma.sqrt(u_ax ** 2 + v_ax ** 2) |
| 420 | |
| 421 | def forward_time(xi, yi): |
| 422 | ds_dt = interpgrid(speed, xi, yi) |
| 423 | if ds_dt == 0: |
| 424 | raise TerminateTrajectory() |
| 425 | dt_ds = 1. / ds_dt |
| 426 | ui = interpgrid(u, xi, yi) |
| 427 | vi = interpgrid(v, xi, yi) |
| 428 | return ui * dt_ds, vi * dt_ds |
| 429 | |
| 430 | def backward_time(xi, yi): |
| 431 | dxi, dyi = forward_time(xi, yi) |
no test coverage detected