passage de la ligne de changement de date avec donnees manquantes

tsg_util/updateTsgStruct.m

@@ -31,11 +31,13 @@ NO_CONTROL = s.code;
 
 % Suppress longitude discontinuity at -180 or 180 crossing 
 % --------------------------------------------------------
-dlon=circshift(tsg.LONX,-1)-tsg.LONX;
+deflon=find(isfinite(tsg.LONX));
+dlon=circshift(tsg.LONX(deflon),-1)-tsg.LONX(deflon);
 dlon=dlon(1:end-1);
-tsg.indcross=find(abs(dlon)>350);
+tsg.indcross=find(abs(dlon)>180);
 if ~isempty(tsg.indcross)
     tsg.loncross=sign(dlon(tsg.indcross))*360;
+    tsg.indcross=deflon(tsg.indcross);
     for i=1:length(tsg.indcross)
         tsg.LONX(tsg.indcross(i)+1:end)=tsg.LONX(tsg.indcross(i)+1:end)-tsg.loncross(i);
     end