diff --git a/docs/src/further_topics/dask_best_practices/dask_parallel_loop.rst b/docs/src/further_topics/dask_best_practices/dask_parallel_loop.rst
index 836503314c..2c19196318 100644
--- a/docs/src/further_topics/dask_best_practices/dask_parallel_loop.rst
+++ b/docs/src/further_topics/dask_best_practices/dask_parallel_loop.rst
@@ -42,7 +42,7 @@ We start out by loading cubes of pressure, temperature, dewpoint temperature and
 We set up the NumPy arrays we will be filling with the output data::
 
     output_arrays = [np.zeros(pressure.shape[0]) for _ in range(6)]
-    cape, cin, lcl, lfc, el, tpw = output_data
+    cape, cin, lcl, lfc, el, tpw = output_arrays
 
 Now we loop over the columns in the data to calculate the soundings::