diff --git a/examples/hybrid/sphere/hadley_circulation.jl b/examples/hybrid/sphere/hadley_circulation.jl
index 441bcb326c..85b2c3571b 100644
--- a/examples/hybrid/sphere/hadley_circulation.jl
+++ b/examples/hybrid/sphere/hadley_circulation.jl
@@ -281,10 +281,10 @@ apply_remap(datafile_rll, datafile_cc, weightfile, ["q1"])
 
 # load remapped data and create statistics for plots
 nt = NCDataset(datafile_rll, "r") do nc
-    lat = nc["lat"][:]
-    z = nc["z"][:]
-    time = nc["time"][:]
-    q1 = nc["q1"][:]
+    lat = Array(nc["lat"])
+    z = Array(nc["z"])
+    time = Array(nc["time"])
+    q1 = Array(nc["q1"])
     (; lat, z, time, q1)
 end
 (; lat, z, time, q1) = nt