[field] Fix resampling via with_values()

phi/field/_field.py

@@ -319,8 +319,11 @@ def closest_values(self, points: Tensor):
             channels = [component.closest_values(points) for component in self.vector.unstack()]
         return math.stack(channels, points.shape['~vector'])
 
-    def with_values(self, values):
+    def with_values(self, values, **sampling_kwargs):
         """ Returns a copy of this field with `values` replaced. """
+        if not isinstance(values, Tensor):
+            from ._resample import sample
+            values = sample(values, self._geometry, self.sampled_at, self._boundary, **sampling_kwargs)
         return Field(self._geometry, values, self._boundary)
 
     def with_boundary(self, boundary):