@@ -2005,7 +2005,7 @@ def advect_rk4(x_g, y_g, u_g, v_g, m_g, x, y, m, nb_step):
2005
2005
x_ref , y_ref = x_g [0 ], y_g [0 ]
2006
2006
x_step , y_step = x_g [1 ] - x_ref , y_g [1 ] - y_ref
2007
2007
is_circular = abs (x_g [- 1 ] % 360 - (x_g [0 ] - x_step ) % 360 ) < 1e-5
2008
- nb_x_ = x_g .size - 2
2008
+ nb_x_ = x_g .size
2009
2009
nb_x = nb_x_ if is_circular else 0
2010
2010
# cache
2011
2011
i_cache , j_cache = - 1000000 , - 1000000
@@ -2108,7 +2108,7 @@ def advect(x_g, y_g, u_g, v_g, m_g, x, y, m, nb_step):
2108
2108
x_ref , y_ref = x_g [0 ], y_g [0 ]
2109
2109
x_step , y_step = x_g [1 ] - x_ref , y_g [1 ] - y_ref
2110
2110
is_circular = abs (x_g [- 1 ] % 360 - (x_g [0 ] - x_step ) % 360 ) < 1e-5
2111
- nb_x_ = x_g .size - 2
2111
+ nb_x_ = x_g .size
2112
2112
nb_x = nb_x_ if is_circular else 0
2113
2113
# Indices which should be never exist
2114
2114
i0_old , j0_old = - 100000 , - 100000
@@ -2472,7 +2472,7 @@ def advect_t(x_g, y_g, u_g0, v_g0, m_g0, u_g1, v_g1, m_g1, x, y, m, weigths, hal
2472
2472
x_ref , y_ref = x_g [0 ], y_g [0 ]
2473
2473
x_step , y_step = x_g [1 ] - x_ref , y_g [1 ] - y_ref
2474
2474
is_circular = abs (x_g [- 1 ] % 360 - (x_g [0 ] - x_step ) % 360 ) < 1e-5
2475
- nb_x_ = x_g .size - 2
2475
+ nb_x_ = x_g .size
2476
2476
nb_x = nb_x_ if is_circular else 0
2477
2477
# Indices which should be never exist
2478
2478
i0_old , j0_old = - 100000 , - 100000
@@ -2598,7 +2598,7 @@ def advect_t_rk4(
2598
2598
x_ref , y_ref = x_g [0 ], y_g [0 ]
2599
2599
x_step , y_step = x_g [1 ] - x_ref , y_g [1 ] - y_ref
2600
2600
is_circular = abs (x_g [- 1 ] % 360 - (x_g [0 ] - x_step ) % 360 ) < 1e-5
2601
- nb_x_ = x_g .size - 2
2601
+ nb_x_ = x_g .size
2602
2602
nb_x = nb_x_ if is_circular else 0
2603
2603
# cache
2604
2604
i_cache , j_cache = - 1000000 , - 1000000
0 commit comments