diff --git a/meshmode/array_context.py b/meshmode/array_context.py
index 6118e9d5530c1ee46aff3ff75ed4d0413698420b..81461b71aa9008e52fcb379532cfc51ea5de4b2c 100644
--- a/meshmode/array_context.py
+++ b/meshmode/array_context.py
@@ -153,6 +153,14 @@ class ArrayContext:
 
     @memoize_method
     def _get_scalar_func_loopy_program(self, name, nargs, naxes):
+        if name == "arctan2":
+            name = "atan2"
+        elif name == "atan2":
+            from warnings import warn
+            warn("'atan2' in ArrayContext.np is deprecated. Use 'arctan2', "
+                    "as in numpy2. This will be disallowed in 2021.",
+                    DeprecationWarning, stacklevel=3)
+
         from pymbolic import var
 
         var_names = ["i%d" % i for i in range(naxes)]