lab_sea/matlab/lookat_exp6.m

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

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

% load model output
area1=readbin([exp1_path 'AREAtave.0000000010.data'],[20 16 1],1);
heff1=readbin([exp1_path 'HEFFtave.0000000010.data'],[20 16 1],1);
area2=readbin([exp2_path 'AREAtave.0000000010.data'],[20 16 1],1);
heff2=readbin([exp2_path 'HEFFtave.0000000010.data'],[20 16 1],1);

% plot comparison figures
clf, subplot(221)
mypcolor(lon,lat,area2'); caxis([0 1]), colorbar
title('Sea-ice concentration, no dynamics')
set(gca,'xticklabel',[])

subplot(222)
mypcolor(lon,lat,area2'-area1'); colorbar
title('Difference with LSR solver')
set(gca,'xticklabel',[])

subplot(223)
mypcolor(lon,lat,heff2'); caxis([0 .3]), colorbar
title('Effective sea-ice thickness (m)')
set(gca,'xticklabel',[])

subplot(224)
mypcolor(lon,lat,heff2'-heff1'); colorbar
title('Difference with LSR solver')
set(gca,'xticklabel',[])
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='../input/exp1/';
5	exp2_path='../input/exp6/';
6
7	% load monthly-mean SMMR-SSM/I data
8	load SSMI
9
10	% load model output
11	area1=readbin([exp1_path 'AREAtave.0000000010.data'],[20 16 1],1);
12	heff1=readbin([exp1_path 'HEFFtave.0000000010.data'],[20 16 1],1);
13	area2=readbin([exp2_path 'AREAtave.0000000010.data'],[20 16 1],1);
14	heff2=readbin([exp2_path 'HEFFtave.0000000010.data'],[20 16 1],1);
15
16	% plot comparison figures
17	clf, subplot(221)
18	mypcolor(lon,lat,area2'); caxis([0 1]), colorbar
19	title('Sea-ice concentration, no dynamics')
20	set(gca,'xticklabel',[])
21
22	subplot(222)
23	mypcolor(lon,lat,area2'-area1'); colorbar
24	title('Difference with LSR solver')
25	set(gca,'xticklabel',[])
26
27	subplot(223)
28	mypcolor(lon,lat,heff2'); caxis([0 .3]), colorbar
29	title('Effective sea-ice thickness (m)')
30	set(gca,'xticklabel',[])
31
32	subplot(224)
33	mypcolor(lon,lat,heff2'-heff1'); colorbar
34	title('Difference with LSR solver')
35	set(gca,'xticklabel',[])