Fix up

Author: pbrubeck

firedrake/interpolation.py

@@ -462,8 +462,6 @@ def __init__(
             # Probably just need to pass freeze_expr to the various
             # interpolators for this to work.
             raise NotImplementedError("freeze_expr not implemented")
-        if access != op2.WRITE:
-            raise NotImplementedError("access other than op2.WRITE not implemented")
         if bcs:
             raise NotImplementedError("bcs not implemented")
         if V.ufl_element().mapping() != "identity":
@@ -476,6 +474,9 @@ def __init__(
             )
         super().__init__(expr, V, subset, freeze_expr, access, bcs, allow_missing_dofs, matfree)
 
+        if self.access != op2.WRITE:
+            raise NotImplementedError("access other than op2.WRITE not implemented")
+
         expr = self.expr_renumbered
         self.arguments = extract_arguments(expr)
         self.nargs = len(self.arguments)