diff --git a/PyMPDATA/scalar_field.py b/PyMPDATA/scalar_field.py
index 972ae11e..d3a64df4 100644
--- a/PyMPDATA/scalar_field.py
+++ b/PyMPDATA/scalar_field.py
@@ -55,7 +55,7 @@ def make_null(n_dims, traversals):
         see `Field._make_null` other properties of the returned field"""
         return Field._make_null(
             ScalarField(
-                np.empty([INVALID_NULL_VALUE] * n_dims),
+                np.empty(shape=[0] * n_dims),
                 halo=0,
                 boundary_conditions=tuple([Constant(np.nan)] * n_dims),
             ),