Changeset 26966

Timestamp:

04/23/22 13:03:52 (3 years ago)

Author:

dmenemen

Message:

M test/MITgcm/tools/lookat_4004.m

File:

• issm/trunk-jpl/test/MITgcm/tools/lookat_4004.m (modified) (2 diffs)

issm/trunk-jpl/test/MITgcm/tools/lookat_4004.m

@@ -24 +24 @@
 for ts=0:5
     figure(ts+7), clf, orient tall, wysiwyg
+    r2=readbin([pn 'R_shelfIce2_' myint2str(ts,10) '.data'],[3 200]);
+    if ts<5
+        r1=readbin([pn 'R_shelfIce1_' myint2str(ts+1,10) '.data'],[3 200]);
+    end
     s =readbin([pn 'S.' myint2str(ts,10) '.data'],[3 200 90]);
     t =readbin([pn 'T.' myint2str(ts,10) '.data'],[3 200 90]);
@@ -30 +34 @@
     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
+    clf,
+    subplot(711), mypcolor(r2); 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