diff --git a/cf/read_write/um/umread.py b/cf/read_write/um/umread.py
index 04ef52d8bf..a54a00a439 100644
--- a/cf/read_write/um/umread.py
+++ b/cf/read_write/um/umread.py
@@ -1970,7 +1970,7 @@ def create_data(self):
                 # ----------------------------------------------------
                 # 1-d partition matrix
                 # ----------------------------------------------------
-                z_axis = _axis[self.z_axis]
+                z_axis = _axis.get(self.z_axis)
                 if nz > 1:
                     pmaxes = [z_axis]
                     data_shape = (nz, LBROW, LBNPT)
diff --git a/docs/source/installation.rst b/docs/source/installation.rst
index 8c7311e214..148ff077f9 100644
--- a/docs/source/installation.rst
+++ b/docs/source/installation.rst
@@ -255,7 +255,7 @@ environments for which these features are not required.
 
 .. rubric:: Convolution filters, derivatives and relative vorticity
 
-* `scipy <https://pypi.org/project/scipy>`_, version 1.1.0 or newer.
+* `scipy <https://pypi.org/project/scipy>`_, version 1.10.0 or newer.
 
 .. rubric:: Subspacing based on N-dimensional construct cells (N > 1)
             containing a given value
diff --git a/docs/source/tutorial.py b/docs/source/tutorial.py
index 786df6de3c..bf1ca00fbd 100644
--- a/docs/source/tutorial.py
+++ b/docs/source/tutorial.py
@@ -758,9 +758,8 @@
 import netCDF4
 nc = netCDF4.Dataset('file.nc', 'r')
 v = nc.variables['ta']
-netcdf_array = cf.NetCDFArray(filename='file.nc', ncvar='ta',
-                               dtype=v.dtype, ndim=v.ndim,
-     		  	       shape=v.shape, size=v.size)
+netcdf_array = cf.NetCDFArray(filename='file.nc', address='ta',
+                               dtype=v.dtype, shape=v.shape)
 data_disk = cf.Data(netcdf_array)
 numpy_array = v[...]
 data_memory = cf.Data(numpy_array)
diff --git a/docs/source/tutorial.rst b/docs/source/tutorial.rst
index 2ca8f7f1cb..4c6bd5008d 100644
--- a/docs/source/tutorial.rst
+++ b/docs/source/tutorial.rst
@@ -4539,9 +4539,8 @@ used to initialise a `cf.Data` instance.
    >>> import netCDF4
    >>> nc = netCDF4.Dataset('file.nc', 'r')
    >>> v = nc.variables['ta']
-   >>> netcdf_array = cf.NetCDFArray(filename='file.nc', ncvar='ta',
-   ...	                               dtype=v.dtype, ndim=v.ndim,
-   ...	     		  	       shape=v.shape, size=v.size)
+   >>> netcdf_array = cf.NetCDFArray(filename='file.nc', address='ta',
+   ...	                               dtype=v.dtype, shape=v.shape)
    >>> data_disk = cf.Data(netcdf_array)