From e0c0b7de18875bf2a4a21839cd0dd3966476385a Mon Sep 17 00:00:00 2001
From: Philipp Holl <philipp.holl@tum.de>
Date: Mon, 25 Nov 2024 22:34:59 +0100
Subject: [PATCH] [physics] Fix MacCormack advection

---
 phi/physics/advect.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/phi/physics/advect.py b/phi/physics/advect.py
index 1101939d8..a535745bc 100644
--- a/phi/physics/advect.py
+++ b/phi/physics/advect.py
@@ -200,7 +200,7 @@ def mac_cormack(field: Field,
     Returns:
         Advected field of type `type(field)`
     """
-    v0 = resample(velocity, to=field).values
+    v0 = sample(velocity, field.geometry, at=field.sampled_at, boundary=field.boundary)
     points_bwd = integrator(field, velocity, -dt, v0=v0)
     points_fwd = integrator(field, velocity, dt, v0=v0)
     # --- forward+backward semi-Lagrangian advection ---