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/exp3a/'; |
5 |
exp2_path='../input/exp3b/'; |
6 |
|
7 |
% load model output |
8 |
uice1=readbin([exp1_path 'UICEtave.0000000001.data'],[20 16 1],1); |
9 |
vice1=readbin([exp1_path 'VICEtave.0000000001.data'],[20 16 1],1); |
10 |
uice2=readbin([exp2_path 'UICEtave.0000000001.data'],[20 16 1],1); |
11 |
vice2=readbin([exp2_path 'VICEtave.0000000001.data'],[20 16 1],1); |
12 |
|
13 |
% plot comparison figures |
14 |
clf, subplot(221), mypcolor(uice1'); colorbar |
15 |
title('UICE using ADI solver') |
16 |
subplot(222), mypcolor(vice1'); colorbar |
17 |
title('VICE using ADI solver') |
18 |
subplot(223), mypcolor(uice2'); colorbar |
19 |
title('UICE using LSR solver') |
20 |
subplot(224), mypcolor(vice2'); colorbar |
21 |
title('VICE using LSR solver') |