(xi, yi)
| 500 | speed = np.ma.sqrt(u_ax ** 2 + v_ax ** 2) |
| 501 | |
| 502 | def forward_time(xi, yi): |
| 503 | if not dmap.grid.within_grid(xi, yi): |
| 504 | raise OutOfBounds |
| 505 | ds_dt = interpgrid(speed, xi, yi) |
| 506 | if ds_dt == 0: |
| 507 | raise TerminateTrajectory() |
| 508 | dt_ds = 1. / ds_dt |
| 509 | ui = interpgrid(u, xi, yi) |
| 510 | vi = interpgrid(v, xi, yi) |
| 511 | return ui * dt_ds, vi * dt_ds |
| 512 | |
| 513 | def backward_time(xi, yi): |
| 514 | dxi, dyi = forward_time(xi, yi) |
no test coverage detected
searching dependent graphs…