35 |
|
|
36 |
clf, subplot(321) |
clf, subplot(321) |
37 |
mypcolor(lon,lat,area1'); caxis([0 1]), colorbar |
mypcolor(lon,lat,area1'); caxis([0 1]), colorbar |
38 |
title('Sea-ice concentration') |
title('Sea-ice concentration, this run') |
39 |
set(gca,'xticklabel',[]) |
set(gca,'xticklabel',[]) |
40 |
|
|
41 |
subplot(322) |
subplot(322) |
42 |
mypcolor(lon,lat,area1'-area'); colorbar |
mypcolor(lon,lat,area'); caxis([0 1]), colorbar |
43 |
title('Difference with checkpoint51f') |
title('Sea-ice concentration, checkpoint52l') |
44 |
set(gca,'xticklabel',[]) |
set(gca,'xticklabel',[]) |
45 |
|
|
46 |
subplot(323) |
subplot(323) |
47 |
mypcolor(lon,lat,heff1'); caxis([0 .3]), colorbar |
mypcolor(lon,lat,heff1'); caxis([0 .6]), colorbar |
48 |
title('Effective sea-ice thickness (m)') |
title('Effective sea-ice thickness (m), this run') |
49 |
set(gca,'xticklabel',[]) |
set(gca,'xticklabel',[]) |
50 |
|
|
51 |
subplot(324) |
subplot(324) |
52 |
mypcolor(lon,lat,heff1'-heff'); colorbar |
mypcolor(lon,lat,heff'); caxis([0 .6]), colorbar |
53 |
title('Difference with checkpoint51f') |
title('Effective sea-ice thickness (m), checkpoint52l') |
54 |
set(gca,'xticklabel',[]) |
set(gca,'xticklabel',[]) |
55 |
|
|
56 |
subplot(325) |
subplot(325) |
57 |
mypcolor(Blon,Blat,icespeed1'); caxis([0 10]), colorbar |
mypcolor(Blon,Blat,icespeed1'); caxis([0 12]), colorbar |
58 |
hold on, myquiver(Blon,Blat,udir1',vdir1','k'); |
hold on, myquiver(Blon,Blat,udir1',vdir1','k'); |
59 |
title('Sea-ice velocity (cm/s)') |
title('Sea-ice velocity (cm/s), this run') |
60 |
|
|
61 |
subplot(326) |
subplot(326) |
62 |
mypcolor(Blon,Blat,icespeed1'-icespeed'); colorbar |
mypcolor(Blon,Blat,icespeed'); caxis([0 12]), colorbar |
63 |
if mmax(abs(udir1-udir)) | mmax(abs(vdir1-vdir)) |
hold on, myquiver(Blon,Blat,udir',vdir','k'); |
64 |
hold on, myquiver(Blon,Blat,udir1'-udir',vdir1'-vdir','k'); |
title('Sea-ice velocity (cm/s), checkpoint52l') |
|
end |
|
|
title('Difference with checkpoint51f') |
|