From a653fc8b3169b6058de4c13d2168c4943a70dbc1 Mon Sep 17 00:00:00 2001
From: Gael Alory <gael.alory@legos.obs-mip.fr>
Date: Wed, 2 Sep 2009 10:35:27 +0000
Subject: [PATCH] correction d'un bug occasionnel a la lecture de la climato

---
 tsg_util/plot_Climatology.m | 12 ++++++++++--
 1 file changed, 10 insertions(+), 2 deletions(-)

diff --git a/tsg_util/plot_Climatology.m b/tsg_util/plot_Climatology.m
index 54c804e..8896ec0 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
-- 
GitLab