diff --git a/bfps/DNS.py b/bfps/DNS.py
index e15e3ee847f420c12c08644320533d30259616c0..8cbe8d9cc3c165054212b363a962d10471013a44 100644
--- a/bfps/DNS.py
+++ b/bfps/DNS.py
@@ -847,6 +847,7 @@ class DNS(_code):
                          min(dst_shape[1], src_file[src_dset_name].shape[1]),
                          min(dst_shape[2], src_file[src_dset_name].shape[2]),
                          3)
+            src_shape = src_file[src_dset_name].shape
             dst_file.create_dataset(
                     dst_dset_name,
                     shape = dst_shape,