lab_sea/matlab/lookat_exp4.m

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

% 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
area1=readbin([exp1_path 'AREAtave.0000000010.data'],[20 16 1],1);
heff1=readbin([exp1_path 'HEFFtave.0000000010.data'],[20 16 1],1);
uice1=readbin([exp1_path 'UICEtave.0000000010.data'],[20 16 1],1);
vice1=readbin([exp1_path 'VICEtave.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);
uice2=readbin([exp2_path 'UICEtave.0000000010.data'],[20 16 1],1);
vice2=readbin([exp2_path 'VICEtave.0000000010.data'],[20 16 1],1);

% plot comparison figures
clf, orient tall, wysiwyg
subplot(421), mypcolor(lon,lat,area2'); colorbar
title('Ice concentration during first 10 hours')
subplot(422), mypcolor(lon,lat,area2'-area1'); colorbar
title('2cpu - 1cpu solution')
subplot(423), mypcolor(lon,lat,heff2'); colorbar
title('Effective ice thickness')
subplot(424), mypcolor(lon,lat,heff2'-heff1'); colorbar
title('2cpu - 1cpu solution')
subplot(425), mypcolor(lon,lat,uice2'); colorbar
title('Uice')
subplot(426), mypcolor(lon,lat,uice2'-uice1'); colorbar
title('2cpu - 1cpu solution')
subplot(427), mypcolor(lon,lat,area2'); colorbar
title('Vice')
subplot(428), mypcolor(lon,lat,vice2'-vice1'); colorbar
title('2cpu - 1cpu solution')
1	heimbach	1.1.2.1	% need to be in verification/lab_sea/matlab directory
2			% and to specify location of exp1 and exp2 output
3			cd ../../../verification/lab_sea/matlab
4			exp1_path='../../../exe/exp1/';
5			exp2_path='../../../exe/exp4/';
6
7			% load monthly-mean SMMR-SSM/I data
8			load SSMI
9
10			% B-grid latitude for uice and vice
11			Blat=lat+1; Blon=lon+1;
12
13			% load model output
14			area1=readbin([exp1_path 'AREAtave.0000000010.data'],[20 16 1],1);
15			heff1=readbin([exp1_path 'HEFFtave.0000000010.data'],[20 16 1],1);
16			uice1=readbin([exp1_path 'UICEtave.0000000010.data'],[20 16 1],1);
17			vice1=readbin([exp1_path 'VICEtave.0000000010.data'],[20 16 1],1);
18			area2=readbin([exp2_path 'AREAtave.0000000010.data'],[20 16 1],1);
19			heff2=readbin([exp2_path 'HEFFtave.0000000010.data'],[20 16 1],1);
20			uice2=readbin([exp2_path 'UICEtave.0000000010.data'],[20 16 1],1);
21			vice2=readbin([exp2_path 'VICEtave.0000000010.data'],[20 16 1],1);
22
23			% plot comparison figures
24			clf, orient tall, wysiwyg
25			subplot(421), mypcolor(lon,lat,area2'); colorbar
26			title('Ice concentration during first 10 hours')
27			subplot(422), mypcolor(lon,lat,area2'-area1'); colorbar
28			title('2cpu - 1cpu solution')
29			subplot(423), mypcolor(lon,lat,heff2'); colorbar
30			title('Effective ice thickness')
31			subplot(424), mypcolor(lon,lat,heff2'-heff1'); colorbar
32			title('2cpu - 1cpu solution')
33			subplot(425), mypcolor(lon,lat,uice2'); colorbar
34			title('Uice')
35			subplot(426), mypcolor(lon,lat,uice2'-uice1'); colorbar
36			title('2cpu - 1cpu solution')
37			subplot(427), mypcolor(lon,lat,area2'); colorbar
38			title('Vice')
39			subplot(428), mypcolor(lon,lat,vice2'-vice1'); colorbar
40			title('2cpu - 1cpu solution')