Diff of /MITgcm_contrib/enderton/Diagnostics/DiagLoadGradsData.m

Parent Directory | Revision Log | Revision Graph | Patch

--- MITgcm_contrib/enderton/Diagnostics/DiagLoadGradsData.m	2005/01/31 19:58:15	1.2
+++ MITgcm_contrib/enderton/Diagnostics/DiagLoadGradsData.m	2005/04/29 17:34:39	1.3
@@ -5,7 +5,7 @@
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %                          Parse table file                               %
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
+ShiftData = 0;
 format = 'NONSEQUENTIAL';
 
 tablfile = [dad,'/',Grads,'.tabl'];
@@ -27,7 +27,7 @@
         inc = str2num(tokens{5});
         xax = [ini:inc:ini+(num-1)*inc];
         if min(xax) >= 0 && max(xax) > 180
-            xax = xax - 180;
+            ShiftData = 1;
         end
         nx = length(xax);
     end
@@ -147,4 +147,10 @@
 end
 data = reshape(data,[nx,ny,nz,nv,nt]);
 data = squeeze(data(:,:,:,ivar,:));
-data( abs((data-undef)/undef) < tol ) = NaN;
\ No newline at end of file
+data( abs((data-undef)/undef) < tol ) = NaN;
+if ShiftData
+    indexWestHemi = xax>=180;
+    indexEastHemi = xax<180;
+    data = cat(1,data(indexWestHemi,:,:),data(indexEastHemi,:,:));
+    xax = cat(2,xax(indexWestHemi)-360,xax(indexEastHemi));
+end
\ No newline at end of file