diff --git a/Gallery/Contours/NCL_conLev_1.py b/Gallery/Contours/NCL_conLev_1.py
index 8dc0cc810..72ea6912c 100644
--- a/Gallery/Contours/NCL_conLev_1.py
+++ b/Gallery/Contours/NCL_conLev_1.py
@@ -33,8 +33,10 @@
                      decode_times=False)
 # Extract slice of the data
 temp = ds.TS.isel(time=43).drop_vars(names=['time'])
-# Convert from Celsius to Kelvin
+
+# Convert from Kelvin to Celsius and update units
 temp.data = temp.data - 273.15
+temp.attrs['units'] = 'C'
 
 # Fix the artifact of not-shown-data around 0 and 360-degree longitudes
 temp = gv.xr_add_cyclic_longitudes(temp, "lon")