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') |