fix: Fix velocity in advection2d.py

seismic-anisotropy · Mar 21, 2024 · 22a24a3 · 22a24a3
1 parent 1176ae2
commit 22a24a3
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/examples/fluidity/advection2d/advection2d.py b/examples/fluidity/advection2d/advection2d.py
@@ -51,7 +51,7 @@ def _interp_vel_grad(int_time):
 @nb.njit(fastmath=True)
 def velocity(coords, time):
     """Return prescribed velocity for 2D simple shear flow."""
-    return np.array([coords[0] * STRAIN_RATE, 0.0])
+    return np.array([coords[1] * STRAIN_RATE, 0.0])
 
 
 @nb.njit(fastmath=True)