diff --git a/tsg_util/plot_Climatology.m b/tsg_util/plot_Climatology.m
index 54c804e952ad352ebca7b409092385cad2dac361..8896ec0310a5791f72f751a7c886638149e47114 100644
--- a/tsg_util/plot_Climatology.m
+++ b/tsg_util/plot_Climatology.m
@@ -91,8 +91,16 @@ for i = 1:2
   for ii=1:length(ind)
     ilat     = find(tsg.levitus.data.WOA_LATX == lat2(ii));
     ilon     = find(tsg.levitus.data.WOA_LONX == lon2(ii));
-    mean(ii) = tsg.levitus.data.(['WOA_MEAN_' para])(time_dim,depth,ilat,ilon);
-    std(ii)  = tsg.levitus.data.(['WOA_STD_'  para])(time_dim,depth,ilat,ilon);
+    if ~isempty(tsg.levitus.data.(['WOA_MEAN_' para])(time_dim,depth,ilat,ilon))
+        mean(ii) = tsg.levitus.data.(['WOA_MEAN_' para])(time_dim,depth,ilat,ilon);
+    else
+        mean(ii) = nan;
+    end
+    if ~isempty(tsg.levitus.data.(['WOA_STD_'  para])(time_dim,depth,ilat,ilon))
+        std(ii)  = tsg.levitus.data.(['WOA_STD_'  para])(time_dim,depth,ilat,ilon);
+    else
+        std(ii)  = nan;
+    end
   end
 
   % Select the axes