source: issm/oecreview/Archive/26740-27031/ISSM-26962-26963.diff

Index: ../trunk-jpl/test/MITgcm/tools/lookat_4004.m
===================================================================
--- ../trunk-jpl/test/MITgcm/tools/lookat_4004.m        (revision 26962)
+++ ../trunk-jpl/test/MITgcm/tools/lookat_4004.m        (revision 26963)
@@ -7,8 +7,8 @@
     subplot(5,1,ts+1), mypcolor(r2-r1); colorbar
 end
 
-figure(2), orient tall, wysiwyg
 for ts=1:5
+    figure(ts+1), clf, orient tall, wysiwyg
     e=readbin([pn 'Eta.' myint2str(ts,10) '.data'],[3 200]);
     p=readbin([pn 'PHL.' myint2str(ts,10) '.data'],[3 200]);
     f=readbin([pn 'SHICE_fwFlux.' myint2str(ts,10) '.data'],[3 200]);
@@ -19,5 +19,21 @@
     subplot(513), mypcolor(f); title('SHICE fwFlux'), colorbar
     subplot(514), mypcolor(h); title('SHICE heatFlux'), colorbar
     subplot(515), mypcolor(1:200,-1:-1:-90,squeeze(ph(2,:,:))'); title('PH'), colorbar
-    pause
 end
+
+for ts=0:5
+    figure(ts+7), clf, orient tall, wysiwyg
+    s =readbin([pn 'S.' myint2str(ts,10) '.data'],[3 200 90]);
+    t =readbin([pn 'T.' myint2str(ts,10) '.data'],[3 200 90]);
+    u =readbin([pn 'U.' myint2str(ts,10) '.data'],[3 200 90]);
+    v =readbin([pn 'V.' myint2str(ts,10) '.data'],[3 200 90]);
+    w =readbin([pn 'W.' myint2str(ts,10) '.data'],[3 200 90]);
+    in=find(~s); s(in)=nan; t(in)=nan;
+    clf, subplot(711), mypcolor(r1); title(['draft @ ts=' int2str(ts)]), colorbar
+    subplot(712), mypcolor(r2-r1); title('draft change'), colorbar
+    subplot(713), pcolorcen(1:200,-1:-1:-90,squeeze(s(2,:,:))'); title('S'), colorbar
+    subplot(714), pcolorcen(1:200,-1:-1:-90,squeeze(t(2,:,:))'); title('T'), colorbar
+    subplot(715), pcolorcen(1:200,-1:-1:-90,squeeze(u(2,:,:))'); title('U'), colorbar
+    subplot(716), pcolorcen(1:200,-1:-1:-90,squeeze(v(2,:,:))'); title('V'), colorbar
+    subplot(717), pcolorcen(1:200,-1:-1:-90,squeeze(w(2,:,:))'); title('W'), colorbar
+end