lab_sea/matlab/lookat_exp5.m

% need to be in verification/lab_sea/matlab directory
% and to specify location of exp1 output
cd ../../../verification/lab_sea/matlab
exp1_path='../../../exe/';

% load monthly-mean SMMR-SSM/I data
load SSMI

% B-grid latitude for uice and vice
Blat=lat+1; Blon=lon+1;

% load model output
area=readbin([exp1_path 'AREAtave.0000010950.data'],[20 16 1],1);
heff=readbin([exp1_path 'HEFFtave.0000010950.data'],[20 16 1],1);
uice=readbin([exp1_path 'UICEtave.0000010950.data'],[20 16 1],1);
vice=readbin([exp1_path 'VICEtave.0000010950.data'],[20 16 1],1);

% compute ice speed and direction
udir=0*uice; vdir=0*vice;
icespeed=100*sqrt(uice.^2+vice.^2);
in=find(icespeed);
udir(in)=uice(in)./icespeed(in);
vdir(in)=vice(in)./icespeed(in);

% plot comparison figures
clf
subplot(221)
mypcolor(lon,lat,SSMI(:,:,4)'/100); caxis([0 1]), colorbar
title('Observed March sea-ice concentration')
subplot(222), mypcolor(lon,lat,area'); caxis([0 1]), colorbar
title('Modeled March sea-ice concentration')
subplot(223), mypcolor(lon,lat,heff'); caxis([0 1]), colorbar
title('Effective thickness in m')
subplot(224)
mypcolor(Blon,Blat,icespeed'); caxis([0 10]), colorbar
hold on,myquiver(Blon,Blat,udir',vdir','k');
title('Sea-ice velocity (cm/s)')
1	% need to be in verification/lab_sea/matlab directory
2	% and to specify location of exp1 output
3	cd ../../../verification/lab_sea/matlab
4	exp1_path='../../../exe/';
5
6	% load monthly-mean SMMR-SSM/I data
7	load SSMI
8
9	% B-grid latitude for uice and vice
10	Blat=lat+1; Blon=lon+1;
11
12	% load model output
13	area=readbin([exp1_path 'AREAtave.0000010950.data'],[20 16 1],1);
14	heff=readbin([exp1_path 'HEFFtave.0000010950.data'],[20 16 1],1);
15	uice=readbin([exp1_path 'UICEtave.0000010950.data'],[20 16 1],1);
16	vice=readbin([exp1_path 'VICEtave.0000010950.data'],[20 16 1],1);
17
18	% compute ice speed and direction
19	udir=0uice; vdir=0vice;
20	icespeed=100*sqrt(uice.^2+vice.^2);
21	in=find(icespeed);
22	udir(in)=uice(in)./icespeed(in);
23	vdir(in)=vice(in)./icespeed(in);
24
25	% plot comparison figures
26	clf
27	subplot(221)
28	mypcolor(lon,lat,SSMI(:,:,4)'/100); caxis([0 1]), colorbar
29	title('Observed March sea-ice concentration')
30	subplot(222), mypcolor(lon,lat,area'); caxis([0 1]), colorbar
31	title('Modeled March sea-ice concentration')
32	subplot(223), mypcolor(lon,lat,heff'); caxis([0 1]), colorbar
33	title('Effective thickness in m')
34	subplot(224)
35	mypcolor(Blon,Blat,icespeed'); caxis([0 10]), colorbar
36	hold on,myquiver(Blon,Blat,udir',vdir','k');
37	title('Sea-ice velocity (cm/s)')