--- 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